The human immune system displays substantial variation between individuals, leading to differences in susceptibility to autoimmune disease. We present single-cell RNA sequencing (scRNA-seq) data from ...1,267,758 peripheral blood mononuclear cells from 982 healthy human subjects. For 14 cell types, we identified 26,597 independent cis-expression quantitative trait loci (eQTLs) and 990 trans-eQTLs, with most showing cell type-specific effects on gene expression. We subsequently show how eQTLs have dynamic allelic effects in B cells that are transitioning from naïve to memory states and demonstrate how commonly segregating alleles lead to interindividual variation in immune function. Finally, using a Mendelian randomization approach, we identify the causal route by which 305 risk loci contribute to autoimmune disease at the cellular level. This work brings together genetic epidemiology with scRNA-seq to uncover drivers of interindividual variation in the immune system.
Myopia (near-sightedness) is an important public health issue. Spending more time outdoors can prevent myopia but the long-term association between this exposure and myopia has not been well ...characterised. We investigated the relationship between time spent outdoors in childhood, adolescence and young adulthood and risk of myopia in young adulthood. The Kidskin Young Adult Myopia Study (KYAMS) was a follow-up of the Kidskin Study, a sun exposure-intervention study of 1776 children aged 6-12 years. Myopia status was assessed in 303 (17.6%) KYAMS participants (aged 25-30 years) and several subjective and objective measures of time spent outdoors were collected in childhood (8-12 years) and adulthood. Index measures of total, childhood and recent time spent outdoors were developed using confirmatory factor analysis. Logistic regression was used to assess the association between a 0.1-unit change in the time outdoor indices and risk of myopia after adjusting for sex, education, outdoor occupation, parental myopia, parental education, ancestry and Kidskin Study intervention group. Spending more time outdoors during childhood was associated with reduced risk of myopia in young adulthood (multivariable odds ratio OR 0.82, 95% confidence interval CI 0.69, 0.98). Spending more time outdoors in later adolescence and young adulthood was associated with reduced risk of late-onset myopia (≥ 15 years of age, multivariable OR 0.79, 95% CI 0.64, 0.98). Spending more time outdoors in both childhood and adolescence was associated with less myopia in young adulthood.
Using latent variables in gene expression data can help correct unobserved confounders and increase statistical power for expression quantitative trait Loci (eQTL) detection. The probabilistic ...estimation of expression residuals (PEER) and principal component analysis (PCA) are widely used methods that can remove unwanted variation and improve eQTL discovery power in bulk RNA-seq analysis. However, their performance has not been evaluated extensively in single-cell eQTL analysis, especially for different cell types. Potential challenges arise due to the structure of single-cell RNA-seq data, including sparsity, skewness, and mean-variance relationship. Here, we show by a series of analyses that PEER and PCA require additional quality control and data transformation steps on the pseudo-bulk matrix to obtain valid latent variables; otherwise, it can result in highly correlated factors (Pearson's correlation r = 0.63 ~ 0.99). Incorporating valid PFs/PCs in the eQTL association model would identify 1.7 ~ 13.3% more eGenes. Sensitivity analysis showed that the pattern of change between the number of eGenes detected and fitted PFs/PCs varied significantly in different cell types. In addition, using highly variable genes to generate latent variables could achieve similar eGenes discovery power as using all genes but save considerable computational resources (~ 6.2-fold faster).
Inflammation and cigarette smoking predispose to macular diseases, and choroidal and retinal thinning. We explored the choroidal and retinal thicknesses in young adults against their 7-year ...C-reactive protein (CRP) level trajectory and pack-years smoked. Participants from the Raine study, a longitudinal cohort study, had serum CRP levels analysed at the 14-, 17-, and 20-year follow-ups. Group-based trajectory modelling was used to classify participants according to their 7-year CRP levels. At the 20-year follow-up (at 18-22 years old), participants completed questionnaires on their smoking history, and underwent optical coherence tomography imaging to obtain their choroidal and retinal thicknesses at the macula. Three CRP trajectories were identified: consistently low CRP levels (78% of sample), increasing (11%), or consistently high (11%). 340 and 1035 participants were included in the choroidal and retinal thickness analyses, respectively. Compared to those in the "Low" trajectory group, participants in the "Increasing" and "High" groups had 14-21 μm thinner choroids at most macular regions. Every additional pack-year smoked was linked with a 0.06-0.10 μm thinner retina at the inner and outer macular rings, suggesting a dose-dependent relationship between smoking and thinner retinas. These associations may suggest that an increased risk of future visual impairment or eye disease associated with these risk factors may be present since young adulthood.
The total number of acquired melanocytic nevi on the skin is strongly correlated with melanoma risk. Here we report a meta-analysis of 11 nevus GWAS from Australia, Netherlands, UK, and USA ...comprising 52,506 individuals. We confirm known loci including MTAP, PLA2G6, and IRF4, and detect novel SNPs in KITLG and a region of 9q32. In a bivariate analysis combining the nevus results with a recent melanoma GWAS meta-analysis (12,874 cases, 23,203 controls), SNPs near GPRC5A, CYP1B1, PPARGC1B, HDAC4, FAM208B, DOCK8, and SYNE2 reached global significance, and other loci, including MIR146A and OBFC1, reached a suggestive level. Overall, we conclude that most nevus genes affect melanoma risk (KITLG an exception), while many melanoma risk loci do not alter nevus count. For example, variants in TERC and OBFC1 affect both traits, but other telomere length maintenance genes seem to affect melanoma risk only. Our findings implicate multiple pathways in nevogenesis.
Changes in retinal thickness are common in various ocular diseases. Transverse magnification due to differing ocular biometrics, in particular axial length, affects measurement of retinal thickness ...in different regions. This study evaluated the effect of axial length and refractive error on measured macular thickness in two community-based cohorts of healthy young adults.
A total of 2160 eyes of 1247 community-based participants (18-30 years; 23.4% myopes, mean axial length = 23.6mm) were included in this analysis. Macular thickness measurements were obtained using a spectral-domain optical coherence tomography (which assumes an axial length of 24.385mm). Using a custom program, retinal thickness data were extracted at the 9 Early Treatment of Diabetic Retinopathy Study (ETDRS) regions with and without correction for transverse magnificent effects, with the corrected measurements adjusting according to the participant's axial length. Linear mixed models were used to analyse the effect of correction and its interaction with axial length or refractive group on retinal thickness.
The raw measures (uncorrected for axial length) underestimated the true retinal thickness at the central macula, while overestimating at most non-central macular regions. There was an axial length by correction interaction effect in all but the nasal regions (all p<0.05). For each 1mm increase in axial length, the central macular thickness is overestimated by 2.7-2.9μm while thicknesses at other regions were underestimated by 0.2-4.1μm. Based on the raw thickness measurements, myopes have thinner retinas than non-myopes at most non-central macular. However, this difference was no longer significant when the corrected data was used.
In a community-based sample, the raw measurements underestimate the retinal thickness at the central macula and overestimate the retinal thickness at non-central regions of the ETDRS grid. The effect of axial length and refractive error on retinal thickness is reduced after correcting for transverse magnification effects resulting from axial length differences.
To characterize and quantify Bruch membrane opening (BMO)-based optic nerve head (ONH) parameters in a large, young and healthy, predominantly white population.
Cross-sectional study and reliability ...analysis.
The ONH of 1344 predominantly white subjects were imaged with spectral-domain optical coherence tomography (SD-OCT). A customized script, coded in Matlab, was used to manually segment and measure multiple BMO-based parameters of the ONH. Measurements were compared to those obtained with confocal scanning laser ophthalmoscopy (Heidelberg Retina Tomograph; HRT). Regression analysis was performed to assess the relationship between BMO parameters and other ocular and demographic variables.
Mean BMO disc and neuroretinal rim (NRR) areas ranged from 0.94 to 4.06 mm2 (mean 1.77 ± 0.38 mm2) and 0.94 to 3.99 mm2 (mean 1.56 ± 0.31 mm2), respectively. When compared to the equivalent HRT measurements, SD-OCT-derived measures differed significantly for all comparable ONH parameters (P < .001). The limits of agreement computed from Bland-Altman plots comparing SD-OCT and HRT measurements showed suboptimal agreement between the techniques. Linear regression analysis demonstrated an effect of ethnicity, axial length, and refractive error on BMO-based parameters.
We have quantified BMO-based parameters in a large cohort of young adults using SD-OCT. These data will be informative in constructing normative profiles for clinical and research purposes in glaucoma diagnosis and management.
SIX1/SIX6 polymorphism has been shown to be associated with glaucoma. Studies have also found that, in older adults, retinal nerve fibre layer (RNFL) thickness is significantly thinned with each copy ...of the risk allele in SIX1/SIX6. However, it is not known whether these genetic variants exert their effects in younger individuals. Comparing a healthy young adult with an older adult cohort (mean age 20 vs 63 years), both of Northern European descent, we found that there was no significant RNFL thinning in each copy of the risk alleles in SIX1/SIX6 in the eyes of younger individuals. The older cohort showed an unexpectedly thicker RNFL in the nasal sector with each copy of the risk allele for both the SIX1 (rs10483727) and SIX6 (rs33912345) variants. In the temporal sector, thinner RNFL was found with each copy of the risk allele in rs33912345 with a decrease trend observed in rs10483727. Our results suggest that SIX1/SIX6 gene variants exert their influence later in adult life.
Purpose
Astigmatism is a common cause of refractive error and is known to vary in prevalence with age. Although the search for genes associated with spherical refractive errors (especially myopia) ...has met with limited success, current efforts to identify genetic variants implicated in astigmatism development have been less rewarding. We aimed to assess the association between astigmatism and age to identify appropriate age cut‐offs for maximizing power in genetic studies of astigmatism.
Methods
We performed a cross‐sectional analysis of right eye astigmatism data from four Australian‐based eye studies comprising 3841 participants aged 5–90 years. Measurements were performed under cycloplegia using an autorefractor, and individuals with a history of cataract, refractive surgery or corneal pathology were excluded from the analysis. In addition to the magnitude and type (against‐the‐rule, with‐the‐rule, and oblique) of astigmatism, we calculated the vector components (J0, J45) and evaluated the association of these outcome measures with age.
Results
The magnitude of refractive astigmatism (RA) remained relatively stable mean ± SD (−0.44 D ± 0.50) until individuals reached the age of 50, thereafter increasing in average magnitude by approximately 1.00 D for those subjects aged 90. In contrast, corneal astigmatism (CA) remained relatively stable from childhood until the age of 80 (−0.76 D ± 0.61). The prevalence of clinically significant RA (≥1.00 D) increased with age and was highest in those aged >70 years 55.1% (47.2–62.7%). Age was significantly associated with RA in adults odds ratio (OR) = 1.04 per 1 year, p < 0.001. A weaker relationship was observed between CA and age (OR = 1.007 per 1 year, p = 0.02).
Conclusions
We have confirmed the previously documented association between RA and age. Our results indicate that most of the observed change occurs after the age of 50, providing a recommended cut‐off for participants in genetic studies of this refractive condition.