A test of Hardy-Weinberg equilibrium in structured populations

The 5' untranslated region variant rs3811050 C/T of the interleukin-38 encoding gene is associated with susceptibility to rheumatoid arthritis in Iraqi women

BACKGROUND

Interleukin (IL)-38, the latest member of the IL-1 cytokine family, is proposed to have a pathogenic role in rheumatoid arthritis (RA). It is encoded by the IL1F10 gene, which harbors single nucleotide polymorphisms (SNPs) that may predict the risk of autoimmune diseases. Among them are 5' untranslated region (UTR) SNPs, which play a key role in post-transcriptional control, but have not been studied in Iraqi RA patients.

METHODS

Two novel IL1F10 5'UTR SNPs (rs3811050 C/T and rs3811051 T/G) were explored in RA and control women (n = 120 and 110, respectively). SNPs were genotyped using TaqMan assay. An ELISA kit was used to measure serum IL-38 concentrations.

RESULTS

A reduced risk of RA was associated with rs3811050 T allele and CT genotype (corrected probability [pc] = 0.01 and < 0.001, respectively), while there was no significant association with rs3811051. Haplotype analysis demonstrated that C-T haplotype was associated with a 1.65-fold greater risk of RA, whereas a reduced risk was linked to T-G haplotype. IL-38 concentrations were higher in patients than in controls (p < 0.001). In addition, IL-38 showed acceptable performance in distinguishing between RA and control women (p < 0.001). When IL-38 concentrations were stratified according to SNP genotypes, no significant differences were found.

CONCLUSIONS

The rs3811050 variant was more likely to affect RA susceptibility in Iraqi women, and the T allele may play a role in reducing disease risk. IL-38 concentrations were elevated in RA patients, but were not affected by the rs3811050 and rs3811051 genotypes.

Factors affecting HBV DNA suppression in chronic hepatitis B patients treated with tenofovir disoproxil fumarate

Background: This study aims to determine the factors affecting HBV DNA suppression in chronic hepatitis B patients with tenofovir disoproxil fumarate (TDF). Methods: A case-control was carried out from October 2021 to August 2022 on 182 chronic hepatitis B patients who had TDF therapy regularly for 24 weeks at H. Adam Malik and USU Hospitals in Medan, Indonesia. The history of the samples was obtained, followed by physical examination, and blood collection. CTLA-4 polymorphism examination was carried out using real-time PCR, while the serum CTLA-4 levels were assessed with ELISA. Results: The CTLA-4 -1661G>A polymorphism, genotype GG+AG, increased 1.52 times risk of not achieving HBV DNA suppression to TDF compared to genotype AA (p=0.041). High CTLA-4 levels increased 2.28 times risk, high HBV DNA levels increased 2.09 times risk, low ALT levels increased 1.95 times risk of not achieving HBV DNA suppression (p= 0.009, 0.026, 0.036, respectively). There was no relationship between gender, age, ethnicity, obesity, baseline AST, HBeAg, genotype, liver fibrosis and HBV DNA suppression after 24 weeks of treatment (p>0.05). Conclusions: The levels of CTLA-4, HBV DNA, ALT, and CTLA-4 -1661G>A polymorphism have a potential relationship with the suppression of HBV DNA in chronic hepatitis B patients with TDF.

Transmission Distortion of MCT1 rs1049434 among Polish Elite Athletes

Background: To date, nearly 300 genetic markers were linked to endurance and power/strength traits. The current study aimed to compare genotype distributions and allele frequencies of the common polymorphisms: MCT1 rs1049434, NRF2 rs12594956, MYBPC3 rs1052373 and HFE rs1799945 in Polish elite athletes versus nonathletes. Methods: The study involved 101 male elite Polish athletes and 41 healthy individuals from the Polish population as a control group. SNP data were extracted from whole-genome sequencing (WGS) performed using the following parameters: paired reads of 150 bps, at least 90 Gb of data per sample with 300 M reads and 30× mean coverage. Results: All the analyzed polymorphisms conformed to Hardy–Weinberg equilibrium (HWE) in athletes and the control group, except the MCT1 rs1049434, where allele T was over-represented in the elite trainers’ group. No significant between-group differences were found for analyzed polymorphisms. Conclusions: The MCT1 rs1049434 transmission distortion might be characteristic of Polish athletes and the effect of strict inclusion criteria. This result and the lack of statistically significant changes in the frequency of other polymorphisms between the groups might result from the small group size.

Robust, flexible, and scalable tests for Hardy-Weinberg equilibrium across diverse ancestries

Traditional Hardy-Weinberg equilibrium (HWE) tests (the χ2 test and the exact test) have long been used as a metric for evaluating genotype quality, as technical artifacts leading to incorrect genotype calls often can be identified as deviations from HWE. However, in data sets composed of individuals from diverse ancestries, HWE can be violated even without genotyping error, complicating the use of HWE testing to assess genotype data quality. In this manuscript, we present the Robust Unified Test for HWE (RUTH) to test for HWE while accounting for population structure and genotype uncertainty, and to evaluate the impact of population heterogeneity and genotype uncertainty on the standard HWE tests and alternative methods using simulated and real sequence data sets. Our results demonstrate that ignoring population structure or genotype uncertainty in HWE tests can inflate false-positive rates by many orders of magnitude. Our evaluations demonstrate different tradeoffs between false positives and statistical power across the methods, with RUTH consistently among the best across all evaluations. RUTH is implemented as a practical and scalable software tool to rapidly perform HWE tests across millions of markers and hundreds of thousands of individuals while supporting standard VCF/BCF formats. RUTH is publicly available at https://www.github.com/statgen/ruth.

Extending Tests of Hardy-Weinberg Equilibrium to Structured Populations

Testing for Hardy-Weinberg equilibrium (HWE) is an important component in almost all analyses of population genetic data. Genetic markers that violate HWE are often treated as special cases; for example, they may be flagged as possible genotyping errors, or they may be investigated more closely for evolutionary signatures of interest. The presence of population structure is one reason why genetic markers may fail a test of HWE. This is problematic because almost all natural populations studied in the modern setting show some degree of structure. Therefore, it is important to be able to detect deviations from HWE for reasons other than structure. To this end, we extend statistical tests of HWE to allow for population structure, which we call a test of "structural HWE." Additionally, our new test allows one to automatically choose tuning parameters and identify accurate models of structure. We demonstrate our approach on several important studies, provide theoretical justification for the test, and present empirical evidence for its utility. We anticipate the proposed test will be useful in a broad range of analyses of genome-wide population genetic data.

Testing for Hardy-Weinberg equilibrium in structured populations using genotype or low-depth next generation sequencing data

Testing for deviations from Hardy-Weinberg equilibrium (HWE) is a common practice for quality control in genetic studies. Variable sites violating HWE may be identified as technical errors in the sequencing or genotyping process, or they may be of particular evolutionary interest. Large-scale genetic studies based on next-generation sequencing (NGS) methods have become more prevalent as cost is decreasing but these methods are still associated with statistical uncertainty. The large-scale studies usually consist of samples from diverse ancestries that make the existence of some degree of population structure almost inevitable. Precautions are therefore needed when analysing these data set, as population structure causes deviations from HWE. Here we propose a method that takes population structure into account in the testing for HWE, such that other factors causing deviations from HWE can be detected. We show the effectiveness of PCAngsd in low-depth NGS data, as well as in genotype data, for both simulated and real data set, where the use of genotype likelihoods enables us to model the uncertainty.

Genetic effects of vascular endothelial growth factor and its receptor 2 on feather maturity in three chicken breeds

1. The goal of the current study was to evaluate the genetic effects of the vascular endothelial growth factor (VEGF) and its receptor (VEGFR-2) on feather maturity in the Qingyuan partridge chicken, Guangxi sanhuang chicken and Princess chicken. 2. Both SSCP-PCR and qPCR were employed to detect the polymorphism and gene expression of the VEGF and VEGFR-2 genes. 3. Four SNPs were identified in the VEGFR-2 gene. Exon10-A69G was associated with feather maturity (P < 0.01). Princess chickens with the genotype EF had higher feather maturity scores (P < 0.01). Higher expression levels of VEGF and VEGFR-2 were detected in the immature feather group of Qingyuan partridge chickens, especially in the skin. 4. The VEGF and VEGFR-2 genes play critical roles in feather maturity. In addition, exon10-A69G and genotype EF in the Princess chicken could potentially be utilised as genetic markers to improve efficiency in breeding.

Gene-Metabolite Interaction in the One Carbon Metabolism Pathway: Predictors of Colorectal Cancer in Multi-Ethnic Families

For personalized healthcare, the purpose of this study was to examine the key genes and metabolites in the one-carbon metabolism (OCM) pathway and their interactions as predictors of colorectal cancer (CRC) in multi-ethnic families. In this proof-of-concept study, we included a total of 30 participants, 15 CRC cases and 15 matched family/friends representing major ethnic groups in southern California. Analytics based on supervised machine learning were applied, with the target variable being specified as cancer, including the ensemble method and generalized regression (GR) prediction. Elastic Net with Akaike's Information Criterion with correction (AICc) and Leave-One-Out cross validation GR methods were used to validate the results for enhanced optimality, prediction, and reproducibility. The results revealed that despite some family members sharing genetic heritage, the CRC group had greater combined gene polymorphism-mutations than the family controls (p < 0.1) for five genes including MTHFR C677T, MTHFR A1298C, MTR A2756G, MTRR A66G, and DHFR 19bp. Blood metabolites including homocysteine (7 µmol/L), methyl-folate (40 nmol/L) with total gene mutations (≥4); age (51 years) and vegetable intake (2 cups), and interactions of gene mutations and methylmalonic acid (MMA) (400 nmol/L) were significant predictors (all p < 0.0001) using the AICc. The results were validated by a 3% misclassification rate, AICc of 26, and >99% area under the receiver operating characteristic curve. These results point to the important roles of blood metabolites as potential markers in the prevention of CRC. Future intervention studies can be designed to target the ways to mitigate the enzyme-metabolite deficiencies in the OCM pathway to prevent cancer.

Meta-Analysis of NOS3 G894T Polymorphisms with Air Pollution on the Risk of Ischemic Heart Disease Worldwide

The purpose of this updated meta-analysis was to investigate the effect of nitric oxide synthase-3 (NOS3) G894T polymorphisms, air pollution and their interaction on ischemic heart disease (IHD) risk across populations worldwide. Recursive partition trees, nonlinear association curve fit and geographic information system maps were incorporated to verify results of conventional pooled analyses for sources of heterogeneity. Results from 61 studies (16,219 cases, 12,222 controls) revealed a significant increased relative risk (RR) of IHD associated with NOS3 894 polymorphisms TT (RR = 1.44) and GT (RR = 1.37). Subgroup analysis revealed that the TT polymorphism genotype had significantly increased risk of IHD in Caucasian, East Asian, South Asian, and Middle Eastern populations (all p < 0.05). It is important to point out that many countries demonstrated an average risk of greater than two, which identifies the NOS3 894 TT polymorphism as a potential causal factor and biological marker of IHD, based on criteria for strong evidence used in international consensus panels. These 10 countries include Ukraine, the United Kingdom, Brazil, Chile, Japan, South Korea, India, Iran, Egypt and Morocco. For these countries with elevated risk (RR > 2) from the NOS3 894 TT polymorphism, meta-predictive analysis demonstrated an increasing trend in air pollution association with increased NOS3 894 polymorphisms. Further studies are needed to explore the complexity of the associations among NOS3 gene polymorphisms per population stratifications within countries, detailed air pollution data for added specificity for geographic location across time, and disease risk.

Lung cancer susceptibility from GSTM1 deletion and air pollution with smoking status: a meta-prediction of worldwide populations

Glutathione S transferase mu 1 (GSTM1) gene has been associated with lung cancer (LC) risk, for GSTM1 enzyme playing a vital role in detoxification pathway and protective against toxic insults. The major objective of this study was to investigate GSTM1 deletion pattern and its association with LC in the world’s population by using meta-prediction techniques. The secondary objective was to examine the effects of air pollution, smoking status, and other factors for gene-environment interactions with GSTM1 deletion and LC risk. We completed a comprehensive search to yield a total of 170 studies (40,296 cases and 48,346 controls) published from 1999 to 2017 for meta-analyses. The results revealed that GSTM1 deletion type was associated with increased risk of LC, while GSTM1 present type provided protective effect for all populations combined worldwide. Subgroup analysis on the rank order of risks from highest to lowest, among racial–ethnic groups, were Chinese, South East Asian, other North Asian, European, and finally American. Additional predictive analyses presented that air pollution played a significant role with increased risks of GSTM1 deletion and LC susceptibility, and the risks increased for smokers with higher levels of air pollution. Based on the findings of meta-predictive analysis, increased air pollution levels and smoking status presented additive effects to the LC risk susceptibilities and GSTM1 gene polymorphisms, for gene-environment interactions. Future studies are needed to examine gene-environment interactions for GSTM1 interacting with environmental factors and dietary interventions to mitigate the toxic effects, for LC prevention.

Gene-environment interactions and predictors of breast cancer in family-based multi-ethnic groups

Breast cancer (BC) is the most common cancer in women worldwide and second leading cause of cancer-related death. Understanding gene-environment interactions could play a critical role for next stage of BC prevention efforts. Hence, the purpose of this study was to examine the key gene-environmental factors affecting the risks of BC in a diverse sample. Five genes in one-carbon metabolism pathway including MTHFR 677, MTHFR 1298, MTR 2756, MTRR 66, and DHFR 19bp together with demographics, lifestyle, and dietary intake factors were examined in association with BC risks. A total of 80 participants (40 BC cases and 40 family/friend controls) in southern California were interviewed and provided salivary samples for genotyping. We presented the first study utilizing both conventional and new analytics including ensemble method and predictive modeling based on smallest errors to predict BC risks. Predictive modeling of Generalized Regression Elastic Net Leave-One-Out demonstrated alcohol use (p = 0.0126) and age (p < 0.0001) as significant predictors; and significant interactions were noted between body mass index (BMI) and alcohol use (p = 0.0027), and between BMI and MTR 2756 polymorphisms (p = 0.0090). Our findings identified the modifiable lifestyle factors in gene-environment interactions that are valuable for BC prevention.

Gene Environment Interactions and Predictors of Colorectal Cancer in Family-Based, Multi-Ethnic Groups

For the personalization of polygenic/omics-based health care, the purpose of this study was to examine the gene-environment interactions and predictors of colorectal cancer (CRC) by including five key genes in the one-carbon metabolism pathways. In this proof-of-concept study, we included a total of 54 families and 108 participants, 54 CRC cases and 54 matched family friends representing four major racial ethnic groups in southern California (White, Asian, Hispanics, and Black). We used three phases of data analytics, including exploratory, family-based analyses adjusting for the dependence within the family for sharing genetic heritage, the ensemble method, and generalized regression models for predictive modeling with a machine learning validation procedure to validate the results for enhanced prediction and reproducibility. The results revealed that despite the family members sharing genetic heritage, the CRC group had greater combined gene polymorphism rates than the family controls (p < 0.05), on MTHFR C677T, MTR A2756G, MTRR A66G, and DHFR 19 bp except MTHFR A1298C. Four racial groups presented different polymorphism rates for four genes (all p < 0.05) except MTHFR A1298C. Following the ensemble method, the most influential factors were identified, and the best predictive models were generated by using the generalized regression models, with Akaike's information criterion and leave-one-out cross validation methods. Body mass index (BMI) and gender were consistent predictors of CRC for both models when individual genes versus total polymorphism counts were used, and alcohol use was interactive with BMI status. Body mass index status was also interactive with both gender and MTHFR C677T gene polymorphism, and the exposure to environmental pollutants was an additional predictor. These results point to the important roles of environmental and modifiable factors in relation to gene-environment interactions in the prevention of CRC.

Meta-Prediction of MTHFR Gene Polymorphisms and Air Pollution on the Risk of Hypertensive Disorders in Pregnancy Worldwide

Hypertensive disorders in pregnancy (HDP) are devastating health hazards for both women and children. Both methylenetetrahydrofolate reductase (MTHFR) gene polymorphisms and air pollution can affect health status and result in increased risk of HDP for women. The major objective of this study was to investigate the effect of MTHFR polymorphisms, air pollution, and their interaction on the risk of HDP by using meta-predictive analytics. We searched various databases comprehensively to access all available studies conducted for various ethnic populations from countries worldwide, from 1997 to 2017. Seventy-one studies with 8064 cases and 13,232 controls for MTHFR C677T and 11 studies with 1425 cases and 1859 controls for MTHFR A1298C were included. MTHFR C677T homozygous TT (risk ratio (RR) = 1.28, p < 0.0001) and CT plus TT (RR = 1.07, p = 0.0002) were the risk genotypes, while wild-type CC played a protective role (RR = 0.94, p = 0.0017) for HDP. The meta-predictive analysis found that the percentage of MTHFR C677T TT plus CT (p = 0.044) and CT (p = 0.043) genotypes in the HDP case group were significantly increased with elevated levels of air pollution worldwide. Additionally, in countries with higher air pollution levels, the pregnant women with wild-type CC MTHFR 677 had a protection effect against HDP (p = 0.014), whereas, the homozygous TT of MTHFR C677T polymorphism was a risk genotype for developing HDP. Air pollution level is an environmental factor interacting with increased MTHFR C677T polymorphisms, impacting the susceptibility of HDP for women.

Population Structure, Diversity and Reproductive Mode of the Grape Phylloxera (Daktulosphaira vitifoliae) across Its Native Range

Grape Phylloxera, Daktulosphaira vitifoliae, is a gall-forming insect that feeds on the leaves and roots of many Vitis species. The roots of the cultivated V. vinifera cultivars and hybrids are highly susceptible to grape phylloxera feeding damage. The native range of this insect covers most of North America, and it is particularly abundant in the eastern and central United States. Phylloxera was introduced from North America to almost all grape-growing regions across five of the temperate zone continents. It devastated vineyards in each of these regions causing large-scale disruptions to grape growers, wine makers and national economies. In order to understand the population diversity of grape phylloxera in its native range, more than 500 samples from 19 States and 34 samples from the introduced range (northern California, Europe and South America) were genotyped with 32 simple sequence repeat markers. STRUCTURE, a model based clustering method identified five populations within these samples. The five populations were confirmed by a neighbor-joining tree and principal coordinate analysis (PCoA). These populations were distinguished by their Vitis species hosts and their geographic locations. Samples collected from California, Europe and South America traced back to phylloxera sampled in the northeastern United States on V. riparia, with some influence from phylloxera collected along the Atlantic Coast and Central Plains on V. vulpina. Reproductive statistics conclusively confirmed that sexual reproduction is common in the native range and is combined with cyclical parthenogenesis. Native grape phylloxera populations were identified to be under Hardy-Weinberg equilibrium. The identification of admixed samples between many of these populations indicates that shared environments facilitate sexual reproduction between different host associated populations to create new genotypes of phylloxera. This study also found that assortative mating might occur across the sympatric range of the V. vulpina west and V. cinerea populations.

Meta-Prediction of MTHFR Gene Polymorphism Mutations and Associated Risk for Colorectal Cancer

The methylenetetrahydrofolate reductase (MTHFR) gene is one of the most investigated of the genes associated with chronic human diseases because of its associations with hyperhomocysteinemia and toxicity. It has been proposed as a prototype gene for the prevention of colorectal cancer (CRC). The major objectives of this meta-analysis were to examine the polymorphism-mutation patterns of MTHFR and their associations with risk for CRC as well as potential contributing factors for mutations and disease risks. This analysis included 33,626 CRC cases and 48,688 controls across 92 studies for MTHFR 677 and 16,367 cases and 24,874 controls across 54 studies for MTHFR 1298, comprising data for various racial and ethnic groups, both genders, and multiple cancer sites. MTHFR 677 homozygous TT genotype was protective (p < .05) for CRC for all included populations; however, with heterogeneity across various racial–ethnic groups and opposing findings, it was a risk genotype for the subgroup of Hispanics (p < .01). Additional countries for which subgroup analyses resulted in 677 TT as a risk genotype included Turkey, Romania, Croatia, Hungary, Portugal, Mexico, Brazil, U.S. Hawai’i, Taiwan, India, and Egypt. Countries with the highest mutation rates and risks for both MTHFR 677 and 1298 genotypes are presented using global maps to visualize the grouping patterns. Meta-predictive analyses revealed that air pollution levels were associated with gene polymorphisms for both genotypes. Future nursing research should be conducted to develop proactive measures to protect populations in cities where air pollution causes more deaths.

Polymorphisms in the FK506 binding protein 5 gene are associated with attention deficit hyperactivity disorder and diurnal cortisol levels

AIM

Previous studies have shown an association between childhood attention deficit hyperactivity disorder (ADHD) and a down-regulated hypothalamus-pituitary-adrenal axis (HPA axis) with low diurnal cortisol levels. Given the role of the FK506 binding protein 5 (FKBP5) as an important regulator of the negative feedback system of the HPA axis, we set out to investigate possible associations between single nucleotide polymorphisms (SNPs) in FKBP5 in relation to ADHD and diurnal cortisol levels.

METHODS

Children with ADHD (n = 81) and healthy comparisons (n = 88) collected saliva four times during a regular school day for radioimmunoassay analysis of cortisol and for genotyping of five SNPs in FKBP5 (rs9296158, rs1360780, rs9470080, rs7748266 and rs9394309).

RESULTS

We found associations between SNP genotypes and ADHD as well as between genotypes and diurnal cortisol levels. One of these SNPs, rs9470080, was significantly associated with both ADHD and lower cortisol levels.

CONCLUSION

This study contributes to previous findings on a down-regulated HPA axis in children with ADHD by demonstrating an association between ADHD, lower cortisol levels and SNPs of the FKBP5-gene. The relevance of these findings for the development and shaping of ADHD symptoms needs to be approached in larger samples, preferably also taking stress reactivity into consideration.

Contributory role of five common polymorphisms of RAGE and APE1 genes in lung cancer among Han Chinese

BACKGROUND

Lung cancer is the leading cause of cancer mortality in China. Given the ubiquitous nature of gene-to-gene interaction in lung carcinogenesis, we sought to evaluate five common polymorphisms from advanced glycosylation end product-specific receptor (RAGE) and apurinic/apyrimidinic endonuclease 1 (APE1) genes in association with lung cancer among Han Chinese.

METHODS AND RESULTS

819 patients with lung cancer and 803 cancer-free controls were recruited from Qiqihar city. Genotypes of five examined polymorphisms (RAGE gene: rs1800625, rs1800624, rs2070600; APE1 gene: rs1760944, rs1130409) were determined by ligase detection reaction method. Data were analyzed by R software and multifactor dimensionality reduction (MDR). Hardy-Weinberg equilibrium was satisfied for all five polymorphisms. Overall differences in the genotype and allele distributions were significant for rs1800625 (Pgenotype<0.0005; Pallele<0.0005), rs2070600 (Pgenotype = 0.005; Pallele = 0.004) and rs1130409 (Pgenotype = 0.009; Pallele = 0.004) polymorphisms. Haplotype C-A-A (alleles in order of rs1800625, rs1800624 and rs2070600) of RAGE gene was overrepresented in patients, and conferred a 2.1-fold increased risk of lung cancer (95% confidence interval: 1.52-2.91), independent of confounding factors. Further application of MDR method to five examined polymorphisms identified the overall best interaction model including rs2070600 and rs1130409 polymorphisms. This model had a maximal testing accuracy of 64.63% and a maximal cross-validation consistency of 9 out of 10 at the significant level of 0.006.

CONCLUSIONS

Our findings demonstrated a potential interactive contribution of RAGE and APE1 genes to the pathogenesis of lung cancer among Han Chinese. Further studies are warranted to confirm or refute these findings.