Association study reveals a susceptibility locus with male pattern baldness in the Han Chinese population.

Introduction: Male pattern baldness (MPB), also known as androgenetic alopecia, represents the most prevalent form of progressive hair loss in humans. It is characterized by a distinctive pattern of hair loss progression from the scalp; however, its underlying mechanism remains elusive and is influenced by hereditary, immune, and environmental factors. Genome-wide association studies (GWASs) have uncovered numerous risk genes/loci among European individuals with MPB. However, the validation of these susceptibility genes/loci within Han Chinese men remains largely unexplored. The aim of this study was to investigate whether the 71 susceptibility loci identified in a recent GWAS among European men also confer risk for MPB in Chinese men. Methods: Forty-seven single nucleotide polymorphisms (SNPs) previously reported in GWASs of MPB were selected and genotyped in independent individuals comprising 499 Han Chinese cases and 1,489 controls using the Sequenom MassArray system. After stringent quality control measures, 25 SNPs were subjected to statistical analyses. Cochran-Armitage trend test was used to evaluate the association between SNPs and disease susceptibility. To address multiple tests, Bonferroni correction was conducted, setting the threshold for statistical significance at a p-value <2 × 10-3 (0.05/25). Results: The rs13405699 SNP located at 2q31.1 exhibited a significant association with MPB in Han Chinese men (p = 4.84 × 10-5, OR = 1.37, 95% CI: 1.18-1.59). Moreover, the difference in rs13405699 genotype distribution between MPB cases and controls was statistically significant (p = 7.00 × 10-5). Genotype-based association analysis suggested that the recessive model provided the best fit for the rs13405699 polymorphism. Conclusion: This study represents the first confirmation of the association between the rs13405699 SNP at 2q31.1 and MPB within the Han Chinese population, thereby enhancing our understanding of the genetic underpinnings of MPB.

MethylGenotyper: Accurate Estimation of SNP Genotypes and Genetic Relatedness from DNA Methylation Data.

Epigenome-wide association studies (EWAS) are susceptible to widespread confounding caused by population structure and genetic relatedness. Nevertheless, kinship estimation is challenging in EWAS without genotyping data. Here, we proposed MethylGenotyper, a method that for the first time enables accurate genotyping at thousands of single nucleotide polymorphisms (SNPs) directly from commercial DNA methylation microarrays. We modeled the intensities of methylation probes near SNPs with a mixture of three beta distributions corresponding to different genotypes and estimated parameters with an expectation-maximization algorithm. We conducted extensive simulations to demonstrate the performance of the method. When applying MethylGenotyper to the Infinium EPIC array data of 4662 Chinese samples, we obtained genotypes at 4319 SNPs with a concordance rate of 98.26%, enabling the identification of 255 pairs of close relatedness. Furthermore, we showed that MethylGenotyper allows for the estimation of both population structure and cryptic relatedness among 702 Australians of diverse ancestry. We also implemented MethylGenotyper in a publicly available R package (https://github.com/Yi-Jiang/MethylGenotyper) to facilitate future large-scale EWAS.

Identification of Key Genes Controlling Sugar and Organic Acid Accumulation in Wampee Fruit (Clausena lansium) via Genome Assembly and Genome-wide Association Analysis.

Wampee (Clausena lansium) is an economically significant subtropical fruit tree widely cultivated in Southern China. To provide high-quality genomic resources for C. lansium, we report a chromosome-level genome sequence for the "JinFeng" cultivar. The 297.1 Mb C. lansium genome contained nine chromosomes with a scaffold N50 of 29.2 Mb and encoded 23,468 protein-coding genes. Selective sweep analysis between sweet and sour C. lansium varieties and genome-wide association analysis identified 14 candidate genes putatively involved in sugar and acid accumulation. ClERF061, encoding an ethylene response factor, and ClSWEET7, encoding a Sugars Will Eventually be Exported Transporters (SWEET) family protein, were proposed as key regulators of the sweet and sour tastes of the wampee fruit. ClERF061 and ClSWEET7 overexpression in tomatoes increased the total sugar and acid content in fruits. ClSWEET7 promoter activation by ClERF061 was confirmed via Nicotiana benthamiana transient expression. Our study provides valuable genomic resources for C. lansium genetics and breeding.

Immune cell traits and causal relationships with cholecystitis: a mendelian randomization analysis.

Cholecystitis, characterized by inflammation of the gallbladder, is intricately linked to immune cells and the cytokines they produce. Despite this association, the specific contributions of immune cells to the onset and progression of cholecystitis remain to be fully understood. To delineate this relationship, we utilized the Mendelian randomization (MR) method to scrutinize the causal connections between 731 immune cell phenotypes and cholecystitis. By conducting MR analysis on 731 immune cell markers from public datasets, this study seeks to understand their potential impact on the risk of cholecystitis. It aims to elucidate the interactions between immune phenotypes and the disease, aiming to lay the groundwork for advancing precision medicine and developing effective treatment strategies for cholecystitis. Taking immune cell phenotypes as the exposure factor and cholecystitis as the outcome event, this study used single nucleotide polymorphisms (SNPs) closely associated with both immune cell phenotypes and cholecystitis as genetic instrumental variables. We conducted a two-sample MR analysis on genome-wide association studies (GWAS) data. Our research thoroughly examined 731 immune cell markers, to determine potential causal relationships with susceptibility to cholecystitis. Sensitivity analyses were performed to ensure the robustness of our findings, excluding the potential impacts of heterogeneity and pleiotropy. To avoid reverse causality, we conducted reverse MR analyses with cholecystitis as the exposure factor and immune cell phenotypes as the outcome event. Among the 731 immune phenotypes, our study identified 21 phenotypes with a causal relationship to cholecystitis (P < 0.05). Of these, eight immune phenotypes exhibited a protective effect against cholecystitis (odds ratio (OR) < 1), while the other 13 immune phenotypes were associated with an increased risk of developing cholecystitis (OR > 1). Additionally, employing the false discovery rate (FDR) method at a significance level of 0.2, no significant causal relationship was found between cholecystitis and immune phenotypes. Our research has uncovered a significant causal relationship between immune cell phenotypes and cholecystitis. This discovery not only enhances our understanding of the role of immune cells in the onset and progression of cholecystitis but also establishes a foundation for developing more precise biomarkers and targeted therapeutic strategies. It provides a scientific basis for more effective and personalized treatments in the future. These findings are expected to substantially improve the quality of life for patients with cholecystitis and mitigate the impact of the disease on patients and their families.

Association of a composite trait for anthropometrics, adiposity and energy expenditure with cardiometabolic diseases: An age-stratified cohort and genetic risk score analysis.

AIM: Various anthropometric measures capture distinct as well as overlapping characteristics of an individual's body composition. To characterize independent body composition measures, we aimed to reduce easily-obtainable individual measures reflecting adiposity, anthropometrics and energy expenditure into fewer independent constructs, and to assess their potential sex- and age-specific relation with cardiometabolic diseases. METHODS: Analyses were performed within European ancestry participants from UK Biobank (N = 418,963, mean age 58.0 years, 56% women). Principal components (PC) analyses were used for the dimension reduction of 11 measures of adiposity, anthropometrics and energy expenditure. PCs were studied in relation to incident type 2 diabetes mellitus (T2D) and coronary artery disease (CAD). Multivariable-adjusted Cox regression analyses, adjusted for confounding factors, were performed in all and stratified by age. Genome-wide association studies were performed in half of the cohort (N = 156,295) to identify genetic variants as instrumental variables. Genetic risk score analyses were performed in the other half of the cohort stratified by age of disease onset (N = 156,295). RESULTS: We identified two PCs, of which PC1 reflected lower overall adiposity (negatively correlated with all adiposity aspects) and PC2 reflected more central adiposity (mainly correlated with higher waist-hip ratio, but with lower total body fat) and increased height, collectively capturing 87.8% of the total variance. Similar to that observed in the multivariable-adjusted regression analyses, we found associations between the PC1 genetic risk score and lower risks of CAD and T2D [CAD cases <50 years, odds ratio: 0.91 (95% confidence interval 0.87, 0.94) per SD; T2D cases <50 years, odds ratio: 0.76 (0.72, 0.81)], which attenuated with higher age (p-values 8.13E-4 and 2.41E-6, respectively). No associations were found for PC2. CONCLUSIONS: The consistently observed weaker associations of the composite traits with cardiometabolic disease suggests the need for age-specific cardiometabolic disease prevention strategies.

Development of genomic and genetic resources facilitating molecular genetic studies on untapped Myanmar rice germplasms.

To counteract the growing population and climate changes, resilient varieties adapted to regional environmental changes are required. Landraces are valuable genetic resources for achieving this goal. Recent advances in sequencing technology have enabled national seed/gene banks to share genomic and genetic information from their collections including landraces, promoting the more efficient utilization of germplasms. In this study, we developed genomic and genetic resources for Myanmar rice germplasms. First, we assembled a diversity panel consisting of 250 accessions representing the genetic diversity of Myanmar indica varieties, including an elite lowland variety, Inn Ma Yebaw (IMY). Our population genetic analyses illustrated that the diversity panel represented Myanmar indica varieties well without any apparent population structure. Second, de novo genome assembly of IMY was conducted. The IMY assembly was constructed by anchoring 2888 contigs, which were assembled from 30× coverage of long reads, into 12 chromosomes. Although many gaps existed in the IMY genome assembly, our quality assessments indicated high completeness in the gene-coding regions, identical to other near-gap-free assemblies. Together with dense variant information, the diversity panel and IMY genome assembly will facilitate deeper genetic research and breeding projects that utilize the untapped Myanmar rice germplasms.

Phenome-wide profiling identifies genotype-phenotype associations in Phelan-McDermid syndrome using family-sourced data from an international registry.

BACKGROUND: Phelan-McDermid syndrome (PMS) is a rare neurodevelopmental disorder caused by 22q13 deletions that include the SHANK3 gene or pathogenic sequence variants in SHANK3. It is characterized by global developmental delay, intellectual disability, speech impairment, autism spectrum disorder, and hypotonia; other variable features include epilepsy, brain and renal malformations, and mild dysmorphic features. Here, we conducted genotype-phenotype correlation analyses using the PMS International Registry, a family-driven registry that compiles clinical data in the form of family-reported outcomes and family-sourced genetic test results. METHODS: Data from the registry were harmonized and integrated into the i2b2/tranSMART clinical and genomics data warehouse. We gathered information from 401 individuals with 22q13 deletions including SHANK3 (n = 350, ranging in size from 10 kb to 9.1 Mb) or pathogenic or likely pathogenic SHANK3 sequence variants (n = 51), and used regression models with deletion size as a potential predictor of clinical outcomes for 328 phenotypes. RESULTS: Our results showed that increased deletion size was significantly associated with delay in gross and fine motor acquisitions, a spectrum of conditions related to poor muscle tone, renal malformations, mild dysmorphic features (e.g., large fleshy hands, sacral dimple, dysplastic toenails, supernumerary teeth), lymphedema, congenital heart defects, and more frequent neuroimaging abnormalities and infections. These findings indicate that genes upstream of SHANK3 also contribute to some of the manifestations of PMS in individuals with larger deletions. We also showed that self-help skills, verbal ability and a range of psychiatric diagnoses (e.g., autism, ADHD, anxiety disorder) were more common among individuals with smaller deletions and SHANK3 variants. LIMITATIONS: Some participants were tested with targeted 22q microarrays rather than genome-wide arrays, and karyotypes were unavailable in many cases, thus precluding the analysis of the effect of other copy number variants or chromosomal rearrangements on the phenotype. CONCLUSIONS: This is the largest reported case series of individuals with PMS. Overall, we demonstrate the feasibility of using data from a family-sourced registry to conduct genotype-phenotype analyses in rare genetic disorders. We replicate and strengthen previous findings, and reveal novel associations between larger 22q13 deletions and congenital heart defects, neuroimaging abnormalities and recurrent infections.

A genome-wide Association study of the Count of Codeine prescriptions.

Opioid prescription records in existing electronic health record (EHR) databases are a potentially useful, high-fidelity data source for opioid use-related risk phenotyping in genetic analyses. Prescriptions for codeine derived from EHR records were used as targeting traits by screening 16 million patient-level medication records. Genome-wide association analyses were then conducted to identify genomic loci and candidate genes associated with different count patterns of codeine prescriptions. Both low- and high-prescription counts were captured by developing 8 types of phenotypes with selected ranges of prescription numbers to reflect potentially different levels of opioid risk severity. We identified one significant locus associated with low-count codeine prescriptions (1, 2 or 3 prescriptions), while up to 7 loci were identified for higher counts (≥ 4, ≥ 5, ≥6, or ≥ 7 prescriptions), with a strong overlap across different thresholds. We identified 9 significant genomic loci with all-count phenotype. Further, using the polygenic risk approach, we identified a significant correlation (Tau = 0.67, p = 0.01) between an externally derived polygenic risk score for opioid use disorder and numbers of codeine prescriptions. As a proof-of-concept study, our research provides a novel and generalizable phenotyping pipeline for the genomic study of opioid-related risk traits.

Evidence of survival bias in the association between APOE-Є4 and age at ischemic stroke onset.

Introduction: Large genome-wide association studies (GWASs) using case-control study designs have now identified tens of loci associated with ischemic stroke (IS). As a complement to these studies, we performed GWAS in a case-only design to identify loci influencing the age at onset (AAO) of ischemic stroke. Methods: Analyses were conducted in a discovery cohort of 10,857 ischemic stroke cases using a linear regression framework. We meta-analyzed all SNPs with p-value <1 x 10-5 in a sexcombined or sex-stratified analysis using summary data from two additional replication cohorts. Results: In the women-only meta-analysis, we detected significant evidence for the association of AAO with rs429358, an exonic variant in apolipoprotein E (APOE) that encodes for the APOE-Є4 allele. Each copy of the rs429358:T>C allele was associated with a 1.29-year earlier stroke AAO (meta p-value = 2.48 x 10-11). This APOE variant has previously been associated with increased mortality and ischemic stroke AAO. We hypothesized that the association with AAO may reflect a survival bias attributable to an age-related decrease in mortality among APOE-Є4 carriers and have no association to stroke AAO per se. A simulation study showed that a variant associated with overall mortality might indeed be detected with an AAO analysis. A variant with a 2-fold increase in mortality risk would lead to an observed effect of AAO that is comparable to what we found. Discussion: In conclusion, we detected a robust association of the APOE locus with stroke AAO and provided simulations to suggest that this association may be unrelated to ischemic stroke per se but related to a general survival bias.

Genetic association study of preterm birth and gestational age in a population-based case-control study in Peru.

BACKGROUND: Preterm birth (PTB) affects ∼15 million pregnancies worldwide. Genetic studies have identified several candidate loci for PTB, but results remain inconclusive and limited to European populations. Thus, we conducted a genome-wide association study (GWAS) of PTB and gestational age at delivery (GA) among 2,212 Peruvian women. METHODS: PTB cases delivered≥20 weeks' but < 37 weeks' gestation, while controls delivered at term (≥37 weeks but <42 weeks). Multivariable regressions were used to identify genetic markers for PTB and GA (∼6 million SNPs), adjusting for maternal age and the first two genetic principal components. In silico functional analysis was conducted among top signals detected with an arbitrary P < 1.0×10-5 . We sought to replicate genetic markers for PTB and GA identified in Europeans, and we developed a genetic risk score for GA based on European markers. RESULTS: Mean GA was 30 ± 4 weeks in PTB cases (N = 933) and 39 ± 1 in the controls (N = 1,279). No associatiosn were identified at genome-wide level. Nominal PTB variants were enriched for biological pathways associated with polyketide, progesterone, steroid hormones, and glycosyl metabolism. Nominal GA variants were enriched in intronic regions and cancer pathways. Variants in WNT4 associated with GA in Europeans were replicated in our study. A genetic risk score was associated with a 2-day longer GA (P = 0.002) in our sample. CONCLUSIONS: This study identified various signals suggestively associated with PTB and GA in pregnant Peruvian women. None of these variants overlapped with signals previously identified in Europeans.