Genotype-phenotype Association Analysis Research Articles

Genotype-phenotype association and functional analysis of hnRNPA1 mutations in amyotrophic lateral sclerosis

Background Pathogenic variants in hnRNPA1 have been reported in amyotrophic lateral sclerosis (ALS) patients. However, studies on hnRNPA1 mutant spectrum and pathogenicity of variants were rare. Methods We performed whole exome sequencing of ALS-associated genes and subsequent verification of rare variants in hnRNPA1 in our ALS patients. The hnRNPA1 mutations reported in literature were reviewed and combined with our results to determine the genotype-phenotype relationship. Functional analysis of the novel variant p.G195A was performed in vitro by transfection of mutant hnRNPA1 into 293T cell. Results Among 207 ALS patients recruited, 3 rare hnRNPA1 variants were identified (mutant frequency 1.45%), including two recurrent mutations (p.P340S and p.G283R), and a novel rare variant p.G195A. In combination with previous reports, there are 27 ALS patients with 15 hnRNPA1 mutations identified. Disease onset age was 47.90 ± 1.52 years with predominant limb onset. The p.P340S mutation caused flail arm syndrome (FAS) in two independent families with extended life expectancy. The newly identified p.G195A mutation, lying at the start of the PrLD (“prion-like” domain)/LCD (low-complexity domain), causes local structural changes in 3D protein prediction. Upon sodium arsenite exposure, mutant hnRNPA1 retained in the nucleus but deficit of cytoplasmic G3BP1-positive stress granule clearance was observed. This is different from the p.P340S mutation which caused both cytoplasmic translocation and stress granule formation. No cytoplasmic TDP-43 translocation was observed. Conclusion Mutations in hnRNPA1 are overall minor in ALS patients. The p.P340S mutation is associated with manifestation of FAS. Mutations in LCD of hnRNPA1 cause stress granule misprocessing.

Exploring the causality between intestinal flora and hyperplastic scars of human based on two-sample Mendelian randomization analysis

Objective: To investigate the causality between intestinal flora and hypertrophic scars (HS) of human. Methods: This study was a study based on two-sample Mendelian randomization (TSMR) analysis. The data on intestinal flora (n=18 473) and HS (n=208 248) of human were obtained from the genome-wide association study database. Genetically variable genes at five levels (phylum, class, order, family, and genus) of known intestinal flora, i.e., single nucleotide polymorphisms (SNPs), were extracted as instrumental variables for linkage disequilibrium (LD) analysis. Human genotype-phenotype association analysis was performed using PhenoScanner V2 database to exclude SNPs unrelated to HS in intestinal flora and analyze whether the selected SNPs were weak instrumental variables. The causal relationship between intestinal flora SNPs and HS was analyzed through four methods of TSMR analysis, namely inverse variance weighted (IVW), MR-Egger regression, weighted median, and weighted mode. Scatter plots of significant results from the four aforementioned analysis methods were plotted to analyze the correlation between intestinal flora SNPs and HS. Both IVW test and MR-Egger regression test were used to assess the heterogeneity of intestinal flora SNPs, MR-Egger regression test and MR-PRESSO outlier test were used to assess the horizontal multiplicity of intestinal flora SNPs, and leave-one-out sensitivity analysis was used to determine whether HS was caused by a single SNP in the intestinal flora. Reverse TSMR analyses were performed for HS SNPs and genus Intestinimonas or genus Ruminococcus2, respectively, to detect whether there was reverse causality between them. Results: A total of 196 known intestinal flora, belonging to 9 phyla, 16 classes, 20 orders, 32 families, and 119 genera, were obtained, and multiple SNPs were obtained from each flora as instrumental variables. LD analysis showed that the SNPs of the intestinal flora were consistent with the hypothesis that genetic variation was strongly associated with exposure factors, except for rs1000888, rs12566247, and rs994794. Human genotype-phenotype association analysis showed that none of the selected SNPs after LD analysis was excluded and there were no weak instrumental variables. IVW, MR-Egger regression, weighted median, and weighted mode of TSMR analysis showed that both genus Intestinimonas and genus Ruminococcus2 were causally associated with HS. Among them, forest plots of IVW and MR-Egger regression analyses also showed that 16 SNPs (the same SNPs number of this genus below) of genus Intestinimonas and 15 SNPs (the same SNPs number of this genus below) of genus Ruminococcus2 were protective factors for HS. Further, IVW analysis showed that genus Intestinimonas SNPs (with odds ratio of 0.62, 95% confidence interval of 0.41-0.93, P<0.05) and genus Ruminococcus2 SNPs (with odds ratio of 0.62, 95% confidence interval of 0.40-0.97, P<0.05) were negatively correlated with the risk of HS. Scatter plots showed that SNPs of genus Intestinimonas and genus Ruminococcus2 were protective factors of HS. Both IVW test and MR-Egger regression test showed that SNPs of genus Intestinimonas (with Q values of 5.73 and 5.76, respectively, P>0.05) and genus Ruminococcus2 (with Q values of 13.67 and 15.61, respectively, P>0.05) were not heterogeneous. MR-Egger regression test showed that the SNPs of genus Intestinimonas and genus Ruminococcus2 had no horizontal multiplicity (with intercepts of 0.01 and 0.06, respectively, P>0.05); MR-PRESSO outlier test showed that the SNPs of genus Intestinimonas and genus Ruminococcus2 had no horizontal multiplicity (P>0.05). Leave-one-out sensitivity analysis showed that no single intestinal flora SNP drove the occurrence of HS. Reverse TSMR analysis showed no reverse causality between HS SNPs and genus Intestinimonas or genus Ruminococcus2 (with odds ratios of 1.01 and 0.99, respectively, 95% confidence intervals of 0.97-1.06 and 0.96-1.04, respectively, P>0.05). Conclusions: There is a causal relationship between intestinal flora and HS of human, in which genus Intestinimonas and genus Ruminococcus2 have a certain effect on inhibiting HS.

Fine mapping of TFL, a major gene regulating fruit length in snake gourd (Trichosanthes anguina L).

Fruit length is a crucial agronomic trait of snake gourd (Trichosanthes anguina L); however, genes associated with fruit length have not been characterised. In this study, F2 snake gourd populations were generated by crossing the inbred lines, S1 and S2 (fruit lengths: 110 and 20cm, respectively). Subsequently, bulk segregant analysis, sequencing, and fine-mapping were performed on the F2 population to identify target genes. Our findings suggest that the fruit length of snake gourd is regulated by a major-effect regulatory gene. Mining of genes regulating fruit length in snake gourd to provide a basis for subsequent selection and breeding of new varieties. Genotype-phenotype association analysis was performed on the segregating F2 population comprising 6,000 plants; the results indicate that the target gene is located on Chr4 (61,846,126-61,865,087bp, 18.9-kb interval), which only carries the annotated candidate gene, Tan0010544 (designated TFL). TFL belongs to the MADS-box family, one of the largest transcription factor families. Sequence analysis revealed a non-synonymous mutation of base C to G at position 202 in the coding sequence of TFL, resulting in the substitution of amino acid Gln to Glu at position 68 in the protein sequence. Subsequently, an InDel marker was developed to aid the marker-assisted selection of TFL. The TFL in the expression parents within the same period was analysed using quantitative real-time PCR; the TFL expression was significantly higher in short fruits than long fruits. Therefore, TFL can be a candidate gene for determining the fruit length in snake gourd. Collectively, these findings improve our understanding of the genetic components associated with fruit length in snake gourds, which could aid the development of enhanced breeding strategies for plant species.

Primary carnitine deficiency: Estimation of prevalence in Chinese population and insights into newborn screening.

Primary carnitine deficiency (PCD) caused by pathogenic variants in the solute carrier family 22 member 5 (SLC22A5) gene is a rare autosomal recessive disease that results in defective fatty acid oxidation. PCD can be detected through tandem mass spectrometry (MS/MS), but transplacental transport of free carnitine from mothers may cause false negatives or positives during newborn screening (NBS). This study aimed to analyze the genetic characteristics of SLC22A5 and estimate the prevalence of PCD in the Chinese population, providing useful information for NBS and genetic counseling. We manually curated SLC22A5 pathogenic or likely pathogenic (P/LP) variants according to the American College of Medical Genetics and Genomics (ACMG) guidelines and identified 128P/LP variants. Based on the China Neonatal Genomes Project (CNGP), the estimated PCD prevalence was 1:17,456, which was higher than that in other populations. The genotype-phenotype association analysis showed that patients carrying homozygous c.760C>T and c.844C>T were more likely to present cardiomyopathy, whereas those carrying homozygous c.1400C>G were more likely to be asymptomatic (all p-values < 0.05). We found that there was no significant difference in initial C0 concentrations between patients and carriers, but there was a significant difference in the second-tier screening of C0 concentration between them (p-value < 0.05). We established a cost-effective variant panel containing 10 high-frequency sites and developed a screening algorithm incorporating gene panels with MS/MS, which could rescue one more patient who was undetected from MS/MS. In conclusion, the prevalence of PCD in the Chinese population is relatively high. The combination of conventional NBS with genetic sequencing is suggested for early diagnosis of PCD.

Identification of genetic associations and functional SNPs of bovine KLF6 gene on milk production traits in Chinese holstein

BackgroundOur previous research identified the Kruppel like factor 6 (KLF6) gene as a prospective candidate for milk production traits in dairy cattle. The expression of KLF6 in the livers of Holstein cows during the peak of lactation was significantly higher than that during the dry and early lactation periods. Notably, it plays an essential role in activating peroxisome proliferator-activated receptor α (PPARα) signaling pathways. The primary aim of this study was to further substantiate whether the KLF6 gene has significant genetic effects on milk traits in dairy cattle.ResultsThrough direct sequencing of PCR products with pooled DNA, we totally identified 12 single nucleotide polymorphisms (SNPs) within the KLF6 gene. The set of SNPs encompasses 7 located in 5′ flanking region, 2 located in exon 2 and 3 located in 3′ untranslated region (UTR). Of these, the g.44601035G > A is a missense mutation that resulting in the replacement of arginine (CGG) with glutamine (CAG), consequently leading to alterations in the secondary structure of the KLF6 protein, as predicted by SOPMA. The remaining 7 regulatory SNPs significantly impacted the transcriptional activity of KLF6 following mutation (P < 0.005), manifesting as changes in transcription factor binding sites. Additionally, 4 SNPs located in both the UTR and exons were predicted to influence the secondary structure of KLF6 mRNA using the RNAfold web server. Furthermore, we performed the genotype-phenotype association analysis using SAS 9.2 which found all the 12 SNPs were significantly correlated to milk yield, fat yield, fat percentage, protein yield and protein percentage within both the first and second lactations (P < 0.0001 ~ 0.0441). Also, with Haploview 4.2 software, we found the 12 SNPs linked closely and formed a haplotype block, which was strongly associated with five milk traits (P < 0.0001 ~ 0.0203).ConclusionsIn summary, our study represented the KLF6 gene has significant impacts on milk yield and composition traits in dairy cattle. Among the identified SNPs, 7 were implicated in modulating milk traits by impacting transcriptional activity, 4 by altering mRNA secondary structure, and 1 by affecting the protein secondary structure of KLF6. These findings provided valuable molecular insights for genomic selection program of dairy cattle.

Putatively adaptive missense variant at Adrenergic Receptor Alpha-1A associated with decreased risk of Chronic Mountain Sickness in Andean highlanders

Andean highlanders living at high altitude (&gt;4200 m) for hundreds of generations exhibit distinct cardiovascular and ventilatory traits in response to environmental stress. Adrenergic Receptor Alpha-1A ( ADRA1A), a gene coding for the Alpha 1-A adrenergic receptor, has been previously shown to be under positive natural selection in the Andeans. ADRA1A has also been demonstrated to have a protective effect against cardiovascular complications in other lowland populations. We hypothesized that putatively adaptive single nucleotide variants (SNVs) of strong functional impact can be found in ADRA1A and associated with cardiovascular traits in the Andean highlanders.We examined data from the Composite of Multiple Signals test in 40 high-coverage whole genomes of Andean highlanders living in Cerro de Pasco, Peru (~4340 m), specifically focusing on SNVs with strong predicted functional impact. Genetic variants of interests were then genotyped and tested for phenotypic associations in an extended Andean cohort (132 males, 46 females). The presence of Chronic Mountain Sickness (CMS), a maladaptive syndrome observed in Andean highlanders characterized by excessive erythrocytosis, and resting heart rate during hypoxia were prioritized for the regression analysis. The genotype-phenotype association analysis adapted a multivariate linear model with sex as a stratum and age as a covariate.A missense variant (rs1048101, A&gt;G) in the second exon of ADRA1A with strong predicted functional impact was tagged by the Composite of Multiple Signals analysis. In Andean males, the rs1048101 A allele was associated with 1.52-fold decreased risk for CMS (p&lt;0.018) and attenuated resting heart rate at room air (~10% inspired O 2 ) and at 21% inspired O 2 (room air: p&lt;0.019, 21% O 2 : p&lt;0.012). In women, the association between ADRA1A rs1048101 and CMS was not significant (p&gt;0.917). In addition, women with more copies of the putatively adaptive A allele demonstrated increased heart rate both at room air and 21% inspired O 2 (room air: p&lt;0.033, 21% O 2 : p&lt;0.044) in contrast to observations in men.In conclusion, the putatively adaptive missense variant rs1048101 at ADRA1A demonstrated a protective effect against Chronic Mountain Sickness and was further associated with decreased resting heart rate in Andean males. The lack of replicability of these associations in females suggests potential sex-dependent effect of ADRA1A. Further experiments functionalizing the variant and validating the associations in large genome cohorts may shed light on the mechanisms of the gene. This project is funded by R01HL145470. This is the full abstract presented at the American Physiology Summit 2023 meeting and is only available in HTML format. There are no additional versions or additional content available for this abstract. Physiology was not involved in the peer review process.

Association analysis between the TLR9 gene polymorphism rs352140 and type 1 diabetes.

To a great extent, genetic factors contribute to the susceptibility to type 1 diabetes (T1D) development, and by triggering immune imbalance, Toll-like receptor (TLR) 9 is involved in the development of T1D. However, there is a lack of evidence supporting a genetic association between polymorphisms in the TLR9 gene and T1D. In total, 1513 individuals, including T1D patients (n=738) and healthy control individuals (n=775), from the Han Chinese population were recruited for an association analysis of the rs352140 polymorphism of the TLR9 gene and T1D. rs352140 was genotyped by MassARRAY. The allele and genotype distributions of rs352140 in the T1D and healthy groups and those in different T1D subgroups were analyzed by the chi-squared test and binary logistic regression model. The chi-square test and Kruskal-Wallis H test were performed to explore the association between genotype and phenotype in T1D patients. The allele and genotype distributions of rs352140 were significantly different in T1D patients and healthy control individuals (p=0.019, p=0.035). Specifically, the T allele and TT genotype of rs352140 conferred a higher risk of T1D (OR=1.194, 95% CI=1.029-1.385, p=0.019, OR=1.535, 95% CI=1.108-2.126, p=0.010). The allele and genotype distributions of rs352140 were not significantly different between childhood-onset and adult-onset T1D and between T1D with a single islet autoantibody and T1D with multiple islet autoantibodies (p=0.603, p=0.743). rs352140 was associated with T1D susceptibility according to the recessive and additive models (p=0.015, p=0.019) but was not associated with T1D susceptibility in the dominant and overdominant models (p=0.117, p=0.928). Moreover, genotype-phenotype association analysis showed that the TT genotype of rs352140 was associated with higher fasting C-peptide levels (p=0.017). In the Han Chinese population, the TLR9 polymorphism rs352140 is associated with T1D and is a risk factor for susceptibility to T1D.

Uncovering the specificity and predictability of tryptophan metabolism in lactic acid bacteria with genomics and metabolomics.

Tryptophan is metabolized by microorganisms into various indole derivatives that have been proven to alleviate diseases and promote human health. Lactic acid bacteria (LAB) are a broad microbial concept, some of which have been developed as probiotics. However, the capacity of most LAB to metabolize tryptophan is unknown. In this study, the aim is to reveal the rule of tryptophan metabolism in LAB by multi-omics. The findings showed that LAB were rich in genes for tryptophan catabolism and that multiple genes were shared among LAB species. Although the number of their homologous sequences was different, they could still form the same metabolic enzyme system. The metabolomic analysis revealed that LAB were capable of producing a variety of metabolites. Strains belonging to the same species can produce the same metabolites and have similar yields. A few strains showed strain-specificity in the production of indole-3-lactic acid (ILA), indole-3-acetic acid, and 3-indolealdehyde (IAld). In the genotype-phenotype association analysis, the metabolites of LAB were found to be highly consistent with the outcomes of gene prediction, particularly ILA, indole-3-propionic acid, and indole-3-pyruvic acid. The overall prediction accuracy was more than 87% on average, which indicated the predictability of tryptophan metabolites of LAB. Additionally, genes influenced the concentration of metabolites. The levels of ILA and IAld were significantly correlated with the numbers of aromatic amino acid aminotransferase and amidase, respectively. The unique indolelactate dehydrogenase in Ligilactobacillus salivarius was the primary factor contributing to its large production of ILA. In summary, we demonstrated the gene distribution and production level of tryptophan metabolism in LAB and explored the correlation between genes and phenotypes. The predictability and specificity of the tryptophan metabolites in LAB were proven. These results provide a novel genomic method for the discovery of LAB with tryptophan metabolism potential and offer experimental data for probiotics that produce specific tryptophan metabolites.

Genome-wide sequencing identifies a thermal-tolerance related synonymous mutation in the mussel, Mytilisepta virgata

The roles of synonymous mutations for adapting to stressful thermal environments are of fundamental biological and ecological interests but poorly understood. To study whether synonymous mutations influence thermal adaptation at specific microhabitats, a genome-wide genotype-phenotype association analysis is carried out in the black mussels Mytilisepta virgata. A synonymous mutation of Ubiquitin-specific Peptidase 15 (MvUSP15) is significantly associated with the physiological upper thermal limit. The individuals carrying GG genotype (the G-type) at the mutant locus possess significantly lower heat tolerance compared to the individuals carrying GA and AA genotypes (the A-type). When heated to sublethal temperature, the G-type exhibit higher inter-individual variations in MvUSP15 expression, especially for the mussels on the sun-exposed microhabitats. Taken together, a synonymous mutation in MvUSP15 can affect the gene expression profile and interact with microhabitat heterogeneity to influence thermal resistance. This integrative study sheds light on the ecological importance of adaptive synonymous mutations as an underappreciated genetic buffer against heat stress and emphasizes the importance of integrative studies at a microhabitat scale for evaluating and predicting the impacts of climate change.

Genotype-phenotype association and biochemical analyses of glucose-6-phosphate dehydrogenase variants: Implications for the hemolytic risk of using 8-aminoquinolines for radical cure.

Background: Plasmodium vivax remains the malaria species posing a major threat to human health worldwide owing to its relapse mechanism. Currently, the only drugs of choice for radical cure are the 8-aminoquinolines (primaquine and tafenoquine), which are capable of killing hypnozoites and thus preventing P. vivax relapse. However, the therapeutic use of primaquine and tafenoquine is restricted because these drugs can cause hemolysis in individuals with glucose-6-phosphate dehydrogenase (G6PD) deficiency. This study aimed to assess and understand the hemolytic risk of using 8-aminoquinolines for radical treatment in a malaria endemic area of Thailand. Methods: The prevalence of G6PD deficiency was determined using a quantitative test in 1,125 individuals. Multiplexed high-resolution meltinging (HRM) assays were developed and applied to detect 12 G6PD mutations. Furthermore, biochemical and structural characterization of G6PD variants was carried out to understand the molecular basis of enzyme deficiency. Results: The prevalence of G6PD deficiency was 6.76% (76/1,125), as assessed by a phenotypic test. Multiplexed HRM assays revealed G6PD Mahidol in 15.04% (77/512) of males and 28.38% (174/613) of females, as well as G6PD Aures in one female. G6PD activity above the 30% cut-off was detected in those carrying G6PD Mahidol, even in hemizygous male individuals. Two variants, G6PD Murcia Oristano and G6PD Songklanagarind + Viangchan, were identified for the first time in Thailand. Biochemical characterization revealed that structural instability is the primary cause of enzyme deficiency in G6PD Aures, G6PD Murcia Oristano, G6PD Songklanagarind + Viangchan, and G6PD Chinese 4 + Viangchan, with double G6PD mutations causing more severe enzyme deficiency. Conclusion: In western Thailand, up to 22% of people may be ineligible for radical cure. Routine qualitative tests may be insufficient for G6PD testing, so quantitative tests should be implemented. G6PD genotyping should also be used to confirm G6PD status, especially in female individuals suspected of having G6PD deficiency. People with double G6PD mutations are more likely to have hemolysis than are those with single G6PD mutations because the double mutations significantly reduce the catalytic activity as well as the structural stability of the protein.

MetaRNN: differentiating rare pathogenic and rare benign missense SNVs and InDels using deep learning

Multiple computational approaches have been developed to improve our understanding of genetic variants. However, their ability to identify rare pathogenic variants from rare benign ones is still lacking. Using context annotations and deep learning methods, we present pathogenicity prediction models, MetaRNN and MetaRNN-indel, to help identify and prioritize rare nonsynonymous single nucleotide variants (nsSNVs) and non-frameshift insertion/deletions (nfINDELs). We use independent test sets to demonstrate that these new models outperform state-of-the-art competitors and achieve a more interpretable score distribution. Importantly, prediction scores from both models are comparable, enabling easy adoption of integrated genotype-phenotype association analysis methods. All pre-computed nsSNV scores are available at http://www.liulab.science/MetaRNN. The stand-alone program is also available at https://github.com/Chang-Li2019/MetaRNN.

Estimation of hereditary fructose intolerance prevalence in the Chinese population

BackgroundHereditary fructose intolerance (HFI) caused by aldolase B reduction or deficiency that results in fructose metabolism disorder. The disease prevalence in the Chinese population is unknown, which impedes the formulation of HFI screening and diagnosis strategies.Materials and methodsBy searching a local cohort (Chinese Children’s Rare Disease Genetic Testing Clinical Collaboration System, CCGT) and public databases (ClinVar and Human Gene Mutation Database) and reviewing HFI-related literature, we manually curated ALDOB pathogenic or likely pathogenic (P/LP) variants according to ACMG guidelines. Allele frequency (AF) information from the local database CCGT and the public databases HuaBiao and gnomAD for ALDOB P/LP variants was used to estimate and the HFI prevalence in the Chinese population and other populations by the Bayesian framework. We collected the genotype and clinical characteristics of HFI patients from the CCGT database and published literature to study genotype–phenotype relationships.ResultIn total, 81 variants of ALDOB were curated as P/LP. The estimated Chinese HFI prevalence was approximately 1/504,678, which was much lower than that for non-Finland European (1/23,147), Finnish in Finland (1/55,539), admixed American (1/132,801) and Ashkenazi Jewish (1/263,150) populations. By analyzing the genetic characteristics of ALDOB in the Chinese population, two variants (A338V, A338G) had significantly higher AFs in the Chinese population than in the non-Finland European population from gnomAD (all P values < 0.05). Five variants (A150P, A175D, N335K, R60*, R304Q) had significantly lower AFs (all P values < 0.1). The genotype–phenotype association analyses were based on 68 reported HFI patients from a literature review and the CCGT database. The results showed that patients carrying homozygous variant sites (especially A150P) were more likely to present nausea, and patients carrying two missense variant sites were more likely to present aversion to sweets and fruit (all P values < 0.05). Our research reveals that some gastrointestinal symptoms seem to be associated with certain genotypes.ConclusionThe prevalence of HFI in the Chinese population is extremely low, and there is no need to add HFI testing to the current newborn screening programs if medical costs are considered. A genetic testing strategy is suggested for early diagnosis of HFI.

Early marker of ocular neurodegeneration in children and adolescents with type 1 diabetes: the contributing role of polymorphisms in mir146a and mir128a genes

BackgroundEarly ocular neurodegenerative signs of diabetic neuropathy (DN) can be found in children and adolescents with type 1 diabetes (T1D). No data are available on the potential role of polymorphisms in miRNAs genes in predisposing T1D subjects to these signs.AimsTo determine whether MIR146A rs2910164 and MIR128A rs11888095 polymorphisms are associated with early retinal and corneal neurodegenerative changes in pediatric patients with T1D.MethodsA total of 140 T1D children/adolescents underwent spectral domain-optical coherence tomography (SD-OCT) and in vivo confocal microscopy (IVCM) with measurement of retinal and corneal nerve fiber parameters. Risk factors for diabetes complications (diabetes duration, blood pressure, HbA1c) were recorded. Genotyping of rs2910164 and rs1188095 SNPs and genotype–phenotype association analysis were performed.ResultsThe C allele of rs2910164 in MIR146A was associated with higher values of IVCM parameters and minimum rim width (MRW) of the peripapillary region of optic nerve head measured in the retina, whereas the T allele of rs1188095 in MIR128A was associated with a significant impairment of them. Multiple regression analysis showed that MIR146A and MIR128A polymorphisms were significantly associated with corneal nerve fiber length (beta = 0.225 and − 0.204, respectively) and other IVCM parameters, independently from age, diabetes duration, HbA1c and systolic blood pressure percentile. Similar results were found for MRW (beta = 0.213 and − 0.286, respectively).ConclusionsThese results provide new insight into the genetic predisposition to DN showing that two polymorphisms in MIR146A and MIR128A genes could significantly contribute to the development of early ocular preclinical signs of DN.

Genetic basis of growth reaction to drought stress differs in contrasting high-latitude treeline ecotones of a widespread conifer.

Climate change is increasing the frequency and intensity of drought events in many boreal forests. Trees are sessile organisms with a long generation time, which makes them vulnerable to fast climate change and hinders fast adaptations. Therefore, it is important to know how forests cope with drought stress and to explore the genetic basis of these reactions. We investigated three natural populations of white spruce (Picea glauca) in Alaska, located at one drought-limited and two cold-limited treelines with a paired plot design of one forest and one treeline plot. We obtained individual increment cores from 458 trees and climate data to assess dendrophenotypes, in particular the growth reaction to drought stress. To explore the genetic basis of these dendrophenotypes, we genotyped the individual trees at 3000 single nucleotide polymorphisms in candidate genes and performed genotype-phenotype association analysis using linear mixed models and Bayesian sparse linear mixed models. Growth reaction to drought stress differed in contrasting treeline populations. Therefore, the populations are likely to be unevenly affected by climate change. We identified 40 genes associated with dendrophenotypic traits that differed among the treeline populations. Most genes were identified in the drought-limited site, indicating comparatively strong selection pressure of drought-tolerant phenotypes. Contrasting patterns of drought-associated genes among sampled sites and in comparison to Canadian populations in a previous study suggest that drought adaptation acts on a local scale. Our results highlight genes that are associated with wood traits which in turn are critical for the establishment and persistence of future forests under climate change.

Applicability of the IrisPlex system for eye color prediction in an admixed population from Argentina.

Eye color prediction based on an individual's genetic information is of interest in the field of forensic genetics. In recent years, researchers have studied different genes and markers associated with this externally visible characteristic and have developed methods for its prediction. The IrisPlex represents a validated tool for homogeneous populations, though its applicability in populations of mixed ancestry is limited, mainly regarding the prediction of intermediate eye colors. With the aim of validating the applicability of this system in an admixed population from Argentina (n=302), we analyzed the six single nucleotide variants used in that multiplex for eye color and four additional SNPs, and evaluated its prediction ability. We also performed a genotype-phenotype association analysis. This system proved to be useful when dealing with the extreme ends of the eye color spectrum (blue and brown) but presented difficulties in determining the intermediate phenotypes (green), which were found in a large proportion of our population. We concluded that these genetic tools should be used with caution in admixed populations and that more studies are required in order to improve the prediction of intermediate phenotypes.

Human adaptation to arsenic in Bolivians living in the Andes.

Humans living in the Andes Mountains have been historically exposed to arsenic from natural sources, including drinking water. Enzymatic methylation of arsenic allows it to be excreted more efficiently by the human body. Adaptation to high-arsenic environments via enhanced methylation and excretion of arsenic was first reported in indigenous women in the Argentinean Andes, but whether adaptation to arsenic is a general phenomenon across native populations from the Andes Mountains remains unclear. Therefore, we evaluated whether adaptation to arsenic has occurred in the Bolivian Andes by studying indigenous groups who belong to the Aymara-Quechua and Uru ethnicities and have lived in the Bolivian Andes for generations. Our population genetics methods, including genome-wide selection scans based on linkage disequilibrium patterns and allele frequency differences, in combination with targeted and whole-genome sequencing and genotype-phenotype association analyses, detected signatures of positive selection near the gene encoding arsenite methyltransferase (AS3MT), the main arsenic methylating enzyme. This was among the strongest selection signals (top 0.5% signals via locus-specific branch length and extended haplotype homozygosity tests) at a genome-wide level in the Bolivian study groups. We found a large haplotype block of 676kb in the AS3MT region and identified candidate functional variants for further analysis. Moreover, our analyses revealed associations between AS3MT variants and the fraction of mono-methylated arsenic in urine and showed that the Bolivian study groups had the highest frequency of alleles associated with more efficient arsenic metabolism reported so far. Our data support the idea that arsenic exposure has been a driver for human adaptation to tolerate arsenic through more efficient arsenic detoxification in different Andean populations.

Upside down: dissecting impulsivity in attention-deficit hyperactivity disorder through genotype-phenotype association analyses

IntroductionBetter evaluation and understanding of the core symptoms have key importance both in clinical practice and the research of attention-deficit hyperactivity disorder (ADHD). One hallmark neurocognitive feature of ADHD is impaired inhibition, related to impulsivity. Given the high heritability of ADHD, the assessment of genetic background of impaired inhibition may contribute to our knowledge about the genetic background of the disorder.ObjectivesIn our study we investigated whether different forms of impulsivity (attentive, motor, and nonplanning) and polymorphisms in genes of the noradrenergic, serotonergic, and dopaminergic neurotransmission, i.e. dopamine transporter-1 (DAT1), cathecoloamin-O-metiltransferase (COMT), and serotonin receptor-1B (HTR1B genes show association.Methods208 aADHD patients diagnosed according to DSM-5 criteria from a clinical sample and 142 individuals from a population sample who screened positive for aADHD were included in the study. DNA samples were genotyped for the HTR-1B gene rs1321041 and the COMT gene rs4680 SNPs, moreover the DAT-1 VNTR polymorphism. Dimensional variables for impulsivity were compared between genotypes with the Generalized Linear Model procedure corrected for sex and age, using the PLINK 1.9 statistical software.ResultsThe 9 repeat polymorphism in DAT1 was associated with the severity of hyperactivity, moreover, all impulsivity factors. The A allele in COMT was associated with hyperactivity and better motor inhibition activity. In carriers of the G allele in HTR1B we detected significantly higher inattention scores and increased reaction time.ConclusionsOur results support the putative role of the investigated genetic polymorphisms in the etiology of impulsivity. Nevertheless, these polymorphisms demonstrate a heterogeneous associations.DisclosureNo significant relationships.

The Genomic and Phenotypic Landscape of Ichthyosis

Ichthyoses are clinically and genetically heterogeneous disorders characterized by scaly skin. Despite decades of investigation identifying pathogenic variants in more than 50 genes, clear genotype-phenotype associations have been difficult to establish. To expand the genotypic and phenotypic spectra of ichthyosis and delineate genotype-phenotype associations. This cohort study recruited an international group of individuals with ichthyosis and describes characteristic and distinguishing features of common genotypes, including genotype-phenotype associations, during a 10-year period from June 2011 to July 2021. Participants of all ages, races, and ethnicities were included and were enrolled worldwide from referral centers and patient advocacy groups. A questionnaire to assess clinical manifestations was completed by those with a genetic diagnosis. Genetic analysis of saliva or blood DNA, a phenotyping questionnaire, and standardized clinical photographs. Descriptive statistics, such as frequency counts, were used to describe the cases in the cohort. Fisher exact tests identified significant genotype-phenotype associations. Results were reported for 1000 unrelated individuals enrolled from around the world (mean [SD] age, 50.0 [34.0] years; 524 [52.4%] were female, 427 [42.7%] were male, and 49 [4.9%] were not classified); 75% were from the US, 12% from Latin America, 4% from Canada, 3% from Europe, 3% from Asia, 2% from Africa, 1% from the Middle East, and 1% from Australia and New Zealand. A total of 266 novel disease-associated variants in 32 genes were identified among 869 kindreds. Of these, 241 (91%) pathogenic variants were found through multiplex amplicon sequencing and 25 (9%) through exome sequencing. Among the 869 participants with a genetic diagnosis, 304 participants (35%) completed the phenotyping questionnaire. Analysis of clinical manifestations in these 304 individuals revealed that pruritus, hypohydrosis, skin pain, eye problems, skin odor, and skin infections were the most prevalent self-reported features. Genotype-phenotype association analysis revealed that the presence of a collodion membrane at birth (odds ratio [OR], 6.7; 95% CI, 3.0-16.7; P < .001), skin odor (OR, 2.8; 95% CI, 1.1-6.8; P = .02), hearing problems (OR, 2.9; 95% CI, 1.6-5.5; P < .001), eye problems (OR, 3.0; 95% CI, 1.5-6.0; P < .001), and alopecia (OR, 4.6; 95% CI, 2.4-9.0; P < .001) were significantly associated with TGM1 variants compared with other ichthyosis genotypes studied. Skin pain (OR, 6.8; 95% CI, 1.6-61.2; P = .002), odor (OR, 5.7; 95% CI, 2.0-19.7; P < .001), and infections (OR, 3.1; 95% CI, 1.4-7.7; P = .03) were significantly associated with KRT10 pathogenic variants compared with disease-associated variants in other genes that cause ichthyosis. Pathogenic variants were identified in 869 (86.9%) participants. Most of the remaining individuals had unique phenotypes, enabling further genetic discovery. This cohort study expands the genotypic and phenotypic spectrum of ichthyosis, establishing associations between clinical manifestations and genotypes. Collectively, the findings may help improve clinical assessment, assist with developing customized management plans, and improve clinical course prognostication.

A Genome-Wide Association Study Identifies Novel Risk Loci for Sarcopenia in a Taiwanese Population.

PurposeA genome-wide association study (GWAS) of sarcopenia unraveled the importance of genetic contribution to decline in muscle. The current study investigated sarcopenia-related single nucleotide polymorphisms (SNPs) in Asian older adults, and further constructed a genotype score that tests the combined effect of these SNPs on risk of sarcopenia.Patients and MethodsNinety-six subjects aged 60 or above were recruited from the database of annual geriatric health examination at Tri-Service General Hospital during 2020. Eligible criteria included: 1) not having severe comorbidities; 2) agreed to join the Taiwan Precision Medicine Initiative project; and 3) having sufficient information of required sarcopenic measurements. Genotype–phenotype association analysis was performed to find SNPs that were significantly associated with each of three sarcopenic indices (low muscle mass, muscle strength, and physical performance). Subsequently, these SNPs comprised a sarcopenia-related genotype score that summed up the number of SNPs carrying unfavorable allele(s).ResultsTwelve SNPs revealed suggestive genome-wide significance with the three sarcopenic indices, and eight of them revealed a relationship with more than one index. Low muscle strength was the item that had the most (eight) related SNPs. Among them, rs10282247 affects cholesterol binding and rs7022373 participates in cellular apoptosis. In addition, higher genotype score demonstrated higher risk of sarcopenia (≥4 points: OR=630.6; 2–3 points: OR=408, p-value<0.001).ConclusionSeveral newly discovered SNPs suggest that genetic contribution plays a part in the pathogenesis of sarcopenia. Further studies are warranted to verify the underlying mechanisms. Moreover, a genotype score provides an estimate of the combined effect of genetic association with sarcopenia, which may modestly improve clinical risk classification.

Genome-wide association study revealed novel candidate gene loci associated with soluble E-selectin levels in a Taiwanese population

Background and aimsIncrease soluble E-selectin (sE-selectin) levels are associated with various inflammation and cardiometabolic disorders. MethodsThis study aimed to investigate the genetic determinants of circulating sE-selectin levels by genome-wide association study (GWAS) in 4,525 Taiwan Biobank (TWB) participants and genotype-phenotype association analysis for sE-selectin level-determining alleles in over 80,000 TWB participants. ResultsBy GWAS, ABO, SELE, and FUT6 gene variants were identified as the determinants of sE-selectin levels, which reach genome-wide significance (maximum p = 3.25 × 10−271, 4.81 × 10−14, and 9.64 × 10−12, respectively). After further adjustment for the lead ABO rs2519093 genotypes, three novel gene loci, EVI5, FER and DMAC1, were associated with sE-selectin levels at p < 5 × 10−7. Three other previously reported gene loci, CELSR2, ST3GAL6-AS1, and HNF1A-AS1, also showed supportive evidence for the association with sE-selectin levels (maximum p < 0.0073). A multivariate analysis revealed age, body mass index, current smoking, hemoglobin A1C, hematocrit, leukocyte and platelet counts, serum alanine aminotransferase, triglycerides, and uric acid levels were independently associated with sE-selectin levels, in which the above ten gene loci contribute to 27.68% of the variance. For genotype-phenotype association analysis, a pleiotropic effect was demonstrated with genome-wide significant association between ABO gene variants and total and low-density-lipoprotein cholesterol levels, leukocyte counts and hematocrit. ConclusionsOur data provide novel insight into the regulation of sE-selectin levels. These results may open new avenues in understanding the critical role of E-selectin on the pathogenesis of inflammatory and cardiometabolic disorders.