Decoding depression by exploring the exposome-genome edge amidst COVID-19 lockdown.

Risk of depression increased in the general population after the COVID-19 pandemic outbreak. By examining the interplay between genetics and individual environmental exposures during the COVID-19 lockdown, we have been able to gain an insight as to why some individuals are more vulnerable to depression, while others are more resilient. This study, conducted on a Spanish cohort of 9218 individuals (COVICAT), includes a comprehensive non-genetic risk analysis, the exposome, complemented by a genomics analysis in a subset of 2442 participants. Depression levels were evaluated using the Hospital Anxiety and Depression Scale. Together with Polygenic Risk Scores (PRS), we introduced a novel score; Poly-Environmental Risk Scores (PERS) for non-genetic risks to estimate the effect of each cumulative score and gene-environment interaction. We found significant positive associations for PERSSoc (Social and Household), PERSLife (Lifestyle and Behaviour), and PERSEnv (Wider Environment and Health) scores across all levels of depression severity, and for PRSB (Broad depression) only for moderate depression (OR 1.2, 95% CI 1.03-1.40). On average OR increased 1.2-fold for PERSEnv and 1.6-fold for PERLife and PERSoc from mild to severe depression level. The complete adjusted model explained 16.9% of the variance. We further observed an interaction between PERSEnv and PRSB showing a potential mitigating effect. In summary, stressors within the social and behavioral domains emerged as the primary drivers of depression risk in this population, unveiling a mitigating interaction effect that should be interpreted with caution.

Biological basis of extensive pleiotropy between blood traits and cancer risk.

BACKGROUND: The immune system has a central role in preventing carcinogenesis. Alteration of systemic immune cell levels may increase cancer risk. However, the extent to which common genetic variation influences blood traits and cancer risk remains largely undetermined. Here, we identify pleiotropic variants and predict their underlying molecular and cellular alterations. METHODS: Multivariate Cox regression was used to evaluate associations between blood traits and cancer diagnosis in cases in the UK Biobank. Shared genetic variants were identified from the summary statistics of the genome-wide association studies of 27 blood traits and 27 cancer types and subtypes, applying the conditional/conjunctional false-discovery rate approach. Analysis of genomic positions, expression quantitative trait loci, enhancers, regulatory marks, functionally defined gene sets, and bulk- and single-cell expression profiles predicted the biological impact of pleiotropic variants. Plasma small RNAs were sequenced to assess association with cancer diagnosis. RESULTS: The study identified 4093 common genetic variants, involving 1248 gene loci, that contributed to blood-cancer pleiotropism. Genomic hotspots of pleiotropism include chromosomal regions 5p15-TERT and 6p21-HLA. Genes whose products are involved in regulating telomere length are found to be enriched in pleiotropic variants. Pleiotropic gene candidates are frequently linked to transcriptional programs that regulate hematopoiesis and define progenitor cell states of immune system development. Perturbation of the myeloid lineage is indicated by pleiotropic associations with defined master regulators and cell alterations. Eosinophil count is inversely associated with cancer risk. A high frequency of pleiotropic associations is also centered on the regulation of small noncoding Y-RNAs. Predicted pleiotropic Y-RNAs show specific regulatory marks and are overabundant in the normal tissue and blood of cancer patients. Analysis of plasma small RNAs in women who developed breast cancer indicates there is an overabundance of Y-RNA preceding neoplasm diagnosis. CONCLUSIONS: This study reveals extensive pleiotropism between blood traits and cancer risk. Pleiotropism is linked to factors and processes involved in hematopoietic development and immune system function, including components of the major histocompatibility complexes, and regulators of telomere length and myeloid lineage. Deregulation of Y-RNAs is also associated with pleiotropism. Overexpression of these elements might indicate increased cancer risk.

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.

SARS-CoV-2 infection, vaccination, and antibody response trajectories in adults: a cohort study in Catalonia.

BACKGROUND: Heterogeneity of the population in relation to infection, COVID-19 vaccination, and host characteristics is likely reflected in the underlying SARS-CoV-2 antibody responses. METHODS: We measured IgM, IgA, and IgG levels against SARS-CoV-2 spike and nucleocapsid antigens in 1076 adults of a cohort study in Catalonia between June and November 2020 and a second time between May and July 2021. Questionnaire data and electronic health records on vaccination and COVID-19 testing were available in both periods. Data on several lifestyle, health-related, and sociodemographic characteristics were also available. RESULTS: Antibody seroreversion occurred in 35.8% of the 64 participants non-vaccinated and infected almost a year ago and was related to asymptomatic infection, age above 60 years, and smoking. Moreover, the analysis on kinetics revealed that among all responses, IgG RBD, IgA RBD, and IgG S2 decreased less within 1 year after infection. Among vaccinated, 2.1% did not present antibodies at the time of testing and approximately 1% had breakthrough infections post-vaccination. In the post-vaccination era, IgM responses and those against nucleoprotein were much less prevalent. In previously infected individuals, vaccination boosted the immune response and there was a slight but statistically significant increase in responses after a 2nd compared to the 1st dose. Infected vaccinated participants had superior antibody levels across time compared to naïve-vaccinated people. mRNA vaccines and, particularly the Spikevax, induced higher antibodies after 1st and 2nd doses compared to Vaxzevria or Janssen COVID-19 vaccines. In multivariable regression analyses, antibody responses after vaccination were predicted by the type of vaccine, infection age, sex, smoking, and mental and cardiovascular diseases. CONCLUSIONS: Our data support that infected people would benefit from vaccination. Results also indicate that hybrid immunity results in superior antibody responses and infection-naïve people would need a booster dose earlier than previously infected people. Mental diseases are associated with less efficient responses to vaccination.

A Polygenic Risk Score Based on a Cardioembolic Stroke Multitrait Analysis Improves a Clinical Prediction Model for This Stroke Subtype.

Background: Occult atrial fibrillation (AF) is one of the major causes of embolic stroke of undetermined source (ESUS). Knowing the underlying etiology of an ESUS will reduce stroke recurrence and/or unnecessary use of anticoagulants. Understanding cardioembolic strokes (CES), whose main cause is AF, will provide tools to select patients who would benefit from anticoagulants among those with ESUS or AF. We aimed to discover novel loci associated with CES and create a polygenetic risk score (PRS) for a more efficient CES risk stratification. Methods: Multitrait analysis of GWAS (MTAG) was performed with MEGASTROKE-CES cohort (n = 362,661) and AF cohort (n = 1,030,836). We considered significant variants and replicated those variants with MTAG p-value < 5 × 10-8 influencing both traits (GWAS-pairwise) with a p-value < 0.05 in the original GWAS and in an independent cohort (n = 9,105). The PRS was created with PRSice-2 and evaluated in the independent cohort. Results: We found and replicated eleven loci associated with CES. Eight were novel loci. Seven of them had been previously associated with AF, namely, CAV1, ESR2, GORAB, IGF1R, NEURL1, WIPF1, and ZEB2. KIAA1755 locus had never been associated with CES/AF, leading its index variant to a missense change (R1045W). The PRS generated has been significantly associated with CES improving discrimination and patient reclassification of a model with age, sex, and hypertension. Conclusion: The loci found significantly associated with CES in the MTAG, together with the creation of a PRS that improves the predictive clinical models of CES, might help guide future clinical trials of anticoagulant therapy in patients with ESUS or AF.

GCAT|Panel, a comprehensive structural variant haplotype map of the Iberian population from high-coverage whole-genome sequencing.

The combined analysis of haplotype panels with phenotype clinical cohorts is a common approach to explore the genetic architecture of human diseases. However, genetic studies are mainly based on single nucleotide variants (SNVs) and small insertions and deletions (indels). Here, we contribute to fill this gap by generating a dense haplotype map focused on the identification, characterization, and phasing of structural variants (SVs). By integrating multiple variant identification methods and Logistic Regression Models (LRMs), we present a catalogue of 35 431 441 variants, including 89 178 SVs (≥50 bp), 30 325 064 SNVs and 5 017 199 indels, across 785 Illumina high coverage (30x) whole-genomes from the Iberian GCAT Cohort, containing a median of 3.52M SNVs, 606 336 indels and 6393 SVs per individual. The haplotype panel is able to impute up to 14 360 728 SNVs/indels and 23 179 SVs, showing a 2.7-fold increase for SVs compared with available genetic variation panels. The value of this panel for SVs analysis is shown through an imputed rare Alu element located in a new locus associated with Mononeuritis of lower limb, a rare neuromuscular disease. This study represents the first deep characterization of genetic variation within the Iberian population and the first operational haplotype panel to systematically include the SVs into genome-wide genetic studies.

Assessment of kinship detection using RNA-seq data.

Analysis of RNA sequencing (RNA-seq) data from related individuals is widely used in clinical and molecular genetics studies. Prediction of kinship from RNA-seq data would be useful for confirming the expected relationships in family based studies and for highlighting samples from related individuals in case-control or population based studies. Currently, reconstruction of pedigrees is largely based on SNPs or microsatellites, obtained from genotyping arrays, whole genome sequencing and whole exome sequencing. Potential problems with using RNA-seq data for kinship detection are the low proportion of the genome that it covers, the highly skewed coverage of exons of different genes depending on expression level and allele-specific expression. In this study we assess the use of RNA-seq data to detect kinship between individuals, through pairwise identity by descent (IBD) estimates. First, we obtained high quality SNPs after successive filters to minimize the effects due to allelic imbalance as well as errors in sequencing, mapping and genotyping. Then, we used these SNPs to calculate pairwise IBD estimates. By analysing both real and simulated RNA-seq data we show that it is possible to identify up to second degree relationships using RNA-seq data of even low to moderate sequencing depth.