Association between genetically predicted leukocyte telomere length and non-scarring alopecia: A two-sample Mendelian randomization study

Background: The most commonly acknowledged non-scarring alopecia are androgenetic alopecia (AGA) and alopecia areata (AA). Previous studies have revealed various risk factors associated with alopecia. However, the relationship between leukocyte telomere length (LTL) and non-scarring alopecia remains unclear.

Methods: A two-sample Mendelian randomization (MR) analysis was performed to evaluate the causality between genetically predicted LTL and the risk of non-scarring alopecia. MR analyses were performed using the inverse variance-weighted (IVW) method and complemented with other MR methods.

Results: The summary statistics of the genome-wide association studies (GWAS) for AGA and AA were obtained from the FinnGen biobank, which included 119,185 and 211,428 individuals, respectively. A total of 126 single nucleotide polymorphisms (SNPs) with genome-wide significance were selected as the instrumental variables for LTL. The MR analyses suggested a causal relationship between LTL and AGA, and the risk of AGA increased by 3.19 times as the genetically predicted LTL was shortened by one standard deviation in log transformed form under the IVW method (OR = 4.19, 95% CI = 1.20–14.61, p = 0.024). The other MR methods also demonstrated a similar trend of the effect of LTL on AGA. There was no causal relationship between LTL and AA (p > 0.05). Sensitivity analyses further demonstrated that the current results were less likely to be affected by confounders and bias.

Conclusion: Our results suggested a potential causal relationship between LTL and AGA, and shortened LTL was associated with an increased risk of AGA.

Introduction

The two most commonly known non-scarring alopecia are androgenetic alopecia (AGA) and alopecia areata (AA). Although known as male pattern baldness (MPB), AGA is one of the most common types of chronic hair loss for both sexes, which affects at least 80% of men and half of women at the age of 70 years, with the incidence increasing with age (1, 2). AGA is characterized by an alteration of the hair growth cycle, which includes a reduced anagen phase duration and an increased telogen phase duration, resulting in gradual transformation of terminal hair into intermediate or vellus hair and, finally, balding. The most common forms of AGA consist of the representative receding of the frontal hairline for men together with the thinning of the hair between the frontal and the vertex scalp, but a maintained frontal hairline for women (2–4).

AA is a chronic hair follicle-specific autoimmune disorder characterized by patches of non-scarring alopecia, which influences approximately 2% of the population (5, 6). Nevertheless, the exact mechanisms remain under investigation, but a common feature is the collapse of the immune privilege of the hair follicles induced by immunological mechanisms (7). Telogen effluvium (TE) is also one of the important subtypes of non-scarring alopecia (8). However, due to the lack of existing summary data from the results of genome-wide association studies (GWAS), we excluded TE from this study.

Complicated genetic susceptibility, oxidative stress injury, environmental factors, lifestyle, and advanced age cooperatively contribute to the occurrence and severity of hair loss, among which, aging may play an important role in its pathological changes. Telomeres are the special nucleoprotein structures located at the ends of linear chromosomes that could protect their integrity (9, 10). The leukocyte telomere length (LTL) is shortened gradually with cellular divisions, and gene mutations or DNA damage sometimes occurs during this process, resulting in alterations in the cellular life span and senescence (11). Therefore, LTL has become widely recognized as a measure of an individual’s biological age and as a biomarker for certain age-related diseases (12). Previous observational studies have revealed some risk factors for alopecia; however, the relationship between aging and non-scarring alopecia has not been well studied. Therefore, we applied the Mendelian randomization (MR) method to determine whether there is a causal relationship between LTL and the two main subtypes of non-scarring alopecia.

MR is a method of instrumental variable (IV) analysis that uses single nucleotide polymorphisms (SNPs) from GWAS as instruments, which are representative of the exposure characteristics, to reveal causal relationships between complex characteristics, provided that the necessary assumptions have been satisfied (13). IVs are not correlated with random error (existence of deviation, such as age, sex, and weight, which need to be matched in observational studies) because SNP allocation follows the Mendelian inheritance rule, suggesting random allocation during the formation of a fertilized ovum. Therefore, based on SNP acting as the IV, the basic principle of MR research is comparable to that of randomized controlled trials (RCTs). This strategy is relatively convenient, cost-effective, and less likely to be confounded by covariables (14–17). Therefore, in this study, we used available summary statistics from open-access GWAS databases to perform a two-sample MR analysis and evaluated the causality between LTL and non-scarring alopecia. The results of the current study may partly explain the underlying mechanism of non-scarring alopecia and may offer support in the future development of prevention and intervention strategies.

Materials and methods

Mendelian randomization study design

MR is a statistical method used to evaluate causality without potential bias caused by confounders (18). Genetic variant is the most important and effective IV in an MR study. Eligible IVs should satisfy three important assumptions according to the MR theory: 1) the IVs are strongly correlated with the interested exposure factor; 2) the IVs follow the rule of random assignment and do not affect the outcome through confounders; and 3) the IVs do not exert influence on the outcome directly, but indirectly through the hypothesized causal pathway of exposure under investigation (19). The detailed flowchart of this MR study is shown in Figure 1A, while the MR assumptions are depicted in Figure 1B.

FIGURE 1

Figure 1 (A) Detailed flowchart of the current two-sample Mendelian randomization study. (B) Mendelian randomization assumptions. GWAS, genome-wide summary association study; LTL, leukocyte telomere length; AGA, androgenetic alopecia; AA, alopecia areata; SNP, single nucleotide polymorphism; MR, Mendelian randomization; MR-Egger, Mendelian randomization-Egger; MR-PRESSO, Mendelian randomization pleiotropy residual sum and outlier.

In this study, we selected SNPs closely associated with LTL as the IVs using summary data from open-access GWAS databases to evaluate the potential causal effect of the exposure (i.e., LTL) on the outcomes (i.e., AGA and AA). Moreover, the data for exposure and outcomes were obtained from different independent samples. Ethical approval was provided in the preliminary studies and was no longer required for this study.

Instrument variable selection

The summary statistics were obtained from the European Network for Genetic and Genomic Epidemiology (ENGAGE), which involved 472,174 European participants in the GWAS project for LTL. Quantitative polymerase chain reaction (PCR) analyses were performed to obtain LTL measurements for these participants from the UK Biobank (UKB) (9). The UKB recruited a prospective cohort of more than 500,000 volunteers aged 40–69 years (20). A total of 197 SNPs showed independent associations with LTL, as identified by the above-mentioned GWAS meta-analysis.

SNPs with genome-wide significance were selected as the IVs of telomere length for further study (p < 5 × 10⁸). Linkage disequilibrium (LD) was used to determine whether the SNPs were genetically linked or not, and the threshold was set as r² < 0.001, with a window size of 10,000 kb. SNPs associated with potential confounding factors of AGA and AA were removed from the final MR analyses according to the PhenoScanner V2 database (http://www.phenoscanner.medschl.cam.ac.uk/), an open-access database for genotype–phenotype associations. The palindromic variants were also excluded from this study, and the directions of the effects of SNPs on LTL and AGA, as well as AA, were harmonized. Eventually, a total of 126 eligible SNPs were included for the final MR analyses. Information on these SNPs is detailed in Supplementary Table S1.

Summary data sources for the subtypes of non-scarring alopecia

The aggregated data of the GWAS for AGA and AA were obtained from the open-access GWAS website (https://gwas.mrcieu.ac.uk/), with 119,185 and 211,428, respectively, European participants from the FinnGen biobank. The mean age at the first event was 47.7 years for AGA and 41.9 years for AA. The data for AGA and AA were obtained from participants from another independent consortium different from those in the GWAS for LTL. However, these SNPs were all derived from GWAS on European ancestry to minimize potential bias caused by population heterogeneity.

Statistical analyses of the MR study

For this study, the inverse variance-weighted (IVW), the Mendelian randomization-Egger (MR-Egger), the weighted median, the penalized weighted median, the simple median, and the maximum likelihood estimation methods were adopted to evaluate the risk association between shortened LTL and AGA, as well as AA (15, 21). Of the above methods, IVW is the most commonly used variant-specific causal estimations in a two-sample MR analysis, especially when all the enrolled IVs were of robust validity. Therefore, we adopted IVW as the primary and most efficient analysis method in this study (22). When heterogeneity was statistically significant in the research studies, a random effects IVW model was applied; otherwise, a fixed effects IVW model was adopted.

Other supplementary methods were also used for validation of the consistency and efficiency of the MR results. The MR-Egger regression is mainly used to evaluate horizontal pleiotropy and may provide causal estimation in the case of weaker IVs (23, 24). The penalized weighted median estimation is a new analytic method modified from the original weighted median method; thus, the IVs substantially contributing to the heterogeneity would be penalized, but the results would still be reliable provided that at least half of the IVs were valid (25). The maximum likelihood method was used to estimate the probability distribution by maximizing the likelihood function with low standard errors (26). Furthermore, the simple median method has also been applied in many other studies (21, 27). Moreover, scatter, forest, and funnel plots were drawn to visualize the results and the efficiency and stability of the MR study.

Sensitivity analyses

Sensitivity analyses were performed, including heterogeneity, pleiotropy, leave-one-out tests, and the Mendelian randomization pleiotropy residual sum and outlier (MR-PRESSO). The heterogeneity of the IVs applied in the IVW method was assessed using Cochran’s Q test. A p-value smaller than 0.05 was considered to suggest statistically significant heterogeneity. And the MR result should be carefully interpreted under this circumstances. Pleiotropy refers to a single locus possibly affecting multiple phenotypes. The MR-Egger intercept analysis was applied to evaluate directional pleiotropy, and the pleiotropy effect could be ignored when p > 0.05 (24). In addition, the MR-PRESSO analysis was performed to assess the influence of outliers. The leave-one-out test was performed by removing the SNPs one by one. When there is no great change to the remaining results from such a removal process (the result lines remained on the same side of the zero point), the causal relationship is stable and reliable.

All statistical analyses were performed using the “Two-Sample MR” and “MRPRESSO” packages in R software (version 4.1.1).

Results

Genetically predicted shorter LTL is associated with increased risk of AGA instead of AA

In the two-sample MR analysis, after removal of the SNPs associated with confounding factors and outcomes, along with removal of the SNPs of the incompatible alleles or those that are palindromic with intermediate allele frequencies, a total of 126 LTL-related SNPs that met the three main assumptions of the MR study were retained for further MR analysis. The scatter plots revealed the effects of the above 126 SNPs on LTL and AGA, as well as on AA (Figures 2A, B). The funnel plots showed the distribution of the effect of a single SNP (Figures 3A, B), while the forest plots visualized the effect of an individual SNP on the estimation of the outcomes (Supplementary Figures S1A, B).

FIGURE 2

Figure 2 Scatter plots for the effects of single-nucleotide polymorphisms (SNPs) on telomere length and non-scarring alopecia. The x-axis represents the effects of each genetic variant on telomere length, and the y-axis represents the effects of each genetic variant on androgenetic alopecia (A) or alopecia areata (B).

FIGURE 3

Figure 3 Funnel plots of leukocyte telomere length genetic variants and non-scarring alopecia. (A, B). Funnel plots for androgenetic alopecia (A) and alopecia areata (B).

The potential causal relationship between LTL and AGA was evaluated using the two-sample MR method. Table 1 shows that there was no clear evidence of a heterogeneity between the LTL-related SNPs, as the MR-Egger test had a Cochran’s Q statistic of 112.0, with a p-value greater than 0.05 (MR-Egger method, p = 0.771) for AGA. Therefore, the IVW method with a fixed effects model was adopted primarily for causal estimation. The MR analyses suggested a causal relationship between LTL and AGA, and the risk of AGA increased by 3.19 times as the genetically predicted LTL was shortened by one standard deviation in log transformed form under the fixed effects IVW method (OR = 4.19, 95% CI = 1.20–14.61, p = 0.024). This causal association was confirmed by maximum likelihood estimation (OR = 4.25, 95% CI = 1.21–14.90, p = 0.023) and the simple median method (OR = 7.14, 95% CI = 1.04–49.29, p = 0.046) (Figure 4). Moreover, the results from the MR-Egger, the weighted median, and the penalized weighted median estimations demonstrated a similar trend of the effect of LTL on AGA, although there was no statistical significance. As mentioned above, the IVW method was adopted as the primary and most efficient method for analysis in this study; therefore, the conclusion that the genetically predicted shorter LTL was associated with an increased risk of AGA was drawn.

TABLE 1

Table 1 Heterogeneity analysis of the leukocyte telomere length-related genetic variants in the androgenetic alopecia and alopecia areata GWAS datasets.

FIGURE 4

Figure 4 Forest plot of the association between genetically predicted telomere length and androgenetic alopecia. OR refers to the change of alopecia risk associated with a one standard deviation (SD) decrease in telomere length. IVW, inverse variance weighted; OR, odds ratio; CI, confidence interval.

We also assessed the relationship between LTL and AA using the same MR methods. There was no significant heterogeneity detected, with a Q statistic of 113.15 (p = 0.747) using the MR-Egger method (Table 1). However, the results of the fixed effects IVW, the MR-Egger, the weighted median, the penalized weighted median, and the simple median methods, as well as the maximum likelihood estimation, suggested that there was no clear causal relationship between genetically predicted LTL and risk of AA, with all p-values greater than 0.05 (Figure 5).

FIGURE 5

Figure 5 Forest plot of the association between genetically predicted telomere length and alopecia areata.

Sensitivity analysis

No clear evidence for significant statistical heterogeneity was found among the LTL-related SNPs on the effects of AGA and AA, as mentioned previously. In addition, no evidence for directional pleiotropy was found according to the MR-Egger intercept analysis (intercept = 0.028, SE = 0.032, p = 0.37), indicating that the exposure–outcome relationship was unlikely to be affected by potential confounding factors through other pathways. The MR-PRESSO analysis suggested no directional pleiotropy and no outliers for the IVs (p-global test = 0.795). The leave-one-out analysis yielded comparable results to that of the primary MR studies and indicated that no single SNP has had significant effects on the results, suggesting the robustness and reliability of the MR studies (Supplementary Figures S2A, B).

Discussion

Hair loss has a significant negative impact on body image, social perceptions, and the quality of life of patients and may result in psychological disorders (7, 28, 29). New research works have revealed that complicated factors may contribute to the etiopathogenesis of hair loss, including genetic susceptibility, advanced age, autoimmunological changes, oxidative stress injury, skin microbiome alterations, and epigenetic factors, among other factors (30). To our knowledge AA has affected approximately 2% of the general population, including all sex- and age-based groups (29), while the incidence of AGA has revealed an increasing trend with advanced age. Therefore, this study was performed to assess whether there is an association between non-scarring alopecia, special types of hair loss, and LTL, a well-recognized marker for biological age, rather than chronological age, which may be involved in the pathogenesis of some age-related disorders. In this study, we applied the two-sample MR method to evaluate a potential causal association, which may minimize the bias resulting from confounding factors and reverse causality. The results suggested that, for the European population, the genetically predicted shortened LTL was correlated with a higher risk of AGA, but not of AA.

AA is an autoimmune disorder characterized by non-scarring hair loss on the scalp or any hair-bearing skin. To date, the specific pathogenesis of AA remains unknown, although it has been mainly recognized to be caused by immune disorders and genetic factors. Research studies have revealed that melanocytes may be the trigger point that leads to the hair follicles being attacked by the immune system through oxidative stress injuries or apoptosis (31). On the other hand, different immune cells cooperatively contribute to the progression of AA (32, 33). The infiltration of T helper cells and cytolytic CD8⁺ cells may lead to the disruption of the hair growth cycle. T cells, natural killer (NK) cells, and plasmacytoid dendritic cells surround the lower part of the hair bulb during the anagen and growth phases, and their immunological activities may cause the collapse of the immune privilege of hair follicles and therefore lead to alopecia (7, 34, 35). This process is different from the natural aging process of hair follicles, which may partly account for the result that no association was found between LTL and AA risk.

AGA, a progressive thinning of scalp hair with specific patterns characterized by the hairline receding at the temples (the widow’s peak) and the hair eventually left on the sides and the back of head (the Hippocratic wreath), is now commonly regarded to be correlated with polygenic susceptibility, increased sensitivity to dihydrotestosterone (DHT) of the hair follicle, and chronic scalp inflammation (36). DHT affects hair growth and leads to hair loss through binding to the receptors located in the oil glands of hair follicles, therefore causing the affected follicles to shrink and shortening their anagen cycles, eventually inhibiting the ability of the hair to grow (37, 38). The conversion of testosterone to DHT is catalyzed by 5-α-reductase, an enzyme stored in the oil glands located in hair follicles. Based on the above theories, one recent important treatment strategy emphasizes the suppression of DHT and 5-α-reductase (39, 40); therefore, injuries of the hair follicles may be alleviated and hair loss may be reduced. Previous epidemiological studies suggested that serum DHT and the estradiol (E2) level are correlated with LTL independently of age in men (41), while the testosterone level is not associated with LTL (42).

Observational studies also suggested that, compared with AGA in men, the onset of female androgenetic alopecia (FAGA) occurs in those of more advanced ages and is particularly more common among those undergoing menopause, and the role of androgens remains uncertain in FAGA (1). However, the incidence increasing with age has been observed in both sexes. Advanced age is accompanied by an accumulation of health conditions and the shortening of chromosomal telomere length, signifying an individual’s biological aging. In this study, the results revealed a causal relationship between shortened LTL and a higher AGA risk, which is in accordance with previous observational results. However, aging may be just one of the underlying mechanisms of alopecia, and more studies are still needed to investigate the mechanisms of its onset. One significant strength of the MR method is that the bias resulting from the effects of confounders are removed from the final results. In this study, no SNPs correlated with the DHT concentration were found; therefore, the relationship between telomere length and the risk of AGA was less likely to be affected by the influence of DHT. We also excluded the SNPs related to baldness, hyperthyroidism, and hypothyroidism in order to avoid confounding effects.

As mentioned in other published MR research studies (17, 21, 43), causal inference is considered to be significant if the following criteria were satisfied: 1) the p-value using the IVW method is statistically significant (p < 0.05); 2) the directions of estimates by the IVW and other methods are all the same (in the current study on AGA, all ORs were greater than 1, prompting the conclusion that the genetically predicted shortened LTL was associated with an increased risk of AGA); and 3) in the sensitivity analyses, the MR-Egger intercept test and the MR-PRESSO global test are not significant (p > 0.05). Sensitivity analysis is one of the most important components of MR research, which is designed to detect whether the results are influenced by confounders and bias. In this study, the above criteria were all satisfied; therefore, the conclusion that there might be a causal relationship between LTL and AGA was drawn.

The major strengths of the current MR study are as follows. Firstly, the eligible IVs were obtained from the largest GWAS study from the European Network for Genetic and Genomic Epidemiology (ENGAGE), which satisfied the three main important principles for MR investigation, allowing for reasonable causal estimations. Secondly, the results were valid and consistent in the different MR methods and remained stable in the sensitivity analyses, suggesting the validity and consistency of the causal results. However, this study has several limitations. Firstly, all enrolled individuals in the preliminary GWAS studies were of European ancestry; thus, whether or not the current results could be extended to other populations remains unknown, which needs further investigation. Secondly, we only investigated the relationship between telomere length and the risk of AGA and AA, and due to the lack of summary data for other types of alopecia, the stratification statistics for the degree of AGA, and the stratification data based on sex, further studies were limited. Thirdly, there are also other factors, such as the limitations of the already known phenotype characteristics and LTL determined by genetic factors and that may be affected by environmental factors and potential epigenetic modifications, that have not been thoroughly investigated yet and therefore might contradict the current causal assumptions in future practice. However, these potential effects should not be ignored, and the current conclusions were drawn on the basis of the existing statistics.

Our study suggests that genetically predicted LTL may have a causal effect on the onset of AGA, but not of AA, and that shortened LTL may be one of the risk factors for AGA. The results may offer new insights into the mechanisms of AGA. However, future high-quality randomized controlled studies based on large sample sizes are still needed to further verify this association.

Conclusion

This MR study suggested that the shortening of the LTL has a negative effect on AGA, while LTL has no association with the risk of AA. Further exploration into the role of LTL on AGA will promote an understanding of the pathogenesis of AGA. A new insight put forward is that preventing the loss of the leukocyte telomere may be a novel target to reduce the risk of AGA.

Data availability statement

The original contributions presented in the study are partly included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.

Author contributions

YCL and LY put forward the conception of this study and wrote the manuscript. YCL, STY and MJL performed the statistical analyses. ZJZ and MYL helped revise this manuscript. XRW and MQL collected the related data. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the Science and Technology Innovation Project of Guangdong Province (no. 2018KJYZ005), the Natural Science Foundation of Guangdong Province (no. 2020A151501107), the Natural Science Foundation of Tibet Autonomous Region (no. XZ2017ZR-ZY021), and the Guangdong Province Key Field R&D Programme Project (no. 2020B1111150001).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2022.1072573/full#supplementary-material

Supplementary Figure 1 | Forest plot of leukocyte telomere length and risk of androgenetic alopecia (A) and alopecia areata (B).

Supplementary Figure 2 | Leave-one-out analyses for the effect of leukocyte telomere length and non-scarring alopecia. (A) Analysis of leukocyte telomere length and androgenetic alopecia; (B) Analysis of leukocyte telomere length and alopecia areata.

Abbreviations

AGA, androgenetic alopecia; AA, alopecia areata; LTL, leukocyte telomere length; MR, Mendelian randomization; IVW, inverse variance weighted; GWAS, genome-wide association study; SNPs, single nucleotide polymorphisms; MPB, male pattern baldness; TE, telogen effluvium; IVs, instrumental variants; ENGAGE, European Network for Genetic and Genomic Epidemiology; MR-Egger, Mendelian randomization-Egger; LD, linkage disequilibrium; MR-PRESSO, Mendelian randomization pleiotropy residual sum and outlier; DHT, dihydrotestosterone; FAGA, female androgenetic alopecia.

References

1. Kelly Y, Blanco A, Tosti A. Androgenetic alopecia: An update of treatment options. Drugs (2016) 76(14):1349–64. doi: 10.1007/s40265-016-0629-5

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Nestor MS, Ablon G, Gade A, Han H, Fischer DL. Treatment options for androgenetic alopecia: Efficacy, side effects, compliance, financial considerations, and ethics. J Cosmet Dermatol (2021) 20(12):3759–81. doi: 10.1111/jocd.14537

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Starace M, Orlando G, Alessandrini A, Piraccini BM. Female androgenetic alopecia: An update on diagnosis and management. Am J Clin Dermatol (2020) 21(1):69–84. doi: 10.1007/s40257-019-00479-x

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Kitagawa T, Matsuda K, Inui S, Takenaka H, Katoh N, Itami S, et al. Keratinocyte growth inhibition through the modification of wnt signaling by androgen in balding dermal papilla cells. J Clin Endocrinol Metab (2009) 94(4):1288–94. doi: 10.1210/jc.2008-1053

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Zhou C, Li X, Wang C, Zhang J. Alopecia areata: an update on etiopathogenesis, diagnosis, and management. Clin Rev Allergy Immunol (2021) 61(3):403–23. doi: 10.1007/s12016-021-08883-0

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Lee HH, Gwillim E, Patel KR, Hua T, Rastogi S, Ibler E, et al. Epidemiology of alopecia areata, ophiasis, totalis, and universalis: A systematic review and meta-analysis. J Am Acad Dermatol (2020) 82(3):675–82. doi: 10.1016/j.jaad.2019.08.032

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Simakou T, Butcher JP, Reid S, Henriquez FL. Alopecia areata: A multifactorial autoimmune condition. J Autoimmun (2019) 98:74–85. doi: 10.1016/j.jaut.2018.12.001

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Alessandrini A, Bruni F, Piraccini BM, Starace M. Common causes of hair loss - clinical manifestations, trichoscopy and therapy. J Eur Acad Dermatol Venereol (2021) 35(3):629–40. doi: 10.1111/jdv.17079

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Codd V, Wang Q, Allara E, Musicha C, Kaptoge S, Stoma S, et al. Polygenic basis and biomedical consequences of telomere length variation. Nat Genet (2021) 53(10):1425–33. doi: 10.1038/s41588-021-00944-6

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Li C, Stoma S, Lotta LA, Warner S, Albrecht E, Allione A, et al. Genome-wide association analysis in humans links nucleotide metabolism to leukocyte telomere length. Am J Hum Genet (2020) 106(3):389–404. doi: 10.1016/j.ajhg.2020.02.006

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Blasco MA. Telomeres and human disease: ageing, cancer and beyond. Nat Rev Genet (2005) 6(8):611–22. doi: 10.1038/nrg1656

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Andreu-Sanchez S, Aubert G, Ripoll-Cladellas A, Henkelman S, Zhernakova DV, Sinha T, et al. Genetic, parental and lifestyle factors influence telomere length. Commun Biol (2022) 5(1):565. doi: 10.1038/s42003-022-03521-7

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Boehm FJ ZX. Statistical methods for mendelian randomization in genome-wide association studies: A review. Comput Struct Biotechnol J (2022) 20:2338–51. doi: 10.1016/j.csbj.2022.05.015

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Hu X, Zhao J, Lin Z, Wang Y, Peng H, Zhao H, et al. Mendelian randomization for causal inference accounting for pleiotropy and sample structure using genome-wide summary statistics. Proc Natl Acad Sci U S A. (2022) 119(28):e2106858119. doi: 10.1073/pnas.2106858119

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Burgess S, Butterworth A, Thompson SG. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol (2013) 37(7):658–65. doi: 10.1002/gepi.21758

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Rothman KJ, Greenland S. Causation and causal inference in epidemiology. Am J Public Health (2005) 95 Suppl 1:S144–50. doi: 10.2105/AJPH.2004.059204

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Zhang Y LZ, Choudhury T, Cornelis MC, Liu W. Habitual coffee intake and risk for nonalcoholic fatty liver disease: a two-sample mendelian randomization study. Eur J Nutr (2021) 60(4):1761–7. doi: 10.1007/s00394-020-02369-z

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Davey Smith G, Hemani G. Mendelian randomization: genetic anchors for causal inference in epidemiological studies. Hum Mol Genet (2014) 23(R1):R89–98. doi: 10.1093/hmg/ddu328

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Didelez V, Sheehan N. Mendelian randomization as an instrumental variable approach to causal inference. Stat Methods Med Res (2007) 16(4):309–30. doi: 10.1177/0962280206077743

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Allen NE SC, Peakman T, Collins R, Biobank UK. UK Biobank data: come and get it. Sci Transl Med (2014) 6(224):224ed4. doi: 10.1126/scitranslmed.3008601

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Xia K, Zhang L, Zhang G, Wang Y, Huang T, Fan D. Leukocyte telomere length and amyotrophic lateral sclerosis: a mendelian randomization study. Orphanet J Rare Dis (2021) 16(1):508. doi: 10.1186/s13023-021-02135-2

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Burgess S, Scott RA, Timpson NJ, Davey Smith G, Thompson SG, Consortium E-I. Using published data in mendelian randomization: a blueprint for efficient identification of causal risk factors. Eur J Epidemiol (2015) 30(7):543–52. doi: 10.1007/s10654-015-0011-z

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Bowden J, Del Greco MF, Minelli C, Davey Smith G, Sheehan NA, Thompson JR. Assessing the suitability of summary data for two-sample mendelian randomization analyses using MR-egger regression: the role of the I2 statistic. Int J Epidemiol (2016) 45(6):1961–74. doi: 10.1093/ije/dyw220

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Burgess S, Thompson SG. Interpreting findings from mendelian randomization using the MR-egger method. Eur J Epidemiol (2017) 32(5):377–89. doi: 10.1007/s10654-017-0255-x

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent estimation in mendelian randomization with some invalid instruments using a weighted median estimator. Genet Epidemiol (2016) 40(4):304–14. doi: 10.1002/gepi.21965

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Milligan BG. Maximum-likelihood estimation of relatedness. Genetics (2003) 163(3):1153–67. doi: 10.1093/genetics/163.3.1153

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Zhang N, Wang Y, Chen Z, Liu D, Tse G, Korantzopoulos P, et al. Circulating vitamin d concentrations and risk of atrial fibrillation: A mendelian randomization study using non-deficient range summary statistics. Front Nutr (2022) 9:842392. doi: 10.3389/fnut.2022.842392

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Saini K, Mysore V. Role of vitamin d in hair loss: A short review. J Cosmet Dermatol (2021) 20(11):3407–14. doi: 10.1111/jocd.14421

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Jadeja SD, Tobin DJ. Autoantigen discovery in the hair loss disorder, alopecia areata: Implication of post-translational modifications. Front Immunol (2022) 13:890027. doi: 10.3389/fimmu.2022.890027

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Strazzulla LC, Wang EHC, Avila L, Lo Sicco K, Brinster N, Christiano AM, et al. Alopecia areata: Disease characteristics, clinical evaluation, and new perspectives on pathogenesis. J Am Acad Dermatol (2018) 78(1):1–12. doi: 10.1016/j.jaad.2017.04.1141

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Xie B, Sun J, Song X. Hair follicle melanocytes initiate autoimmunity in alopecia areata: a trigger point. Clin Rev Allergy Immunol (2022) 63:(3) 417–30. doi: 10.1007/s12016-022-08954-w

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Gilhar A, Laufer-Britva R, Keren A, Paus R. Frontiers in alopecia areata pathobiology research. J Allergy Clin Immunol (2019) 144(6):1478–89. doi: 10.1016/j.jaci.2019.08.035

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Trueb RM, Dias M. Alopecia areata: a comprehensive review of pathogenesis and management. Clin Rev Allergy Immunol (2018) 54(1):68–87. doi: 10.1007/s12016-017-8620-9

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Pratt CH, King LE Jr., Messenger AG, Christiano AM, Sundberg JP. Alopecia areata. Nat Rev Dis Primers (2017) 3:17011. doi: 10.1038/nrdp.2017.11

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Abou Rahal J, Kurban M, Kibbi AG, Abbas O. Plasmacytoid dendritic cells in alopecia areata: missing link? J Eur Acad Dermatol Venereol (2016) 30(1):119–23. doi: 10.1111/jdv.12932

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Pirastu N, Joshi PK, de Vries PS, Cornelis MC, McKeigue PM, Keum N, et al. GWAS for male-pattern baldness identifies 71 susceptibility loci explaining 38% of the risk. Nat Commun (2017) 8(1):1584. doi: 10.1038/s41467-017-01490-8

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Fu D, Huang J, Li K, Chen Y, He Y, Sun Y, et al. Dihydrotestosterone-induced hair regrowth inhibition by activating androgen receptor in C57BL6 mice simulates androgenetic alopecia. BioMed Pharmacother (2021) 137:111247. doi: 10.1016/j.biopha.2021.111247

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Jung YH, Chae CW, Choi GE, Shin HC, Lim JR, Chang HS, et al. Cyanidin 3-o-arabinoside suppresses DHT-induced dermal papilla cell senescence by modulating p38-dependent ER-mitochondria contacts. J BioMed Sci (2022) 29(1):17. doi: 10.1186/s12929-022-00800-7

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Gupta AK, Venkataraman M, Talukder M, Bamimore MA. Relative efficacy of minoxidil and the 5-alpha reductase inhibitors in androgenetic alopecia treatment of Male patients: A network meta-analysis. JAMA Dermatol (2022) 158(3):266–74. doi: 10.1001/jamadermatol.2021.5743

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Kaufman KD. Androgens and alopecia. Mol Cell Endocrinol (2002) 198(1-2):89–95. doi: 10.1016/s0303-7207(02)00372-6

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Yeap BB, Knuiman MW, Divitini ML, Hui J, Arscott GM, Handelsman DJ, et al. Epidemiological and mendelian randomization studies of dihydrotestosterone and estradiol and leukocyte telomere length in men. J Clin Endocrinol Metab (2016) 101(3):1299–306. doi: 10.1210/jc.2015-4139

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Yeap BB, Hui J, Knuiman MW, SAP C, Ken KYH, Flicker L, et al. Associations of plasma IGF1, IGFBP3 and estradiol with leucocyte telomere length, a marker of biological age, in men. Eur J Endocrinol (2020) 182(1):23–33. doi: 10.1530/EJE-19-0638

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Chen X, Hong X, Gao W, Luo S, Cai J, Liu G, et al. Causal relationship between physical activity, leisure sedentary behaviors and COVID-19 risk: a mendelian randomization study. J Transl Med (2022) 20(1):216. doi: 10.1186/s12967-022-03407-6

PubMed Abstract | CrossRef Full Text | Google Scholar