Differences In Allele Frequencies Research Articles

Identification of genetic variants by using a whole exome sequencing approach in Taiwanese patients with Brugada Syndrome

Abstract Background Brugada Syndrome (BrS) is a life-threatening arrhythmia disease primarily affecting young men. Many genetic studies have been conducted mainly in Caucasian BrS patients to identify potential genetic drivers. SCN5A mutations are less common in Asian BrS patients (around 10%) compared to Caucasian patients (20-30%). This means the pathogenic process of more than 70% of BrS patients remains unclear. With the advancement in NGS technologies, we can conduct genome-wide screening of DNA variants in BrS patients efficiently. Purpose Since less than 30% of the BrS patients have known genetic drivers for the disease, we aimed to identify more DNA variants in BrS patients by using the whole exome sequencing (WES) approach. Furthermore, we evaluated the landscape of DNA variants in the two groups of BrS patients according to whether they have SCN5A mutations or not. Methods A total of 201 BrS patients were recruited for the WES experiments. After obtaining the raw sequencing data, we performed the quality check procedures using the multiQC and Trimmomatic softwares (Figure 1). For reads passing the quality check, the BWA algorithm was used for the alignment. Next, the samtools, bamtools, and Picard methods were used to remove duplicates and perform the sorting. Subsequently, those sorted sequencing files were used to call the DNA variants by using the GATK approach and annotate them using the ANNOVAR method. The 201 BrS patients were classified into two different groups based on their mutation status of SCN5A, presentation of sudden cardiac death (SCD), and SCD along with syncope. The gene-based algorithm, SKAT, was used to analyze the summarized effects of DNA variants located in a single gene. For each DNA variant, Fisher's exact test was used to compare the allele frequency in the BrS patients under their groupings. Results The WES approach in the 201 BrS patients revealed a total of 205,454 variants, which were further categorized into two groups, common variants (51,867) and rare variants (153,587), based on their allele frequencies (5%). The SKAT algorithm was performed in the 30 genes that have been previously reported to be associated with BrS. No significant differences were observed in any gene in BrS patients for the common variants. Intriguingly, KCNJ8, SCN1B, ZFPM2, MAPRE2, and KCNE5 showed significant differences in allele frequencies between the SCN5A+ and SCN5A- patients (Figure 2). Fisher's exact test further revealed that three DNA variants (chr19:35524980, chr12:21926440, chr6:126210346) were significantly enriched in SCN5A+ patients. The 5 DNA variants located in ZFPM2 showed a linkage disequilibrium block and were enriched in SCN5A+ patients. Conclusions The WES approach in BrS patients revealed DNA variants that were reported in previous studies. Several DNA variants were significantly enriched in SCN5A+ patients, suggesting that the genome-wide mutation landscape was different in SCN5A+ and SCN5A- BrS patients.Figure 1Figure 2

Association of polymorphism of apolipoprotein e gene with acute myocardial infarction and its short term clinical outcome

Abstract Background/Introduction Myocardial infarction (MI) is a multifactorial disease influenced by environmental and genetic factors. Apolipoprotein E (ApoE) is a protein incorporated into serum lipoproteins directing their catabolism via binding to receptors. It is a structural component of both chylomicrons and very low density lipoprotein (VLDL) remnants. ApoE is thought to play a central role in atherosclerosis by participating in overall plasma cholesterol homeostasis. In addition to other classical risk factors, positive family history of cardiovascular disease (CVD) has been found to be significantly and independently associated with MI in epidemiological studies (1). Besides traditional cardiovascular risk factors it is likely that additional genetic factors like ApoE gene polymorphism may interact with environmental factors to determine overall cardiovascular risk of an individual to develop MI. Purpose The purpose of our study was to determine the risk of association of variants of ApoE gene with Acute Myocardial Infarction. Moreover, we also wanted to study the impact of ApoE gene polymorphism in short term major adverse cardiac events (MACE) in patients with MI not willing for coronary intervention. Methods In this case control study a total of 200 consecutively admitted AMI patients admitted to the coronary care unit of our hospital were enrolled to the study along with 200 age and gender matched controls admitted during the same time period . ApoE gene polymorphism was detected by DNA extraction, DNA amplification by PCR in a DNA thermal cycler followed by restriction enzyme digestion(2). Chi square test was used to investigate the difference in genotype and allele frequencies in casend controls and also between observed and expected frequencies in the realm of Hardy-Weinberg equilibrium. Odds Ratio (OR) with 95 % confidence interval and p values were calculated while analysing the association of different genotypes of ApoE genetic polymorphism with the disease and MACE. Logistic regression analysis was carried out amongst the groups under consideration to find the unadjusted OR and adjusted OR of the genotypes in relation to the disease.p<0.05 was considered as significant Analysis was done in SPSS software. Results Besides the traditional risk factors Apo E4 allele (OR 2.298, CI 1.305-4.044 ) and Apo E3/E4 genotype (OR 2.116, CI 1.008-4.443, p=0.048) showed strong of association with the disease. There was a higher population of Apo E3/E4 genotypes and Apo E2/E3 genotypes with MACE events (25.6% and 28.6%), while the common ApoE3/E3 genotypes had a MACE event of 10.5 %. However, the difference was not significant (p=0.071, Χ 2 =8.63 df =4). Similarly, ApoE4 allele showed a non significantly higher MACE events (p=0.22, Χ 2 =3.01 df =4) Conclusion ApoE4 allele and ApoE3/E4 genotype are strongly and independently associated with AMI.RISK FACTORS AND RISK OF MIGenotype and MACE

That's Not a Hybrid: How to Distinguish Patterns ofAdmixture and Isolation By Distance.

Describing naturally occurring genetic variation is a fundamental goal of molecular phylogeography and population genetics. Popular methods for this task include STRUCTURE, a model-based algorithm that assigns individuals to genetic clusters, and principal component analysis (PCA), a parameter-free method. The ability of STRUCTURE to infer mixed ancestry makes it popular for documenting natural hybridisation, which is of considerable interest to evolutionary biologists, given that such systems provide a window into the speciation process. Yet, STRUCTURE can produce misleading results when its underlying assumptions are violated, like when genetic variation is distributed continuously across geographic space. To test the ability of STRUCTURE and PCA to accurately distinguish admixture from continuous variation, we use forward-time simulations to generate population genetic data under three demographic scenarios: two involving admixture and one with isolation by distance (IBD). STRUCTURE and PCA alone cannot distinguish admixture from IBD, but complementing these analyses with triangle plots, which visualise hybrid index against interclass heterozygosity, provides more accurate inference of demographic history, especially in cases of recent admixture. We demonstrate that triangle plots are robust to missing data, while STRUCTURE and PCA are not, and show that setting a low allele frequency difference threshold for ancestry-informative marker (AIM) identification can accurately characterise the relationship between hybrid index and interclass heterozygosity across demographic histories of admixture and range expansion. While STRUCTURE and PCA provide useful summaries of genetic variation, results should be paired with triangle plots before admixture is inferred.

Investigating Leptin Gene Variants and Methylation Status in Relation to Breastfeeding and Preventing Obesity

Objective: We investigated whether the results of leptin gene (LEP) 2548G/A (rs7799039) and leptin receptor gene (LEPR) 668 A/G (rs1137101) variants, as well as the methylation analysis of CpG regions at nucleotides −31 and −51 of the LEP gene, showed any differences between breastfed and non-breastfed children in this study. Materials and Methods: The cross-sectional study included 100 children aged 2–5 years who were attending nursery and kindergarten and had been accepted to the Department of General Paediatrics. Infants who were exclusively breastfed for the first six months after birth constituted the study group, and those who were not only breastfeed constituted the control group. Methylation percentages at CpG islands of the LEP gene were compared between exclusively breastfed and non-exclusively breastfed infants, and the statistical significance was analyzed by looking for changes in LEP −31 and −51 nt methylation and LEP 2548G/A ve LEPR 668 A/G variants. Results: Both groups were compared by feeding, and the association of LEPR and LEP gene polymorphisms and −51 nt and −31 nt methylations were analyzed. There were no significant differences between the groups regarding genotype and allele frequency for the LEPR 668 A/G, LEP 2548 G/A gene variant, −31 nt methylation, and −51 nt methylation status. Similarly, there was no significant difference in genotype and allele frequency for the LEPR 668 A/G gene variant in terms of duration of exclusive breastfeeding, total breastfeeding, body mass index, family obesity, and satiety status. However, maternal support from family elders and physical activity increased the 51 nt methylation, but this methylation was not significantly affected by BMI, age, or satiety status. Conclusions: Maternal support from family elders and physical activity were associated with increased 51 nt methylation, but this methylation was not significantly affected by BMI, age, or satiety status. However, there are not enough studies in this area to reach a definitive conclusion, and further research is needed.

Insights into trait-association of selection signatures and adaptive eQTL in indigenous African cattle

BackgroundAfrican cattle represent a unique resource of genetic diversity in response to adaptation to numerous environmental challenges. Characterising the genetic landscape of indigenous African cattle and identifying genomic regions and genes of functional importance can contribute to targeted breeding and tackle the loss of genetic diversity. However, pinpointing the adaptive variant and determining underlying functional mechanisms of adaptation remains challenging.ResultsIn this study, we use selection signatures from whole-genome sequence data of eight indigenous African cattle breeds in combination with gene expression and quantitative trait loci (QTL) databases to characterise genomic targets of artificial selection and environmental adaptation and to identify the underlying functional candidate genes. In general, the trait-association analyses of selection signatures suggest the innate and adaptive immune system and production traits as important selection targets. For example, a large genomic region, with selection signatures identified for all breeds except N’Dama, was located on BTA27, including multiple defensin DEFB coding-genes. Out of 22 analysed tissues, genes under putative selection were significantly enriched for those overexpressed in adipose tissue, blood, lung, testis and uterus. Our results further suggest that cis-eQTL are themselves selection targets; for most tissues, we found a positive correlation between allele frequency differences and cis-eQTL effect size, suggesting that positive selection acts directly on regulatory variants.ConclusionsBy combining selection signatures with information on gene expression and QTL, we were able to reveal compelling candidate selection targets that did not stand out from selection signature results alone (e.g. GIMAP8 for tick resistance and NDUFS3 for heat adaptation). Insights from this study will help to inform breeding and maintain diversity of locally adapted, and hence important, breeds.

Genomic association of SNPs rs4077582 of CYP11A1 and rs700519 of CYP19A1 genes with polycystic ovarian syndrome.

Polycystic Ovarian Syndrome (PCOS) is a complex endocrine disorder that affects women of reproductive age. Several candidate genes have been shown to be associated with PCOS. Previous studies have shown that variations in CYP11A1 and CYP19A1 genes are associated with hormonal dysregulation associated with PCOS in different ethnic populations. This study aims to investigate the genomic association between SNPs rs4077582 of CYP11A1 and rs700519 of CYP19A1 and the development of PCOS in Pakistani population. A total of 280 subjects were recruited for the study, including 142 PCOS cases diagnosed based on Rotterdam criteria and 138 age-matched controls. The anthropometric, hormonal and biochemical parameters of all subjects were analyzed. Genomic DNA was extracted and genotyping of the selected SNPs was performed using Sanger sequencing. Further, we also examined the genotypic-phenotypic correlation analysis for various clinical and biochemical parameters for SNP rs4077582 of CYP11A1. We found significant differences in allele frequency (OR = 0.42, 95% CI = 0.30-0.60, χ2 = 16.3693, p = 0.000052) and genotypic frequency (χ2 = 26.4376, p = 0.00001) between PCOS women and controls for SNP rs4077582 of CYP11A1. Genotype-phenotype correlation analysis showed a significant difference in FAI (p = 0.005), testosterone (p = 0.001), androstenedione (p = 0.005) and urea (p = 0.049) levels between the three genotypes. No association between SNP rs700519 of CYP19A1 and PCOS was observed. Our results suggest the role of SNP rs4077582 of CYP11A1 gene in the clinical manifestation of PCOS in Pakistani women.

A Dominance Hypothesis Argument for Historical Genetic Gains and the Fixation of Heterosis in Octoploid Strawberry.

Heterosis was the catalyst for the domestication of cultivated strawberry (Fragaria × ananassa), an interspecific hybrid species that originated in the 1700s. The hybrid origin was discovered because the phenotypes of spontaneous hybrids transgressed those of their parent species. The transgressions included fruit yield increases and other genetic gains in the twentieth century that sparked the global expansion of strawberry production. The importance of heterosis to the agricultural success of the hybrid species, however, has remained a mystery. Here we show that heterosis has disappeared (become fixed) among improved hybrids within a population (the California population) that has been under long-term selection for increased fruit yield, weight, and firmness. We found that the highest yielding hybrids are among the most highly inbred (59-79%), which seems counterintuitive for a highly heterozygous, outbreeder carrying heavy genetic loads. Although faint remnants of heterosis were discovered, the between-parent allele frequency differences and dispersed favorable dominant alleles necessary for heterosis have decreased nearly genome-wide within the California population. Conversely, heterosis was prevalent and significant among wide hybrids, especially for fruit count, a significant driver of genetic gains for fruit yield. We attributed the disappearance (fixation) of heterosis within the California population to increased homozygosity of favorable dominant alleles and inbreeding associated with selection, random genetic drift, and selective sweeps. Despite historical inbreeding, the highest yielding hybrids reported to-date are estimated to be heterozygous for 20,370-44,280 of 97,000-108,000 genes in the octoploid genome, the equivalent of an entire diploid genome or more.

Genotype distribution and allele frequency of thioester-containing protein 1(Tep1) and its effect on development of Plasmodium oocyst in populations of Anopheles arabiensis in Ethiopia.

Thioester-containing protein 1 (TEP1) is a crucial component of mosquitoes' natural resistance to parasites. To effectively combat malaria, there is a need to better understand how TEP1 polymorphism affects phenotypic traits during infections. Therefore, the purpose of this study was to determine the Tep1 genotype frequency in malaria vector populations from south-western Ethiopia and investigate its effect on Plasmodium oocyst development in Anopheles arabiensis populations. Using standard dippers, Anopheles mosquito larvae were collected from aquatic habitats in Asendabo, Arjo Dedessa, and Gambella in 2019 and 2020. Collected larvae were reared to adults and identified morphologically. Female An. gambiae s.l. were allowed to feed on infected blood containing the same number of gametocytes obtained from P. falciparum and P. vivax gametocyte-positive individuals using indirect membrane feeding methods. Polymerase Chain Reaction (PCR) was used to identify An. gambiae s.l. sibling species. Three hundred thirty An. gambiae s.l. were genotyped using Restricted Fragment Length Polymorphism (RFLP) PCR and sub samples were sequenced to validate the TEP1 genotyping. Among the 330 samples genotyped, two TEP1 alleles, TEP1*S1 (82% frequency) and TEP1*R1 (18% frequency), were identified. Three equivalent genotypes, TEP1*S1/S1, TEP1*R1/R1, and TEP1*S1/R1, had mean frequencies of 65.15%, 2.12%, and 32.73%, respectively. The nucleotide diversity was ranging from 0.36554 to 0. 46751 while haplotype diversity ranged from 0.48871 to 0.63161, across all loci. All sample sites had positive Tajima's D and Fu's Fs values. There was a significant difference in the TEP1 allele frequency and genotype frequency among mosquito populations (p < 0.05), except populations of Anopheles arabiensis from Asendabo and Gambella (p > 0.05). In addition, mosquitoes with the TEP1 *RR genotype were susceptible and produced fewer Plasmodium oocysts than mosquitoes with the TEP1 *SR and TEP1 *SS genotypes. The alleles identified in populations of An. arabiensis were TEP1*R1 and TEP1*S1. There was no significant variation in TEP1*R1 allele frequency between the high and low transmission areas. Furthermore, An. arabiensis carrying the TEP1*R1 allele was susceptible to Plasmodium infection. Further studies on vector-parasite interactions, particularly on the TEP1 gene, are required for vector control techniques.

Heterogeneous genetic architectures of prostate cancer susceptibility in sub-Saharan Africa.

Men of African descent have the highest prostate cancer incidence and mortality rates, yet the genetic basis of prostate cancer in African men has been understudied. We used genomic data from 3,963 cases and 3,509 controls from Ghana, Nigeria, Senegal, South Africa and Uganda to infer ancestry-specific genetic architectures and fine-map disease associations. Fifteen independent associations at 8q24.21, 6q22.1 and 11q13.3 reached genome-wide significance, including four new associations. Intriguingly, multiple lead associations are private alleles, a pattern arising from recent mutations and the out-of-Africa bottleneck. These African-specific alleles contribute to haplotypes with odds ratios above 2.4. We found that the genetic architecture of prostate cancer differs across Africa, with effect size differences contributing more to this heterogeneity than allele frequency differences. Population genetic analyses reveal that African prostate cancer associations are largely governed by neutral evolution. Collectively, our findings emphasize the utility of conducting genetic studies that use diverse populations.

Cross-ancestry analysis of brain QTLs enhances interpretation of schizophrenia genome-wide association studies.

Research on brain expression quantitative trait loci (eQTLs) has illuminated the genetic underpinnings of schizophrenia (SCZ). Yet most of these studies have been centered on European populations, leading to a constrained understanding of population diversities and disease risks. To address this gap, we examined genotype and RNA-seq data from African Americans (AA, n= 158), Europeans (EUR, n= 408), and East Asians (EAS, n= 217). When comparing eQTLs between EUR and non-EUR populations, we observed concordant patterns of genetic regulatory effect, particularly in terms of the effect sizes of the eQTLs. However, 343,737 cis-eQTLs linked to 1,276 genes and 198,769 SNPs were found to be specific to non-EUR populations. Over 90% of observed population differences in eQTLs could be traced back to differences in allele frequency. Furthermore, 35% of these eQTLs were notably rare in the EUR population. Integrating brain eQTLs with SCZ signals from diverse populations, we observed a higher disease heritability enrichment of brain eQTLs in matched populations compared to mismatched ones. Prioritization analysis identified five risk genes (SFXN2, VPS37B, DENR, FTCDNL1, and NT5DC2) and three potential regulatory variants in known risk genes (CNNM2, MTRFR, and MPHOSPH9) that were missed in the EUR dataset. Our findings underscore that increasing genetic ancestral diversity is more efficient for power improvement than merely increasing the sample size within single-ancestry eQTLs datasets. Such a strategy will not only improve our understanding of the biological underpinnings of population structures but also pave the way for the identification of risk genes in SCZ.

Advancing pharmacogenomic research in US Hmong populations: prevalence of key single nucleotide variations in California Hmong.

In collaboration with the Minnesota Hmong community, we have previously discovered significant differences in allele frequencies for key Single Nucleotide Variations (SNVs) within Very Important Pharmacogenes (VIPs) between Hmong and East Asians. Recognizing the potential clinical implications of these observed differences, we sought to validate these observations in a Hmong cohort residing in California, the state with the largest Hmong population in the US. Robust validation of these differences would affect motivation for clinicians treating individuals who identify as Hmong to consider pharmacogenomic (PGx) testing as a means to improve clinical decision making when using therapeutic agents in this unique population. Guided by California Hmong community leaders and utilizing the basic approach of community-based participatory research, demographic, clinical information and a buccal swab was obtained from Hmong adults residing in California. A commercial PGx testing panel was performed on these samples and specific allele frequencies of interest were compared between California and Minnesota Hmong. Allele frequency differences between California Hmong, East Asians and Europeans, were also compared. Return-of-PGx-results and presentations of group data were made to members of the Hmong along with PGx educational sessions to help interpret the observations. In 118 California Hmong who completed the study, the allele frequencies for SNV's were similar to previous Minnesota Hmong results. Furthermore, out of the 18 SNVs that were not previously reported in Hmong, allele frequencies were statistically different in 38% (7/18) of SNVs comparing California Hmong to East Asians, and in 77.8% (14/18) SNVs comparing California Hmong to Europeans. These results validate the original study's findings that Hmong people living in different US locations have similar allele frequencies for key PGx genes. Further, for many of these PGx genes, their allele frequencies are significantly different compared to either East Asians or Europeans. Clinicians should consider these important differences when prescribing medications for people who identify as Hmong.

Abstract C051: Impact of population pharmacogenomics on cisplatin-induced neurotoxicities

Abstract Background: Studies have shown that taxane-treated breast cancer patients with African (AFR) ancestry experience worse peripheral sensory neuropathy (PSN), while African American vincristine-treated leukemia patients have less PSN than Europeans (EUR). We hypothesized that ancestral disparities may exist for cisplatin-induced neurotoxicities and may relate to differences in allele frequency and resultant gene expression of SNPs associated with these toxicities. Methods: In our Platinum Study of 1680 genotyped testicular cancer survivors (TCS), multidimensional scaling scores were anchored by the 1000 Genomes population to classify geographic ancestry. PSN (EORTC-CIPN20) and other toxicities plus clinical and lifestyle variables were examined for associations between geographic ancestry and phenotypes using logistic regression and Wilcoxon rank sum test. Genotyped SNPs with Fst &amp;gt; 0.25 in 1000 Genomes were identified and evaluated with GTEx data to find expression or splicing quantitative trait loci in nerve/brain tissue or for cisplatin-associated candidate genes, across all tissues. LDmatrix was used to find linkage disequilibrium. If a block of SNPs had R2 &amp;gt; 0.75, one independent SNP was kept for association analysis. Logistic regression with age at questionnaire completion and 10 genetic principal components as covariates was used to calculate the association between genotype and toxicity. Results: In an analysis of patients who received 400-450 mg/m2 of cisplatin including EUR (n=681), Asian (ASN) axis (n=44) and AFR (n=13), the odds ratio (OR) for any neuropathy vs. none in the AFR vs. EUR was 7.8 (P=0.049) and in the AFR vs. ASN axis (including East Asian 1000 Genomes populations) was 10.0 (P=0.034). TCS with AFR ancestry had significantly more vertigo. There were 19,874 SNPs for analysis from the genome-wide approach and 118 SNPs from the candidate gene approach. While associations between genotype and toxicity did not reach Bonferroni threshold, some SNPs had p-values in the 10-5 range. For PSN, one SNP had increased frequency of a risk allele in the AFR population while the other had increased frequency of a protective allele, while for vertigo, four SNPs had increased frequency of risk alleles in the AFR population. Of interest was rs34904346, in which the TT vs. AA genotype had an OR of 1.9 (p=2x10-5) for any vs. no PSN, with the T allele having 90% frequency in AFR ancestry patients versus ∼60% in other populations. This SNP and three other independent SNPs associated with PSN incidence (p∼10-3) are eQTLs for RNF24 in nerve tissue, which interacts with membrane proteins that regulate intracellular calcium levels. Conclusions: We show preliminary evidence for the role of differing risk allele frequencies across populations in explaining disparities in cisplatin-induced PSN and vertigo, likely in combination with structural and environmental factors. Genotyping risk alleles could customize dosing, treatment, and post-treatment monitoring to benefit all cisplatin-treated patients. Citation Format: Swetha Nakshatri, Paul C. Dinh, Darren R. Feldman, Robert J. Hamilton, Chunkit Fung, David J. Vaughn, Christian Kollmannsberger, Robert Huddart, Lawrence H. Einhorn, Nancy J. Cox, Lois B. Travis, Eileen Dolan. Impact of population pharmacogenomics on cisplatin-induced neurotoxicities [abstract]. In: Proceedings of the 17th AACR Conference on the Science of Cancer Health Disparities in Racial/Ethnic Minorities and the Medically Underserved; 2024 Sep 21-24; Los Angeles, CA. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2024;33(9 Suppl):Abstract nr C051.

Ancestry-Related Differences in Allele and Genotype Frequencies of EGF A61G Polymorphism in the Cuban Population

Background: The polymorphism rs4444903 of the Epidermal Growth Factor gene (EGF A61G) causes differences in the EGF serum levels. It has become a biomarker for genetic susceptibility to cancer and a pharmacogenomic marker for therapies involving the EGF/EGF-receptor pathway. Objective: The present study aimed to characterize the allele and genotype frequencies of the rs4444903 in a Cuban sample and its relationship to a specific genetic ancestry. Methods: A cross-sectional study was carried out. Genomic data was collected from a dense genome-wide genotyping array analysis of 948 Cubans from all provinces. The allele and genotype frequencies of the rs4444903 were calculated. Analysis of ancestryrelated allelic/genotypic differences was performed. Results: The frequencies for both alleles were found to be very similar (0.52 for G vs. 0.48 for A allele), and genotype frequencies were 24.3%, 47.9%, and 27.8% for AA, AG, and GG, respectively. Greater differences were found between Cuban provinces, with frequencies for the G allele ranging from 0.38 in Artemisa to 0.69 in Guantánamo and for the GG genotype from 14.29% in Mayabeque to 50.88% in Guantánamo. An increased African-ancestry proportion was related to a higher probability of carrying G allele and GG genotype, with a significant (p=0.0038, q=0.024) African-ancestry-enrichment pattern. Conclusion: African ancestry seems to contribute to an increase in the EGF61*G allele in Cubans. Geographic patterns in admixture proportions for African and European ancestry are a determinant factor in the allelic and genotypic frequency differences between Cuban provinces. Such differences should be observed when designing association studies and implementing therapeutic approaches based on the EGF/EGF receptor pathway in Cuba.

Does HLA explain the high incidence of childhood-onset type 1 diabetes in the Canary Islands? The role of Asp57 DQB1 molecules

The Canary Islands inhabitants, a recently admixed population with significant North African genetic influence, has the highest incidence of childhood-onset type 1 diabetes (T1D) in Spain and one of the highest in Europe. HLA accounts for half of the genetic risk of T1D.AimsTo characterize the classical HLA-DRB1 and HLA-DQB1 alleles in children from Gran Canaria with and without T1D.MethodsWe analyzed classic HLA-DRB1 and HLA-DQB1 alleles in childhood-onset T1D patients (n = 309) and control children without T1D (n = 222) from the island of Gran Canaria. We also analyzed the presence or absence of aspartic acid at position 57 in the HLA-DQB1 gene and arginine at position 52 in the HLA-DQA1 gene. Genotyping of classical HLA-DQB1 and HLA-DRB1 alleles was performed at two-digit resolution using Luminex technology. The chi-square test (or Fisher's exact test) and odds ratio (OR) were computed to assess differences in allele and genotype frequencies between patients and controls. Logistic regression analysis was also used.ResultsMean age at diagnosis of T1D was 7.4 ± 3.6 years (46% female). Mean age of the controls was 7.6 ± 1.1 years (55% female). DRB1*03 (OR = 4.2; p = 2.13–13), DRB1*04 (OR = 6.6; p ≤ 2.00–16), DRB1* 07 (OR = 0.37; p = 9.73–06), DRB1*11 (OR = 0.17; p = 6.72–09), DRB1*12, DRB1*13 (OR = 0.38; p = 1.21–05), DRB1*14 (OR = 0.0; p = 0.0024), DRB1*15 (OR = 0.13; p = 7.78–07) and DRB1*16 (OR = 0.21; p = 0.003) exhibited significant differences in frequency between groups. Among the DQB1* alleles, DQB1*02 (OR: 2.3; p = 5.13–06), DQB1*03 (OR = 1.7; p = 1.89–03), DQB1*05 (OR = 0.64; p = 0.027) and DQB1*06 (OR = 0.19; p = 6.25–14) exhibited significant differences. A total of 58% of the studied HLA-DQB1 genes in our control population lacked aspartic acid at position 57.ConclusionsIn this population, the overall distributions of the HLA-DRB1 and HLA-DQB1 alleles are similar to those in other European populations. However, the frequency of the non-Asp-57 HLA-DQB1 molecules is greater than that in other populations with a lower incidence of T1D. Based on genetic, historical and epidemiological data, we propose that a common genetic background might help explain the elevated pediatric T1D incidence in the Canary Islands, North-Africa and middle eastern countries.

Selection with two alleles of X-linkage and its application to the fitness component analysis of OdsH in Drosophila.

In organisms with the XY sex-determination system, there is an imbalance in the inheritance and transmission of the X chromosome between males and females. Unlike an autosomal allele, an X-linked recessive allele in a female will have phenotypic effects on its male counterpart. Thus, genes located on the X chromosome are of particular interest to researchers in molecular evolution and genetics. Here we present a model for selection with two alleles of X-linkage to understand fitness components associated with genes on the X chromosome. We apply this model to the fitness analysis of an X-linked gene, OdsH (16D), in the fruit fly Drosophila melanogaster. The function of OdsH is involved in sperm production and the gene is rapidly evolving under positive selection. Using site-directed gene targeting, we generated functional and defective OdsH variants tagged with the eye-color marker gene white. We compare the allele frequency changes of the two OdsH variants, each directly competing against a wild-type OdsH allele in concurrent but separate experimental populations. After 20 generations, the two genetically modified OdsH variants displayed a 40% difference in allele frequencies, with the functional OdsH variant demonstrating an advantage over the defective variant. Using maximum likelihood estimation, we determined the fitness components associated with the OdsH alleles in males and females. Our analysis revealed functional aspects of the fitness determinants associated with OdsH, and that sex-specific fertility and viability consequences both contribute to selection on an X-linked gene.

Haplotype-Aware Detection of SERPINA1 Variants by Nanopore Sequencing

Alpha-1 antitrypsin (AAT) is an acute-phase reactant with immunomodulatory properties that mainly inhibits neutrophil elastase. Low serum levels cause AAT deficiency (AATD), an underdiagnosed condition that predisposes to pulmonary and hepatic diseases. The SERPINA1 gene, which encodes AAT, contains more than 500 variants. PI*Z and PI*S alleles are the most diagnosed causes of AATD, but the role of the SERPINA1 haplotypes in AAT function remains unknown. SERPINA1 gene was PCR amplified from 94 asthma patients, using primers with tails for indexing. Sequencing libraries were loaded into a MinION-Mk1C, and MinKNOW was used for basecalling and demultiplexing. Nanofilt and Minimap2 were used for filtering and mapping/alignment. Variant calling/phasing were performed with PEPPER-Margin-DeepVariant. SERPINA1 gene was 100% covered for all samples, with a minimum sequencing depth of 500X. 75 single nucleotide variants and 4 indels were detected, with 45 and 2 of them highly polymorphic (MAF>0.1), respectively. Nine of the SNVs showed differences in allele frequencies when compared with the overall Spanish population. More than 90% of heterozygous SNVs were phased, yielding 91 and 58 different haplotypes for each SERPINA1 amplified region. Haplotype-based Linkage Disequilibrium (LD) analysis suggests that a recombination hotspot could generate variation in the SERPINA1 gene. The proposed workflow enables haplotype-aware genotyping of the SERPINA1 gene by nanopore sequencing, which will allow the development of novel AATD diagnostic strategies.

Molecular characteristics of NPM1 mutations in AML patients using targeted NGS

ABSTRACT Background Acute myeloid leukemia (AML) is a clonal hematopoietic stem cell malignancy characterized by the accumulation of myeloid progenitor cells with arrested differentiation, leading to the suppression of normal hematopoiesis. The disease exhibits significant heterogeneity, with less than 30% five-year overall survival. Nucleophosmin 1 (NPM1) is the most commonly mutated gene in adult AML, found in approximately 25–35% of patients, and is recognized as a distinct subtype by the WHO. Methods The mutational status of NPM1 was assessed using the Oncomine Myeloid research panel on the Ion5S NGS platform (Thermo Fisher Scientific, USA). Results In our AML cohort, NPM1 was the most frequently mutated gene (75%). Indel mutations in the NPM1 gene, frequently resulting in frameshift consequences, show a statistically significant difference in variant allele frequency (VAF) compared to other mutation types. Specifically, NPM1-mutated cases exhibit a VAF that is positively correlated with initial bone marrow (BM) blast counts. Univariate logistic regression analysis reveals that a high VAF is significantly associated with non-complete remission (Non-CR), while a low VAF is significantly associated with complete remission (CR), indicating a statistically significant difference in outcomes (p = 0.012). Conclusion Targeted NGS technology enables the simultaneous analysis of multiple AML-relevant genes, offering a comprehensive mutational profile that supports risk stratification, diagnosis, and personalized treatment. Our findings highlight the significant impact of NPM1 mutations on patient outcomes, revealing their molecular characteristics and their correlation with different clinical parameters.

QTL-Seq identified a genomic region on chromosome 1 for soil-salinity tolerance in F2 progeny of Thai salt-tolerant rice donor line "Jao Khao".

Owing to advances in high-throughput genome sequencing, QTL-Seq mapping of salt tolerance traits is a major platform for identifying soil-salinity tolerance QTLs to accelerate marker-assisted selection for salt-tolerant rice varieties. We performed QTL-BSA-Seq in the seedling stage of rice from a genetic cross of the extreme salt-sensitive variety, IR29, and "Jao Khao" (JK), a Thai salt-tolerant variety. A total of 462 F2 progeny grown in soil and treated with 160 mM NaCl were used as the QTL mapping population. Two high- and low-bulk sets, based on cell membrane stability (CMS) and tiller number at the recovery stage (TN), were equally sampled. The genomes of each pool were sequenced, and statistical significance of QTL was calculated using QTLseq and G prime (G') analysis, which is based on calculating the allele frequency differences or Δ(SNP index). Both methods detected the overlapping interval region, wherein CMS-bulk was mapped at two loci in the 38.41-38.85 Mb region with 336 SNPs on chromosome 1 (qCMS1) and the 26.13-26.80 Mb region with 1,011 SNPs on chromosome 3 (qCMS3); the Δ(SNP index) peaks were -0.2709 and 0.3127, respectively. TN-bulk was mapped at only one locus in the overlapping 38.26-38.95 Mb region on chromosome 1 with 575 SNPs (qTN1) and a Δ(SNP index) peak of -0.3544. These identified QTLs in two different genetic backgrounds of segregating populations derived from JK were validated. The results confirmed the colocalization of the qCMS1 and qTN1 traits on chromosome 1. Based on the CMS trait, qCMS1/qTN1 stably expressed 6%-18% of the phenotypic variance in the two validation populations, while qCMS1/qTN1 accounted for 16%-20% of the phenotypic variance in one validation population based on the TN trait. The findings confirm that the CMS and TN traits are tightly linked to the long arm of chromosome 1 rather than to chromosome 3. The validated qCMS-TN1 QTL can be used for gene/QTL pyramiding in marker-assisted selection to expedite breeding for salt resistance in rice at the seedling stage.

The Quantitative Genetics of Human Disease: 2 Polygenic Risk Scores

In this the second of an anticipated four papers, we examine polygenic risk scores from a quantitative genetics perspective. In its most simplistic form, a polygenic risk score (PRS) analysis involves estimating the genetic effects of alleles in one study and then using those estimates to predict phenotype in another sample of individuals. Almost since the first application of these types of analyses it has been noted that PRSs often give unexpected and difficult-to-interpret results, particularly when applying effect-size estimates taken from individuals with ancestry very different than those to whom it is applied (applying PRSs across differing populations). To understand these seemingly perplexing observations, we deconstruct the effects of applying valid statistical estimates taken from one population to another when the two populations have differing allele frequencies at the sites contributing effect, when alleles with effects in one population are absent from the other, and finally when there is differing linkage disequilibrium (LD) patterns in the two populations. It will be shown that many of the seemingly most confusing results in the field are natural consequences of these factors. Given our best current understanding of human demographic history, most of the patterns seen in PRS analysis can be predicted as resulting from systematic differences in allele frequency and LD. Put the other way around, the most challenging and confusing results seen in cross population application of PRSs are likely to be the result of allele frequency and LD differences, not differences in the genetic effects of individual alleles. PRS analysis is an important tool both for understanding the genetic basis of complex phenotypes and, potentially, for identifying individuals at risk of developing disease before such disease manifests. As such it has the potential to be among the most important analysis frameworks in human genetics. Nevertheless, when a PRS is trained in people with one ancestry and then applied to people with another, the PRS’s behavior is often unpredictable, and sometimes is seemingly perverse. PRS distributions are often nearly non-overlapping between individuals with differing ancestry, i.e., odds ratios for unaffected people with one ancestry might be vastly larger than affected individuals from another. The correlation between a PRS and known phenotype might differ substantially, and sometimes the correlation is higher among people with ancestry different than the one used to create the PRS. Naively, one might conclude from these observations that the genetic basis of traits differs substantially among people of differing ancestry, and that the behavior of a PRS is difficult to predict when applied to new study populations. Differing definitions of genetic effect sizes are discussed, and key observations are made. It is shown that when populations differ in allele frequency, a locus affecting phenotype could have equal differences in allelic (additive) effects or equal additive variances, but not both. They cannot have equal additive effects, equal allelic penetrances, or equal odds ratios. PRS is defined, and its moments are derived. The effect of differing allele frequency and LD patterns is described. Perplexing PRS observations are discussed in light of theory and human demographic history. Suggestions for best practices for PRS construction are made. The most confusing results seen in cross population application of PRSs are often the predictable result of allele frequency and LD differences. There is relatively little evidence for systematic differences in the genetic basis of disease in individuals of differing ancestry, other than that which results from environmental, allele frequency, and LD differences.

Genome-wide SNP assessment of contemporary European red deer genetic structure highlights the distinction of peripheral populations and the main admixture zones in Europe.

Genome-wide technologies open up new possibilities to clarify questions on genetic structure and phylogeographic history of taxa previously studied with microsatellite loci and mitochondrial sequences. Here, we used 736 individual red deer (Cervus elaphus) samples genotyped at 35,701 single nucleotide polymorphism loci (SNPs) to assess the population structure of the species throughout Europe. The results identified 28 populations, with higher degrees of genetic distinction in peripheral compared to mainland populations. Iberian red deer show high genetic differentiation, with lineages in Western and Central Iberia maintaining their distinctiveness, which supports separate refugial ranges within Iberia along with little recent connection between Iberian and the remaining Western European populations. The Norwegian population exhibited the lowest variability and the largest allele frequency differences from mainland European populations, compatible with a history of bottlenecks and drift during post-glacial colonization from southern refugia. Scottish populations showed high genetic distance from the mainland but high levels of diversity. Hybrid zones were found between Eastern and Western European lineages in Central Europe as well as in the Pyrenees, where red deer from France are in close contact with Iberian red deer. Anthropogenic restocking has promoted the Pyrenean contact zone, admixture events in populations on the Isle of Rum and in the Netherlands, and at least partly the admixture of the two main lineages in central-eastern Europe. Our analysis enabled detailed resolution of population structure of a large mammal widely distributed throughout Europe and contributes to resolving the evolutionary history, which can also inform conservation and management policies.