Rare and low-frequency variants in human common diseases and other complex traits

Aberrant gene expression yet undiminished retinal ganglion cell genesis in iPSC-derived models of optic nerve hypoplasia

BACKGROUND

Optic nerve hypoplasia (ONH), the leading congenital cause of permanent blindness, is characterized by a retinal ganglion cell (RGC) deficit at birth. Multifactorial developmental events are hypothesized to underlie ONH and its frequently associated neurologic and endocrine abnormalities; however, environmental influences are unclear and genetic underpinnings are unexplored. This work investigates the genetic contribution to ONH RGC production and gene expression using patient induced pluripotent stem cell (iPSC)-derived retinal organoids (ROs).

MATERIALS AND METHODS

iPSCs produced from ONH patients and controls were differentiated to ROs. RGC genesis was assessed using immunofluorescence and flow cytometry. Flow-sorted BRN3+ cells were collected for RNA extraction for RNA-Sequencing. Differential gene expression was assessed using DESeq2 and edgeR. PANTHER was employed to identify statistically over-represented ontologies among the differentially expressed genes (DEGs). DEGs of high interest to ONH were distinguished by assessing function, mutational constraint, and prior identification in ONH, autism and neurodevelopmental disorder (NDD) studies.

RESULTS

RGC genesis and survival were similar in ONH and control ROs. Differential expression of 70 genes was identified in both DESeq2 and edgeR analyses, representing a ~ 4-fold higher percentage of DEGs than in randomized study participants. DEGs showed trends towards over-representation of validated NDD genes and ONH exome variant genes. Among the DEGs, RAPGEF4 and DMD had the greatest number of disease-relevant features.

CONCLUSIONS

ONH genetic background was not associated with impaired RGC genesis but was associated with DEGs exhibiting disease contribution potential. This constitutes some of the first evidence of a genetic contribution to ONH.

An Update on the Genetics of IgA Nephropathy

Immunoglobulin A (IgA) nephropathy (IgAN), the most common form of glomerulonephritis, is one of the leading causes of end-stage kidney disease (ESKD). It is widely believed that genetic factors play a significant role in the development of IgAN. Previous studies of IgAN have provided important insights to unravel the genetic architecture of IgAN and its potential pathogenic mechanisms. The genome-wide association studies (GWASs) together have identified over 30 risk loci for IgAN, which emphasizes the importance of IgA production and regulation in the pathogenesis of IgAN. Follow-up fine-mapping studies help to elucidate the candidate causal variant and the potential pathogenic molecular pathway and provide new potential therapeutic targets. With the rapid development of next-generation sequencing technologies, linkage studies based on whole-genome sequencing (WGS)/whole-exome sequencing (WES) also identify rare variants associated with IgAN, accounting for some of the missing heritability. The complexity of pathogenesis and phenotypic variability may be better understood by integrating genetics, epigenetics, and environment. We have compiled a review summarizing the latest advancements in genetic studies on IgAN. We similarly summarized relevant studies examining the involvement of epigenetics in the pathogenesis of IgAN. Future directions and challenges in this field are also proposed.

Balancing selection on genomic deletion polymorphisms in humans

A key question in biology is why genomic variation persists in a population for extended periods. Recent studies have identified examples of genomic deletions that have remained polymorphic in the human lineage for hundreds of millennia, ostensibly owing to balancing selection. Nevertheless, genome-wide investigation of ancient and possibly adaptive deletions remains an imperative exercise. Here, we demonstrate an excess of polymorphisms in present-day humans that predate the modern human-Neanderthal split (ancient polymorphisms), which cannot be explained solely by selectively neutral scenarios. We analyze the adaptive mechanisms that underlie this excess in deletion polymorphisms. Using a previously published measure of balancing selection, we show that this excess of ancient deletions is largely owing to balancing selection. Based on the absence of signatures of overdominance, we conclude that it is a rare mode of balancing selection among ancient deletions. Instead, more complex scenarios involving spatially and temporally variable selective pressures are likely more common mechanisms. Our results suggest that balancing selection resulted in ancient deletions harboring disproportionately more exonic variants with GWAS (genome-wide association studies) associations. We further found that ancient deletions are significantly enriched for traits related to metabolism and immunity. As a by-product of our analysis, we show that deletions are, on average, more deleterious than single nucleotide variants. We can now argue that not only is a vast majority of common variants shared among human populations, but a considerable portion of biologically relevant variants has been segregating among our ancestors for hundreds of thousands, if not millions, of years.

Genetic Variants and Functional Analyses of the ATG16L1 Gene Promoter in Acute Myocardial Infarction

Background

Acute myocardial infarction (AMI), a common complex disease caused by an interaction between genetic and environmental factors, is a serious type of coronary artery disease and is also a leading cause of death worldwide. Autophagy-related 16-like 1 (ATG16L1) is a key regulatory factor of autophagy and plays an important role in induced autophagy. In the cardiovascular system, autophagy is essential to preserve the homeostasis and function of the heart and blood vessels. No studies have hitherto examined the association between AMI and ATG16L1 gene promoter.

Methods

We conducted a case-control study, using polymerase chain reaction and sequencing techniques, dual luciferase reporter assay, and electrophoretic mobility shift assay, to analyze genetic and functional variation in the ATG16L1 gene promoter between AMI and controls. A variety of statistical analyses were used to analyze the allele and genotype frequencies and the relationship between single-nucleotide polymorphisms (SNPs) and AMI.

Results

In all, 10 SNPs and two DNA-sequence variants (DSVs) were identified in 688 subjects, and three ATG16L1 gene promoter mutations [g.233250693 T > C (rs185213911), g.233250946 G > A (rs568956599), g.233251133 C > G (rs1301744254)] that were identified in AMI patients significantly altered the transcriptional activity of ATG16L1 gene promoter in HEH2, HEK-293, and H9c2 cells (P < 0.05). Further electrophoretic mobility shift assays indicated that the SNPs affected the binding of transcription factors (P < 0.01).

Conclusion

ATG16L1 gene promoter mutations in AMI patients may affect the binding of transcription factors and change the transcriptional activity of the ATG16L1 gene, changing the level of autophagy and contributing to the occurrence and development of AMI as rare and low-frequency risk factors.

Emerging roles of rare and low-frequency genetic variants in type 1 diabetes mellitus

Type 1 diabetes mellitus (T1DM) is defined as an autoimmune disorder and has enormous complexity and heterogeneity. Although its precise pathogenic mechanisms are obscure, this disease is widely acknowledged to be precipitated by environmental factors in individuals with genetic susceptibility. To date, the known susceptibility loci, which have mostly been identified by genome-wide association studies, can explain 80%–85% of the heritability of T1DM. Researchers believe that at least a part of its missing genetic component is caused by undetected rare and low-frequency variants. Most common variants have only small to modest effect sizes, which increases the difficulty of dissecting their functions and restricts their potential clinical application. Intriguingly, many studies have indicated that rare and low-frequency variants have larger effect sizes and play more significant roles in susceptibility to common diseases, including T1DM, than common variants do. Therefore, better recognition of rare and low-frequency variants is beneficial for revealing the genetic architecture of T1DM and for providing new and potent therapeutic targets for this disease. Here, we will discuss existing challenges as well as the great significance of this field and review current knowledge of the contributions of rare and low-frequency variants to T1DM.

Molecular genetic study on GATA5 gene promoter in acute myocardial infarction

Background

Acute myocardial infarction (AMI) is a severe type of coronary artery disease, caused by coronary occlusion and followed by cardiac ischaemia. GATA binding protein 5 (GATA5) is an important member of GATA family and plays an important role in vascular inflammation, endothelial function, oxidative stress and cell metabolism. Previous studies have shown that the DNA sequence variants (DSVs) in GATA4 and GATA6 promoter can increase susceptibility to AMI. In this study, we explored the relationship between GATA5 promoter and AMI for the first time, hoping to provide a new genetic basis for understanding the pathogenesis of AMI.

Methods

GATA5 promoter was sequenced in 683 individuals (332 AMI patients and 351 controls). The transcriptional activity of the GATA5 promoter with or without DSVs in HEK-293 cells, H9c2 cells and primary neonatal rat cardiomyocytes were examined by Promega Dual-Luciferase® Reporter Assay system. Electrophoretic mobility shift assay (EMSA) was performed to explore whether the DSVs interfered with the binding of transcription factors (TFs).

Results

Nine mutations have been found in GATA5 promoter, eight of them evidently altered the transcriptional activity of the GATA5 promoter, five of them disrupted the binding of TFs (such as farnesoid X receptor). Furthermore, haplotype AT (across rs80197101 and rs77067995) is a dangerous haplotype of AMI. Genotype GA and allele A of rs80197101 and genotype CT and allele T of rs77067995 are the risk factors of AMI.

Conclusions

DSVs in GATA5 promoter can increase susceptibility to AMI. But the mechanism remains to be verified in vivo.

Targeted resequencing reveals rare variants enrichment in multiple sclerosis susceptibility genes

Although genome-wide association studies have identified a number of common variants associated with multiple sclerosis (MS) susceptibility, little is known about the relevance of rare variants. Here, we aimed to explore the role of rare variants in 14 MS risk genes (FCRL1, RGS1, TIMMDC1, HHEX, CXCR5, LTBR, TSFM, GALC, TRAF3, STAT3, TNFSF14, IFI30, CD40, and CYP24A1) by targeted resequencing in an Iberian population of 524 MS cases and 546 healthy controls. Four rare variants-enriched regions within CYP24A1, FCRL1, RGS1, and TRAF3 were identified as significantly associated with MS. Functional studies revealed significantly decreased regulator of G protein signaling 1 (RGS1) gene expression levels in peripheral blood mononuclear cells from MS patients with RGS1 rare variants compared to noncarriers, whereas no significant differences in gene expression were observed for CYP24A1, FCRL1, and TRAF3 between rare variants carriers and noncarriers. Immunophenotyping showed significant decrease in RGS1 expression in peripheral blood B lymphocytes from MS patients with RGS1 rare variants relative to noncarriers. Lastly, peripheral blood mononuclear cell from MS patients carrying RGS1 rare variants showed significantly lower induction of RGS1 gene expression by interferon-β compared to MS patients lacking RGS1 variants. The presence of rare variants in RGS1 reinforce the ideas of high genetic heterogeneity and a role of rare variants in MS pathogenesis.

Rare Variants of Putative Candidate Genes Associated With Sporadic Meniere's Disease in East Asian Population

Objectives: The cause of Meniere's disease (MD) is unclear but likely involves genetic and environmental factors. The aim of this study was to investigate the genetic basis underlying MD by screening putative candidate genes for MD.

Methods: Sixty-eight patients who met the diagnostic criteria for MD of the Barany Society were included. We performed targeted gene sequencing using next generation sequencing (NGS) panel composed of 45 MD-associated genes. We identified the rare variants causing non-synonymous amino acid changes, stop codons, and insertions/deletions in the coding regions, and excluded the common variants with minor allele frequency >0.01 in public databases. The pathogenicity of the identified variants was analyzed by various predictive tools and protein structural modeling.

Results: The average read depth for the targeted regions was 1446.3-fold, and 99.4% of the targeted regions were covered by 20 or more reads, achieving the high quality of the sequencing. After variant filtering, annotation, and interpretation, we identified a total of 15 rare heterozygous variants in 12 (17.6%) sporadic patients. Among them, four variants were detected in familial MD genes (DTNA, FAM136A, DPT), and the remaining 11 in MD-associated genes (PTPN22, NFKB1, CXCL10, TLR2, MTHFR, SLC44A2, NOS3, NOTCH2). Three patients had the variants in two or more genes. All variants were not detected in our healthy controls (n = 100). No significant differences were observed between patients with and without a genetic variant in terms of sex, mean age of onset, bilaterality, the type of MD, and hearing threshold at diagnosis.

Conclusions: Our study identified rare variants of putative candidate genes in some of MD patients. The genes were related to the formation of inner ear structures, the immune-associated process, or systemic hemostasis derangement, suggesting the multiple genetic predispositions in the development of MD.

Genetic and Functional Variants Analysis of the GATA6 Gene Promoter in Acute Myocardial Infarction

Background: Acute myocardial infarction (AMI) which is a specific type of coronary artery disease (CAD), is caused by the combination of genetic factors and acquired environment. Although some common genetic variations have been recorded to contribute to the development of CAD and AMI, more genetic factors and potential molecular mechanisms remain largely unknown. The GATA6 gene is expressed in the heart during embryogenesis and is also detected in vascular smooth muscle cells (VSMCs), different human primary endothelial cells (ECs), and vascular ECs in mice. To date, no studies have directly linked GATA6 gene with regulation of the CAD.

Methods: In this study, we used a case-control study to investigate and analyze the genetic variations and functional variations of the GATA6 gene promoter region in AMI patients and controls. A variety of statistical analysis methods were utilized to analyze the association of single nucleotide polymorphisms (SNPs) with AMI. Functional analysis of DNA sequence variants (DSVs) was performed using a dual luciferase reporter assay. In vitro, electrophoretic mobility shift assay (EMSA) was selected to examine DNA-protein interactions.

Results: A total of 705 subjects were enrolled in the study. Ten DSVs were found in AMI patients (n = 352) and controls (n = 353), including seven SNPs. One novel heterozygous DSV, (g.22168409 A > G), and two SNPs, [g.22168362 C > A(rs1416421760) and g.22168521 G > T(rs1445501474)], were reported in three AMI patients, which were not found in controls. The relevant statistical analysis, including allele and genotype frequencies between AMI patients and controls, five genetic models, linkage disequilibrium (LD) and haplotype analysis, and SNP–SNP interactions, suggested no statistical significance (P > 0.05). The transcriptional activity of GATA6 gene promoter was significantly increased by the DSV (g.22168409 A > G) and SNP [g.22168362 C > A(rs1416421760)]. The EMSA revealed that the DSV (g.22168409 A > G) and SNP [g.22168362 C > A(rs1416421760)] evidently influenced the binding of transcription factors.

Conclusion: In conclusion, the DSV (g.22168409 A > G) and SNP [g.22168362 C > A(rs1416421760)] may increase GATA6 levels in both HEK-293 and H9c2 cell lines by affecting the binding of transcription factors. Whether the two variants identified in the GATA6 gene promoter can promote the development and progression of human AMI by altering GATA6 levels still requires further studies to verify.

Functional variants of the ATG7 gene promoter in acute myocardial infarction

BACKGROUND

Coronary artery disease including acute myocardial infarction (AMI) is mainly caused by atherosclerosis, an inflammatory and metabolic disease. Autophagy has been demonstrated to play critical roles in lipid metabolism and inflammation. Altered autophagic activity has been reported in AMI patients. However, molecular basis for dysfunctional autophagy in AMI remains unexplained.

METHODS

In this study, the promoter of the ATG7 gene, encoding a core protein for autophagy, was genetically and functionally analyzed in large cohorts of AMI patients (n = 355) and ethnic-matched healthy controls (n = 363). Related molecular mechanisms were also explored.

RESULTS

A total of 19 DNA sequence variants (DSVs) including single-nucleotide polymorphisms (SNPs) were found in the ATG7 gene promoter. Two novel DSVs and five SNPs were only identified in AMI patients group. These DSVs and SNPs, except one SNP, significantly altered the transcriptional activity of the ATG7 gene promoter in both HEK-293 and H9c2 cells (p < 0.05). Further electrophoretic mobility shift assay revealed that the DSVs and SNPs evidently affected the binding of transcription factors.

CONCLUSIONS

ATG7 gene DSVs and SNPs identified in AMI patients may alter the transcriptional activity of the ATG7 gene promoter and change ATG7 level, contributing to the AMI development as a rare risk factor.

A novel SNP-set analytical method without distinguishing common variants or rare variants in genome-wide association study

Single nucleotide polymorphism (SNP)-set analysis in genome-wide association studies (GWASs) has become a hot topic. Most existing SNP-set analystic methods are designed and work well according to the different natures of common or rare variants and associated diseases. But the information that the disease associated variants are common or rare cannot be gained in advance. Therefore, in this research, we proposed a new and powerful weighted function method without distinguishing common or rare variants to select tagging SNP-set. We applied our selection method to sequence kernel association test (SKAT) and compared the power with some existing methods. The simulation results showed that our method has higher power not only than SKAT in un-weighted case, but also than SKAT in other weighted functions. Moreover, the power is improved significantly when the minor allele frequency (MAF) of causal SNP is relatively small.

Functional variants in the LC3B gene promoter in acute myocardial infarction

Acute myocardial infarction (AMI) is a common disease mainly caused by atherosclerosis, for which genetic causes remain largely unknown. Recently, low frequency and rare genetic variants have been proposed as risk factors. Autophagy has been involved in many cellular processes, such as lipid metabolism and inflammation, and implicated in human diseases, including cardiovascular diseases. In previous studies, we have reported reduced levels of LC3B, a core protein and a marker for autophagy, in AMI patients. In this study, the LC3B gene promoter was genetically and functionally analyzed in large cohorts of AMI patients (n = 383) and healthy controls (n = 390). A total of 25 DNA sequence variants (DSVs) including SNPs were found. Seven DSVs and three SNPs were only identified in AMI patients. All the DSVs and SNPs (except one) significantly decreased the transcriptional activity of the LC3B gene promoter in both HEK-293 and H9c2 cells (P < 0.05). Further electrophoretic mobility shift assay suggested that the DSVs affected the binding of transcription factors. In contrast, the DSVs and SNPs found only in controls or in both AMI patients and controls did not significantly affected LC3B gene promoter activity (P > 0.05). Therefore, our data suggested that the DSVs identified in AMI patients may change LC3B level by affecting the transcriptional activity of LC3B gene promoter, contributing to the AMI development. Upregulation of the LC3B gene expression may provide a novel and potential therapy for AMI patients.

Estimating the mutational load for cardiovascular diseases in Pakistani population

The deleterious genetic variants contributing to certain diseases may differ in terms of number and allele frequency from population to population depending on their evolutionary background. Here, we prioritize the deleterious variants from Pakistani population in manually curated gene list already reported to be associated with common, Mendelian, and congenital cardiovascular diseases (CVDs) using the genome/exome sequencing data of Pakistani individuals publically available in 1000 Genomes Project (PJL), and Exome Aggregation Consortium (ExAC) South Asia. By applying a set of tools such as Combined Annotation Dependent Depletion (CADD), ANNOVAR, and Variant Effect Predictor (VEP), we highlighted 561 potentially detrimental variants from PJL data, and 7374 variants from ExAC South Asian data. Likewise, filtration from ClinVar for CVDs revealed 03 pathogenic and 02 likely pathogenic variants from PJL and 112 pathogenic and 42 likely pathogenic variants from ExAC South Asians. The comparison of derived allele frequencies (DAF) revealed many of these prioritized variants having two fold and higher DAF in Pakistani individuals than in other populations. The highest number of deleterious variants contributing to common CVDs in descending order includes hypertension, atherosclerosis, heart failure, aneurysm, and coronary heart disease, and for Mendelian and congenital CVDs cardiomyopathies, cardiac arrhythmias, and atrioventricular septal defects.

Functional genetic variants within the SIRT2 gene promoter in acute myocardial infarction

Coronary artery disease (CAD), including acute myocardial infarction (AMI) is the complication of atherosclerosis. Recently, genome-wide association studies have identified a large number of CAD-related genetic variants. However, only 10% of CAD cases could be explained. Low frequent and rare genetic variants have been recently proposed to be main causes for CAD. SIRT2 is a member of sirtuin family, NAD(+)-dependent class III deacetylases. SIRT2 is involved in genomic stability, metabolism, inflammation, oxidative stress and autophagy, as well as in platelet function. Thus, we hypothesized that genetic variants in SIRT2 gene may contribute to AMI. In this study, SIRT2 gene promoter was analyzed in large cohorts of AMI patients (n = 375) and ethnic-matched controls (n = 377). Three novel heterozygous DSVs (g.38900888_91delTAAA, g.38900270A>G and g.38899853C>T) were identified in three AMI patients, but in none of controls. These DSVs significantly altered the transcriptional activity of the SIRT2 gene promoter (P<0.05) in both HEK-293 and H9c2 cells. Five novel heterozygous DSVS (g.38900562C>T, g.38900413A>C, g.38900030G>A, g.38899925A>C and g.38899852C>T) were only found in controls, which did not significantly affected SIRT2 gene promoter activity (P>0.05). In addition, four novel heterozygous DSVs and five SNPs were found in both AMI patients and control with similar frequencies (P>0.05), two SNPs of which were examined and did not affect SIRT2 gene promoter activity (P>0.05). Taken together, the DSVs identified in AMI patients may change SIRT2 level by affecting the transcriptional activity of SIRT2 gene promoter, contributing to the AMI development as a rare risk factor.

Otosclerosis Associated with a De Novo Mutation -832G > A in the TGFB1 Gene Promoter Causing a Decreased Expression Level

Otosclerosis (OTSC) is defined by abnormal bone remodeling in the otic capsule of middle ear which leads to conductive hearing loss. In our previous study, we have identified a de novo heterozygous mutation -832G > A in the promoter of TGFB1 in an otosclerosis patient. In the present study, we progressively screened this mutation in a cohort of 254 cases and 262 controls. The family members of the patient positive for -832G > A variation were also screened and found inheritance of this variation only to her daughter. Interestingly, this variation is associated with a decreased level of the TGFB1 transcript in the patient compared to her parents and controls. In silico analysis of this mutation predicted the altered binding of two transcription factors v-Myb and MZF1 in the mutated promoter sequence. Further, functional analysis of this mutation using in vitro luciferase and electrophoretic mobility shift assays revealed that this variation is associated with decreased gene expression. In conclusion, this study established the fact that TGFB1 mutation -832G > A altered the TGFB1 promoter activity, which could affect the susceptibility to otosclerosis development. Further, systemic analysis of TGFB1 gene sequence and expression analysis of this gene might reveal its precise role in the pathogenesis of otosclerosis.

Application of computational methods in genetic study of inflammatory bowel disease

Genetic factors play an important role in the etiology of inflammatory bowel disease (IBD). The launch of genome-wide association study (GWAS) represents a landmark in the genetic study of human complex disease. Concurrently, computational methods have undergone rapid development during the past a few years, which led to the identification of numerous disease susceptibility loci. IBD is one of the successful examples of GWAS and related analyses. A total of 163 genetic loci and multiple signaling pathways have been identified to be associated with IBD. Pleiotropic effects were found for many of these loci; and risk prediction models were built based on a broad spectrum of genetic variants. Important gene-gene, gene-environment interactions and key contributions of gut microbiome are being discovered. Here we will review the different types of analyses that have been applied to IBD genetic study, discuss the computational methods for each type of analysis, and summarize the discoveries made in IBD research with the application of these methods.

Assessing the Power of Exome Chips

Genotyping chips for rare and low-frequent variants have recently gained popularity with the introduction of exome chips, but the utility of these chips remains unclear. These chips were designed using exome sequencing data from mainly American-European individuals, enriched for a narrow set of common diseases. In addition, it is well-known that the statistical power of detecting associations with rare and low-frequent variants is much lower compared to studies exclusively involving common variants. We developed a simulation program adaptable to any exome chip design to empirically evaluate the power of the exome chips. We implemented the main properties of the Illumina HumanExome BeadChip array. The simulated data sets were used to assess the power of exome chip based studies for varying effect sizes and causal variant scenarios. We applied two widely-used statistical approaches for rare and low-frequency variants, which collapse the variants into genetic regions or genes. Under optimal conditions, we found that a sample size between 20,000 to 30,000 individuals were needed in order to detect modest effect sizes (0.5% < PAR > 1%) with 80% power. For small effect sizes (PAR <0.5%), 60,000–100,000 individuals were needed in the presence of non-causal variants. In conclusion, we found that at least tens of thousands of individuals are necessary to detect modest effects under optimal conditions. In addition, when using rare variant chips on cohorts or diseases they were not originally designed for, the identification of associated variants or genes will be even more challenging.

Analysis of Low Frequency Protein Truncating Stop-Codon Variants and Fasting Concentration of Growth Hormone

Background

The genetic background of Growth Hormone (GH) secretion is not well understood. Mutations giving rise to a stop codon have a high likelihood of affecting protein function.

Objectives

To analyze likely functional stop codon mutations that are associated with fasting plasma concentration of Growth Hormone.

Methods

We analyzed stop codon mutations in 5451 individuals in the Malmö Diet and Cancer study by genotyping the Illumina Exome Chip. To enrich for stop codon mutations with likely functional effects on protein function, we focused on those disrupting >80% of the predicted amino acid sequence, which were carried by ≥10 individuals. Such mutations were related to GH concentration, measured with a high sensitivity assay (hs-GH) and, if nominally significant, to GH related phenotypes, using linear regression analysis.

Results

Two stop codon mutations were associated with the fasting concentration of hs-GH. rs121909305 (NP_005370.1:p.R93*) [Minor Allele Frequency (MAF) = 0.8%] in the Myosin 1A gene (MYO1A) was associated with a 0.36 (95%CI, 0.04 to 0.54; p=0.02) increment of the standardized value of the natural logarithm of hs-GH per 1 minor allele and rs35699176 (NP_067040.1:p.Q100*) in the Zink Finger protein 77 gene (ZNF77) (MAF = 4.8%) was associated with a 0.12 (95%CI, 0.02 to 0.22; p = 0.02) increase of hs-GH. The mutated high hs-GH associated allele of MYO1A was related to lower BMI (β-coefficient, -0.22; p = 0.05), waist (β-coefficient, -0.22; p = 0.04), body fat percentage (β-coefficient, -0.23; p = 0.03) and with higher HDL (β-coefficient, 0.23; p = 0.04). The ZNF77 stop codon was associated with height (β-coefficient, 0.11; p = 0.02) but not with cardiometabolic risk factors.

Conclusion

We here suggest that a stop codon of MYO1A, disrupting 91% of the predicted amino acid sequence, is associated with higher hs-GH and GH-related traits suggesting that MYO1A is involved in GH metabolism and possibly body fat distribution. However, our results are preliminary and need replication in independent populations.