Genome-wide study of hair colour in UK Biobank explains most of the SNP heritability

Genome wide association study and meta-analysis identified multiple new risk loci for freckles in 4813 Chinese individuals

Freckle is a prevalent pigmentary dermatosis with an obvious hereditary component. Dozens of freckles risk loci have been discovered through research on multiple traits or other diseases, rather than as an independent trait. To discover novel variants associated with freckles, we performed GWAS and meta-analysis in 4813 Chinese individuals. We conducted GWAS and meta-analysis of two cohorts: 197 patients and 1603 controls (Cohort I), and 336 patients and 2677 controls (Cohort II), both from China. Then we performed linkage disequilibrium (LD) analysis, eQTL study, and enrichment analysis with association results for functional implications. Finally, we discovered 59 new SNPs and 13 novel susceptibility genes associated with freckles (Pmeta <5 × 10-8), which has enriched the genetic research on freckles.

Genetic influence on scar vascularity after burn injury in individuals of European ancestry: A prospective cohort study

After burn injury there is considerable variation in scar outcome, partially due to genetic factors. Scar vascularity is one characteristic that varies between individuals, and this study aimed to identify genetic variants contributing to different scar vascularity outcomes. An exome-wide array association study and gene pathway analysis was performed on a prospective cohort of 665 patients of European ancestry treated for burn injury, using their scar vascularity (SV) sub-score, part of the modified Vancouver Scar Scale (mVSS), as an outcome measure. DNA was genotyped using the Infinium HumanCoreExome-24 BeadChip, imputed to the Haplotype Reference Consortium panel. Associations between genetic variants (single nucleotide polymorphisms) and SV were estimated using an additive genetic model adjusting for sex, age, % total body surface area and number of surgical procedures, utilising linear and multinomial logistic regression. No individual genetic variants achieved the cut-off threshold for significance. Gene sets were also analysed using the Functional Mapping and Annotation (FUMA) platform, in which biological processes indirectly related to angiogenesis were significantly represented. This study suggests that SNPs in genes associated with angiogenesis may influence SV, but further studies with larger datasets are essential to validate these findings.

The Skin-Brain Axis: From UV and Pigmentation to Behaviour Modulation

The skin-brain axis has been suggested to play a role in several pathophysiological conditions, including opioid addiction, Parkinson's disease and many others. Recent evidence suggests that pathways regulating skin pigmentation may directly and indirectly regulate behaviour. Conversely, CNS-driven neural and hormonal responses have been demonstrated to regulate pigmentation, e.g., under stress. Additionally, due to the shared neuroectodermal origins of the melanocytes and neurons in the CNS, certain CNS diseases may be linked to pigmentation-related changes due to common regulators, e.g., MC1R variations. Furthermore, the HPA analogue of the skin connects skin pigmentation to the endocrine system, thereby allowing the skin to index possible hormonal abnormalities visibly. In this review, insight is provided into skin pigment production and neuromelanin synthesis in the brain and recent findings are summarised on how signalling pathways in the skin, with a particular focus on pigmentation, are interconnected with the central nervous system. Thus, this review may supply a better understanding of the mechanism of several skin-brain associations in health and disease.

Melanocortin-1 receptor (MC1R): a review for dermatologists

Melanocortin-1 receptor (MC1R) and its variants have a pivotal role in melanin synthesis. However, MC1R has been associated to non-pigmentary pathways related to DNA-repair activities and inflammation. The aim of this review is to provide an up-to-date overview about the role of MC1R in the skin. Specifically, after summarizing the current knowledge about MC1R structure and polymorphisms, we report data concerning the correlation between MC1R, phenotypic traits, skin aging, other diseases and skin cancers and their risk assessment through genetic testing.

Red-Haired People's Altered Responsiveness to Pain, Analgesics, and Hypnotics: Myth or Fact?-A Narrative Review

Red hair has been linked to altered sensitivity to pain, analgesics, and hypnotics. This alteration may be impacted by variants in the melanocortin-1 receptor (MC1R) gene, which are mainly found in redheads. The aim of this narrative review was to explore and present the current state of knowledge on red hair and its plausible associations with altered responsiveness to pain, analgesics, and hypnotics. Structured searches in the PubMed, CINAHL Complete, and Scopus electronic databases were conducted. Evidence suggests that women with red hair have an increased sensitivity to pain. Conversely, data also indicate a higher pain tolerance in homozygous carriers of MC1R variant alleles. Varied responses to analgesia have been reported, with both increased analgesic responsiveness in homozygous carriers of MC1R variant alleles and less analgesia in redheads. Data indicate an increased need for hypnotics in redheads. However, failed attempts to find statistical associations between red hair and altered responsiveness to hypnotics are also evident. Even though there seems to be an association between red hair and an altered responsiveness to pain, analgesics, and/or hypnotics, the results of this narrative review are inconclusive. Further research studies with larger populations and MC1R testing are needed.

Epistasis regulates genetic control of cardiac hypertrophy

The combinatorial effect of genetic variants is often assumed to be additive. Although genetic variation can clearly interact non-additively, methods to uncover epistatic relationships remain in their infancy. We develop low-signal signed iterative random forests to elucidate the complex genetic architecture of cardiac hypertrophy. We derive deep learning-based estimates of left ventricular mass from the cardiac MRI scans of 29,661 individuals enrolled in the UK Biobank. We report epistatic genetic variation including variants close to CCDC141 , IGF1R , TTN , and TNKS. Several loci where variants were deemed insignificant in univariate genome-wide association analyses are identified. Functional genomic and integrative enrichment analyses reveal a complex gene regulatory network in which genes mapped from these loci share biological processes and myogenic regulatory factors. Through a network analysis of transcriptomic data from 313 explanted human hearts, we found strong gene co-expression correlations between these statistical epistasis contributors in healthy hearts and a significant connectivity decrease in failing hearts. We assess causality of epistatic effects via RNA silencing of gene-gene interactions in human induced pluripotent stem cell-derived cardiomyocytes. Finally, single-cell morphology analysis using a novel high-throughput microfluidic system shows that cardiomyocyte hypertrophy is non-additively modifiable by specific pairwise interactions between CCDC141 and both TTN and IGF1R . Our results expand the scope of genetic regulation of cardiac structure to epistasis.

Biallelic loss-of-function variants of SLC12A9 cause lysosome dysfunction and a syndromic neurodevelopmental disorder

PURPOSE

Pathogenic variants of FIG4 generate enlarged lysosomes and neurological and developmental disorders. To identify additional genes regulating lysosomal volume, we carried out a genome-wide activation screen to detect suppression of enlarged lysosomes in FIG4-/- cells.

METHODS

The CRISPR-a gene activation screen utilized sgRNAs from the promoters of protein-coding genes. Fluorescence-activated cell sorting separated cells with correction of the enlarged lysosomes from uncorrected cells. Patient variants of SLC12A9 were identified by exome or genome sequencing and studied by segregation analysis and clinical characterization.

RESULTS

Overexpression of SLC12A9, a solute co-transporter, corrected lysosomal swelling in FIG4-/- cells. SLC12A9 (NP_064631.2) colocalized with LAMP2 at the lysosome membrane. Biallelic variants of SLC12A9 were identified in 3 unrelated probands with neurodevelopmental disorders. Common features included intellectual disability, skeletal and brain structural abnormalities, congenital heart defects, and hypopigmented hair. Patient 1 was homozygous for nonsense variant p.(Arg615∗), patient 2 was compound heterozygous for p.(Ser109Lysfs∗20) and a large deletion, and proband 3 was compound heterozygous for p.(Glu290Glyfs∗36) and p.(Asn552Lys). Fibroblasts from proband 1 contained enlarged lysosomes that were corrected by wild-type SLC12A9 cDNA. Patient variant p.(Asn552Lys) failed to correct the lysosomal defect.

CONCLUSION

Impaired function of SLC12A9 results in enlarged lysosomes and a recessive disorder with a recognizable neurodevelopmental phenotype.

iGWAS: Image-based genome-wide association of self-supervised deep phenotyping of retina fundus images

Existing imaging genetics studies have been mostly limited in scope by using imaging-derived phenotypes defined by human experts. Here, leveraging new breakthroughs in self-supervised deep representation learning, we propose a new approach, image-based genome-wide association study (iGWAS), for identifying genetic factors associated with phenotypes discovered from medical images using contrastive learning. Using retinal fundus photos, our model extracts a 128-dimensional vector representing features of the retina as phenotypes. After training the model on 40,000 images from the EyePACS dataset, we generated phenotypes from 130,329 images of 65,629 British White participants in the UK Biobank. We conducted GWAS on these phenotypes and identified 14 loci with genome-wide significance (p<5×10-8 and intersection of hits from left and right eyes). We also did GWAS on the retina color, the average color of the center region of the retinal fundus photos. The GWAS of retina colors identified 34 loci, 7 are overlapping with GWAS of raw image phenotype. Our results establish the feasibility of this new framework of genomic study based on self-supervised phenotyping of medical images.

Learning epistatic polygenic phenotypes with Boolean interactions

Detecting epistatic drivers of human phenotypes is a considerable challenge. Traditional approaches use regression to sequentially test multiplicative interaction terms involving pairs of genetic variants. For higher-order interactions and genome-wide large-scale data, this strategy is computationally intractable. Moreover, multiplicative terms used in regression modeling may not capture the form of biological interactions. Building on the Predictability, Computability, Stability (PCS) framework, we introduce the epiTree pipeline to extract higher-order interactions from genomic data using tree-based models. The epiTree pipeline first selects a set of variants derived from tissue-specific estimates of gene expression. Next, it uses iterative random forests (iRF) to search training data for candidate Boolean interactions (pairwise and higher-order). We derive significance tests for interactions, based on a stabilized likelihood ratio test, by simulating Boolean tree-structured null (no epistasis) and alternative (epistasis) distributions on hold-out test data. Finally, our pipeline computes PCS epistasis p-values that probabilisticly quantify improvement in prediction accuracy via bootstrap sampling on the test set. We validate the epiTree pipeline in two case studies using data from the UK Biobank: predicting red hair and multiple sclerosis (MS). In the case of predicting red hair, epiTree recovers known epistatic interactions surrounding MC1R and novel interactions, representing non-linearities not captured by logistic regression models. In the case of predicting MS, a more complex phenotype than red hair, epiTree rankings prioritize novel interactions surrounding HLA-DRB1, a variant previously associated with MS in several populations. Taken together, these results highlight the potential for epiTree rankings to help reduce the design space for follow up experiments.

Multi-ancestry tandem repeat association study of hair colour using exome-wide sequencing

Hair colour variation is influenced by hundreds of positions across the human genome but this genetic contribution has only been narrowly explored. Genome-wide association studies identified single nucleotide polymorphisms (SNPs) influencing hair colour but the biology underlying these associations is challenging to interpret. We report 16 tandem repeats (TRs) with effects on different models of hair colour plus two TRs associated with hair colour in diverse ancestry groups. Several of these TRs expand or contract amino acid coding regions of their localized protein such that structure, and by extension function, may be altered. We also demonstrate that independent of SNP variation, these TRs can be used to great an additive polygenic score that predicts darker hair colour. This work adds to the growing body of evidence regarding TR influence on human traits with relatively large and independent effects relative to surrounding SNP variation.

Integrative functional genomic analyses identify genetic variants influencing skin pigmentation in Africans

Skin color is highly variable in Africans, yet little is known about the underlying molecular mechanism. Here we applied massively parallel reporter assays to screen 1,157 candidate variants influencing skin pigmentation in Africans and identified 165 single-nucleotide polymorphisms showing differential regulatory activities between alleles. We combine Hi-C, genome editing and melanin assays to identify regulatory elements for MFSD12, HMG20B, OCA2, MITF, LEF1, TRPS1, BLOC1S6 and CYB561A3 that impact melanin levels in vitro and modulate human skin color. We found that independent mutations in an OCA2 enhancer contribute to the evolution of human skin color diversity and detect signals of local adaptation at enhancers of MITF, LEF1 and TRPS1, which may contribute to the light skin color of Khoesan-speaking populations from Southern Africa. Additionally, we identified CYB561A3 as a novel pigmentation regulator that impacts genes involved in oxidative phosphorylation and melanogenesis. These results provide insights into the mechanisms underlying human skin color diversity and adaptive evolution.

Natural hair color and skin cancers: A two-sample Mendelian randomization study

Previous observational studies have indicated an association between hair color and the risk of melanoma and keratinocyte skin cancer (KSC); however, different hair colors show inconsistent effects on skin cancers. Here, we conducted a two-sample Mendelian randomization (MR) study to evaluate the causal relationship between natural hair color and skin cancers by using 211 single nucleotide polymorphisms as genetic instruments from a genome-wide meta-analysis of 360,270 individuals of European ancestry. Light hair colors (red, blonde, and light brown) were associated with high levels of cutaneous melanoma (CM) and KSC (CM-inverse variance weighted [IVW] odds ratio [OR]-red: 1.034, 95% confidence interval [CI]: 1.025-1.044, P < 0.001; OR-blonde: 1.008, 95% CI: 1.003-1.014, P = 0.003; OR-light brown: 1.006, 95% CI: 1.002-1.011, P = 0.009; KSC-IVW OR-red: 1.078, 95% CI: 1.053-1.103, P < 0.001; OR-blonde: 1.024, 95% CI: 1.009-1.040, P = 0.002; OR-light brown: 1.018, 95% CI: 1.004-1.033, P = 0.01). However, dark brown hair showed an inverse causal relationship with skin cancers (CM IVW OR: 0.987, 95% CI: 0.984-0.990, P < 0.001; KSC IVW OR: 0.979, 95% CI: 0.970-0.988, P < 0.001). Black hair was associated with a decreased risk of KSC (IVW OR: 0.954, 95% CI: 0.913-0.997, P = 0.036) but showed no causal relationship with CM. The present study provides strong MR evidence of a causal association between hair color and skin cancer. Secondary MR analyses enhances result robustness by replicating findings, exploring gender-specific effects, and providing a more comprehensive understanding of the complex relationship between hair color and skin cancers. More large-scale MR studies or randomized controlled trials are required to further investigate the mechanisms of the association between hair color and skin cancers.

Linear hair growth rates in preschool children

BACKGROUND

Human scalp hair is a validated bio-substrate for monitoring various exposures in childhood including contextual stressors, environmental toxins, prescription or non-prescription drugs. Linear hair growth rates (HGR) are required to accurately interpret hair biomarker concentrations.

METHODS

We measured HGR in a prospective cohort of preschool children (N = 266) aged 9-72 months and assessed demographic factors, anthropometrics, and hair protein content (HPC). We examined HGR differences by age, sex, race, height, hair pigment, and season, and used univariable and multivariable linear regression models to identify HGR-related factors.

RESULTS

Infants below 1 year (288 ± 61 μm/day) had slower HGR than children aged 2-5 years (p = 0.0073). Dark-haired children (352 ± 52 μm/day) had higher HGR than light-haired children (325 ± 50 μm/day; p = 0.0019). Asian subjects had the highest HGR overall (p = 0.016). Younger children had higher HPC (p = 0.0014) and their HPC-adjusted HGRs were slower than older children (p = 0.0073). Age, height, hair pigmentation, and HPC were related to HGR in multivariable regression models.

CONCLUSIONS

We identified age, height, hair pigment, and hair protein concentration as significant determinants of linear HGRs. These findings help explain the known hair biomarker differences between children and adults and aid accurate interpretation of hair biomarker results in preschool children.

IMPACT

Discovery of hair biomarkers in the past few decades has transformed scientific disciplines like toxicology, pharmacology, epidemiology, forensics, healthcare, and developmental psychology. Identifying determinants of hair growth in children is essential for accurate interpretation of hair biomarker results in pediatric clinical studies. Childhood hair growth rates define the time-periods of biomarker incorporation into growing hair, essential for interpreting the biomarkers associated with environmental exposures and the mind-brain-body connectome. Our study describes age-, sex-, and height-based distributions of linear hair growth rates and provides determinants of linear hair growth rates in a large population of children. Age, height, hair pigmentation, and hair protein content are determinants of hair growth rates and should be accounted for in child hair biomarkers studies. Our findings on hair protein content and linear hair growth rates may provide physiological explanations for differences in hair growth rates and biomarkers in preschool children as compared to adults.

Genetic landscape of forensic DNA phenotyping markers among Mediterranean populations

Forensic DNA Phenotyping can reveal the appearance of an unknown individual by predicting the External Visible Characteristics (EVC) from DNA obtained at the crime scene. Our aim is to characterize the genetic landscape of Human identification markers responsible for EVC among Mediterranean populations compared to other worldwide groups. We conducted an exhaustive search for genes involved in EVC variation. Then, variants located on these genes were extracted from public genotypic data of Mediterranean, American, African and East Asiatic populations. The genetic landscape of these Human identification markers, their allelic distribution and admixture analyses, were determined using plink, R and ADMIXTURE softwares. Our results showed that the Mediterranean populations appear close to the Mexican populations and distinguished from sub Saharan African populations living in the USA and from East Asiatic populations. We highlighted a total of 103454 common variants shared between the studied populations and among them, 25 common variants associated with EVC. Interestingly, genotype frequencies results showed that the rs17646946, rs13016869, rs977588, rs1805008 and rs2240751 variants located respectively in the TCHH, PRKCE, OCA2, MC1R and MFSD12 genes are significantly different between the Mediterranean and Asiatic populations. The genotype frequencies of the variants rs977589 and rs7179994 located in the OCA2 gene, and of rs12913832 and rs2240751 located respectively in HERC2 and MFSD12 genes are significantly different between the Mediterranean and American populations. Our work generates a large number of EVC variants that could be a valuable resource for future studies in the forensic field.

Signatures of adaptation in myopia-related genes on the sunlight exposure hypothesis

BACKGROUND

Myopia is a common eye disorder that results from gene-environment interactions. The prevalence of myopia varies across populations, and exposure to bright sunlight may prevent its development. We hypothesize that local adaptation to light environments during human migration played a role in shaping the genetic basis of myopia, and we aim to investigate how the environment influences the genetic basis of myopia.

METHOD

We utilized the whole-genome variant data of the 1000 Genomes Project for analysis. We searched myopia-associated loci that were under selection in Europeans using population branch statistics and the number of segregating sites by length statistics. The outliers of these statistics were enriched in the Kyoto Encyclopedia of Genes and Genomes pathways and the gene ontology biological process terms in searching for pathways that were under selection. We applied Bayesian inference to estimate the correlation between environmental factors and allele frequencies of the selected loci and performed causal inference of myopia using two-sample Mendelian randomization analysis.

RESULTS

We detected signatures of adaptation in vision and light perception pathways, supporting our hypothesis of sunlight adaptation. We discovered a strong correlation between latitude and allele frequencies in genes that are under significant selection, and we found pleiotropic effects of pigmentation or circadian rhythm genes on myopia, indicating that sunlight exposure influences the genetic diversity of myopia.

CONCLUSIONS

Myopia genes involved in light perception showed signs of selection. Local adaptation during human migration shaped the genetic basis of myopia and may have influenced its global prevalence distribution.

Genetics of yellow-orange color variation in a pair of sympatric sulphur butterflies

Continuous color polymorphisms can serve as a tractable model for the genetic and developmental architecture of traits. Here we investigated continuous color variation in Colias eurytheme and Colias philodice, two species of sulphur butterflies that hybridize in sympatry. Using quantitative trait locus (QTL) analysis and high-throughput color quantification, we found two interacting large-effect loci affecting orange-to-yellow chromaticity. Knockouts of red Malpighian tubules (red), likely involved in endosomal maturation, result in depigmented wing scales. Additionally, the transcription factor bric-a-brac can act as a modulator of orange pigmentation. We also describe the QTL architecture of other continuously varying traits, together supporting a large-X effect model where the genetic control of species-defining traits is enriched on sex chromosomes. This study sheds light on the range of possible genetic architectures that can underpin a continuously varying trait and illustrates the power of using automated measurement to score phenotypes that are not always conspicuous to the human eye.

Addressing the ethical and societal challenges posed by genome-wide association studies of behavioral and brain-related traits

Genome-wide association studies have led to the identification of robust statistical associations of genetic variants with numerous brain-related traits, including neurological and psychiatric conditions, and psychological and behavioral measures. These results may provide insight into the biology underlying these traits and may facilitate clinically useful predictions. However, these results also carry the risk of harm, including possible negative effects of inaccurate predictions, violations of privacy, stigma and genomic discrimination, raising serious ethical and legal implications. Here, we discuss ethical concerns surrounding the results of genome-wide association studies for individuals, society and researchers. Given the success of genome-wide association studies and the increasing availability of nonclinical genomic prediction technologies, better laws and guidelines are urgently needed to regulate the storage, processing and responsible use of genetic data. Also, researchers should be aware of possible misuse of their results, and we provide guidance to help avoid such negative impacts on individuals and society.

Fine-mapping of retinal vascular complexity loci identifies Notch regulation as a shared mechanism with myocardial infarction outcomes

There is increasing evidence that the complexity of the retinal vasculature measured as fractal dimension, Df, might offer earlier insights into the progression of coronary artery disease (CAD) before traditional biomarkers can be detected. This association could be partly explained by a common genetic basis; however, the genetic component of Df is poorly understood. We present a genome-wide association study (GWAS) of 38,000 individuals with white British ancestry from the UK Biobank aimed to comprehensively study the genetic component of Df and analyse its relationship with CAD. We replicated 5 Df loci and found 4 additional loci with suggestive significance (P < 1e-05) to contribute to Df variation, which previously were reported in retinal tortuosity and complexity, hypertension, and CAD studies. Significant negative genetic correlation estimates support the inverse relationship between Df and CAD, and between Df and myocardial infarction (MI), one of CAD's fatal outcomes. Fine-mapping of Df loci revealed Notch signalling regulatory variants supporting a shared mechanism with MI outcomes. We developed a predictive model for MI incident cases, recorded over a 10-year period following clinical and ophthalmic evaluation, combining clinical information, Df, and a CAD polygenic risk score. Internal cross-validation demonstrated a considerable improvement in the area under the curve (AUC) of our predictive model (AUC = 0.770 ± 0.001) when comparing with an established risk model, SCORE, (AUC = 0.741 ± 0.002) and extensions thereof leveraging the PRS (AUC = 0.728 ± 0.001). This evidences that Df provides risk information beyond demographic, lifestyle, and genetic risk factors. Our findings shed new light on the genetic basis of Df, unveiling a common control with MI, and highlighting the benefits of its application in individualised MI risk prediction.

Genome-wide analysis identifies genetic effects on reproductive success and ongoing natural selection at the FADS locus

Identifying genetic determinants of reproductive success may highlight mechanisms underlying fertility and identify alleles under present-day selection. Using data in 785,604 individuals of European ancestry, we identified 43 genomic loci associated with either number of children ever born (NEB) or childlessness. These loci span diverse aspects of reproductive biology, including puberty timing, age at first birth, sex hormone regulation, endometriosis and age at menopause. Missense variants in ARHGAP27 were associated with higher NEB but shorter reproductive lifespan, suggesting a trade-off at this locus between reproductive ageing and intensity. Other genes implicated by coding variants include PIK3IP1, ZFP82 and LRP4, and our results suggest a new role for the melanocortin 1 receptor (MC1R) in reproductive biology. As NEB is one component of evolutionary fitness, our identified associations indicate loci under present-day natural selection. Integration with data from historical selection scans highlighted an allele in the FADS1/2 gene locus that has been under selection for thousands of years and remains so today. Collectively, our findings demonstrate that a broad range of biological mechanisms contribute to reproductive success.

Genome-Wide Scan of Wool Production Traits in Akkaraman Sheep

The objective of this study was to uncover the genetic background of wool quality, a production trait, by estimating genomic heritability and implementing GWAS in Akkaraman sheep. The wool characteristics measured included fibre diameter (FD) and staple length (SL) at the age of 8 months and yearling fibre diameter (YFD), yearling staple length (YSL) and yearling greasy fleece weight (YGFW) at 18 months of age. Animals were genotyped using the Axiom 50 K Ovine Genotyping Array. Maximum likelihood estimations of a linear mixed model (LMM) were used to estimate genomic heritability, where GWAS was conducted following a score test of each trait. Genomic heritability estimates for the traits ranged between 0.22 and 0.63, indicating that phenotypes have a moderate range of heritability. One genome- and six chromosome-wide significant SNPs were associated with the wool traits in Akkaraman lambs. Accordingly, TRIM2, MND1, TLR2, RNF175, CEP290, TMTC3, RERE, SLC45A1, SOX2, MORN1, SKI, FAAP20, PRKCZ, GABRD, CFAP74, CALML6 and TMEM52 genes as well as nine uncharacterized regions (LOC101118971, LOC105609137, LOC105603067, LOC101122892, LOC106991694, LOC106991467, LOC106991455, LOC105616534 and LOC105609719) were defined as plausible candidates. The findings of this study shed light on the genetics of wool quality and yield for the Akkaraman breed and suggests targets for breeders during systematic breeding programmes.

A Comprehensive Review of Mammalian Pigmentation: Paving the Way for Innovative Hair Colour-Changing Cosmetics

The natural colour of hair shafts is formed at the bulb of hair follicles, and it is coupled to the hair growth cycle. Three critical processes must happen for efficient pigmentation: (1) melanosome biogenesis in neural crest-derived melanocytes, (2) the biochemical synthesis of melanins (melanogenesis) inside melanosomes, and (3) the transfer of melanin granules to surrounding pre-cortical keratinocytes for their incorporation into nascent hair fibres. All these steps are under complex genetic control. The array of natural hair colour shades are ascribed to polymorphisms in several pigmentary genes. A myriad of factors acting via autocrine, paracrine, and endocrine mechanisms also contributes for hair colour diversity. Given the enormous social and cosmetic importance attributed to hair colour, hair dyeing is today a common practice. Nonetheless, the adverse effects of the long-term usage of such cosmetic procedures demand the development of new methods for colour change. In this context, case reports of hair lightening, darkening and repigmentation as a side-effect of the therapeutic usage of many drugs substantiate the possibility to tune hair colour by interfering with the biology of follicular pigmentary units. By scrutinizing mammalian pigmentation, this review pinpoints key targetable processes for the development of innovative cosmetics that can safely change the hair colour from the inside out.

Skin Phototype and Disease: A Comprehensive Genetic Approach to Pigmentary Traits Pleiotropy Using PRS in the GCAT Cohort

Human pigmentation has largely been associated with different disease prevalence among populations, but most of these studies are observational and inconclusive. Known to be genetically determined, pigmentary traits have largely been studied by Genome-Wide Association Study (GWAS), mostly in Caucasian ancestry cohorts from North Europe, identifying robustly, several loci involved in many of the pigmentary traits. Here, we conduct a detailed analysis by GWAS and Polygenic Risk Score (PRS) of 13 pigmentary-related traits in a South European cohort of Caucasian ancestry (n = 20,000). We observed fair phototype strongly associated with non-melanoma skin cancer and other dermatoses and confirmed by PRS-approach the shared genetic basis with skin and eye diseases, such as melanoma (OR = 0.95), non-melanoma skin cancer (OR = 0.93), basal cell carcinoma (OR = 0.97) and darker phototype with vitiligo (OR = 1.02), cataracts (OR = 1.04). Detailed genetic analyses revealed 37 risk loci associated with 10 out of 13 analyzed traits, and 16 genes significantly associated with at least two pigmentary traits. Some of them have been widely reported, such as MC1R, HERC2, OCA2, TYR, TYRP1, SLC45A2, and some novel candidate genes C1QTNF3, LINC02876, and C1QTNF3-AMACR have not been reported in the GWAS Catalog, with regulatory potential. These results highlight the importance of the assess phototype as a genetic proxy of skin functionality and disease when evaluating open mixed populations.

Investigating the morphology and genetics of scalp and facial hair characteristics for phenotype prediction

Microscopic traits and ultrastructure of hair such as cross-sectional shape, pigmentation, curvature, and internal structure help determine the level of variations between and across human populations. Apart from cosmetics and anthropological applications, such as determining species, somatic origin (body area), and biogeographic ancestry, the evidential value of hair has increased with rapid progression in the area of forensic DNA phenotyping (FDP). Individuals differ in the features of their scalp hair (greying, shape, colour, balding, thickness, and density) and facial hair (eyebrow thickness, monobrow, and beard thickness) features. Scalp and facial hair characteristics are genetically controlled and lead to visible inter-individual variations within and among populations of various ethnic origins. Hence, these characteristics can be exploited and made more inclusive in FDP, thereby leading to more comprehensive, accurate, and robust prediction models for forensic purposes. The present article focuses on understanding the genetics of scalp and facial hair characteristics with the goal to develop a more inclusive approach to better understand hair biology by integrating hair microscopy with genetics for genotype-phenotype correlation research.

A multi-phenotype analysis reveals 19 susceptibility loci for basal cell carcinoma and 15 for squamous cell carcinoma

Basal cell carcinoma and squamous cell carcinoma are the most common skin cancers, and have genetic overlap with melanoma, pigmentation traits, autoimmune diseases, and blood biochemistry biomarkers. In this multi-trait genetic analysis of over 300,000 participants from Europe, Australia and the United States, we reveal 78 risk loci for basal cell carcinoma (19 previously unknown and replicated) and 69 for squamous cell carcinoma (15 previously unknown and replicated). The previously unknown risk loci are implicated in cancer development and progression (e.g. CDKL1), pigmentation (e.g. TPCN2), cardiometabolic (e.g. FADS2), and immune-regulatory pathways for innate immunity (e.g. IFIH1), and HIV-1 viral load modulation (e.g. CCR5). We also report an optimised polygenic risk score for effective risk stratification for keratinocyte cancer in the Canadian Longitudinal Study of Aging (794 cases and 18139 controls), which could facilitate skin cancer surveillance e.g. in high risk subpopulations such as transplantees.

Association of TYR SNP rs1042602 with Melanoma Risk and Prognosis

Cutaneous melanoma is the most aggressive of skin tumors. In order to discover new biomarkers that could help us improve prognostic prediction in melanoma patients, we have searched for germline DNA variants associated with melanoma progression. Thus, after exome sequencing of a set of melanoma patients and healthy control individuals, we identified rs1042602, an SNP within TYR, as a good candidate. After genotyping rs1042602 in 1025 patients and 773 healthy donors, we found that the rs1042602-A allele was significantly associated with susceptibility to melanoma (CATT test: p = 0.0035). Interestingly, we also observed significant differences between patients with good and bad prognosis (5 years of follow-up) (n = 664) (CATT test for all samples p = 0.0384 and for men alone p = 0.0054). Disease-free-survival (DFS) analyses also showed that patients with the A allele had shorter DFS periods. In men, the association remained significant even in a multivariate Cox Proportional-hazards model, which was adjusted for age at diagnosis, Breslow thickness, ulceration and melanoma subtype (HR 0.4; 95% confidence interval (CI) 0.20-0.83; p = 0.0139). Based on our results, we propose that rs1042602-A is a risk allele for melanoma, which also seems to be responsible for a poorer prognosis of the disease, particularly in men.

Aberrant expression of agouti signaling protein (ASIP) as a cause of monogenic severe childhood obesity

Here we report a heterozygous tandem duplication at the ASIP (agouti signaling protein) gene locus causing ubiquitous, ectopic ASIP expression in a female patient with extreme childhood obesity. The mutation places ASIP under control of the ubiquitously active itchy E3 ubiquitin protein ligase promoter, driving the generation of ASIP in patient-derived native and induced pluripotent stem cells for all germ layers and hypothalamic-like neurons. The patient's phenotype of early-onset obesity, overgrowth, red hair and hyperinsulinemia is concordant with that of mutant mice ubiquitously expressing the homolog nonagouti. ASIP represses melanocyte-stimulating hormone-mediated activation as a melanocortin receptor antagonist, which might affect eating behavior, energy expenditure, adipocyte differentiation and pigmentation, as observed in the index patient. As the type of mutation escapes standard genetic screening algorithms, we rescreened the Leipzig Childhood Obesity cohort of 1,745 patients and identified four additional patients with the identical mutation, ectopic ASIP expression and a similar phenotype. Taken together, our data indicate that ubiquitous ectopic ASIP expression is likely a monogenic cause of human obesity.

Positive association between research competitiveness of Chinese academic hospitals and the scale of their biobanks: A national survey

Biobanks are important research infrastructure developed rapidly by Chinese hospitals. The objective of this study is to investigate the association between the comprehensive research competitiveness of hospitals and the development of hospital biobanks. In 2018, we conducted a national survey among Chinese biobank managers and directors. An online questionnaire was used to collect data of biobank characteristics. Of the 70 academic hospital biobanks responded to our survey, 49 of their hospitals were listed in the Science and Technology Evaluation Metrics (STEM) and 46 of their hospitals were listed in the Fudan Hospital Rankings, respectively, in 2018. Hospital scores from the STEM and Fudan Hospital Rankings were identified from their official websites. Multivariate linear regression analyses were used to assess the associations of STEM scores and Fudan Hospital Rankings with the scale of biobanks. The overall STEM score, Scientific and Technological Output, and Academic Impact in hospitals with large-scale biobanks were 48.35%, 55.16%, and 58.65% higher than those with small-scale biobanks, respectively. The scale of biobanks was positively associated with STEM score (β = 0.367, p = 0.009), Scientific and Technological Output (β = 0.441, p = 0.001), and Academic Impact (β = 0.304, p = 0.044) after adjustment for potential confounders. For Fudan Hospital Rankings, the comprehensive score and sustainable development ability score were higher in hospitals with large-scale biobanks. Further analyses showed that the scale of the biobanks was positively associated with a higher comprehensive score (β = 0.313, p = 0.037) and a sustainable development ability score (β = 0.463, p < 0.001). The scale of hospital biobanks was positively associated with the research competitiveness of Chinese hospitals.

Genome-wide association study of pain sensitivity assessed by questionnaire and the cold pressor test

ABSTRACT

We deployed an online pain sensitivity questionnaire (PSQ) and an at-home version of the cold pressor test (CPT) in a large genotyped cohort. We performed genome-wide association studies on the PSQ score (25,321 participants) and CPT duration (6853). We identified one new genome-wide significant locus associated with the PSQ score, which was located in the TSSC1 (also known as EIPR1 ) gene (rs58194899, OR = 0.950 [0.933-0.967], P -value = 1.9 × 10 -8 ). Although high pain sensitivity measured by both PSQ and CPT was associated with individual history of chronic and acute pains, genetic correlation analyses surprisingly suggested an opposite direction: PSQ score was inversely genetically correlated with neck and shoulder pain ( rg = -0.71), rheumatoid arthritis (-0.68), and osteoarthritis (-0.38), and with known risk factors, such as the length of working week (-0.65), smoking (-0.36), or extreme BMI (-0.23). Gene-based analysis followed by pathway analysis showed that genome-wide association studies results were enriched for genes expressed in the brain and involved in neuronal development and glutamatergic synapse signaling pathways. Finally, we confirmed that females with red hair were more sensitive to pain and found that genetic variation in the MC1R gene was associated with an increase in self-perceived pain sensitivity as assessed by the PSQ.

CNest: A novel copy number association discovery method uncovers 862 new associations from 200,629 whole-exome sequence datasets in the UK Biobank

Copy number variation (CNV) is known to influence human traits, having a rich history of research into common and rare genetic disease, and although CNV is accepted as an important class of genomic variation, progress on copy-number-based genome-wide association studies (GWASs) from next-generation sequencing (NGS) data has been limited. Here we present a novel method for large-scale copy number analysis from NGS data generating robust copy number estimates and allowing copy number GWASs (CN-GWASs) to be performed genome-wide in discovery mode. We provide a detailed analysis in the UK Biobank resource and a specifically designed software package. We use these methods to perform CN-GWAS analysis across 78 human traits, discovering over 800 genetic associations that are likely to contribute strongly to trait distributions. Finally, we compare CNV and SNP association signals across the same traits and samples, defining specific CNV association classes.

Genome-Wide Association Study for Screening and Identifying Potential Shin Color Loci in Ducks

Shin color diversity is a widespread phenomenon in birds. In this study, ducks were assessed to identify candidate genes for yellow, black, and spotted tibiae. For this purpose, we performed whole-genome resequencing of an F2 population consisting of 275 ducks crossed between Runzhou crested-white ducks and Cherry Valley ducks. We obtained 12.6 Mb of single nucleotide polymorphism (SNP) data, and the three shin colors were subsequently genotyped. Genome-wide association studies (GWASs) were performed to identify candidate and potential SNPs for the three shin colors. According to the results, 2947 and 3451 significant SNPs were associated with black and yellow shins, respectively, and six potential SNPs were associated with spotted shins. Based on the SNP annotations, the MITF, EDNRB2, POU family members, and the SLC superfamily were the candidate genes regulating pigmentation. In addition, the isoforms of EDNRB2, TYR, TYRP1, and MITF-M were significantly different between the black and yellow tibiae. MITF and EDNRB2 may have synergistic roles in the regulation of melanin synthesis, and their mutations may lead to phenotypic differences in the melanin deposition between individuals. This study provides new insights into the genetic factors that may influence tibia color diversity in birds.

Genome-Wide Analysis Identifies Candidate Genes Encoding Beak Color of Duck

Beak color diversity is a broadly occurring phenomenon in birds. Here, we used ducks to identify candidate genes for yellow, black, and spotted beaks. For this, an F2 population consisting of 275 ducks was genotyped using whole genome resequencing containing 12.6 M single-nucleotide polymorphisms (SNPs) and three beak colors. Genome-wide association studies (GWAS) was used to identify the candidate and potential SNPs for three beak colors in ducks (yellow, spotted, and black). The results showed that 2753 significant SNPs were associated with black beaks, 7462 with yellow, and 17 potential SNPs with spotted beaks. Based on SNP annotation, MITF, EDNRB2, members of the POU family, and the SLC superfamily were the candidate genes regulating pigmentation. Meanwhile, isoforms MITF-M and EDNRB2 were significantly different between black and yellow beaks. MITF and EDNRB2 likely play a synergistic role in the regulation of melanin synthesis, and their mutations contribute to phenotypic differences in beak melanin deposition among individuals. This study provides new insights into genetic factors that may influence the diversity of beak color.

Genome-Wide Association Study Identifies Two Common Loci Associated with Pigment Dispersion Syndrome/Pigmentary Glaucoma and Implicates Myopia in its Development

PURPOSE

To identify genetic variants associated with pigment dispersion syndrome (PDS) and pigmentary glaucoma (PG) in unrelated patients and to further understand the genetic and potentially causal relationships between PDS and associated risk factors.

DESIGN

A 2-stage genome-wide association meta-analysis with replication and subsequent in silico analyses including Mendelian randomization.

PARTICIPANTS

A total of 574 cases with PG or PDS and 52 627 controls of European descent.

METHODS

Genome-wide association analyses were performed in 4 cohorts and meta-analyzed in 3 stages: (1) a discovery meta-analysis was performed in 3 cohorts, (2) replication was performed in the fourth cohort, and (3) all 4 cohorts were meta-analyzed to increase statistical power. Two-sample Mendelian randomization was used to determine whether refractive error and intraocular pressure exert causal effects over PDS.

MAIN OUTCOME MEASURES

The association of genetic variants with PDS and whether myopia exerts causal effects over PDS.

RESULTS

Significant association was present at 2 novel loci for PDS/PG. These loci and follow-up analyses implicate the genes gamma secretase activator protein (GSAP) (lead single nucleotide polymorphism [SNP]: rs9641220, P = 6.0×10-10) and glutamate metabotropic receptor 5 (GRM5)/TYR (lead SNP: rs661177, P = 3.9×10-9) as important factors in disease risk. Mendelian randomization showed significant evidence that negative refractive error (myopia) exerts a direct causal effect over PDS (P = 8.86×10-7).

CONCLUSIONS

Common SNPs relating to the GSAP and GRM5/TYR genes are associated risk factors for the development of PDS and PG. Although myopia is a known risk factor, this study uses genetic data to demonstrate that myopia is, in part, a cause of PDS and PG.

The Skin We Live in: Pigmentation Traits and Tanning Behaviour in British Young Adults, an Observational and Genetically-Informed Study

Skin cancer incidence has been increasing worldwide, representing a particularly high burden for populations of European ancestry. Outdoor and indoor tanning using ultraviolet (UV) radiation devices are major risk factors for skin cancer. While tanning behaviours can be modified by targeted interventions to reduce skin cancer rates, there is insufficient evidence on the motivations for tanning preferences and their relationship with pigmentation phenotypes. The present observational and genetically-informed study investigates motives for tanning and the role that pigmentation phenotypes play on outdoor and indoor tanning behaviour in British young adults. This study included 3722 participants from the Avon Longitudinal Study of Parents and Children in South West England, with data on pigmentation features, tanning ability and preferences, and SNP genotypes. Liking to tan and outdoor tanning were strongly influenced by pigmentary traits and tanning ability. However, the association of these phenotypes with UV indoor tanning was weaker. Our results provide evidence to support the implementation of skin cancer preventative interventions that consider individual biological characteristics and motives for undergoing outdoor and indoor tanning.

A genome-wide epistatic network underlies the molecular architecture of continuous color variation of body extremities

Deciphering the molecular architecture of coat coloration for a better understanding of the biological mechanisms underlying pigmentation still remains a challenge. We took advantage of a rabbit French experimental population in which both a pattern and a gradient of coloration from white to brown segregated within the himalayan phenotype. The whole experimental design was genotyped using the high density Affymetrix® AxiomOrcun™ SNP Array and phenotyped into 6 different groups ordered from the lighter to the darker. Genome-wide association analyses pinpointed an oligogenic determinism, under recessive and additive inheritance, involving genes already known in melanogenesis (ASIP, KIT, MC1R, TYR), and likely processed pseudogenes linked to ribosomal function, RPS20 and RPS14. We also identified (i) gene-gene interactions through ASIP:MC1R affecting light cream/beige phenotypes while KIT:RPS responsible of dark chocolate/brown colors and (ii) a genome-wide epistatic network involving several others coloration genes such as POT1 or HPS5. Finally, we determined the recessive inheritance of the English spotting phenotype likely involving a copy number variation affecting at least the end of the coding sequence of the KIT gene. Our analyses of coloration as a continuous trait allowed us to go beyond much of the established knowledge through the detection of additional genes and gene-gene interactions that may contribute to the molecular architecture of the coloration phenotype.

Genome-wide association study of actinic keratosis identifies new susceptibility loci implicated in pigmentation and immune regulation pathways

Actinic keratosis (AK) is a common precancerous cutaneous neoplasm that arises on chronically sun-exposed skin. AK susceptibility has a moderate genetic component, and although a few susceptibility loci have been identified, including IRF4, TYR, and MC1R, additional loci have yet to be discovered. We conducted a genome-wide association study of AK in non-Hispanic white participants of the Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort (n = 63,110, discovery cohort), with validation in the Mass-General Brigham (MGB) Biobank cohort (n = 29,130). We identified eleven loci (P < 5 × 10-8), including seven novel loci, of which four novel loci were validated. In a meta-analysis (GERA + MGB), one additional novel locus, TRPS1, was identified. Genes within the identified loci are implicated in pigmentation (SLC45A2, IRF4, BNC2, TYR, DEF8, RALY, HERC2, and TRPS1), immune regulation (FOXP1 and HLA-DQA1), and cell signaling and tissue remodeling (MMP24) pathways. Our findings provide novel insight into the genetics and pathogenesis of AK susceptibility.

Cis-epistasis at the LPA locus and risk of cardiovascular diseases

AIMS

Coronary artery disease (CAD) has a strong genetic predisposition. However, despite substantial discoveries made by genome-wide association studies (GWAS), a large proportion of heritability awaits identification. Non-additive genetic effects might be responsible for part of the unaccounted genetic variance. Here, we attempted a proof-of-concept study to identify non-additive genetic effects, namely epistatic interactions, associated with CAD.

METHODS AND RESULTS

We tested for epistatic interactions in 10 CAD case-control studies and UK Biobank with focus on 8068 SNPs at 56 loci with known associations with CAD risk. We identified a SNP pair located in cis at the LPA locus, rs1800769 and rs9458001, to be jointly associated with risk for CAD [odds ratio (OR) = 1.37, P = 1.07 × 10-11], peripheral arterial disease (OR = 1.22, P = 2.32 × 10-4), aortic stenosis (OR = 1.47, P = 6.95 × 10-7), hepatic lipoprotein(a) (Lp(a)) transcript levels (beta = 0.39, P = 1.41 × 10-8), and Lp(a) serum levels (beta = 0.58, P = 8.7 × 10-32), while individual SNPs displayed no association. Further exploration of the LPA locus revealed a strong dependency of these associations on a rare variant, rs140570886, that was previously associated with Lp(a) levels. We confirmed increased CAD risk for heterozygous (relative OR = 1.46, P = 9.97 × 10-32) and individuals homozygous for the minor allele (relative OR = 1.77, P = 0.09) of rs140570886. Using forward model selection, we also show that epistatic interactions between rs140570886, rs9458001, and rs1800769 modulate the effects of the rs140570886 risk allele.

CONCLUSIONS

These results demonstrate the feasibility of a large-scale knowledge-based epistasis scan and provide rare evidence of an epistatic interaction in a complex human disease. We were directed to a variant (rs140570886) influencing risk through additive genetic as well as epistatic effects. In summary, this study provides deeper insights into the genetic architecture of a locus important for cardiovascular diseases.

Citrus-Gene interaction and melanoma risk in the UK Biobank

High citrus consumption may increase melanoma risk; however, little is known about the biological mechanisms of this association, or whether it is modified by genetic variants. We conducted a genome-wide analysis of gene-citrus consumption interactions on melanoma risk among 1563 melanoma cases and 193 296 controls from the UK Biobank. Both the 2-degrees-of-freedom (df) joint test of genetic main effect and gene-environment (G-E) interaction and the standard 1-df G-E interaction test were performed. Three index SNPs (lowest P-value SNP among highly correlated variants [r²  > .6]) were identified from among the 365 genome-wide significant 2-df test results (rs183783391 on chromosome 3 [MITF], rs869329 on chromosome 9 [MTAP] and rs11446223 on chromosome 16 [DEF8]). Although all three were statistically significant for the 2-df test (4.25e-08, 1.98e-10 and 4.93e-13, respectively), none showed evidence of interaction according to the 1-df test (P = .73, .24 and .12, respectively). Eight nonindex, 2-df test significant SNPs on chromosome 16 were significant (P < .05) according to the 1-df test, providing evidence of citrus-gene interaction. Seven of these SNPs were mapped to AFG3L1P (rs199600347, rs111822773, rs113178244, rs3803683, rs73283867, rs78800020, rs73283871), and one SNP was mapped to GAS8 (rs74583214). We identified several genetic loci that may elucidate the association between citrus consumption and melanoma risk. Further studies are needed to confirm these findings.

Predicting Physical Appearance from DNA Data-Towards Genomic Solutions

The idea of forensic DNA intelligence is to extract from genomic data any information that can help guide the investigation. The clues to the externally visible phenotype are of particular practical importance. The high heritability of the physical phenotype suggests that genetic data can be easily predicted, but this has only become possible with less polygenic traits. The forensic community has developed DNA-based predictive tools by employing a limited number of the most important markers analysed with targeted massive parallel sequencing. The complexity of the genetics of many other appearance phenotypes requires big data coupled with sophisticated machine learning methods to develop accurate genomic predictors. A significant challenge in developing universal genomic predictive methods will be the collection of sufficiently large data sets. These should be created using whole-genome sequencing technology to enable the identification of rare DNA variants implicated in phenotype determination. It is worth noting that the correctness of the forensic sketch generated from the DNA data depends on the inclusion of an age factor. This, however, can be predicted by analysing epigenetic data. An important limitation preventing whole-genome approaches from being commonly used in forensics is the slow progress in the development and implementation of high-throughput, low DNA input sequencing technologies. The example of palaeoanthropology suggests that such methods may possibly be developed in forensics.

Bayesian model comparison for rare-variant association studies

Whole-genome sequencing studies applied to large populations or biobanks with extensive phenotyping raise new analytic challenges. The need to consider many variants at a locus or group of genes simultaneously and the potential to study many correlated phenotypes with shared genetic architecture provide opportunities for discovery not addressed by the traditional one variant, one phenotype association study. Here, we introduce a Bayesian model comparison approach called MRP (multiple rare variants and phenotypes) for rare-variant association studies that considers correlation, scale, and direction of genetic effects across a group of genetic variants, phenotypes, and studies, requiring only summary statistic data. We apply our method to exome sequencing data (n = 184,698) across 2,019 traits from the UK Biobank, aggregating signals in genes. MRP demonstrates an ability to recover signals such as associations between PCSK9 and LDL cholesterol levels. We additionally find MRP effective in conducting meta-analyses in exome data. Non-biomarker findings include associations between MC1R and red hair color and skin color, IL17RA and monocyte count, and IQGAP2 and mean platelet volume. Finally, we apply MRP in a multi-phenotype setting; after clustering the 35 biomarker phenotypes based on genetic correlation estimates, we find that joint analysis of these phenotypes results in substantial power gains for gene-trait associations, such as in TNFRSF13B in one of the clusters containing diabetes- and lipid-related traits. Overall, we show that the MRP model comparison approach improves upon useful features from widely used meta-analysis approaches for rare-variant association analyses and prioritizes protective modifiers of disease risk.

Macular thickness varies with age-related macular degeneration genetic risk variants in the UK Biobank cohort

To evaluate the influence AMD risk genomic variants have on macular thickness in the normal population. UK Biobank participants with no significant ocular history were included using the UK Biobank Resource (project 2112). Spectral-domain optical coherence tomography (SD-OCT) images were taken and segmented to define retinal layers. The influence of AMD risk single-nucleotide polymorphisms (SNP) on retinal layer thickness was analysed. AMD risk associated SNPs were strongly associated with outer-retinal layer thickness. The inner-segment outer segment (ISOS)-retinal pigment epithelium (RPE) thickness measurement, representing photoreceptor outer segments was most significantly associated with the cumulative polygenic risk score, composed of 33 AMD-associated variants, resulting in a decreased thickness (p = 1.37 × 10-67). Gene-gene interactions involving the NPLOC4-TSPAN10 SNP rs6565597 were associated with significant changes in outer retinal thickness. Thickness of outer retinal layers is highly associated with the presence of risk AMD SNPs. Specifically, the ISOS-RPE measurement. Changes to ISOS-RPE thickness are seen in clinically normal individuals with AMD risk SNPs suggesting structural changes occur at the macula prior to the onset of disease symptoms or overt clinical signs.

A large Canadian cohort provides insights into the genetic architecture of human hair colour

Hair colour is a polygenic phenotype that results from differences in the amount and ratio of melanins located in the hair bulb. Genome-wide association studies (GWAS) have identified many loci involved in the pigmentation pathway affecting hair colour. However, most of the associated loci overlap non-protein coding regions and many of the molecular mechanisms underlying pigmentation variation are still not understood. Here, we conduct GWAS meta-analyses of hair colour in a Canadian cohort of 12,741 individuals of European ancestry. By performing fine-mapping analyses we identify candidate causal variants in pigmentation loci associated with blonde, red and brown hair colour. Additionally, we observe colocalization of several GWAS hits with expression and methylation quantitative trait loci (QTLs) of cultured melanocytes. Finally, transcriptome-wide association studies (TWAS) further nominate the expression of EDNRB and CDK10 as significantly associated with hair colour. Our results provide insights on the mechanisms regulating pigmentation biology in humans.

Melanocytes in black-boned chicken have immune contribution under infectious bursal disease virus infection

In black-boned chicken, melanocytes are widely distributed in their inner organs. However, the roles of these cells are not fully elucidated. In this study, we used 3-wk-old female Silky Fowl to investigate the functions of melanocytes under infection with infectious bursal disease virus (IBDV). We found the melanocytes in the bursa of Fabricius involved in IBDV infection shown as abundant melanin were transported into the nodule and lamina propria where obvious apoptotic cells and higher expression of BAX were detected. Genes related to the toll-like receptor (TLR) signaling pathway were highly detected by quantitative PCR, including TLR1, TLR3, TLR4, TLR15, myeloid differential protein-88, interferon-α, and interferon-β. We then isolated and infected primary melanocytes with IBDV in vitro and found that higher expressions of immune genes were detected at 24 and 48 h after infection; the upregulated innate and adaptive immune genes were involved in the pathogenesis of IBDV infection, including TLR3, TLR7, interleukin 15 (IL15), IL18, IL1rap, CD7, BG2, ERAP1, and SLA2. These changes in gene expression were highly associated with microtubule-based movement, antigen processing and presentation, defense against viruses, and innate immune responses. Our results indicated that the widely distributed melanocytes in Silky Fowl could migrate to play important innate immune roles during virus infection.

Cell-type-specific meQTLs extend melanoma GWAS annotation beyond eQTLs and inform melanocyte gene-regulatory mechanisms

Although expression quantitative trait loci (eQTLs) have been powerful in identifying susceptibility genes from genome-wide association study (GWAS) findings, most trait-associated loci are not explained by eQTLs alone. Alternative QTLs, including DNA methylation QTLs (meQTLs), are emerging, but cell-type-specific meQTLs using cells of disease origin have been lacking. Here, we established an meQTL dataset by using primary melanocytes from 106 individuals and identified 1,497,502 significant cis-meQTLs. Multi-QTL colocalization with meQTLs, eQTLs, and mRNA splice-junction QTLs from the same individuals together with imputed methylome-wide and transcriptome-wide association studies identified candidate susceptibility genes at 63% of melanoma GWAS loci. Among the three molecular QTLs, meQTLs were the single largest contributor. To compare melanocyte meQTLs with those from malignant melanomas, we performed meQTL analysis on skin cutaneous melanomas from The Cancer Genome Atlas (n = 444). A substantial proportion of meQTL probes (45.9%) in primary melanocytes is preserved in melanomas, while a smaller fraction of eQTL genes is preserved (12.7%). Integration of melanocyte multi-QTLs and melanoma meQTLs identified candidate susceptibility genes at 72% of melanoma GWAS loci. Beyond GWAS annotation, meQTL-eQTL colocalization in melanocytes suggested that 841 unique genes potentially share a causal variant with a nearby methylation probe in melanocytes. Finally, melanocyte trans-meQTLs identified a hotspot for rs12203592, a cis-eQTL of a transcription factor, IRF4, with 131 candidate target CpGs. Motif enrichment and IRF4 ChIP-seq analysis demonstrated that these target CpGs are enriched in IRF4 binding sites, suggesting an IRF4-mediated regulatory network. Our study highlights the utility of cell-type-specific meQTLs.

Investigating DNA methylation as a potential mediator between pigmentation genes, pigmentary traits and skin cancer

Pigmentation characteristics are well-known risk factors for skin cancer. Polymorphisms in pigmentation genes have been associated with these traits and with the risk of malignancy. However, the functional relationship between genetic variation and disease is still unclear. This study aims to assess whether pigmentation SNPs are associated with pigmentary traits and skin cancer via DNA methylation (DNAm). Using a meta-GWAS of whole-blood DNAm from 36 European cohorts (N = 27,750; the Genetics of DNA Methylation Consortium, GoDMC), we found that 19 out of 27 SNPs in 10 pigmentation genes were associated with 391 DNAm sites across 30 genomic regions. We examined the effect of 25 selected DNAm sites on pigmentation traits, sun exposure phenotypes and skin cancer and on gene expression in whole blood. We uncovered an association of DNAm site cg07402062 with red hair in the Avon Longitudinal Study of Parents and Children (ALSPAC). We also found that the expression of ASIP and CDK10 was associated with hair colour, melanoma and basal cell carcinoma. Our results indicate that DNAm and expression of pigmentation genes may play a role as potential mediators of the relationship between genetic variants, pigmentation phenotypes and skin cancer and thus deserve further scrutiny.

MC1R Functions, Expression, and Implications for Targeted Therapy

The G protein-coupled MC1R is expressed in melanocytes and has a pivotal role in human skin pigmentation, with reduced function in human genetic variants exhibiting a red hair phenotype and increased melanoma predisposition. Beyond its role in pigmentation, MC1R is increasingly recognized as promoting UV-induced DNA damage repair. Consequently, there is mounting interest in targeting MC1R for therapeutic benefit. However, whether MC1R expression is restricted to melanocytes or is more widely expressed remains a matter of debate. In this paper, we review MC1R function and highlight that unbiased analysis suggests that its expression is restricted to melanocytes, granulocytes, and the brain.

Behind the Scene: Exploiting MC1R in Skin Cancer Risk and Prevention

Melanoma and non-melanoma skin cancers (NMSCs) are the most frequent cancers of the skin in white populations. An increased risk in the development of skin cancers has been associated with the combination of several environmental factors (i.e., ultraviolet exposure) and genetic background, including melanocortin-1 receptor (MC1R) status. In the last few years, advances in the diagnosis of skin cancers provided a great impact on clinical practice. Despite these advances, NMSCs are still the most common malignancy in humans and melanoma still shows a rising incidence and a poor prognosis when diagnosed at an advanced stage. Efforts are required to underlie the genetic and clinical heterogeneity of melanoma and NMSCs, leading to an optimization of the management of affected patients. The clinical implications of the impact of germline MC1R variants in melanoma and NMSCs' risk, together with the additional risk conferred by somatic mutations in other peculiar genes, as well as the role of MC1R screening in skin cancers' prevention will be addressed in the current review.

A custom capture sequence approach for oculocutaneous albinism identifies structural variant alleles at the OCA2 locus

Oculocutaneous albinism (OCA) is a heritable disorder of pigment production that manifests as hypopigmentation and altered eye development. Exon sequencing of known OCA genes is unsuccessful in producing a complete molecular diagnosis for a significant number of affected individuals. We sequenced the DNA of individuals with OCA using short-read custom capture sequencing that targeted coding, intronic, and noncoding regulatory regions of known OCA genes, and genome-wide association study-associated pigmentation loci. We identified an OCA2 complex structural variant (CxSV), defined by a 143 kb inverted segment reintroduced in intron 1, upstream of the native location. The corresponding CxSV junctions were observed in 11/390 probands screened. The 143 kb CxSV presents in one family as a copy number variant duplication for the 143 kb region. In the remaining 10/11 families, the 143 kb CxSV acquired an additional 184 kb deletion across the same region, restoring exons 3-19 of OCA2 to a copy-number neutral state. Allele-associated haplotype analysis found rare SNVs rs374519281 and rs139696407 are linked with the 143 kb CxSV in both OCA2 alleles. For individuals in which customary molecular evaluation does not reveal a biallelic OCA diagnosis, we recommend preliminary screening for these haplotype-associated rare variants, followed by junction-specific validation for the OCA2 143 kb CxSV.

Exploiting within-breed variability in the autochthonous Reggiana breed identified several candidate genes affecting pigmentation-related traits, stature and udder defects in cattle

Autochthonous cattle breeds constitute important reservoirs of genetic diversity. Reggiana is an Italian local cattle breed reared in the north of Italy for the production of a mono‐breed Parmigiano–Reggiano cheese. Reggiana cattle usually have a classical solid red coat colour and pale muzzle. As part of the strategies designed for the sustainable conservation of this genetic resource, we investigated at the genome‐wise level the within‐breed detected variability of three pigmentation‐related traits (intensity of red coat colour, based on three classes – light/diluted, normal and dark; spotted patterns/piebaldism that sometime emerge in the breed; muzzle colour – pink/pale, grey and black), stature, presence/absence and number of supernumerary teats and teat length. A total of 1776 Reggiana cattle (about two‐thirds of the extant breed population) were genotyped with the GeneSeek GGP Bovine 150k SNP array and single‐marker and haplotype‐based GWASs were carried out. The results indicated that two main groups of genetic factors affect the intensity of red coat colour: darkening genes (including EDN3 and a few other genes) and diluting genes (including PMEL and a few other genes). Muzzle colour was mainly determined by MC1R gene markers. Piebaldism was mainly associated with KIT gene markers. Stature was associated with BTA6 markers upstream of the NCAPG–LCORL genes. Teat defects were associated with TBX3/TBX5, MCC and LGR5 genes. Overall, the identified genomic regions not only can be directly used in selection plans in the Reggiana breed, but also contribute to clarifying the genetic mechanisms involved in determining exterior traits in cattle.