A genome-wide study of Hardy-Weinberg equilibrium with next generation sequence data

Association between vitamin D metabolism gene polymorphisms and schizophrenia

Schizophrenia (SZ) is a multifactorial and neurodegenerative disorder that results from the interaction between genetic and environmental factors. Notably, hundreds of single nucleotide polymorphisms (SNPs) are associated with the susceptibility to SZ. Vitamin D (VD) plays an essential role in regulating several genes important for maintaining brain function and health. To the best of the authors' knowledge, no studies have yet been conducted on the association between the VD pathway and patients with SZ. Therefore, the present study aimed to assess the potential association between eight SNPs in genes related to the VD pathway, including CYP2R1, CYP27B1, CYP24A1 and VDR among patients with SZ. A case-control study was conducted, involving a total of 400 blood samples drawn from 200 patients and 200 healthy controls. Genomic DNA was extracted and variants were genotyped using the tetra-amplification refractory mutation system-polymerase chain reaction method. The present study revealed statistically significant differences between patients with SZ and controls regarding the genotypes and allele distributions of three SNPs [CYP2R1 (rs10741657), CYP27B1 (rs10877012) and CYP24A1 (rs6013897) (P<0.0001)]. The AA genotype of rs10741657 was identified to be associated with SZ (P<0.0001) and the frequency of the A allele was higher in patients with SZ (P<0.0001) compared with the control group. Similarly, the TT genotype of rs10877012 was revealed to be associated with SZ (P<0.0001) and the T allele was more frequent in patients with SZ (P<0.0001) than in the control group. Moreover, the AA genotype of rs6013897 was revealed to be associated with SZ (P<0.0001), although no significant difference was detected between the two groups regarding the A allele (P=0.055). VDR (rs2228570, rs1544410, rs731236 and rs7975232) and CYP27B1 (rs4646536) gene polymorphisms did not exhibit a significant association with SZ. While the studied SNPs revealed promising discriminatory capacity between patients with SZ and controls, the rs10741657 SNP exhibited the most optimal area under the curve value at 0.615. A logistic model was applied considering only the significant SNPs and VD levels, which revealed that rs6013897 (T/A) and VD may have protective effects (0.267, P<0.001; 0.888, P<0.001, respectively). Moreover, a low serum VD level was highly prevalent in patients with SZ compared with the controls. Based on this finding, an association between serum 25(OH)D and SZ could be demonstrated. The present study revealed that CYP2R1 (rs10741657), CYP27B1 (rs10877012) and CYP24A1 (rs6013897) gene SNPs may be associated with SZ susceptibility.

Next-generation data filtering in the genomics era

Genomic data are ubiquitous across disciplines, from agriculture to biodiversity, ecology, evolution and human health. However, these datasets often contain noise or errors and are missing information that can affect the accuracy and reliability of subsequent computational analyses and conclusions. A key step in genomic data analysis is filtering - removing sequencing bases, reads, genetic variants and/or individuals from a dataset - to improve data quality for downstream analyses. Researchers are confronted with a multitude of choices when filtering genomic data; they must choose which filters to apply and select appropriate thresholds. To help usher in the next generation of genomic data filtering, we review and suggest best practices to improve the implementation, reproducibility and reporting standards for filter types and thresholds commonly applied to genomic datasets. We focus mainly on filters for minor allele frequency, missing data per individual or per locus, linkage disequilibrium and Hardy-Weinberg deviations. Using simulated and empirical datasets, we illustrate the large effects of different filtering thresholds on common population genetics statistics, such as Tajima's D value, population differentiation (FST), nucleotide diversity (π) and effective population size (Ne).

Association of OX40L gene polymorphism with multiple sclerosis in Iranians

Introduction

The exact etiology of multiple sclerosis is unknown but recent studies indicated a link between tumor necrosis factor superfamily member 4 and the disease. Polymorphisms located in the regulatory region of the gene may affect its phenotype. Hence, we aimed to investigate the association of promoter polymorphisms of the gene with multiple sclerosis and also to estimate the frequency of haplotypes in the patients and healthy subjects.

Methods

Two hundred age- and sex-matched subjects including 100 patients and 100 healthy subjects were investigated in the study. Genotype and allele distributions of rs3850641, rs1234313, and rs10912580 polymorphisms in the promoter region of the gene were investigated by polymerase chain reaction-restriction fragment length polymorphism. In addition, haplotype frequencies estimation and linkage disequilibrium analysis were performed by SNPStats web tool.

Results

The distribution of AA, AG and GG genotypes of rs3850641 was significantly different between the patient and healthy groups (P = 0.009). In addition, frequencies of A and G alleles of rs3850641 were different between the groups (P < 0.001). Also the distribution of rs3850641 genotypes was different between the women of the both groups (P = 0.007). Our analysis revealed that rs3850641 AG (Odds ratio = 0.393, 95 % confidence interval = 0.170-0.907, P = 0.029) and GG (Odds ratio = 0.373, 95 % confidence interval = 0.168-0.830, P = 0.016) genotypes were associated with decreased risk of the disease. However, rs1234313 genotype and allele distributions were not different between the groups. The distribution of rs10912580polymorphism. AA, AG, and GG genotypes was significantly different between the groups (P = 0.007). rs10912580 AG genotype was associated with low risk of the disease (Odds ratio = 0.252, 95 % confidence interval = 0.102-0.623, P = 0.003). The distribution of haplotypes was statistically different between the patient and healthy groups (P < 0.001). A-G-A was the most frequent haplotype among the patients and the estimated frequency was higher than that of the control group (0.5527 versus 0.3739).

Conclusion

The distribution of rs3850641 and rs10912580 genotypes was different between the patients and healthy subjects. Moreover, rs3850641 AG and GG genotypes and also rs10912580 AG genotype were associated with low risk of the disease in Iranians. Further studies with large groups are recommended to show whether genotype variation in the patients could alter the response to treatment or not.

Towards pharmacogenomics-guided tuberculosis (TB) therapy: N-acetyltransferase-2 genotypes among TB-infected Kenyans of mixed ethnicity

BACKGROUND

Though persons of African descent have one of the widest genetic variability, genetic polymorphisms of drug-metabolising enzymes such as N-Acetyltransferase-2 (NAT2) are understudied. This study aimed to identify prevalent NAT2 single nucleotide polymorphisms (SNPs) and infer their potential effects on enzyme function among Kenyan volunteers with tuberculosis (TB) infection. Genotypic distribution at each SNP and non-random association of alleles were evaluated by testing for Hardy-Weinberg Equilibrium (HWE) and Linkage Disequilibrium (LD).

METHODS

We isolated genomic DNA from cryopreserved Peripheral Blood Mononuclear Cells of 79 volunteers. We amplified the protein-coding region of the NAT2 gene by polymerase chain reaction (PCR) and sequenced PCR products using the Sanger sequencing method. Sequencing reads were mapped and aligned to the NAT2 reference using the Geneious software (Auckland, New Zealand). Statistical analyses were performed using RStudio version 4.3.2 (2023.09.1 + 494).

RESULTS

The most frequent haplotype was the wild type NAT2*4 (37%). Five genetic variants: 282C > T (NAT2*13), 341 T > C (NAT2*5), 803A > G (NAT2*12), 590G > A (NAT2*6) and 481C > T (NAT2*11) were observed with allele frequencies of 29%, 18%, 6%, 6%, and 4% respectively. According to the bimodal distribution of acetylation activity, the predicted phenotype was 76% rapid (mainly consisting of the wildtype NAT2*4 and the NAT2*13A variant). A higher proportion of rapid acetylators were female, 72% vs 28% male (p = 0.022, odds ratio [OR] 3.48, 95% confidence interval [CI] 1.21 to 10.48). All variants were in HWE. NAT2 341 T > C was in strong complete LD with the 590G > A variant (D' = 1.0, r2 = - 0.39) but not complete LD with the 282C > T variant (D' = 0.94, r2 = - 0.54).

CONCLUSION

The rapid acetylation haplotypes predominated. Despite the LD observed, none of the SNPs could be termed tag SNP. This study adds to the genetic characterisation data of African populations at NAT2, which may be useful for developing relevant pharmacogenomic tools for TB therapy. To support optimised, pharmacogenomics-guided TB therapy, we recommend genotype-phenotype studies, including studies designed to explore gender-associated differences.

Pharmacogenetic analysis of inter-ethnic variability in the uptake transporter SLCO1B1 gene in Colombian, Mozambican, and Portuguese populations

BACKGROUND

Statin-induced myopathy is reported to be associated with the solute carrier organic anion transporter family member 1B1 gene single nucleotide polymorphism, c.521 T > C. There is no epidemiologic data on this gene polymorphism in several countries. Therefore, this study aimed at assessing the genotype and allele frequencies of the gene variant in three countries.

METHODS

This study involved healthy individuals from Colombia, Mozambique, and Portugal. Genomic DNA was isolated from blood samples using the Qiamp DNA Extraction Kit (Qiagen). The isolated DNA was genotyped using novel Polymerase Chain Reaction-Restriction Fragment Length Polymorphism. Microstat and GraphPad QuickCal software were used for the Chi-square test and the evaluation of Hardy-Weinberg equilibrium respectively.

RESULTS

A total of 181 individuals' blood samples were analyzed. Overall, the TT (74.0%) genotype was the highest and the CC (7.8%) was the lowest. Country wise genotypic frequencies were Colombia 47(70.2%) TT, 12(17.9%) TC and 8(11.9%) CC; Mozambique 47(88.7%) TT, 5(9.4%) TC, and 1(1.9%) CC; and Portugal 40(65.6%) TT, 16(26.2%) TC, and 5(8.2%) CC. The reference (T) allele was highest among Mozambicans (93.4%) compared to Colombians (79.1%) and Portuguese (78.7%). Mozambicans showed statistically significant genotypic and allelic frequency differences compared to Colombians (p < 0.01) and Portuguese (p < 0.01).

CONCLUSIONS

Overall and country-wise, CC genotype was less frequent and it is relatively high for Colombians and Portuguese populations. This finding may imply statins risk-benefit variability associated with CC genotype among these populations that needs further understanding.

The genetic variant rs55986091 HLA-DQB1 is associated with a protective effect against cervical cancer

Introduction

Cervical cancer (CC) is a prevalent malignancy affecting women globally. The primary causative factor of CC is the high-risk oncogenic human papillomavirus (HR-HPV). However, it is noteworthy that not all women infected with HR-HPV develop cancer, indicating the potential involvement of genetic predisposition in the development of CC. This study aims to identify genetic risks and their distribution in groups of women with different epidemiological features of HR-HPV.

Materials and methods

A comparison was conducted among four groups of women, comprising 218 HPV-negative women, 120 HPV-positive women, 191 women diagnosed with cervical intraepithelial neoplasia (CIN) grade 2 or 3, and 124 women diagnosed with CC. The analysis focused on four single nucleotide polymorphisms (SNPs): rs55986091 in HLA-DQB1, rs138446575 in TTC34, rs1048943 in CYP1A1, and rs2910164 in miRNA-146a.

Results

The rs55986091-A allele exhibited a protective effect within the "CC" group when compared to the "HPV-Negative" group (OR = 0.4, 95% CI= 0.25-0.65) using a log-additive model. Additionally, similar protective effects were observed in the "CIN 2/3" group compared to the "HPV-Negative" group (OR = 0.47, 95% CI = 0.28-0.79).

Conclusion

The data obtained emphasize the importance of developing PCR-based diagnostic kits for the identification of SNP alleles, particularly for rs55986091, among HR-HPV-positive women within the Russian population.

Estimates of Genomic Heritability and the Marker-Derived Gene for Re(Production) Traits in Xinggao Sheep

Xinggao sheep are a breed of Chinese domestic sheep that are adapted to the extremely cold climatic features of the Hinggan League in China. The economically vital reproductive trait of ewes (litter size, LS) and productive traits of lambs (birth weight, BWT; weaning weight, WWT; and average daily gain, ADG) are expressed in females and later in life after most of the selection decisions have been made. This study estimated the genetic parameters for four traits to explore the genetic mechanisms underlying the variation, and we performed genome-wide association study (GWAS) tests on a small sample size to identify novel marker trait associations (MTAs) associated with prolificacy and growth. We detected two suggestive significant single-nucleotide polymorphisms (SNPs) associated with LS and eight significant SNPs for BWT, WWT, and ADG. These candidate loci and genes also provide valuable information for further fine-mapping of QTLs and improvement of reproductive and productive traits in sheep.

Association of GSTTI, M1 and Polymorphism in GSTPI with Chronic Periodontal Disease in a Pakistani Population

OBJECTIVE

Chronic periodontal disease (CP) is a multifactorial infectious and inflammatory disease that occurs due to the challenge between the immune response of the host and specific periodontal bacteria, and that can lead to tooth loss due to damage inflicted to the supporting tissue. The current study investigates the genotypes of the GSTM1 and GSTT1 genes, along with the allelic frequency of the single nucleotide polymorphism [SNP; rs1695] in the GSTP1 gene and correlates them individually or in various combinations with the incidence of CP.

METHODS

A total of 203 clinically confirmed CP patients and 201 control subjects were enrolled from Multan and Dera Ghazi Khan Districts in Pakistan from April to July 2022. Multiplex Polymerase Chain Reaction (PCR) and tetra-primer amplification refractory mutation system-polymerase chain reaction (T-ARMS-PCR) approaches were applied to determine the genotypes of the studied GSTs. The association of rs1695 in GSTP1 with CP was studied both individually and in various combinations with GSTM1 and T1.

RESULTS

The absence of GSTM1, the presence of GSTT1 and the presence of the mutant allele (G) at rs1695 in GSTP1 were found to be significantly associated with CP. Patients aged between 10 and 30 years were more affected by CP.

CONCLUSION

Our results indicate that the genotypes of the analyzed GSTs affect the levels of protection from oxidative stress and may therefore influence the disease progression in CP.

Why does the X chromosome lag behind autosomes in GWAS findings?

The X-chromosome is among the largest human chromosomes. It differs from autosomes by a number of important features including hemizygosity in males, an almost complete inactivation of one copy in females, and unique patterns of recombination. We used data from the Catalog of Published Genome Wide Association Studies to compare densities of the GWAS-detected SNPs on the X-chromosome and autosomes. The density of GWAS-detected SNPs on the X-chromosome is 6-fold lower compared to the density of the GWAS-detected SNPs on autosomes. Differences between the X-chromosome and autosomes cannot be explained by differences in the overall SNP density, lower X-chromosome coverage by genotyping platforms or low call rate of X-chromosomal SNPs. Similar differences in the density of GWAS-detected SNPs were found in female-only GWASs (e.g. ovarian cancer GWASs). We hypothesized that the lower density of GWAS-detected SNPs on the X-chromosome compared to autosomes is not a result of a methodological bias, e.g. differences in coverage or call rates, but has a real underlying biological reason-a lower density of functional SNPs on the X-chromosome versus autosomes. This hypothesis is supported by the observation that (i) the overall SNP density of X-chromosome is lower compared to the SNP density on autosomes and that (ii) the density of genic SNPs on the X-chromosome is lower compared to autosomes while densities of intergenic SNPs are similar.

Development and characterization of 148 SNP markers in the caribbean symmetrical brain coral Pseudodiploria strigosa

As coral reefs decline, restoring lost coral populations has been hampered due to the paucity of genetic information available for many coral species. In the Caribbean, the symmetrical brain coral Pseudodiploria strigosa, a prominent framework builder of the reef, has experienced an accelerated decline due to stony coral tissue loss disease (SCTLD). Colonies of P. strigosa gene-banked in response to the threat of SCTLD in Florida, USA, were sampled to develop 148 novel single nucleotide polymorphism (SNP) markers using genotyping-by-sequencing (GBS). The observed heterozygosity (Ho) and expected heterozygosity (He) ranged from 0.008 to 0.888 and 0.018 to 0.500, respectively. Deviations from Hardy–Weinberg equilibrium within populations, measured by the inbreeding coefficient index (Fis), ranged from − 0.799 to 0.923. In total, 96 SNPs were found to deviate significantly from Hardy–Weinberg (p < 0.05). These SNPs can be used for genetic population analysis to assist management and restoration of P. strigosa.

Next Generation Sequencing of Genotype Variants and Genetic Association between Heat Shock Proteins HSPA1B Single Nucleotide Polymorphism at the g.31829044 Locus and Heat Tolerance: A Pilot Quasi-Experimental Study

Heat tolerance and exertional heat stroke (EHS) are rare health conditions that have been described and characterised but have never been genetically solved. Knowledge of the role of single nucleotide polymorphisms (SNPs) in heat shock proteins (HSPs) genes and their associations with heat tolerance and EHS is limited. This pilot study aimed to identify SNP in HSPA1B, HSP90AA2 and DNAJA1 genes and their associations with heat tolerance and EHS history in a quasi-experimental design. Participants comprised Australian Defence Force members (ADF) who had a history of EHS and the general population. Genomic DNA samples were extracted from the venous blood samples of 48 participants, sequenced and analysed for SNP. Forty-four per cent (44%) of the participants were heat intolerant, and 29% had a history of EHS. Among participants with a history of EHS, there was an association between heat tolerance and HSPA1B SNP at the g.31829044 locus. However, there were no associations between HSPA1B and HSP90AA2 SNP and heat tolerance. All participants had the same distribution for the DNAJA1 SNP. In conclusion, the findings indicate an association between the HSPA1B genetic variant at the g.31829044 locus and heat tolerance among ADF participants with a history of EHS. Further research with a larger number of military participants will shed more light on the associations between HSP genes and heat tolerance.

Detection and Co-occurrence of kdr (F1534C and S989P) Mutations in Multiple Insecticides Resistant Aedes aegypti (Diptera: Culicidae) in Nigeria

The outbreak of yellow fever virus transmitted by Aedes aegypti has been of major concern in Nigeria, this mosquito also transmits several other arboviruses globally. The control of many vectors of mosquito-borne diseases relies heavily on the use of insecticides. Therefore, constant monitoring of insecticide resistance status and associated mechanisms is crucial within the local mosquito population. Here, we determined the resistance profile of adult Ae. aegypti from Majidun and Oke Ota communities, Ikorodu Local Government Area of Lagos State, Nigeria to different classes of insecticides using WHO procedures. The resistant phenotypes of Ae. aegypti were screened for the presence of kdr mutations F1534C, S989P, and V1016G, which have been implicated in insecticide resistance in yellow fever vectors. A high level of resistance to DDT and pyrethroids was recorded in Ae. aegypti in this study, although possible resistance to deltamethrin, one of the pyrethroids was reported in one of the locations. Resistance to bendiocarb was recorded in the Majidun community while Ae. aegypti in both locations were susceptible to malathion. The presence of F1534C mutation associated with DDT and deltamethrin resistance in Ae. aegypti population, and the presence of S989P mutation were detected singly and in co-occurrence with F1534C for the first time in Africa, while V1016G mutation was not detected in this study. The role of these mutations in resistance phenotype expressed in Ae. aegypti in larger populations needs to be established.

Landscape connectivity and genetic structure in a mainstem and a tributary stonefly (Plecoptera) species using a novel reference genome

Understanding how environmental variation influences population genetic structure can help predict how environmental change influences population connectivity, genetic diversity, and evolutionary potential. We used riverscape genomics modelling to investigate how climatic and habitat variables relate to patterns of genetic variation in two stonefly species, one from mainstem river habitats (Sweltsa coloradensis) and one from tributaries (Sweltsa fidelis) in 40 sites in northwest Montana, USA. We produced a draft genome assembly for S. coloradensis (N50 = 0.251 Mbp, BUSCO &gt; 95% using “insecta_ob9” reference genes). We genotyped 1930 SNPs in 372 individuals for S. coloradensis and 520 SNPs in 153 individuals for S. fidelis. We found higher genetic diversity for S. coloradensis compared to S. fidelis, but nearly identical genetic differentiation among sites within each species (both had global loci median FST = 0.000), despite differences in stream network location. For landscape genomics and testing for selection, we produced a less stringently filtered data set (3454 and 1070 SNPs for S. coloradensis and S. fidelis, respectively). Environmental variables (mean summer precipitation, slope, aspect, mean June stream temperature, land cover type) were correlated with 19 putative adaptive loci for S. coloradensis. but there was only one putative adaptive locus for S. fidelis (correlated with aspect). Interestingly, we also detected potential hybridization between multiple Sweltsa species which has never been previously detected. Studies like ours, that test for adaptive variation in multiple related species are needed to help assess landscape connectivity and the vulnerability of populations and communities to environmental change.

Major sex differences in allele frequencies for X chromosomal variants in both the 1000 Genomes Project and gnomAD

An unexpectedly high proportion of SNPs on the X chromosome in the 1000 Genomes Project phase 3 data were identified with significant sex differences in minor allele frequencies (sdMAF). sdMAF persisted for many of these SNPs in the recently released high coverage whole genome sequence of the 1000 Genomes Project that was aligned to GRCh38, and it was consistent between the five super-populations. Among the 245,825 common (MAF>5%) biallelic X-chromosomal SNPs in the phase 3 data presumed to be of high quality, 2,039 have genome-wide significant sdMAF (p-value <5e-8). sdMAF varied by location: non-pseudo-autosomal region (NPR) = 0.83%, pseudo-autosomal regions (PAR1) = 0.29%, PAR2 = 13.1%, and X-transposed region (XTR)/PAR3 = 0.85% of SNPs had sdMAF, and they were clustered at the NPR-PAR boundaries, among others. sdMAF at the NPR-PAR boundaries are biologically expected due to sex-linkage, but have generally been ignored in association studies. For comparison, similar analyses found only 6, 1 and 0 SNPs with significant sdMAF on chromosomes 1, 7 and 22, respectively. Similar sdMAF results for the X chromosome were obtained from the high coverage whole genome sequence data from gnomAD V 3.1.2 for both the non-Finnish European and African/African American samples. Future X chromosome analyses need to take sdMAF into account.

Molecular Epidemiology of Na+-Taurocholate Cotransporting Polypeptide Deficiency in Guangdong Province, China: A Pilot Study by Screening for Four Prevalent Variants of the Causative Gene SLC10A1

Na⁺-taurocholate cotransporting polypeptide deficiency (NTCPD) is an autosomal recessive disorder arising from biallelic SLC10A1 mutations. As a newly-described inborn error of bile acid metabolism, the epidemiology of this condition remains largely unclear in Chinese population so far. In this study, a total of 2,828 peripheral blood samples were collected from 12 cities in Guangdong, a province with the largest population in China, and the four prevalent SLC10A1 variants c.800C > T (p.Ser267Phe), c.263T > C (p.Ile88Thr), c.595A > C (p.Ser199Arg) and c.665T > C (p.Leu222Ser) were screened for by using polymerase chain reaction (PCR)- restriction fragment length polymorphism (RFLP). As a result, 663 mutated SLC10A1 alleles were detected, and the mutated allele frequency was calculated to be 11.72% (663/5,656), with a carrier frequency 20.69% (1/5) and a theoretical morbidity rate 1.37% (1/73) of NTCPD in Guangdong province. The variant c.800C > T (p.Ser267Phe) exhibited highest allele frequency among the four prevalent variants (χ² = 1501.27, p < 0.0001) as well as higher allele frequency in the peripheral region than that within the Pearl River Delta (χ² = 4.834, p < 0.05). The results suggested that NTCPD might be a disorder rather common in Guangdong province. The findings depicted the molecular epidemiologic features of NTCPD, providing preliminary but significant laboratory evidences for the subsequent NTCPD diagnosis and management in Guangdong population.

How HLA diversity is apportioned: influence of selection and relevance to transplantation

In his 1972 paper ‘The apportionment of human diversity’, Lewontin showed that, when averaged over loci, genetic diversity is predominantly attributable to differences among individuals within populations. However, selection can alter the apportionment of diversity of specific genes or genomic regions. We examine genetic diversity at the human leucocyte antigen (HLA) loci, located within the major histocompatibility complex (MHC) region. HLA genes code for proteins that are critical to adaptive immunity and are well-documented targets of balancing selection. The single-nucleotide polymorphisms (SNPs) within HLA genes show strong signatures of balancing selection on large timescales and are broadly shared among populations, displaying low F_ST values. However, when we analyse haplotypes defined by these SNPs (which define ‘HLA alleles’), we find marked differences in frequencies between geographic regions. These differences are not reflected in the F_ST values because of the extreme polymorphism at HLA loci, illustrating challenges in interpreting F_ST. Differences in the frequency of HLA alleles among geographic regions are relevant to bone-marrow transplantation, which requires genetic identity at HLA loci between patient and donor. We discuss the case of Brazil's bone marrow registry, where a deficit of enrolled volunteers with African ancestry reduces the chance of finding donors for individuals with an MHC region of African ancestry.

This article is part of the theme issue ‘Celebrating 50 years since Lewontin's apportionment of human diversity’.

Association of Angio-LncRNAs MIAT rs1061540/MALAT1 rs3200401 Molecular Variants with Gensini Score in Coronary Artery Disease Patients Undergoing Angiography

Long non-coding RNAs (lncRNAs) have emerged as essential biomolecules with variable diagnostic and/or prognostic utility in several diseases, including coronary artery disease (CAD). We aimed for the first time to investigate the potential association of five angiogenesis-related lncRNAs (PUNISHER, SENCR, MIAT, MALAT1, and GATA6-AS) variants with CAD susceptibility and/or severity. TaqMan Real-Time genotyping for PUNISHER rs12318065A/C, SENCR rs12420823C/T, MIAT rs1061540C/T, MALAT1 rs3200401T/C, and GATA6-AS1 rs73390820A/G were run on the extracted genomic DNA from 100 unrelated patients with stable CAD undergoing diagnostic coronary angiography and from 100 controls. After adjusting covariates, the studied variants showed no association with disease susceptibility; however, MIAT*T/T genotype was associated with a more severe Gensini score. In contrast, MALAT1*T/C heterozygosity was associated with a lower score. The lipid profile, and to a lesser extent smoking status, male sex, weight, hypertension, and MALAT1 (T > C) (negative correlation), explained the variance between patients/control groups via a principal component analysis. Incorporating the principal components into a logistic regression model to predict CAD yielded a 0.92 AUC. In conclusion: MIAT rs1061540 and MALAT1 rs3200401 variants were associated with CAD severity and Gensini score in the present sample of the Egyptian population. Further large multi-center and functional analyses are needed to confirm the results and identify the underlying molecular mechanisms.

Is there an association of genetic polymorphisms of the catechol-O-methyltransferase gene (rs165656 and rs174675) and the 5-hydroxytryptamine receptor 2A gene (rs4941573 and rs6313) with sleep bruxism in individuals with obstructive sleep apnea?

OBJECTIVE

To evaluate the association of single-nucleotide polymorphisms within the catechol-O-methyltransferase and 5-hydroxytryptamine receptor 2A genes with sleep bruxism in individuals diagnosed with obstructive sleep apnea.

DESIGN

Sixty-nine individuals with suspected sleep-related problems were evaluated by polysomnography, following the recommendations of the American Academy of Sleep Medicine. Deoxyribonucleic acid (DNA) samples were collected only from 48 of the study participants because of missing polysomnographic data. DNA samples were collected and two single-nucleotide polymorphisms in the 5-hydroxytryptamine receptor 2A encoding HTR2A gene (rs4941573 and rs6313) and two in the catechol-O-methyltransferase gene (rs165656 and rs174675) were selected to be genotyped using real-time polymerase chain reaction. The association between sleep bruxism and genetic polymorphisms was investigated by recessive and dominant models. Association analyses were performed using a 95% confidence interval and the level of statistical significance was p < 0.05.

RESULTS

From the 69 study participants, 48 were included in the polymorphism analysis and sleep bruxism was present in 35.4%. No significant differences were observed in the dominant and recessive models (p > 0.05). Haplotype and diplotype analyses revealed the predicted four haplotypes and two diplotypes were not associated with sleep bruxism.

CONCLUSION

Polymorphisms rs174675 and rs165656 in the catechol-O-methyltransferase gene and rs4941573 and rs6313 in the 5-hydroxytryptamine receptor 2A gene were not significantly associated with sleep bruxism in individuals with obstructive sleep apnea.

Next Generation Sequencing of Single Nucleotide Polymorphic DNA-Markers in Selecting for Intramuscular Fat, Fat Melting Point, Omega-3 Long-Chain Polyunsaturated Fatty Acids and Meat Eating Quality in Tattykeel Australian White MARGRA Lamb

Meat quality data can only be obtained after slaughter when selection decisions about the live animal are already too late. Carcass estimated breeding values present major precision problems due to low accuracy, and by the time an informed decision on the genetic merit for meat quality is made, the animal is already dead. We report for the first time, a targeted next-generation sequencing (NGS) of single nucleotide polymorphisms (SNP) of lipid metabolism genes in Tattykeel Australian White (TAW) sheep of the MARGRA lamb brand, utilizing an innovative and minimally invasive muscle biopsy sampling technique for directly quantifying the genetic worth of live lambs for health-beneficial omega-3 long-chain polyunsaturated fatty acids (n-3 LC-PUFA), intramuscular fat (IMF), and fat melting point (FMP). NGS of stearoyl-CoA desaturase (SCD), fatty acid binding protein-4 (FABP4), and fatty acid synthase (FASN) genes identified functional SNP with unique DNA marker signatures for TAW genetics. The SCD g.23881050T>C locus was significantly associated with IMF, C22:6n-3, and C22:5n-3; FASN g.12323864A>G locus with FMP, C18:3n-3, C18:1n-9, C18:0, C16:0, MUFA, and FABP4 g.62829478A>T locus with IMF. These add new knowledge, precision, and reliability in directly making early and informed decisions on live sheep selection and breeding for health-beneficial n-3 LC-PUFA, FMP, IMF and superior meat-eating quality at the farmgate level. The findings provide evidence that significant associations exist between SNP of lipid metabolism genes and n-3 LC-PUFA, IMF, and FMP, thus underpinning potential marker-assisted selection for meat-eating quality traits in TAW lambs.

Local adaptation and archaic introgression shape global diversity at human structural variant loci

Large genomic insertions and deletions are a potent source of functional variation, but are challenging to resolve with short-read sequencing, limiting knowledge of the role of such structural variants (SVs) in human evolution. Here, we used a graph-based method to genotype long-read-discovered SVs in short-read data from diverse human genomes. We then applied an admixture-aware method to identify 220 SVs exhibiting extreme patterns of frequency differentiation - a signature of local adaptation. The top two variants traced to the immunoglobulin heavy chain locus, tagging a haplotype that swept to near fixation in certain southeast Asian populations, but is rare in other global populations. Further investigation revealed evidence that the haplotype traces to gene flow from Neanderthals, corroborating the role of immune-related genes as prominent targets of adaptive introgression. Our study demonstrates how recent technical advances can help resolve signatures of key evolutionary events that remained obscured within technically challenging regions of the genome.

Revealing modifier variations characterizations for elucidating the genetic basis of human phenotypic variations

Epistatic interactions complicate the identification of variants involved in phenotypic effect. In-depth knowledge in modifiers and in pathogenic variants would benefit the mechanistic studies on the genetic basis of complex traits. We systematically compared the modifier variants which have evidence of modifier effect with the pathogenic variants from multiple angles. Our study found that genomic loci of modifier variations differ from pathogenic loci in many aspects, such as population genetics statistics, epigenetic features, evolutionary characteristics and functional properties of the variations. Genes containing modifier variation(s) exhibit higher probability of being haploinsufficient and higher probability of recessive disease causation, and they are relatively more important in network communication. Furthermore, we reinforced that co-expression analysis is an effective methodology to predict functional associations between modifier genes and their potential target genes. In many aspects, we detected statistically significant differences between modifier variants/genes and pathogenic variants/genes, and investigated relationships between modifiers and their potential targets. Our results offer some actionable insights that may provide appropriate guidelines to clinical genetics and researchers to elucidate the molecular mechanism underlying the human phenotypic variation.

Deciphering the association of intronic single nucleotide polymorphisms of crystallin gene family with congenital cataract

Purpose

Introns play an important role in gene regulation and expression. Single nucleotide polymorphisms (SNPs) in introns have the potential to cause disease and alter the genotype-phenotype association. Hence, this study aimed to decipher the association of SNPs in the introns of the crystallin gene in congenital cataracts.

Methods

SNPs in the introns of crystallin gene family - CRYAA (rs3788059), CRYAB (rs2070894), CRYBA4 (rs2071861), and CRYBB2 (rs5752083, rs5996863) - were genotyped in 248 participants consisting of 141 congenital cataracts and 107 healthy controls by allele-specific oligonucleotide polymerase chain reaction method. Around 10% of samples for each SNPs were sequenced to confirm the genotypes. The allele, genotype, and haplotype frequency were evaluated by the SHEsis online tool.

Results

Using dominant model, the "A" allele of rs3788059 was found to have an increased risk toward congenital cataract development whereas the "G" allele was found to be protective (AA + AG vs. GG; odds ratio [OR] 95% confidence interval [CI] = 3.73 [1.71, 8.15], P = 0.0009). The "A" allele of both rs2070894 (AA + AG vs. GG; OR [95% CI] = 0.49 [0.29, 0.84], P = 0.012) and rs5752083 (AA + AC vs. CC; OR [95% CI] = 0.25 [0.08, 0.76], P = 0.016) were suggested to have a protective role by the dominant model. The A-C-T haplotype (rs2071861, rs5752083, and rs5996863) was found to be a significant risk factor for the development of congenital cataract.

Conclusion

Intronic SNPs in crystallin genes may play a role in the predisposition toward congenital cataract. However, the present findings need to be replicated in a large cohort with more number of samples.

Molecular epidemiologic study of citrin deficiency by screening for four reported pathogenic SLC25A13 variants in the Shaanxi and Guangdong provinces, China

BACKGROUND

Citrin deficiency (CD) is an autosomal recessive disease resulting from biallelic mutations of the SLC25A13 gene. This study aimed to investigate the molecular epidemiological features of CD in the Guangdong and Shaanxi provinces of China.

METHODS

A total of 3,409 peripheral blood samples from Guangdong and 2,746 such samples from Shaanxi province were collected. Four prevalent SLC25A13 mutations NG_012247.2 (NM_014251.3): c.852_855del, c.1638_1660dup, c.615+5G>A, and c.1751-5_1751-4ins(2684) were screened by using the conventional polymerase chain reaction (PCR)/PCR-restriction fragment length polymorphism and newly-developed multiplex PCR methods, respectively. The mutated SLC25A13 allele frequencies, carrier frequencies, and CD morbidity rates were calculated and then compared with the Chi-square and Fisher's exact tests.

RESULTS

The mutations were detected in 68 out of 6,818 SLC25A13 alleles in Guangdong and 29 out of 5,492 alleles in the Shaanxi population. The carrier frequencies were subsequently calculated to be 1/51 and 1/95, while the CD morbidity rates were 1/10,053 and 1/35,865, in the 2 populations, respectively. When compared with the Shaanxi population, Guangdong exhibited a higher frequency of mutated SLC25A13 allele (68/6,818 vs. 29/5,492, χ²=8.570, P=0.003) in general, with higher c.852_855del (54/6,818 vs. 13/5,492, χ²=17.328, P=0.000) but lower c.1751-5_1751 -4ins(2684) (2/6,818 vs. 9/5,492, P=0.015) allele frequencies. The distribution of c.615+5G>A and c.1638_1660dup between the 2 provinces, as well as all 4 prevalent mutations among different geographic regions within the 2 provinces, did not differed significantly.

CONCLUSIONS

Our findings depicted the CD molecular epidemiological features in Guangdong and Shaanxi populations, providing preliminary but significant laboratory evidences for the subsequent CD diagnosis and management in the 2 provinces of mainland China.

A novel (ATC)n microsatellite locus is associated with litter size in an indigenous Chinese pig

Simple sequence repeats (SSRs) are an important part of the genome and have become powerful auxiliary DNA markers in animal breeding using marker-assisted selection (MAS). Based on previous sequencing data of Qinghai Bamei pigs, a total of three novel candidate SSR loci were analysed in this study. Time-of-flight mass spectrometry (TOF-MS) was used for SSR genotyping, and association analyses between SSRs and the litter size of Qinghai Bamei sows was also performed. The results of genotyping showed that the (ATC)_n -P1, (AC)_n -P2 and (AC)_n -P3 loci had 2, 3 and 18 genotypes, respectively; 2, 3 and 8 alleles were also identified at these loci. Except for the (AC)_n -P2 locus, the polymorphism information content (PIC) values of other loci were greater than 0.25. Association analyses indicated that only the (ATC)_n -P1 locus was significantly associated with the litter size of Qinghai Bamei sows (p = .047). Compared to 189-/189- genotype, individuals with the 189-/195- genotype had the senior litter size, which was 9.04 ± 0.21. Our results enrich the data on SSRs in Qinghai Bamei pigs and indicate that (ATC)_n -P1 is a candidate locus for MAS in the pig industry.

Expression quantitative trait loci of genes predicting outcome are associated with survival of multiple myeloma patients

Gene expression profiling can be used for predicting survival in multiple myeloma (MM) and identifying patients who will benefit from particular types of therapy. Some germline single nucleotide polymorphisms (SNPs) act as expression quantitative trait loci (eQTLs) showing strong associations with gene expression levels. We performed an association study to test whether eQTLs of genes reported to be associated with prognosis of MM patients are directly associated with measures of adverse outcome. Using the genotype-tissue expression portal, we identified a total of 16 candidate genes with at least one eQTL SNP associated with their expression with P < 10^-7 either in EBV-transformed B-lymphocytes or whole blood. We genotyped the resulting 22 SNPs in 1327 MM cases from the International Multiple Myeloma rESEarch (IMMEnSE) consortium and examined their association with overall survival (OS) and progression-free survival (PFS), adjusting for age, sex, country of origin and disease stage. Three polymorphisms in two genes (TBRG4-rs1992292, TBRG4-rs2287535 and ENTPD1-rs2153913) showed associations with OS at P < .05, with the former two also associated with PFS. The associations of two polymorphisms in TBRG4 with OS were replicated in 1277 MM cases from the International Lymphoma Epidemiology (InterLymph) Consortium. A meta-analysis of the data from IMMEnSE and InterLymph (2579 cases) showed that TBRG4-rs1992292 is associated with OS (hazard ratio = 1.14, 95% confidence interval 1.04-1.26, P = .007). In conclusion, we found biologically a plausible association between a SNP in TBRG4 and OS of MM patients.

Association of NFκβ, TNFα, IL-6, IL-1β, and LPL Polymorphisms with Type 2 Diabetes Mellitus and Biochemical Parameters in a Mexican Population

Chronic low-grade inflammation is strongly related to the etiology of diabetes mellitus type 2 (T2DM), and the expression of inflammatory cytokines may be modulated by polymorphisms located in the regulatory regions of the NFκβ, IL-1β, IL-6, TNFα, and LPL genes. We considered it particularly important to investigate the relationship of gene polymorphisms involved in chronic inflammation with the risk of T2DM or uncontrolled biochemical parameters.

METHODS

We included 199 individuals with a T2DM diagnosis and 213 individuals without a T2DM diagnosis. Restriction fragment length polymorphism (RFLP) analyses were used to assess polymorphisms.

RESULTS

We found a risk association between T2DM and uncontrolled biochemical parameters in a Mexican population for the genotypes del/del of NFκβ, -174 and -572 of IL-6, C/C of IL-1β, -308 and -238 of TNFα, and T/T of LPL. In subjects without diabetes (controls), we found an association between the G/C genotype of the -572 polymorphism and the G/C and C/C genotypes of the -597 polymorphism of IL-6 with the risk of glucose levels > 131 mg/dL. Genotype C/C of polymorphism -174 of the IL-6 gene was associated with high triglyceride levels, and levels > 5.8% of HbA1c were associated with the G/A genotype of TNFα -308.

CONCLUSION

Here, we describe for the first time the relationship of T2DM risk and uncontrolled biochemical parameters with polymorphisms in the NFκβ, IL-6, TNFα, IL-1β, and LPL genes in a Mexican population. We also showed that for the population included in this study, there is an additive effect of the polymorphisms of the studied genes that considerably increases the risk of developing T2DM.We also showed that there are interactions between genes related to chronic inflammation that affect the risk of T2DM.

No support for the hereditarian hypothesis of the Black-White achievement gap using polygenic scores and tests for divergent selection

OBJECTIVES

Debate about the cause of IQ score gaps between Black and White populations has persisted within genetics, anthropology, and psychology. Recently, authors claimed polygenic scores provide evidence that a significant portion of differences in cognitive performance between Black and White populations are caused by genetic differences due to natural selection, the "hereditarian hypothesis." This study aims to show conceptual and methodological flaws of past studies supporting the hereditarian hypothesis.

MATERIALS AND METHODS

Polygenic scores for educational attainment were constructed for African and European samples of the 1000 Genomes Project. Evidence for selection was evaluated using an excess variance test. Education associated variants were further evaluated for signals of selection by testing for excess genetic differentiation (F_st ). Expected mean difference in IQ for populations was calculated under a neutral evolutionary scenario and contrasted to hereditarian claims.

RESULTS

Tests for selection using polygenic scores failed to find evidence of natural selection when the less biased within-family GWAS effect sizes were used. Tests for selection using F_st values did not find evidence of natural selection. Expected mean difference in IQ was substantially smaller than postulated by hereditarians, even under unrealistic assumptions that overestimate genetic contribution.

CONCLUSION

Given these results, hereditarian claims are not supported in the least. Cognitive performance does not appear to have been under diversifying selection in Europeans and Africans. In the absence of diversifying selection, the best case estimate for genetic contributions to group differences in cognitive performance is substantially smaller than hereditarians claim and is consistent with genetic differences contributing little to the Black-White gap.

Accurate, scalable cohort variant calls using DeepVariant and GLnexus

Motivation

Population-scale sequenced cohorts are foundational resources for genetic analyses, but processing raw reads into analysis-ready cohort-level variants remains challenging.

Results

We introduce an open-source cohort-calling method that uses the highly accurate caller DeepVariant and scalable merging tool GLnexus. Using callset quality metrics based on variant recall and precision in benchmark samples and Mendelian consistency in father-mother-child trios, we optimize the method across a range of cohort sizes, sequencing methods and sequencing depths. The resulting callsets show consistent quality improvements over those generated using existing best practices with reduced cost. We further evaluate our pipeline in the deeply sequenced 1000 Genomes Project (1KGP) samples and show superior callset quality metrics and imputation reference panel performance compared to an independently generated GATK Best Practices pipeline.

Availability and implementation

We publicly release the 1KGP individual-level variant calls and cohort callset (https://console.cloud.google.com/storage/browser/brain-genomics-public/research/cohort/1KGP) to foster additional development and evaluation of cohort merging methods as well as broad studies of genetic variation. Both DeepVariant (https://github.com/google/deepvariant) and GLnexus (https://github.com/dnanexus-rnd/GLnexus) are open-source, and the optimized GLnexus setup discovered in this study is also integrated into GLnexus public releases v1.2.2 and later.

Supplementary information

Supplementary data are available at Bioinformatics online.

Culture-free genome-wide locus sequence typing (GLST) provides new perspectives on Trypanosoma cruzi dispersal and infection complexity

Analysis of genetic polymorphism is a powerful tool for epidemiological surveillance and research. Powerful inference from pathogen genetic variation, however, is often restrained by limited access to representative target DNA, especially in the study of obligate parasitic species for which ex vivo culture is resource-intensive or bias-prone. Modern sequence capture methods enable pathogen genetic variation to be analyzed directly from host/vector material but are often too complex and expensive for resource-poor settings where infectious diseases prevail. This study proposes a simple, cost-effective 'genome-wide locus sequence typing' (GLST) tool based on massive parallel amplification of information hotspots throughout the target pathogen genome. The multiplexed polymerase chain reaction amplifies hundreds of different, user-defined genetic targets in a single reaction tube, and subsequent agarose gel-based clean-up and barcoding completes library preparation at under 4 USD per sample. Our study generates a flexible GLST primer panel design workflow for Trypanosoma cruzi, the parasitic agent of Chagas disease. We successfully apply our 203-target GLST panel to direct, culture-free metagenomic extracts from triatomine vectors containing a minimum of 3.69 pg/μl T. cruzi DNA and further elaborate on method performance by sequencing GLST libraries from T. cruzi reference clones representing discrete typing units (DTUs) TcI, TcIII, TcIV, TcV and TcVI. The 780 SNP sites we identify in the sample set repeatably distinguish parasites infecting sympatric vectors and detect correlations between genetic and geographic distances at regional (< 150 km) as well as continental scales. The markers also clearly separate TcI, TcIII, TcIV and TcV + TcVI and appear to distinguish multiclonal infections within TcI. We discuss the advantages, limitations and prospects of our method across a spectrum of epidemiological research.

Genome-Wide Association Study of Metamizole-Induced Agranulocytosis in European Populations

Agranulocytosis is a rare yet severe idiosyncratic adverse drug reaction to metamizole, an analgesic widely used in countries such as Switzerland and Germany. Notably, an underlying mechanism has not yet been fully elucidated and no predictive factors are known to identify at-risk patients. With the aim to identify genetic susceptibility variants to metamizole-induced agranulocytosis (MIA) and neutropenia (MIN), we conducted a retrospective multi-center collaboration including cases and controls from three European populations. Association analyses were performed using genome-wide genotyping data from a Swiss cohort (45 cases, 191 controls) followed by replication in two independent European cohorts (41 cases, 273 controls) and a joint discovery meta-analysis. No genome-wide significant associations (p < 1 × 10⁻⁷) were observed in the Swiss cohort or in the joint meta-analysis, and no candidate genes suggesting an immune-mediated mechanism were identified. In the joint meta-analysis of MIA cases across all cohorts, two candidate loci on chromosome 9 were identified, rs55898176 (OR = 4.01, 95%CI: 2.41–6.68, p = 1.01 × 10⁻⁷) and rs4427239 (OR = 5.47, 95%CI: 2.81–10.65, p = 5.75 × 10⁻⁷), of which the latter is located in the SVEP1 gene previously implicated in hematopoiesis. This first genome-wide association study for MIA identified suggestive associations with biological plausibility that may be used as a stepping-stone for post-GWAS analyses to gain further insight into the mechanism underlying MIA.

Genetic Diversity and Genetic Structure of the Wild Tsushima Leopard Cat from Genome-Wide Analysis

Simple Summary

The Tsushima leopard cat, Prionailurus bengalensis euptilurus, is a small regional population of the Amur leopard cat and is only found on Tsushima Island in Japan. A breeding program will require adequate information on parentage, kinship, and inbreeding for this population. Hence, there is an urgent need to develop this information in order to conserve the population and its genetic diversity. We performed GRAS-Di analysis to investigate the genetic diversity and genetic structure of the Tsushima leopard cat. We identified between 133 and 158 single-nucleotide polymorphism (SNP) markers in three different genotyping methods. These SNP markers can be used in identification of individuals and parentage. In addition, structure analysis using these markers demonstrated the similar genetic composition of the samples from 48 Tsushima leopard cats, and indicated Tsushima leopard cats have no subpopulations. We have provided genetic markers that are useful for conservation of the Tsushima leopard cat, such as individual identification and parentage. Moreover, we have also clarified units for conservation of the Tsushima leopard cat population from structure analysis.

Abstract

The Tsushima leopard cat (Prionailurus bengalensis euptilurus) lives on Tsushima Island in Japan and is a regional population of the Amur leopard cat; it is threatened with extinction. Its genetic management is important because of the small population. We used genotyping by random amplicon sequencing-direct (GRAS-Di) to develop a draft genome and explore single-nucleotide polymorphism (SNP) markers. The SNPs were analyzed using three genotyping methods (mapping de novo, to the Tsushima leopard cat draft genome, and to the domestic cat genome). We examined the genetic diversity and genetic structure of the Tsushima leopard cat. The genome size was approximately 2.435 Gb. The number of SNPs identified was 133–158. The power of these markers was sufficient for individual and parentage identifications. These SNPs can provide useful information about the life of the Tsushima leopard cat and the pairings and for the introduction of founders to conserve genetic diversity with ex situ conservation. We identified that there are no subpopulations of the Tsushima leopard cat. The identifying units will allow for a concentration of efforts for conservation. SNPs can be applied to the analysis of the leopard cat in other regions, making them useful for comparisons among populations and conservation in other small populations.

Hardy-Weinberg Equilibrium in the Large Scale Genomic Sequencing Era

Hardy-Weinberg Equilibrium (HWE) is used to estimate the number of homozygous and heterozygous variant carriers based on its allele frequency in populations that are not evolving. Deviations from HWE in large population databases have been used to detect genotyping errors, which can result in extreme heterozygote excess (HetExc). However, HetExc might also be a sign of natural selection since recessive disease causing variants should occur less frequently in a homozygous state in the population, but may reach high allele frequency in a heterozygous state, especially if they are advantageous. We developed a filtering strategy to detect these variants and applied it on genome data from 137,842 individuals. The main limitations of this approach were quality of genotype calls and insufficient population sizes, whereas population structure and inbreeding can reduce sensitivity, but not precision, in certain populations. Nevertheless, we identified 161 HetExc variants in 149 genes, most of which were specific to African/African American populations (∼79.5%). Although the majority of them were not associated with known diseases, or were classified as clinically "benign," they were enriched in genes associated with autosomal recessive diseases. The resulting dataset also contained two known recessive disease causing variants with evidence of heterozygote advantage in the sickle-cell anemia (HBB) and cystic fibrosis (CFTR). Finally, we provide supporting in silico evidence of a novel heterozygote advantageous variant in the chromodomain helicase DNA binding protein 6 gene (CHD6; involved in influenza virus replication). We anticipate that our approach will aid the detection of rare recessive disease causing variants in the future.

Systematic genetic analysis of the MHC region reveals mechanistic underpinnings of HLA type associations with disease

The MHC region is highly associated with autoimmune and infectious diseases. Here we conduct an in-depth interrogation of associations between genetic variation, gene expression and disease. We create a comprehensive map of regulatory variation in the MHC region using WGS from 419 individuals to call eight-digit HLA types and RNA-seq data from matched iPSCs. Building on this regulatory map, we explored GWAS signals for 4083 traits, detecting colocalization for 180 disease loci with eQTLs. We show that eQTL analyses taking HLA type haplotypes into account have substantially greater power compared with only using single variants. We examined the association between the 8.1 ancestral haplotype and delayed colonization in Cystic Fibrosis, postulating that downregulation of RNF5 expression is the likely causal mechanism. Our study provides insights into the genetic architecture of the MHC region and pinpoints disease associations that are due to differential expression of HLA genes and non-HLA genes.

Genetic Imbalance Is Associated With Functional Outcome After Ischemic Stroke

Background and Purpose- We sought to explore the effect of genetic imbalance on functional outcome after ischemic stroke (IS). Methods- Copy number variation was identified in high-density single-nucleotide polymorphism microarray data of IS patients from the CADISP (Cervical Artery Dissection and Ischemic Stroke Patients) and SiGN (Stroke Genetics Network)/GISCOME (Genetics of Ischaemic Stroke Functional Outcome) networks. Genetic imbalance, defined as total number of protein-coding genes affected by copy number variations in an individual, was compared between patients with favorable (modified Rankin Scale score of 0-2) and unfavorable (modified Rankin Scale score of ≥3) outcome after 3 months. Subgroup analyses were confined to patients with imbalance affecting ohnologs-a class of dose-sensitive genes, or to those with imbalance not affecting ohnologs. The association of imbalance with outcome was analyzed by logistic regression analysis, adjusted for age, sex, stroke subtype, stroke severity, and ancestry. Results- The study sample comprised 816 CADISP patients (age 44.2±10.3 years) and 2498 SiGN/GISCOME patients (age 67.7±14.2 years). Outcome was unfavorable in 122 CADISP and 889 SiGN/GISCOME patients. Multivariate logistic regression analysis revealed that increased genetic imbalance was associated with less favorable outcome in both samples (CADISP: P=0.0007; odds ratio=0.89; 95% CI, 0.82-0.95 and SiGN/GISCOME: P=0.0036; odds ratio=0.94; 95% CI, 0.91-0.98). The association was independent of age, sex, stroke severity on admission, stroke subtype, and ancestry. On subgroup analysis, imbalance affecting ohnologs was associated with outcome (CADISP: odds ratio=0.88; 95% CI, 0.80-0.95 and SiGN/GISCOME: odds ratio=0.93; 95% CI, 0.89-0.98) whereas imbalance without ohnologs lacked such an association. Conclusions- Increased genetic imbalance was associated with poorer functional outcome after IS in both study populations. Subgroup analysis revealed that this association was driven by presence of ohnologs in the respective copy number variations, suggesting a causal role of the deleterious effects of genetic imbalance.

Programmed cell death 1 gene polymorphism as a possible risk for systemic lupus erythematosus in Egyptian females

Systemic lupus erythematosus (SLE) is a multisystem autoimmune disease with a suggested genetic basis. The newly identified human programmed cell death 1 gene could be associated with SLE susceptibility. We aimed to investigate the association between programmed cell death 1 polymorphism (PD1.3G/A (rs11568821) and PD1.5C/T (rs2227981)) with the risk of SLE in the Egyptian female population. This retrospective case-control study included 150 Egyptian females; 70 patients diagnosed to have SLE and 80 age-matched healthy controls. The two single nucleotide polymorphisms of the pdcd1 gene were genotyped by allelic discrimination through TaqMan real-time polymerase chain reaction. The PD1.3GG genotype and G allele as well as the PD1.5CC genotype were significantly more frequent in SLE patients (67.1%; p = 0.023, 82.1%; p = 0.0021, 62.9%; p = 0.0287 respectively). The GC haplotype was the most common haplotype among SLE patients (70.77%) with a reported significant linkage disequilibrium between the two studied polymorphisms (p = 0.0041). Although most of the studies showed significant association of SLE with the minor alleles, we reported a significant association between the dominant genotypes (PD1.3GG and PD1.5CC) as well as the major G allele with the risk of SLE among Egyptian females.

Identification of HLA-DRB1*04:10 allele as risk allele for Japanese moyamoya disease and its association with autoimmune thyroid disease: A case-control study

Background and purpose

Moyamoya disease (MMD) is a progressive cerebrovascular disease with unknown etiology. Growing evidence suggest its involvement of autoimmune and genetic mechanisms in the pathogenesis of MMD. This study aims to clarify the association between HLA allele and MMD.

Methods

Case-control study: the DNA of 136 MMD patients in Japan was extracted and the genotype of human leukocyte antigen (HLA) from this DNA was determined by super-high-resolution single-molecule sequence-based typing using next-generation sequencing. Next, the frequency of each HLA allele (HLA-A, HLA-B, HLA-C, HLA-DRB1, HLA-DQB1, and HLA-DPB1) was compared with those in the Japanese control database. In addition, haplotype estimation was performed using the expectation maximization algorithm.

Results

The frequencies of the HLA-DRB1*04:10 allele (4.77% vs. 1.47% in the control group; P = 1.7 × 10⁻³; odds ratio [OR] = 3.35) and of the HLA-DRB1*04:10–HLA-DQB1*04:02 haplotype (haplotype frequency 4.41% vs. 1.35% in the control group; P = 2.0 × 10⁻³; OR = 3.37) significantly increased. The frequency of thyroid diseases, such as Graves’ disease and Hashimoto thyroiditis, increased in HLA-DRB1*04:10-positive MMD patients compared with that in HLA-DRB1*04:10-negative MMD patients.

Conclusions

HLA-DRB1*04:10 is a risk allele and HLA-DRB1*04:10–HLA-DQB1*04:02 a risk haplotype for MMD. In addition, HLA-DRB1*04:10 is associated with thyroid disease in MMD patients.

Bayesian model selection for the study of Hardy-Weinberg proportions and homogeneity of gender allele frequencies

Standard statistical tests for Hardy-Weinberg equilibrium assume the equality of allele frequencies in the sexes, whereas tests for the equality of allele frequencies in the sexes assume Hardy-Weinberg equilibrium. This produces a circularity in the testing of genetic variants, which has recently been resolved with new frequentist likelihood and exact procedures. In this paper, we tackle the same problem by posing it as a Bayesian model comparison problem. We formulate an exhaustive set of ten alternative scenarios for biallelic genetic variants. Using Dirichlet and Beta priors for genotype and allele frequencies, we derive marginal likelihoods for all scenarios, and select the most likely scenario using the posterior probabilities that each of these scenarios is the one in place. Different from the usual frequentist testing approach, the Bayesian approach allows one to compare any number of models, and not just two at a time, and the models compared do not have to be nested. We illustrate our Bayesian approach with genetic data from the 1,000 genomes project and through a simulation study.

A test for deviations from expected genotype frequencies on the X chromosome for sex-biased admixed populations

Genome-wide scans for deviations from expected genotype frequencies, as determined by the Hardy-Weinberg equilibrium (HWE), are commonly applied to detect genotyping errors and deviations from random mating. In contrast to the autosomes, genotype frequencies on the X chromosome do not reach HWE within a single generation. Instead, if allele frequencies in males and females initially differ, they oscillate for a few generations toward equilibrium. Allele frequency differences between the sexes are expected in populations that have experienced recent sex-biased admixture, namely, their male and female founders differed in ancestry. Sex-biased admixture does not allow testing for HWE on X, because deviations are naturally expected, even under random mating (post admixture) and error-free genotyping. In this paper, we develop a likelihood ratio test and a χ² test to detect deviations from expected genotype frequencies on X, beyond natural deviations due to sex-biased admixture. We demonstrate by simulations that our tests are powerful for detecting deviations due to non-random mating, while at the same time they do not reject the null under historical sex-biased admixture and random mating thereafter. We also demonstrate that when applied to 1000 Genomes project populations, our likelihood ratio test rejects fewer SNPs than other tests, but we describe limitations in the interpretation of the results.

Inflammasome Genes' Polymorphisms in Egyptian Chronic Hepatitis C Patients: Influence on Vulnerability to Infection and Response to Treatment

Chronic inflammation is a pivotal contributor to the liver damage mediated by hepatitis C virus (HCV). The NOD-like receptor, pyrin domain-containing 3 (NLRP3) inflammasome is activated by HCV in both hepatocytes and Kupffer cells. The aim of our study was to investigate the association of nine single-nucleotide polymorphisms in four inflammasome genes (NLRP3, CARD8, IL-1β, and IL-18) with the susceptibility to HCV infection and outcome of interferon treatment in 201 Egyptian chronic hepatitis C patients and 95 healthy controls. The genotyping was conducted using TaqMan predesigned SNP assay. In the comparative analysis, the CC genotype of the NLRP3 rs1539019 was found to be associated with the lower risk to chronic HCV infection (OR: 0.33, 95% CI: 0.17-0.62). This association was also found for the CA genotype and the A allele of the NLRP3 rs35829419 (OR: 0.18 and 0.22, respectively), in addition to the GG genotype and G allele of IL-18 rs1946518 (OR: 0.55 and 0.61, respectively). In contrast, the AA genotype of the IL-1β rs1143629 was significantly more frequent in HCV patients (OR: 1.7, 95% CI: 1-2.86). Notably, the frequency of the AA genotype of NLRP3 rs1539019 was significantly higher in patients with lack of response (NR) to the interferon treatment (OR: 1.95, 95% CI: 1-3.7). A similar association was found for both the CC genotype and C allele of the NLRP3 rs35829419 (OR: 2.78 and 2.73, respectively) and for the TT genotype and T allele of CARD8 rs2043211 (OR: 2.64 and 1.54, respectively). Yet, the IL-1β (rs1143629, rs1143634) and IL-18 (rs187238, rs1946518) polymorphisms did not show any significant association with response to interferon treatment. In conclusion, this study reports, for the first time, the association of genetic variations in NLRP3 with hepatitis C susceptibility and response to treatment in Egyptian patients. However, further large-scale studies are recommended to confirm our findings.

Development and Validation of 58K SNP-Array and High-Density Linkage Map in Nile Tilapia (O. niloticus)

Despite being the second most important aquaculture species in the world accounting for 7.4% of global production in 2015, tilapia aquaculture has lacked genomic tools like SNP-arrays and high-density linkage maps to improve selection accuracy and accelerate genetic progress. In this paper, we describe the development of a genotyping array containing more than 58,000 SNPs for Nile tilapia (Oreochromis niloticus). SNPs were identified from whole genome resequencing of 32 individuals from the commercial population of the Genomar strain, and were selected for the SNP-array based on polymorphic information content and physical distribution across the genome using the Orenil1.1 genome assembly as reference sequence. SNP-performance was evaluated by genotyping 4991 individuals, including 689 offspring belonging to 41 full-sib families, which revealed high-quality genotype data for 43,588 SNPs. A preliminary genetic linkage map was constructed using Lepmap2 which in turn was integrated with information from the O_niloticus_UMD1 genome assembly to produce an integrated physical and genetic linkage map comprising 40,186 SNPs distributed across 22 linkage groups (LGs). Around one-third of the LGs showed a different recombination rate between sexes, with the female being greater than the male map by a factor of 1.2 (1632.9 to 1359.6 cM, respectively), with most LGs displaying a sigmoid recombination profile. Finally, the sex-determining locus was mapped to position 40.53 cM on LG23, in the vicinity of the anti-Müllerian hormone (amh) gene. These new resources has the potential to greatly influence and improve the genetic gain when applying genomic selection and surpass the difficulties of efficient selection for invasively measured traits in Nile tilapia.

Big data in forensic genetics

The potential and difficulties of the application of genome wide data in forensics are analyzed. We argue that, besides statistical, computational, ethical, economic and technical validation problems, the state of the art of population genetics theory is insufficient to deal with the forensic use of this type of data. In order to keep the current standards of quantifying and reporting genetic evidence, namely in kinship analyses and identification, substantial improvement in the theoretical framework should be reached, since to obtain genome-wide results is to provide the experts with data that they cannot quantify the corresponding evidentiary value. Therefore, while a satisfactory, generalized theoretical and biostatistical modelling is not achieved, it may well be wiser to improve the already established approaches to a limited, pre-defined number of validated genetic markers, amenable to a consensual handling and reporting. Whole genome population analyses will prove extremely useful in selecting the best suited and most efficient of those markers.

What has GWAS done for HLA and disease associations?

The major histocompatibility complex (MHC) is located in chromosome 6p21 and contains crucial regulators of immune response, including human leucocyte antigen (HLA) genes, alongside other genes with nonimmunological roles. More recently, a repertoire of noncoding RNA genes, including expressed pseudogenes, has also been identified. The MHC is the most gene dense and most polymorphic part of the human genome. The region exhibits haplotype-specific linkage disequilibrium patterns, contains the strongest cis- and trans-eQTLs/meQTLs in the genome and is known as a hot spot for disease associations. Another layer of complexity is provided to the region by the extreme structural variation and copy number variations. While the HLA-B gene has the highest number of alleles, the HLA-DR/DQ subregion is structurally most variable and shows the highest number of disease associations. Reliance on a single reference sequence has complicated the design, execution and analysis of GWAS for the MHC region and not infrequently, the MHC region has even been excluded from the analysis of GWAS data. Here, we contrast features of the MHC region with the rest of the genome and highlight its complexities, including its functional polymorphisms beyond those determined by single nucleotide polymorphisms or single amino acid residues. One of the several issues with customary GWAS analysis is that it does not address this additional layer of polymorphisms unique to the MHC region. We highlight alternative approaches that may assist with the analysis of GWAS data from the MHC region and unravel associations with all functional polymorphisms beyond single SNPs. We suggest that despite already showing the highest number of disease associations, the true extent of the involvement of the MHC region in disease genetics may not have been uncovered.

Hidden genomic MHC disparity between HLA-matched sibling pairs in hematopoietic stem cell transplantation

Matching classical HLA alleles between donor and recipient is an important factor in avoiding adverse immunological effects in HSCT. Siblings with no differences in HLA alleles, either due to identical-by-state or identical-by-descent status, are considered to be optimal donors. We carried out a retrospective genomic sequence and SNP analysis of 336 fully HLA-A, -B, -DRB1 matched and 14 partially HLA-matched sibling HSCT pairs to determine the level of undetected mismatching within the MHC segment as well as to map their recombination sites. The genomic sequence of 34 genes locating in the MHC region revealed allelic mismatching at 1 to 8 additional genes in partially HLA-matched pairs. Also, fully matched pairs were found to have mismatching either at HLA-DPB1 or at non-HLA region within the MHC segment. Altogether, 3.9% of fully HLA-matched HSCT pairs had large genomic mismatching in the MHC segment. Recombination sites mapped to certain restricted locations. The number of mismatched nucleotides correlated with the risk of GvHD supporting the central role of full HLA matching in HSCT. High-density genome analysis revealed that fully HLA-matched siblings may not have identical MHC segments and even single allelic mismatching at any classical HLA gene often implies larger genomic differences along MHC.

Multi-allelic exact tests for Hardy-Weinberg equilibrium that account for gender

Statistical tests for Hardy–Weinberg equilibrium are important elementary tools in genetic data analysis. X‐chromosomal variants have long been tested by applying autosomal test procedures to females only, and gender is usually not considered when testing autosomal variants for equilibrium. Recently, we proposed specific X‐chromosomal exact test procedures for bi‐allelic variants that include the hemizygous males, as well as autosomal tests that consider gender. In this study, we present the extension of the previous work for variants with multiple alleles. A full enumeration algorithm is used for the exact calculations of tri‐allelic variants. For variants with many alternate alleles, we use a permutation test. Some empirical examples with data from the 1,000 genomes project are discussed.

Departure from Hardy Weinberg Equilibrium and Genotyping Error

Objective: Departure from Hardy Weinberg Equilibrium (HWE) may occur due to a variety of causes, including purifying selection, inbreeding, population substructure, copy number variation or genotyping error. We searched for specific characteristics of HWE-departure due to genotyping error. Methods: Genotypes of a random set of genetic variants were obtained from the Exome Aggregation Consortium (ExAC) database. Variants with <80% successful genotypes or with minor allele frequency (MAF) <1% were excluded. HWE-departure (d-HWE) was considered significant at p < 10E-05 and classified as d-HWE with loss of heterozygosity (LoH d-HWE) or d-HWE with excess heterozygosity (gain of heterozygosity: GoH d-HWE). Missing genotypes, variant type (single nucleotide polymorphism (SNP) vs. insertion/deletion); MAF, standard deviation (SD) of MAF across populations (MAF-SD) and copy number variation were evaluated for association with HWE-departure. Results: The study sample comprised 3,204 genotype distributions. HWE-departure was observed in 134 variants: LoH d-HWE in 41 (1.3%), GoH d-HWE in 93 (2.9%) variants. LoH d-HWE was more likely in variants located within deletion polymorphisms (p < 0.001) and in variants with higher MAF-SD (p = 0.0077). GoH d-HWE was associated with low genotyping rate, with variants of insertion/deletion type and with high MAF (all at p < 0.001). In a sub-sample of 2,196 variants with genotyping rate >98%, LoH d-HWE was found in 29 (1.3%) variants, but no GoH d-HWE was detected. The findings of the non-random distribution of HWE-violating SNPs along the chromosome, the association with common deletion polymorphisms and indel-variant type, and the finding of excess heterozygotes in genomic regions that are prone to cross-hybridization were confirmed in a large sample of short variants from the 1,000 Genomes Project. Conclusions: We differentiated between two types of HWE-departure. GoH d-HWE was suggestive for genotyping error. LoH d-HWE, on the contrary, pointed to natural variabilities such as population substructure or common deletion polymorphisms.