Genomic analysis uncovers novel candidate genes related to adaptation to tropical climates and milk production traits in native goats.

BACKGROUND: Since domestication, both evolutionary forces and human selection have played crucial roles in producing adaptive and economic traits, resulting in animal breeds that have been selected for specific climates and different breeding goals. Pakistani goat breeds have acquired genomic adaptations to their native climate conditions, such as tropical and hot climates. In this study, using next-generation sequencing data, we aimed to assess the signatures of positive selection in three native Pakistani goats, known as milk production breeds, that have been well adapted to their local climate. RESULTS: To explore the genomic relationship between studied goat populations and their population structure, whole genome sequence data from native goat populations in Pakistan (n = 26) was merged with available worldwide goat genomic data (n = 184), resulting in a total dataset of 210 individuals. The results showed a high genetic correlation between Pakistani goats and samples from North-East Asia. Across all populations analyzed, a higher linkage disequilibrium (LD) level (- 0.59) was found in the Pakistani goat group at a genomic distance of 1 Kb. Our findings from admixture analysis (K = 5 and K = 6) showed no evidence of shared genomic ancestry between Pakistani goats and other goat populations from Asia. The results from genomic selection analysis revealed several candidate genes related to adaptation to tropical/hot climates (such as; KITLG, HSPB9, HSP70, HSPA12B, and HSPA12B) and milk production related-traits (such as IGFBP3, LPL, LEPR, TSHR, and ACACA) in Pakistani native goat breeds. CONCLUSIONS: The results from this study shed light on the structural variation in the DNA of the three native Pakistani goat breeds. Several candidate genes were discovered for adaptation to tropical/hot climates, immune responses, and milk production traits. The identified genes could be exploited in goat breeding programs to select efficient breeds for tropical/hot climate regions.

Genome-wide association study on blood pressure traits in the Iranian population suggests ZBED9 as a new locus for hypertension.

High blood pressure is the heritable risk factor for cardiovascular and kidney diseases. Genome-wide association studies(GWAS) on blood pressure traits increase our understanding of its underlying genetic basis. However, a large proportion of GWAS was conducted in Europeans, and some roadblocks deprive other populations to benefit from their results. Iranians population with a high degree of genomic specificity has not been represented in international databases to date, so to fill the gap, we explored the effects of 652,919 genomic variants on Systolic Blood Pressure (SBP), Diastolic Blood Pressure (DBP), and Hypertension (HTN) in 7694 Iranian adults aged 18 and over from Tehran Cardiometabolic Genetic Study (TCGS). We identified consistent signals on ZBED9 associated with HTN in the genome-wide borderline threshold after adjusting for different sets of environmental predictors. Moreover, strong signals on ABHD17C and suggestive signals on FBN1 were detected for DBP and SBP, respectively, while these signals were not consistent in different GWA analysis. Our finding on ZBED9 was confirmed for all BP traits by linkage analysis in an independent sample. We found significant associations with similar direction of effects and allele frequency of genetic variants on ZBED9 with DBP (genome-wide threshold) and HTN (nominal threshold) in GWAS summary data of UK Biobank. Although there is no strong evidence to support the function of ZBED9 in blood pressure regulation, it provides new insight into the pleiotropic effects of hypertension and other cardiovascular diseases.

Familial genetic and environmental risk profile and high blood pressure event: a prospective cohort of cardio-metabolic and genetic study.

BACKGROUND AND AIMS: High blood pressure is the heritable risk factor for cardiovascular diseases. We investigated whether the presence of familial genetic and environmental risk factors are associated with increased risk of high blood pressure. METHODS: A total of 4,559 individuals from 401 families were included in this study. Familial aggregation analysis was carried out on systolic blood pressure (SBP), diastolic blood pressure (DBP), body mass index (BMI) and waist circumference (WC), and heritability was estimated for SBP and DBP. The association between familial risk factors and blood pressure traits including, incidence of hypertension, SBP and DBP was estimated separately using regression-based two-level Haseman-Elston (HE) method, with individual and familial BMI and WC as environmental exposures and familial genetic profile of known variants as genetic risk factors in 210 index families (≥2 hypertensive cases). Models were adjusted for the two nested sets of covariates. RESULTS: During a follow-up of 15 years, the SBP, DBP, BMI and WC were highly correlated in inter class of mother-offspring and intraclass of sister-sister with heritability of 30 and 25% for DBP and SBP, respectively. Among index families, those whose members with higher familial BMI or WC had significantly increased risk of hypertension and consistent, strong signals of rs2493134 (AGT) linked with SBP and DBP, rs976683 (NLGN1) linked with SBP and HTN, and epistasis of rs2021783 (TNXB) and known genetic variants linked with all blood pressure traits. CONCLUSIONS: Findings from this study show that familial genetic and environmental risk profile increase risk for high blood pressure beyond the effect of the individuals' own risk factors.

Prediction of biological age and evaluation of genome-wide dynamic methylomic changes throughout human aging.

The use of DNA methylation signatures to predict chronological age and aging rate is of interest in many fields, including disease prevention and treatment, forensics, and anti-aging medicine. Although a large number of methylation markers are significantly associated with age, most age-prediction methods use a few markers selected based on either previously published studies or datasets containing methylation information. Here, we implemented reproducing kernel Hilbert spaces (RKHS) regression and a ridge regression model in a Bayesian framework that utilized phenotypic and methylation profiles simultaneously to predict chronological age. We used over 450,000 CpG sites from the whole blood of a large cohort of 4409 human individuals with a range of 10-101 years of age. Models were fitted using adjusted and un-adjusted methylation measurements for cell heterogeneity. Un-adjusted methylation scores delivered a significantly higher prediction accuracy than adjusted methylation data, with a correlation between age and predicted age of 0.98 and a root mean square error (RMSE) of 3.54 years in un-adjusted data, and 0.90 (correlation) and 7.16 (RMSE) years in adjusted data. Reducing the number of predictors (CpG sites) through subset selection improved predictive power with a correlation of 0.98 and an RMSE of 2.98 years in the RKHS model. We found distinct global methylation patterns, with a significant increase in the proportion of methylated cytosines in CpG islands and a decreased proportion in other CpG types, including CpG shore, shelf, and open sea (P < 5e-06). Epigenetic drift seemed to be a widespread phenomenon as more than 97% of the age-associated methylation sites had heteroscedasticity. Apparent methylomic aging rate (AMAR) had a sex-specific pattern, with an increase in AMAR in females with age related to males.

Integration of single nucleotide variants and whole-genome DNA methylation profiles for classification of rheumatoid arthritis cases from controls.

This study evaluated the use of multiomics data for classification accuracy of rheumatoid arthritis (RA). Three approaches were used and compared in terms of prediction accuracy: (1) whole-genome prediction (WGP) using SNP marker information only, (2) whole-methylome prediction (WMP) using methylation profiles only, and (3) whole-genome/methylome prediction (WGMP) with combining both omics layers. The number of SNP and of methylation sites varied in each scenario, with either 1, 10, or 50% of these preselected based on four approaches: randomly, evenly spaced, lowest p value (genome-wide association or epigenome-wide association study), and estimated effect size using a Bayesian ridge regression (BRR) model. To remove effects of high levels of pairwise linkage disequilibrium (LD), SNPs were also preselected with an LD-pruning method. Five Bayesian regression models were studied for classification, including BRR, Bayes-A, Bayes-B, Bayes-C, and the Bayesian LASSO. Adjusting methylation profiles for cellular heterogeneity within whole blood samples had a detrimental effect on the classification ability of the models. Overall, WGMP using Bayes-B model has the best performance. In particular, selecting SNPs based on LD-pruning with 1% of the methylation sites selected based on BRR included in the model, and fitting the most significant SNP as a fixed effect was the best method for predicting disease risk with a classification accuracy of 0.975. Our results showed that multiomics data can be used to effectively predict the risk of RA and identify cases in early stages to prevent or alter disease progression via appropriate interventions.

Including Phenotypic Causal Networks in Genome-Wide Association Studies Using Mixed Effects Structural Equation Models.

Network based statistical models accounting for putative causal relationships among multiple phenotypes can be used to infer single-nucleotide polymorphism (SNP) effect which transmitting through a given causal path in genome-wide association studies (GWAS). In GWAS with multiple phenotypes, reconstructing underlying causal structures among traits and SNPs using a single statistical framework is essential for understanding the entirety of genotype-phenotype maps. A structural equation model (SEM) can be used for such purposes. We applied SEM to GWAS (SEM-GWAS) in chickens, taking into account putative causal relationships among breast meat (BM), body weight (BW), hen-house production (HHP), and SNPs. We assessed the performance of SEM-GWAS by comparing the model results with those obtained from traditional multi-trait association analyses (MTM-GWAS). Three different putative causal path diagrams were inferred from highest posterior density (HPD) intervals of 0.75, 0.85, and 0.95 using the inductive causation algorithm. A positive path coefficient was estimated for BM → BW, and negative values were obtained for BM → HHP and BW → HHP in all implemented scenarios. Further, the application of SEM-GWAS enabled the decomposition of SNP effects into direct, indirect, and total effects, identifying whether a SNP effect is acting directly or indirectly on a given trait. In contrast, MTM-GWAS only captured overall genetic effects on traits, which is equivalent to combining the direct and indirect SNP effects from SEM-GWAS. Although MTM-GWAS and SEM-GWAS use the similar probabilistic models, we provide evidence that SEM-GWAS captures complex relationships in terms of causal meaning and mediation and delivers a more comprehensive understanding of SNP effects compared to MTM-GWAS. Our results showed that SEM-GWAS provides important insight regarding the mechanism by which identified SNPs control traits by partitioning them into direct, indirect, and total SNP effects.

In vitro maturation of ovine oocyte in a modified granulosa cells co-culture system and alpha-tocopherol supplementation: effects on nuclear maturation and cleavage.

This study was designed to investigate the effects of α-tocopherol and granulosa cells monolayer on nuclear maturation and cleavage rates of ovine cumulus-oocyte complexes (COCs). The COCs (n = 2814) were matured in maturation medium supplemented with various concentration of α-tocopherol (0, 5, 10, 15 µg/ml), oocytes were incubated at 39 °C with 5 % CO2 for 24 h in three culture systems: (a) maturation medium (MM; n = 884), (b) co-cultured with granulosa cells (CG; n = 982) and (c) co-cultured with granulosa cells and cells were further cultured in MM for 12 h (CG + 12hMM; n = 948). Our results showed that α-tocopherol had no effect on GVBD and MII as compared to control group, but when α-tocopherol added to maturation medium the rate of cleavage decreased. This indicates interaction of above mentioned factors in any of the treatments showed no significant differences on the rate of maturation and cleavage stages (MII, GVBD and cleavage) (p > 0.05). The oocytes co-cultured with granulosa cells for 24 h had beneficial effects on cleavage rate. The maximum MII and cleavage rates were achieved when oocytes had extra 12 h culture in the maturation medium without granulosa cells. Results also showed our modified co-culture system (CG + 12hMM), improved rates of MII and the cleavage in comparison with other studied maturation systems.