Genetics of facial telangiectasia in the Rotterdam Study: a genome-wide association study and candidate gene approach.

BACKGROUND: The severity of facial telangiectasia or red veins is associated with many lifestyle factors. However, the genetic predisposition remains unclear. OBJECTIVES: We performed a genome-wide association study (GWAS) on facial telangiectasia in the Rotterdam Study (RS) and tested for replication in two independent cohorts. Additionally, a candidate gene approach with known pigmentation genes was performed. METHODS: Facial telangiectasia were extracted from standardized facial photographs (collected from 2010-2013) of 2842 northwestern European participants (median age 66.9, 56.8% female) from the RS. Our GWAS top hits (P-value <10-6 ) were tested for replication in 460 elderly women of the SALIA cohort and in 576 additional men and women of the RS. Associations of top single nucleotide polymorphisms (SNPs) with expression quantitative trait loci (eQTL) in various tissues were reviewed (GTEx database) alongside phenotype associations in the UK biobank database. SNP-based associations between known pigmentation genes and facial telangiectasia were tested. Conditional analysis on skin colour was additionally performed. RESULTS: Our most significant GWAS signal was rs4417318 (P-value 5.38*10-7 ), an intergenic SNP on chromosome 12 mapping to the SLC16A7 gene. Other suggestive SNPs tagged genes ZNF211, ZSCAN4, ICOS and KCNN3; SNP eQTLs and phenotype associations tagged links to the vascular system. However, the top signals did not pass significance in the two replication cohorts. The pigmentation genes KIAA0930, SLCA45A2 and MC1R, were significantly associated with telangiectasia in a candidate gene approach but not independently of skin colour. CONCLUSION: In this GWAS on telangiectasia in a northwestern European population, no genome-wide significant SNPs were found, although suggestive signals indicate genes involved in the vascular system might be involved in telangiectasia. Significantly associated pigmentation genes underline the link between skin colour and telangiectasia.

Predicting keratinocyte carcinoma in patients with actinic keratosis: development and internal validation of a multivariable risk-prediction model.

BACKGROUND: Patients with actinic keratosis (AK) are at increased risk for developing keratinocyte carcinoma (KC) but predictive factors and their risk rates are unknown. OBJECTIVES: To develop and internally validate a prediction model to calculate the absolute risk of a first KC in patients with AK. METHODS: The risk-prediction model was based on the prospective population-based Rotterdam Study cohort. We hereto analysed the data of participants with at least one AK lesion at cohort baseline using a multivariable Cox proportional hazards model and included 13 a priori defined candidate predictor variables considering phenotypic, genetic and lifestyle risk factors. KCs were identified by linkage of the data with the Dutch Pathology Registry. RESULTS: Of the 1169 AK participants at baseline, 176 (15·1%) developed a KC after a median follow-up of 1·8 years. The final model with significant predictors was obtained after backward stepwise selection and comprised the presence of four to nine AKs [hazard ratio (HR) 1·68, 95% confidence interval (CI) 1·17-2·42], 10 or more AKs (HR 2·44, 95% CI 1·65-3·61), AK localization on the upper extremities (HR 0·75, 95% CI 0·52-1·08) or elsewhere except the head (HR 1·40, 95% CI 0·98-2·01) and coffee consumption (HR 0·92, 95% CI 0·84-1·01). Evaluation of the discriminative ability of the model showed a bootstrap validated concordance index (c-index) of 0·60. CONCLUSIONS: We showed that the risk of KC in patients with AK can be calculated with the use of four easily assessable predictor variables. Given the c-index, extension of the model with additional, currently unknown predictor variables is desirable. Linked Comment: Kim et al. Br J Dermatol 2020; 183:415-416.

Principal component analysis of seven skin-ageing features identifies three main types of skin ageing.

BACKGROUND: The underlying phenotypic correlations between wrinkles, pigmented spots (PS), telangiectasia and other related facial-ageing subphenotypes are not well understood. OBJECTIVES: To analyse the underlying phenotypic correlation structure between seven features for facial ageing: global wrinkling, perceived age (PA), Griffiths photodamage grading, PS, telangiectasia, actinic keratosis (AK) and keratinocyte cancer (KC). METHODS: This was a cross-sectional study. Facial photographs and a full-body skin examination were used. We used principal component analysis (PCA) to derive principal components (PCs) of common variation between the features. We performed multivariable linear regressions between age, sex, body mass index, smoking and ultraviolet radiation exposure and the PC scores derived from PCA. We also tested the association between the main PC scores and 140 single-nucleotide polymorphisms (SNPs) previously associated with skin-ageing phenotypes. RESULTS: We analysed data from 1790 individuals with complete data on seven features of skin ageing. Three main PCs explained 73% of the total variance of the ageing phenotypes: a hypertrophic/wrinkling component (PC1: global wrinkling, PA and Griffiths grading), an atrophic/skin colour component (PC2: PS and telangiectasia) and a cancerous component (PC3: AK and KC). The associations between lifestyle and host factors differed per PC. The strength of SNP associations also differed per component with the most SNP associations found with the atrophic component [e.g. the IRF4 SNP (rs12203592); P-value = 1·84 × 10-22 ]. CONCLUSIONS: Using a hypothesis-free approach, we identified three major underlying phenotypes associated with extrinsic ageing. Associations between determinants for skin ageing differed in magnitude and direction per component. What's already known about this topic? Facial ageing is a complex phenotype consisting of different features including wrinkles, pigmented changes, telangiectasia and cancerous-related growths; it is not clear how these phenotypes are related to each other and to other phenotypes. A few studies have described two main clinical phenotypes for photoageing, namely hypertrophic ageing and atrophic ageing, which have been based solely on the clinical assessment of photoageing characteristics. What does this study add? We are the first to use epidemiology data to identify three main components associated with photoageing, namely a hypertrophic component (global wrinkling; perceived age; Griffiths grading) and atrophic component (pigmented spots; telangiectasia) and a cancer component (actinic keratosis; keratinocyte cancer). Association analysis showed different effects and direction of environmental determinants and genetic associations with the three components, with the most significant gene variants associations found for the atrophic component.

Epidemiology and determinants of facial telangiectasia: a cross-sectional study.

BACKGROUND: Telangiectasia or red veins are one of the prominent features of facial skin ageing. To date, there are few studies investigating the determinants of telangiectasia. OBJECTIVES: We investigated lifestyle and physiological factors associated with facial telangiectasia in a large prospective Dutch cohort study. METHODS: Telangiectasia was quantified digitally from standardized facial photographs of 2842 North European participants (56.8% female, median age 66.9) from the Rotterdam Study, collected in 2010-2013. Effect estimates from multivariable linear regressions are presented as the percentage difference in the mean value of telangiectasia area per unit increase of a determinant (%Δ) with corresponding 95% CI. RESULTS: Significant determinants were older age [1.7%Δ per year (95% CI 1.4, 2.0)], female sex [18.3%Δ (95% CI 13.2, 23.6)], smoking [current versus never 38.4%Δ (95% CI 30.3, 47.0); former versus never 11.6%Δ (95% CI 6.6, 16.9)], a high susceptibility to sunburn [10.2%Δ (95% CI 5.4, 15.3)] and light skin colour [pale versus white-to-olive 31.4%Δ (95% CI 19.7, 44.1]; white vs. white-to-olive 9.2%Δ (95% CI 2.8, 16.0)]. CONCLUSIONS: In this large cohort study, we confirmed known and described new determinants of facial telangiectasia.

Nasal and skin microbiomes are associated with disease severity in paediatric atopic dermatitis.

BACKGROUND: Alterations of the skin microbiome have been associated with atopic dermatitis (AD) and its severity. The nasal microbiome in relation to AD severity is less well studied. OBJECTIVES: We aimed to characterize the nasal and skin microbiomes in children with AD in relation to disease severity. In addition, we explored the differences and correlations between the nasal and skin communities. METHODS: We characterized the microbial composition of 90 nasal and 108 lesional skin samples cross-sectionally from patients with AD, using 16S-rRNA sequencing. In addition, a quantitative polymerase chain reaction was performed for Staphylococcus aureus and Staphylococcus epidermidis on the skin samples, and AD severity was estimated using the self-administered Eczema Area and Severity Index. RESULTS: We found an association between the microbial composition and AD severity in both the nose and skin samples (R2  = 2·6%; P = 0·017 and R2  = 7·0%; P = 0·004), strongly driven by staphylococci. However, other species also contributed, such as Moraxella in the nose. Skin lesions were positive for S. aureus in 50% of the children, and the presence and the load of S. aureus were not associated with AD severity. Although the nose and skin harbour distinct microbial communities (n = 48 paired samples; P < 0·001), we found that correlations exist between species in the nose and (other) species on the skin. CONCLUSIONS: Our results indicate that both the nasal and the skin microbiomes are associated with AD severity in children and that, next to staphylococci, other species contribute to this association.

Non-genetic and genetic predictors of a superficial first basal cell carcinoma.

BACKGROUND: Several observational studies have suggested differences in the risk factor profile between patients with superficial basal cell carcinomas (BCCs) and non-superficial BCCs. OBJECTIVE: To test the reproducibility of previous study findings and to find new genetic and non-genetic predictors for patients with a superficial first BCC. METHODS: A total of 14.628 participants of northwestern European descent aged 45 years or older from a prospective population-based cohort study (Rotterdam Study) were linked with the Dutch Pathology Registry (PALGA) of whom 1528 were identified as BCC patients. After exclusion, 948 eligible BCC patients remained for further non-genetic analyses and 1014 for genetic analyses. We included 11 phenotypic, environmental and tumour-specific characteristics, and 20 candidate single nucleotide polymorphisms (SNP) as potential predictors for patients with a superficial first BCC. We performed binary logistic multivariable regression analyses. RESULTS: We found that patients with a superficial first BCC were significantly younger, almost two times more often female and 12-18 times more likely to have their BCC on the trunk or extremities than patients with a non-superficial first BCC. One SNP (rs12203592), mapped to IRF4, looked promising (OR 1.83, 95% CI 1.13-2.97, P-value <0.05), but after adjustment for multiple testing, no significant differences in genetic make-up between superficial BCC and non-superficial BCC patients were found. CONCLUSION: We conclude that patients with a superficial BCC differ from non-superficial BCC patients with respect to environmental factors (tumour localization as a proxy for UVR exposure) and phenotypic characteristics (age and sex), but we found no difference in genotype. As superficial BCC patients develop their first BCCs at a younger age, they could be at higher lifetime risk for subsequent skin cancers and therefore be an important group for secondary prevention.

Prevalence and determinants of seborrhoeic dermatitis in a middle-aged and elderly population: the Rotterdam Study.

BACKGROUND: Seborrhoeic dermatitis is a chronic relapsing inflammatory skin disease with unclear pathophysiological mechanisms. OBJECTIVES: To establish which lifestyle and physiological determinants are associated with seborrhoeic dermatitis. METHODS: Seborrhoeic dermatitis was diagnosed by a trained physician during a full-body skin examination within the Rotterdam Study, a prospective population-based cohort study in middle-aged and elderly people. The current design is a comparative cross-sectional study embedded in the Rotterdam Study. Potential factors were identified from the literature and analysed in a multivariable logistic regression, including: age, sex, obesity, skin colour, stress, depression, education level, hypertension, climate, xerosis cutis, alcohol and tobacco use. RESULTS: Of the 5498 participants, 788 participants were diagnosed with seborrhoeic dermatitis (14·3%). We found associations between seborrhoeic dermatitis and male sex [adjusted odds ratio (OR) 2·09, 95% confidence interval (CI) 1·77-2·47], darker skin (adjusted OR 0·39, 95% CI 0·22-0·69), season (summer vs. winter: adjusted OR 0·63, 95% CI 0·48-0·82) and generalized xerosis cutis (adjusted OR 1·41, 95% CI 1·11-1·80). CONCLUSIONS: Seborrhoeic dermatitis is one of the most common inflammatory dermatoses in middle-aged and elderly individuals, especially during winter. Men, and people with a light and dry skin were most likely to have seborrhoeic dermatitis.

IgG response against Staphylococcus aureus is associated with severe atopic dermatitis in children.

BACKGROUND: An altered immune response against Staphylococcus aureus might contribute to inflammation and barrier damage in atopic dermatitis (AD). OBJECTIVES: To profile IgG antibodies against 55 S. aureus antigens in sera of children with mild-to-severe AD and to evaluate the association between IgG levels and disease severity. METHODS: In this cross-sectional study, we included children with AD from two interventional study cohorts, the Shared Medical Appointment (SMA) cohort (n = 131) and the older DAVOS cohort (n = 76). AD severity was assessed using the Self-Administered Eczema Area and Severity Index (SA-EASI) and levels of thymus and activation-regulated chemokine (TARC) in serum. IgG antibody levels against 55 S. aureus antigens were quantified simultaneously using a Luminex assay. Pair-wise correlations were calculated between the 55 IgG levels using the Spearman rank correlation test. Linear regression analysis was performed to test for associations between 55 IgG levels and SA-EASI and TARC, adjusting for age, sex and S. aureus colonization. RESULTS: In the SMA cohort, 16 antigens were associated with SA-EASI and 12 with TARC (10 overlapping antigens; P-values 0·001-0·044). The associated IgG antibodies targeted mainly secreted proteins with immunomodulatory functions. In the DAVOS study, IgG levels against only four and one S. aureus antigen(s) were associated with SA-EASI and TARC, respectively (no overlap). CONCLUSIONS: In young children, severity of AD is associated with an IgG response directed against S. aureus antigens with mainly immunomodulatory functions. These findings encourage further evaluation of the role of S. aureus in the pathogenesis of AD.

Incidence and relative survival of melanoma in children and adolescents in the Netherlands, 1989-2013.

BACKGROUND: Melanoma is rare in the first two decades of life. Trends in incidence differ across countries. OBJECTIVE: To describe incidence and relative survival of children and adolescents with melanoma in the Netherlands for children (0 through 11 years) and adolescents (12 through 19 years) separately. We hypothesized that adolescent melanoma increased in contrast to childhood melanoma, possibly due to a difference in cancer biology and sun exposure patterns. METHODS: Data on all patients of 0-19 years diagnosed between 1989 and 2013 with histologically confirmed cutaneous invasive melanoma were retrieved from the Netherlands Cancer Registry (NCR). Incidence trends were analysed with Joinpoint regression. Relative survival analysis was performed. RESULTS: Between 1989 and 2013, 80 children and 544 adolescents with melanoma were registered in the NCR. Median age at diagnosis was 17 years (IQR 15-18); the female-to-male ratio was 1.7 : 1 Statistically significant incidence trends were found in the older age group (12-19 years): an increasing incidence since 1991 [annual percentage change (APC) 3.2%, 95%CI 1.3-5.1] followed by a decrease from 2005 to 2013 (APC -4.9%, 95%CI -9.6-0.0). No incidence trends for childhood melanoma were observed (APC 0.3%, 95% CI -3.0-3.8). Relative survival at 1, 5 and 10 years was 98% (95% CI 97-99), 94% (95% CI 92-96) and 90% (95% CI 87-92), respectively. Survival was worse in males and higher Breslow thickness. CONCLUSIONS: Melanoma is very rare under the age of 12 with stable incidence rates. In comparison with childhood melanoma, melanomas in adolescents are more common with a decreasing trend in the past decade. Male sex and increasing Breslow thickness are associated with worse survival in paediatric melanoma patients.

Energetic trade-offs: Implications for selection between two bivalve prey species by a carnivorous muricid gastropod.

Active predators obtain energy and nutrients from prey through complex processes in which the energy gained must exceed the energy invested in finding and ingesting the prey. In addition, the amount of energy available will vary with the prey that are selected for consumption. The muricid gastropod Acanthina monodon inhabits rocky shores, where it routinely feeds on the mytilids Semimytilus algosus and Perumytilus purpuratus. In this study, S. algosus was highly preferred by the predator (over 90% were eaten) versus P. purpuratus (only 9% were eaten) when offered a mixed diet. The energetic cost of attacking one S. algosus individual was 91 J bivalve-1 while for P. purpuratus it was slightly higher: 95 J bivalve-1. Also, whereas A. monodon required on average 19 h to consume S. algosus, successful attacks on P. purpuratus required about 32% more time (25 h). In addition, a longer resting time was needed by the predator after preying on P. purpuratus before it initiated another attack. Moreover, the active metabolic costs associated with successfully attacking the prey increased 3.2 times over the basal metabolic costs when attacking S. algosus, but only by 2.5 times when attacking P. purpuratus. The calculations associated with preying on each species showed that the energetic gain per unit time likely accounts for the predator's preference for attacking S. algosus, even though predation on both species provided net energy gains for the predator. However, as S. algosus occurs seasonally at our study site, P. purpuratus would probably also be consumed due to its constant availability throughout the whole year.

Sex-specific genetic effects influence variation in body composition.

AIMS/HYPOTHESIS: Despite well-known sex differences in body composition it is not known whether sex-specific genetic or environmental effects contribute to these differences. METHODS: We assessed body composition in 2,506 individuals, from a young Dutch genetic isolate participating in the Erasmus Rucphen Family study, by dual-energy X-ray absorptiometry and anthropometry. We used variance decomposition procedures to partition variation of body composition into genetic and environmental components common to both sexes and to men and women separately and calculated the correlation between genetic components in men and women. RESULTS: After accounting for age, sex and inbreeding, heritability ranged from 0.39 for fat mass index to 0.84 for height. We found sex-specific genetic effects for fat percentage (fat%), lean mass, lean mass index (LMI) and fat distribution, but not for BMI and height. Genetic correlations between sexes were significantly different from 1 for fat%, lean mass, LMI, android fat, android:gynoid fat ratio and WHR, indicating that there are sex-specific genes contributing to variation of these traits. Genetic variance was significantly higher in women for the waist, hip and thigh circumference and WHR, implying that genes account for more variance of fat distribution in women than in men. Environmental variance was significantly higher in men for the android:gynoid fat ratio. CONCLUSIONS/INTERPRETATION: Sex-specific genetic effects underlie sexual dimorphism in several body composition traits. The findings are relevant for studies on the relationship of body composition with common diseases like cardiovascular disease and type 2 diabetes and for genetic association studies.

Biochemical responses and physiological status in the crab Hemigrapsus crenulatus (Crustacea, Varunidae) from high anthropogenically-impacted estuary (Lenga, south-central Chile).

Estuarine environmental assessment by sub-individual responses is important in order to understand contaminant effects and to find suitable estuarine biomonitor species. Our study aimed to analyze oxidative stress responses, including glutathione-S-transferase (GST) activity, total antioxidant capacity (ACAP) and lipid peroxidation levels (TBARS) in estuarine crabs Hemigrapsus crenulatus from a high anthropogenically-impacted estuary (Lenga) compared to low and non-polluted estuaries (Tubul and Raqui), in a seasonal scale (winter-summer), tissue specific (hepatopancreas and gills) and sex related responses. Results showed that hepatopancreas in male crabs better reflected inter-estuary differences. Morpho-condition traits as Cephalothorax hepatopancreas index (CHI) could be used as an indicator of physiological status of estuarine crabs. Discriminant analysis also showed that GST and TBARS levels in summer are more suitable endpoints for establishing differences between polluted and non-polluted sites. These results suggest the importance of seasonality, target tissue, sex and physiological status of brachyuran crabs for estuarine biomonitoring assessment.

The effect of genetic drift in a young genetically isolated population.

The genetic make-up of genetically isolated populations may differ from a general population as a result of genetic drift and founder effects. We assessed the extent of this deviation in a recently isolated population located in the southwest of the Netherlands and studied as part of the Genetic Research in Isolated Population (GRIP) program. A gene-dropping experiment was performed in a large pedigree from this isolate, assuming different initial frequencies in the population founders came from. Allelic frequencies in the last generations of this pedigree were estimated. Simulation analysis showed large fluctuations, as measured by variation coefficient and sufficient loss probability, when initial frequencies were lower than or equal to 1%. For initial frequencies larger than 1% the fluctuations were small. We also analyzed mean heterozygosity and allele diversity of 592 markers in a random sample from the GRIP population. The results were compared with a general population (CEPH sample), old large isolate (Icelandic sample) and the small-sized population of Talana (Sardinia). GRIP mean heterozygosity and mean number of alleles were significantly lower as compared with CEPH and Iceland, but much higher when compared with the Talana population. We also concluded that the findings from the GRIP population for common variants (>1%) are likely to be extendable to other young isolates in Europe as well as to outbred populations.