Evaluation of genetic variants associated with the risk of thiopurine-related pancreatitis: a case control study from ENEIDA registry

BACKGROUND

Risk factors for developing pancreatitis due to thiopurines in patients with inflammatory bowel disease (IBD) are not clearly identified.

AIM

Our aim was to evaluate the predictive pharmacogenetic risk of pancreatitis in IBD patients treated with thiopurines.

METHODS

We conducted an observational pharmacogenetic study of acute pancreatitis events in a cohort study of IBD patients treated with thiopurines from the prospectively maintained ENEIDA registry biobank of GETECCU. Samples were obtained and the CASR, CEL, CFTR, CDLN2, CTRC, SPINK1, CPA1, and PRSS1 genes, selected based on their known association with pancreatitis, were fully sequenced.

RESULTS

Ninety-five cases and 105 controls were enrolled, 57% were women. Median age at pancreatitis diagnosis was 39 years. We identified 81 benign variants (50 in cases and 67 in controls) and a total of 35 distinct rare pathogenic and unknown significance variants (10 in CEL, 21 in CFTR, 1 in CDLN2, and 3 in CPA1). None of the cases or controls carried pancreatitis-predisposing variants within the CASR, CPA1, PRSS1, and SPINK1 genes, nor a pathogenic CFTR mutation. Four different variants of unknown significance were detected in the CDLN and CPA1 genes; one of them was in the CDLN gene in a single patient with pancreatitis, and 3 in the CPA1 gene in 5 controls. After the analysis of the variants detected, no significant differences were observed between cases and controls.

CONCLUSION

In patients with IBD, genes known to cause pancreatitis seem not to be involved in thiopurine-related pancreatitis onset.

Towards precision medicine of long-acting aripiprazole through population pharmacokinetic modelling

Population pharmacokinetic (popPK) models constitute a valuable tool for characterizing the pharmacokinetic properties of once-monthly long-acting injectable aripiprazole (LAI aripiprazole) and quantifying the sources of variability in drug exposure. Our aim is to develop a popPK model of both aripiprazole and its metabolite dehydro-aripiprazole in patients treated with LAI aripiprazole, and to personalize the dosing regimen of aripiprazole across different sub-groups of patients. This is a prospective study investigating the pharmacokinetics of LAI aripiprazole. A total of 93 patients were included, 21 for model development and 71 for external model evaluation. A one-compartment model with linear absorption and elimination adequately described both aripiprazole and dehydro-aripiprazole concentrations. The weight of the patients has been shown to be the factor that most influences the absorption. However, the metabolizing phenotype for CYP2D6 and the concomitant treatment with strong inhibitors of this cytochrome have been shown to be the covariates that most influence total drug exposure. This is the first popPK model developed for LAI aripiprazole that includes aripiprazole and its main active metabolite, dehydroaripiprazole. It provides a personalized dosage recommendation that maximizes the probability of achieving optimal therapeutic concentrations and minimizes the difficulties associated with trial-and-error therapeutic strategies carried out in clinical practice.

Clonal chromosomal mosaicism and loss of chromosome Y in elderly men increase vulnerability for SARS-CoV-2

The pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2, COVID-19) had an estimated overall case fatality ratio of 1.38% (pre-vaccination), being 53% higher in males and increasing exponentially with age. Among 9578 individuals diagnosed with COVID-19 in the SCOURGE study, we found 133 cases (1.42%) with detectable clonal mosaicism for chromosome alterations (mCA) and 226 males (5.08%) with acquired loss of chromosome Y (LOY). Individuals with clonal mosaic events (mCA and/or LOY) showed a 54% increase in the risk of COVID-19 lethality. LOY is associated with transcriptomic biomarkers of immune dysfunction, pro-coagulation activity and cardiovascular risk. Interferon-induced genes involved in the initial immune response to SARS-CoV-2 are also down-regulated in LOY. Thus, mCA and LOY underlie at least part of the sex-biased severity and mortality of COVID-19 in aging patients. Given its potential therapeutic and prognostic relevance, evaluation of clonal mosaicism should be implemented as biomarker of COVID-19 severity in elderly people.

An Insertion Within SIRPβ1 Shows a Dual Effect Over Alzheimer's Disease Cognitive Decline Altering the Microglial Response

Background

Microglial dysfunction plays a causative role in Alzheimer's disease (AD) pathogenesis. Here we focus on a germline insertion/deletion variant mapping SIRPβ1, a surface receptor that triggers amyloid-β(Aβ) phagocytosis via TYROBP.

Objective

To analyze the impact of this copy-number variant in SIRPβ1 expression and how it affects AD molecular etiology.

Methods

Copy-number variant proxy rs2209313 was evaluated in GERALD and GR@ACE longitudinal series. Hippocampal specimens of genotyped AD patients were also examined. SIRPβ1 isoform-specific phagocytosis assays were performed in HEK393T cells.

Results

The insertion alters the SIRPβ1 protein isoform landscape compromising its ability to bind oligomeric Aβ and its affinity for TYROBP. SIRPβ1 Dup/Dup patients with mild cognitive impairment show an increased cerebrospinal fluid t-Tau/Aβ ratio (p = 0.018) and a higher risk to develop AD (OR = 1.678, p = 0.018). MRIs showed that Dup/Dup patients exhibited a worse initial response to AD. At the moment of diagnosis, all patients showed equivalent Mini-Mental State Examination scores. However, AD patients with the duplication had less hippocampal degeneration (p < 0.001) and fewer white matter hyperintensities. In contrast, longitudinal studies indicate that patients bearing the duplication allele show a slower cognitive decline (p = 0.013). Transcriptional analysis also shows that the SIRPβ1 duplication allele correlates with higher TREM2 expression and an increased microglial activation.

Conclusions

The SIRPβ1 internal duplication has opposite effects over MCI-to-Dementia conversion risk and AD progression, affecting microglial response to Aβ. Given the pharmacological approaches focused on the TREM2-TYROBP axis, we believe that SIRPβ1 structural variant might be considered as a potential modulator of this causative pathway.

Expanding the clinical and molecular spectrum of FOXG1- and ZBTB18-associated neurodevelopmental disorders

Introduction The zinc finger BTB domain-containing protein ZBTB18 binds to FOXG1 to form a transcriptional repressive complex involved in neuronal differentiation. Disruption of the components of this complex results in chromosome 1q43-q44 deletion syndrome/intellectual developmental disorder 22 or in FOXG1 syndrome. Case presentation This study reports on five patients with cognitive and behavioral impairment, seizures, microcephaly, and/or congenital brain abnormalities. Whole exome sequencing identified deleterious ZBTB18 variants in three patients and deleterious FOXG1 variants in the remaining patients. We have detected a missense variant within the BTB domain of ZBTB18 in two affected monozygotic twins. In addition, we observed agenesis of the septum pellucidum in a missense FOXG1 carrier with a severe FOXG1 syndrome. Conclusion Although the ZBTB18 zinc finger domains harbor the majority of known deleterious variants, we report a novel de novo rare missense variant within the BTB domain. The agenesis of the septum pellucidum observed in a missense FOXG1 carrier could be considered as a novel clinical feature associated with FOXG1 syndrome. The severe FOXG1 syndrome in this patient contrasts with the milder phenotype expected for a missense. Genetic or environmental factors may explain this phenotypic variability in FOXG1 syndrome.

HLA-A*11:01 and HLA-C*04:01 are associated with severe COVID-19

We analyzed the association between HLA polymorphisms and susceptibility to SARS-CoV-2 infection and disease severity. Genotyping data from a total of 9373 COVID-19-positive cases from the Spanish Coalition to Unlock Research on Host Genetics on COVID-19 (SCOURGE) consortium and 5943 population controls were included in the study. We found an association of the alleles HLA-B*14:02 and HLA-C*08:02 with a lower risk to COVID-19 infection (p = 0.006, OR = 0.84, 95% CI = [0.75-0.95], p = 0.024, OR = 0.86, 95% CI = [0.78-0.95], respectively). We also found the alleles HLA-A*11:01 and HLA-C*04:01 associated with disease severity (p = 0.033, OR = 1.16, 95% CI = [1.04-1.31], p = 0.045, OR = 1.14, 95% CI = [1.05-1.25], respectively). These results suggest that an effective presentation of viral peptides by HLA class I alleles involve a faster infection clearance, decreasing the susceptibility and severity of COVID-19.

A Genome-Wide Association Study of Small Cell Lung Cancer

INTRODUCTION

Small cell lung cancer (SCLC) comprises 10-15% of all lung cancer cases and is the most aggressive histological type. Survival is poor and the molecular landscape of this disease is extraordinarily complex. The objective of this paper was to perform a Genome-Wide Association Study (GWAS) of this disease using a case-control study specifically designed for small cell lung cancer (SCLC).

METHODS

Incident cases were consecutively recruited from 8 hospitals from different regions of Spain. Controls were recruited from the same hospitals using a frequency sampling based on age and sex distribution of cases. Biological samples were obtained along with detailed information on cases and controls lifestyle, including tobacco and radon exposure.

RESULTS

We included 271 SCLC cases and 557 controls. We found evidence (p-values<10-5) of an association in the complete dataset for several loci, while MAP4 showed a significant association in the gene-based analysis. Pathway analysis suggested that ATR, ATRIP, MCM4, MCM5, ORC4, RPA3 and CDC25A genes have a role on the onset of SCLC.

CONCLUSION

This study provides biological evidence for pathways related to SCLC, offering novel loci for further research.

HLA alleles: important pieces to the COVID-19 puzzle

Research on human leukocyte antigen (HLA) molecules in coronavirus disease 2019 (COVID-19) raised high expectations but has yielded limited results. Augusto et al.'s recent study in Nature unveils a strong association of HLA-B*15:01 with asymptomatic COVID-19, representing an important contribution to genetics in COVID-19.

Prenatal environmental exposures associated with sex differences in childhood obesity and neurodevelopment

BACKGROUND

Obesity and neurodevelopmental delay are complex traits that often co-occur and differ between boys and girls. Prenatal exposures are believed to influence children's obesity, but it is unknown whether exposures of pregnant mothers can confer a different risk of obesity between sexes, and whether they can affect neurodevelopment.

METHODS

We analyzed data from 1044 children from the HELIX project, comprising 93 exposures during pregnancy, and clinical, neuropsychological, and methylation data during childhood (5-11 years). Using exposome-wide interaction analyses, we identified prenatal exposures with the highest sexual dimorphism in obesity risk, which were used to create a multiexposure profile. We applied causal random forest to classify individuals into two environments: E1 and E0. E1 consists of a combination of exposure levels where girls have significantly less risk of obesity than boys, as compared to E0, which consists of the remaining combination of exposure levels. We investigated whether the association between sex and neurodevelopmental delay also differed between E0 and E1. We used methylation data to perform an epigenome-wide association study between the environments to see the effect of belonging to E1 or E0 at the molecular level.

RESULTS

We observed that E1 was defined by the combination of low dairy consumption, non-smokers' cotinine levels in blood, low facility richness, and the presence of green spaces during pregnancy (OR_interaction = 0.070, P = 2.59 × 10^-5). E1 was also associated with a lower risk of neurodevelopmental delay in girls, based on neuropsychological tests of non-verbal intelligence (OR_interaction = 0.42, P = 0.047) and working memory (OR_interaction = 0.31, P = 0.02). In line with this, several neurodevelopmental functions were enriched in significant differentially methylated probes between E1 and E0.

CONCLUSIONS

The risk of obesity can be different for boys and girls in certain prenatal environments. We identified an environment combining four exposure levels that protect girls from obesity and neurodevelopment delay. The combination of single exposures into multiexposure profiles using causal inference can help determine populations at risk.

Childhood exposure to non-persistent endocrine disrupting chemicals and multi-omic profiles: A panel study

BACKGROUND

Individuals are exposed to environmental pollutants with endocrine disrupting activity (endocrine disruptors, EDCs) and the early stages of life are particularly susceptible to these exposures. Previous studies have focused on identifying molecular signatures associated with EDCs, but none have used repeated sampling strategy and integrated multiple omics. We aimed to identify multi-omic signatures associated with childhood exposure to non-persistent EDCs.

METHODS

We used data from the HELIX Child Panel Study, which included 156 children aged 6 to 11. Children were followed for one week, in two time periods. Twenty-two non-persistent EDCs (10 phthalate, 7 phenol, and 5 organophosphate pesticide metabolites) were measured in two weekly pools of 15 urine samples each. Multi-omic profiles (methylome, serum and urinary metabolome, proteome) were measured in blood and in a pool urine samples. We developed visit-specific Gaussian Graphical Models based on pairwise partial correlations. The visit-specific networks were then merged to identify reproducible associations. Independent biological evidence was systematically sought to confirm some of these associations and assess their potential health implications.

RESULTS

950 reproducible associations were found among which 23 were direct associations between EDCs and omics. For 9 of them, we were able to find corroborating evidence from previous literature: DEP - serotonin, OXBE - cg27466129, OXBE - dimethylamine, triclosan - leptin, triclosan - serotonin, MBzP - Neu5AC, MEHP - cg20080548, oh-MiNP - kynurenine, oxo-MiNP - 5-oxoproline. We used these associations to explore possible mechanisms between EDCs and health outcomes, and found links to health outcomes for 3 analytes: serotonin and kynurenine in relation to neuro-behavioural development, and leptin in relation to obesity and insulin resistance.

CONCLUSIONS

This multi-omics network analysis at two time points identified biologically relevant molecular signatures related to non-persistent EDC exposure in childhood, suggesting pathways related to neurological and metabolic outcomes.

Colorectal cancer. Genetic variants in BMP signaling pathway and ancestry in the Mexican population

INTRODUCTION

Genetic variants related to bone morphogenetic proteins (BMP2, BMP4, GREM1, SMAD7) signaling pathway have been associated with colorectal cancer, mainly in Caucasian populations.

OBJECTIVE

To describe the association of variants in members of the BMP signaling pathway in a Mexican population, characterized by its indigenous American and Caucasian ancestry.

METHODS

Genotyping of 1,000 colorectal cancer cases and 1,043 control individuals recruited in Mexico City, Monterrey, and Torreón was carried out using the Sequenom platform. Associations between colorectal cancer and variants were studied with univariate and multivariate analyses.

RESULTS

Variants rs4444235, rs12953717 and rs4939827 replicated the association with the neoplasm (p ≤ 0.05). Caucasian ancestry showed association with the tumor.

CONCLUSIONS

The study replicated the associations between colorectal cancer and SMAD7 and BMP4 variants, with an association being observed with the Caucasian component of the ethnic mix.

Detection of the Copy Number Variants of Genes in Patients with Familial Cardiac Diseases by Massively Parallel Sequencing

INTRODUCTION

The implication of copy number variations in familial heart disease is known, although in-depth knowledge is lacking; hence, more studies are needed to further our understanding. Massively parallel sequencing, thanks to its recent surge in use, is emerging as a valid tool for the detection of this type of variant, through the use of appropriate software.

METHODS

We conducted a study with 182 patients diagnosed with mendelian cardiovascular diseases who underwent sequencing using a cardiac gene panel and then a specific calling process for copy number variations (CNVs) with ExomeDepth software, which provides us with a Bayes factor (BF), a score of the probability that a CNV detected is true.

RESULTS

After a rigorous CNV prioritization process, we confirmed the variants obtained by MLPA or SNP-based array, finding three real CNVs in five individuals in the MYH11, FBN1 and PDMI7 genes.

CONCLUSION

The confirmed CNVs present in all cases BF values > 60, thus establishing a threshold to consider real CNVs in the calling process carried out by ExomeDepth on our gene panel.

Multi-omics signatures of the human early life exposome

Environmental exposures during early life play a critical role in life-course health, yet the molecular phenotypes underlying environmental effects on health are poorly understood. In the Human Early Life Exposome (HELIX) project, a multi-centre cohort of 1301 mother-child pairs, we associate individual exposomes consisting of >100 chemical, outdoor, social and lifestyle exposures assessed in pregnancy and childhood, with multi-omics profiles (methylome, transcriptome, proteins and metabolites) in childhood. We identify 1170 associations, 249 in pregnancy and 921 in childhood, which reveal potential biological responses and sources of exposure. Pregnancy exposures, including maternal smoking, cadmium and molybdenum, are predominantly associated with child DNA methylation changes. In contrast, childhood exposures are associated with features across all omics layers, most frequently the serum metabolome, revealing signatures for diet, toxic chemical compounds, essential trace elements, and weather conditions, among others. Our comprehensive and unique resource of all associations ( https://helixomics.isglobal.org/ ) will serve to guide future investigation into the biological imprints of the early life exposome.

A genome-wide association study of survival in patients with sepsis

BACKGROUND

Sepsis is a severe systemic inflammatory response to infections that is accompanied by organ dysfunction and has a high mortality rate in adult intensive care units. Most genetic studies have identified gene variants associated with development and outcomes of sepsis focusing on biological candidates. We conducted the first genome-wide association study (GWAS) of 28-day survival in adult patients with sepsis.

METHODS

This study was conducted in two stages. The first stage was performed on 687 European sepsis patients from the GEN-SEP network and 7.5 million imputed variants. Association testing was conducted with Cox regression models, adjusting by sex, age, and the main principal components of genetic variation. A second stage focusing on the prioritized genetic variants was performed on 2,063 ICU sepsis patients (1362 European Americans and 701 African-Americans) from the MESSI study. A meta-analysis of results from the two stages was conducted and significance was established at p < 5.0 × 10-8. Whole-blood transcriptomic, functional annotations, and sensitivity analyses were evaluated on the identified genes and variants.

FINDINGS

We identified three independent low-frequency variants associated with reduced 28-day sepsis survival, including a missense variant in SAMD9 (hazard ratio [95% confidence interval] = 1.64 [1.37-6.78], p = 4.92 × 10-8). SAMD9 encodes a possible mediator of the inflammatory response to tissue injury.

INTERPRETATION

We performed the first GWAS of 28-day sepsis survival and identified novel variants associated with reduced survival. Larger sample size studies are needed to better assess the genetic effects in sepsis survival and to validate the findings.

Gravedad de COVID-19 atribuible a obesidad según IMC y CUN-BAE

Introducción

La obesidad es considerada un factor de riesgo en casos graves de la COVID-19, habiendo sido analizada mediante el índice de masa corporal (IMC), estimador que no correlaciona adecuadamente con el porcentaje de grasa corporal (GC). El objetivo de este estudio ha sido analizar la fracción atribuible poblacional a la GC en formas graves de COVID-19 atendiendo al IMC y al CUN-BAE.

Material y métodos

Estudio multicéntrico observacional de prevalencia. Se recogió información sociodemográfica, antecedentes personales, IMC y CUN-BAE, de casos positivos SARS-CoV-2, de las provincias de León y La Rioja. Mediante modelos de regresión logística se calcularon odds ratio con sus respectivos intervalos de confianza del 95% ajustando por edad y antecedentes personales, así como la fracción atribuible poblacional a la GC.

Resultados

Participaron 785 pacientes, 123 (15,7%) fueron graves. Se detectaron como factores de riesgo la edad, la obesidad (tanto por IMC como por CUN-BAE) y los antecedentes personales. Un 51,6% de casos graves podrían ser atribuidos a un exceso de IMC y un 61,4% a exceso de GC estimada según CUN-BAE, observándose una mayor infraestimación del riesgo en mujeres.

Conclusiones

El exceso de GC es un factor de riesgo para formas graves de la COVID-19 junto con la edad avanzada y la presencia de enfermedades cardiovasculares, respiratorias crónicas u oncohematológicas. El IMC infraestima el riesgo, especialmente en mujeres, siendo el CUN-BAE el predictor seleccionado por su mejor estimación del porcentaje de GC.

Novel genes and sex differences in COVID-19 severity

Here we describe the results of a genome-wide study conducted in 11 939 COVID-19 positive cases with an extensive clinical information that were recruited from 34 hospitals across Spain (SCOURGE consortium). In sex-disaggregated genome-wide association studies for COVID-19 hospitalization, genome-wide significance (p < 5x10-8) was crossed for variants in 3p21.31 and 21q22.11 loci only among males (p = 1.3x10-22 and p = 8.1x10-12, respectively), and for variants in 9q21.32 near TLE1 only among females (p = 4.4x10-8). In a second phase, results were combined with an independent Spanish cohort (1598 COVID-19 cases and 1068 population controls), revealing in the overall analysis two novel risk loci in 9p13.3 and 19q13.12, with fine-mapping prioritized variants functionally associated with AQP3 (p = 2.7x10-8) and ARHGAP33 (p = 1.3x10-8), respectively. The meta-analysis of both phases with four European studies stratified by sex from the Host Genetics Initiative confirmed the association of the 3p21.31 and 21q22.11 loci predominantly in males and replicated a recently reported variant in 11p13 (ELF5, p = 4.1x10-8). Six of the COVID-19 HGI discovered loci were replicated and an HGI-based genetic risk score predicted the severity strata in SCOURGE. We also found more SNP-heritability and larger heritability differences by age (<60 or ≥ 60 years) among males than among females. Parallel genome-wide screening of inbreeding depression in SCOURGE also showed an effect of homozygosity in COVID-19 hospitalization and severity and this effect was stronger among older males. In summary, new candidate genes for COVID-19 severity and evidence supporting genetic disparities among sexes are provided.

The early-life exposome modulates the effect of polymorphic inversions on DNA methylation

Polymorphic genomic inversions are chromosomal variants with intrinsic variability that play important roles in evolution, environmental adaptation, and complex traits. We investigated the DNA methylation patterns of three common human inversions, at 8p23.1, 16p11.2, and 17q21.31 in 1,009 blood samples from children from the Human Early Life Exposome (HELIX) project and in 39 prenatal heart tissue samples. We found inversion-state specific methylation patterns within and nearby flanking each inversion region in both datasets. Additionally, numerous inversion-exposure interactions on methylation levels were identified from early-life exposome data comprising 64 exposures. For instance, children homozygous at inv-8p23.1 and higher meat intake were more susceptible to TDH hypermethylation (P = 3.8 × 10⁻²²); being the inversion, exposure, and gene known risk factors for adult obesity. Inv-8p23.1 associated hypermethylation of GATA4 was also detected across numerous exposures. Our data suggests that the pleiotropic influence of inversions during development and lifetime could be substantially mediated by allele-specific methylation patterns which can be modulated by the exposome.

Analysis of the relationship between presence of common DNA sequence inversions and DNA methylation patterns suggests a role for environmental exposures (such as food intake) in mediating inversion state-specific methylation patterns.

Identification of Sodium Transients Through NaV1.5 Channels as Regulators of Differentiation in Immortalized Dorsal Root Ganglia Neurons

Neuronal differentiation is a complex process through which newborn neurons acquire the morphology of mature neurons and become excitable. We employed a combination of functional and transcriptomic approaches to deconvolute and identify key regulators of the differentiation process of a DRG neuron-derived cell line, and we focused our study on the NaV1.5 ion channel (encoded by Scn5a) as a channel involved in the acquisition of DRG neuronal features. Overexpression of Scn5a enhances the acquisition of neuronal phenotypic features and increases the KCl-elicited hyperexcitability response in a DRG-derived cell line. Moreover, pharmacologic inhibition of the NaV1.5 channel during differentiation hinders the acquisition of phenotypic features of neuronal cells and the hyperexcitability increase in response to changes in the extracellular medium ionic composition. Taken together, these data highlight the relevance of sodium transients in regulating the neuronal differentiation process in a DRG neuron-derived cell line.

New insights into the genetic etiology of Alzheimer's disease and related dementias

Characterization of the genetic landscape of Alzheimer's disease (AD) and related dementias (ADD) provides a unique opportunity for a better understanding of the associated pathophysiological processes. We performed a two-stage genome-wide association study totaling 111,326 clinically diagnosed/'proxy' AD cases and 677,663 controls. We found 75 risk loci, of which 42 were new at the time of analysis. Pathway enrichment analyses confirmed the involvement of amyloid/tau pathways and highlighted microglia implication. Gene prioritization in the new loci identified 31 genes that were suggestive of new genetically associated processes, including the tumor necrosis factor alpha pathway through the linear ubiquitin chain assembly complex. We also built a new genetic risk score associated with the risk of future AD/dementia or progression from mild cognitive impairment to AD/dementia. The improvement in prediction led to a 1.6- to 1.9-fold increase in AD risk from the lowest to the highest decile, in addition to effects of age and the APOE ε4 allele.

Utility of Circulating Cell-Free DNA in Assessing Microsatellite Instability and Loss of Heterozygosity in Breast Cancer Using Human Identification Approach

The diagnostic and prognostic utility of circulating cell-free DNA (cfDNA) in breast cancer (BC) patients was recently reported. Here, we investigated the use of cfDNA to examine microsatellite instability (MSI) and loss of heterozygosity (LOH) for early BC diagnosis. cfDNA and genomic DNA from 41 female BC patients and 40 healthy controls were quantified using NanoDrop spectrophotometry and real-time PCR. The stability of genomic and cfDNA was assessed using a high-resolution AmpFlSTR MiniFiler human identification kit. Significant increases in cfDNA plasma concentrations were observed in BC patients compared to controls. The genotype distribution of the eight autosomal short tandem repeat (STR) loci D7S820, D13S317, D21S11, D2S1338, D18S51, D16S539, FGA, and CSF1PO were in Hardy–Weinberg equilibrium. Significant differences in the allele frequencies of D7S820 allele-8, D21S11 allele-29, allele-30.2, allele-32.2, and CSF1PO allele-11 were seen between BC patients and controls. LOH and MSI were detected in 36.6% of the cfDNA of patients compared to genomic DNA. This study highlights the utility of plasma-derived cfDNA for earlier, less invasive, and cost-effective cancer diagnosis and molecular stratification. It also highlights the potential value of cfDNA in molecular profiling and biomarkers discovery in precision and forensic medicine.

Trio-based exome sequencing reveals a high rate of the de novo variants in intellectual disability

Intellectual disability (ID), a neurodevelopmental disorder affecting 1-3% of the general population, is characterized by limitations in both intellectual function and adaptive skills. The high number of conditions associated with ID underlines its heterogeneous origin and reveals the difficulty of obtaining a rapid and accurate genetic diagnosis. However, the Next Generation Sequencing, and the whole exome sequencing (WES) in particular, has boosted the diagnosis rate associated with ID. In this study, WES performed on 244 trios of patients clinically diagnosed with isolated or syndromic ID and their respective unaffected parents has allowed the identification of the underlying genetic basis of ID in 64 patients, yielding a diagnosis rate of 25.2%. Our results suggest that trio-based WES facilitates ID's genetic diagnosis, particularly in patients who have been extensively waiting for a definitive molecular diagnosis. Moreover, genotypic information from parents provided by trio-based WES enabled the detection of a high percentage (61.5%) of de novo variants inside our cohort. Establishing a quick genetic diagnosis of ID would allow early intervention and better clinical management, thus improving the quality of life of these patients and their families.

Identification of autosomal cis expression quantitative trait methylation (cis eQTMs) in children's blood

Background

The identification of expression quantitative trait methylation (eQTMs), defined as associations between DNA methylation levels and gene expression, might help the biological interpretation of epigenome-wide association studies (EWAS). We aimed to identify autosomal cis eQTMs in children's blood, using data from 832 children of the Human Early Life Exposome (HELIX) project.

Methods

Blood DNA methylation and gene expression were measured with the Illumina 450K and the Affymetrix HTA v2 arrays, respectively. The relationship between methylation levels and expression of nearby genes (1 Mb window centered at the transcription start site, TSS) was assessed by fitting 13.6 M linear regressions adjusting for sex, age, cohort, and blood cell composition.

Results

We identified 39,749 blood autosomal cis eQTMs, representing 21,966 unique CpGs (eCpGs, 5.7% of total CpGs) and 8,886 unique transcript clusters (eGenes, 15.3% of total transcript clusters, equivalent to genes). In 87.9% of these cis eQTMs, the eCpG was located at <250 kb from eGene's TSS; and 58.8% of all eQTMs showed an inverse relationship between the methylation and expression levels. Only around half of the autosomal cis-eQTMs eGenes could be captured through annotation of the eCpG to the closest gene. eCpGs had less measurement error and were enriched for active blood regulatory regions and for CpGs reported to be associated with environmental exposures or phenotypic traits. In 40.4% of the eQTMs, the CpG and the eGene were both associated with at least one genetic variant. The overlap of autosomal cis eQTMs in children's blood with those described in adults was small (13.8%), and age-shared cis eQTMs tended to be proximal to the TSS and enriched for genetic variants.

Conclusions

This catalogue of autosomal cis eQTMs in children's blood can help the biological interpretation of EWAS findings and is publicly available at https://helixomics.isglobal.org/ and at Dryad (doi:10.5061/dryad.fxpnvx0t0).

Funding

The study has received funding from the European Community's Seventh Framework Programme (FP7/2007-206) under grant agreement no 308333 (HELIX project); the H2020-EU.3.1.2. - Preventing Disease Programme under grant agreement no 874583 (ATHLETE project); from the European Union's Horizon 2020 research and innovation programme under grant agreement no 733206 (LIFECYCLE project), and from the European Joint Programming Initiative "A Healthy Diet for a Healthy Life" (JPI HDHL and Instituto de Salud Carlos III) under the grant agreement no AC18/00006 (NutriPROGRAM project). The genotyping was supported by the projects PI17/01225 and PI17/01935, funded by the Instituto de Salud Carlos III and co-funded by European Union (ERDF, "A way to make Europe") and the Centro Nacional de Genotipado-CEGEN (PRB2-ISCIII). BiB received core infrastructure funding from the Wellcome Trust (WT101597MA) and a joint grant from the UK Medical Research Council (MRC) and Economic and Social Science Research Council (ESRC) (MR/N024397/1). INMA data collections were supported by grants from the Instituto de Salud Carlos III, CIBERESP, and the Generalitat de Catalunya-CIRIT. KANC was funded by the grant of the Lithuanian Agency for Science Innovation and Technology (6-04-2014_31V-66). The Norwegian Mother, Father and Child Cohort Study is supported by the Norwegian Ministry of Health and Care Services and the Ministry of Education and Research. The Rhea project was financially supported by European projects (EU FP6-2003-Food-3-NewGeneris, EU FP6. STREP Hiwate, EU FP7 ENV.2007.1.2.2.2. Project No 211250 Escape, EU FP7-2008-ENV-1.2.1.4 Envirogenomarkers, EU FP7-HEALTH-2009- single stage CHICOS, EU FP7 ENV.2008.1.2.1.6. Proposal No 226285 ENRIECO, EU- FP7- HEALTH-2012 Proposal No 308333 HELIX), and the Greek Ministry of Health (Program of Prevention of obesity and neurodevelopmental disorders in preschool children, in Heraklion district, Crete, Greece: 2011-2014; "Rhea Plus": Primary Prevention Program of Environmental Risk Factors for Reproductive Health, and Child Health: 2012-15). We acknowledge support from the Spanish Ministry of Science and Innovation through the "Centro de Excelencia Severo Ochoa 2019-2023" Program (CEX2018-000806-S), and support from the Generalitat de Catalunya through the CERCA Program. MV-U and CR-A were supported by a FI fellowship from the Catalan Government (FI-DGR 2015 and #016FI_B 00272). MC received funding from Instituto Carlos III (Ministry of Economy and Competitiveness) (CD12/00563 and MS16/00128).

CIBERER: Spanish National Network for Research on Rare Diseases: a highly productive collaborative initiative

CIBER (Center for Biomedical Network Research; Centro de Investigación Biomédica En Red) is a public national consortium created in 2006 under the umbrella of the Spanish National Institute of Health Carlos III (ISCIII). This innovative research structure comprises 11 different specific areas dedicated to the main public health priorities in the National Health System. CIBERER, the thematic area of CIBER focused on rare diseases (RDs) currently consists of 75 research groups belonging to universities, research centers, and hospitals of the entire country. CIBERER's mission is to be a center prioritizing and favoring collaboration and cooperation between biomedical and clinical research groups, with special emphasis on the aspects of genetic, molecular, biochemical, and cellular research of RDs. This research is the basis for providing new tools for the diagnosis and therapy of low‐prevalence diseases, in line with the International Rare Diseases Research Consortium (IRDiRC) objectives, thus favoring translational research between the scientific environment of the laboratory and the clinical setting of health centers. In this article, we intend to review CIBERER's 15‐year journey and summarize the main results obtained in terms of internationalization, scientific production, contributions toward the discovery of new therapies and novel genes associated to diseases, cooperation with patients' associations and many other topics related to RD research.

Meta-analysis of epigenome-wide association studies in newborns and children show widespread sex differences in blood DNA methylation

BACKGROUND

Among children, sex-specific differences in disease prevalence, age of onset, and susceptibility have been observed in health conditions including asthma, immune response, metabolic health, some pediatric and adult cancers, and psychiatric disorders. Epigenetic modifications such as DNA methylation may play a role in the sexual differences observed in diseases and other physiological traits.

METHODS

We performed a meta-analysis of the association of sex and cord blood DNA methylation at over 450,000 CpG sites in 8438 newborns from 17 cohorts participating in the Pregnancy And Childhood Epigenetics (PACE) Consortium. We also examined associations of child sex with DNA methylation in older children ages 5.5-10 years from 8 cohorts (n = 4268).

RESULTS

In newborn blood, sex was associated at Bonferroni level significance with differences in DNA methylation at 46,979 autosomal CpG sites (p < 1.3 × 10-7) after adjusting for white blood cell proportions and batch. Most of those sites had lower methylation levels in males than in females. Of the differentially methylated CpG sites identified in newborn blood, 68% (31,727) met look-up level significance (p < 1.1 × 10-6) in older children and had methylation differences in the same direction.

CONCLUSIONS

This is a large-scale meta-analysis examining sex differences in DNA methylation in newborns and older children. Expanding upon previous studies, we replicated previous findings and identified additional autosomal sites with sex-specific differences in DNA methylation. Differentially methylated sites were enriched in genes involved in cancer, psychiatric disorders, and cardiovascular phenotypes.

Genome-wide meta-analysis for Alzheimer's disease cerebrospinal fluid biomarkers

Amyloid-beta 42 (Aβ42) and phosphorylated tau (pTau) levels in cerebrospinal fluid (CSF) reflect core features of the pathogenesis of Alzheimer's disease (AD) more directly than clinical diagnosis. Initiated by the European Alzheimer & Dementia Biobank (EADB), the largest collaborative effort on genetics underlying CSF biomarkers was established, including 31 cohorts with a total of 13,116 individuals (discovery n = 8074; replication n = 5042 individuals). Besides the APOE locus, novel associations with two other well-established AD risk loci were observed; CR1 was shown a locus for Aβ42 and BIN1 for pTau. GMNC and C16orf95 were further identified as loci for pTau, of which the latter is novel. Clustering methods exploring the influence of all known AD risk loci on the CSF protein levels, revealed 4 biological categories suggesting multiple Aβ42 and pTau related biological pathways involved in the etiology of AD. In functional follow-up analyses, GMNC and C16orf95 both associated with lateral ventricular volume, implying an overlap in genetic etiology for tau levels and brain ventricular volume.

Genomic Characterization of Host Factors Related to SARS-CoV-2 Infection in People with Dementia and Control Populations: The GR@ACE/DEGESCO Study

Emerging studies have suggested several chromosomal regions as potential host genetic factors involved in the susceptibility to SARS-CoV-2 infection and disease outcome. We nested a COVID-19 genome-wide association study using the GR@ACE/DEGESCO study, searching for susceptibility factors associated with COVID-19 disease. To this end, we compared 221 COVID-19 confirmed cases with 17,035 individuals in whom the COVID-19 disease status was unknown. Then, we performed a meta-analysis with the publicly available data from the COVID-19 Host Genetics Initiative. Because the APOE locus has been suggested as a potential modifier of COVID-19 disease, we added sensitivity analyses stratifying by dementia status or by disease severity. We confirmed the existence of the 3p21.31 region (LZTFL1, SLC6A20) implicated in the susceptibility to SARS-CoV-2 infection and TYK2 gene might be involved in COVID-19 severity. Nevertheless, no statistically significant association was observed in the COVID-19 fatal outcome or in the stratified analyses (dementia-only and non-dementia strata) for the APOE locus not supporting its involvement in SARS-CoV-2 pathobiology or COVID-19 prognosis.

xDEEP-MSI: Explainable Bias-Rejecting Microsatellite Instability Deep Learning System in Colorectal Cancer

The prediction of microsatellite instability (MSI) using deep learning (DL) techniques could have significant benefits, including reducing cost and increasing MSI testing of colorectal cancer (CRC) patients. Nonetheless, batch effects or systematic biases are not well characterized in digital histology models and lead to overoptimistic estimates of model performance. Methods to not only palliate but to directly abrogate biases are needed. We present a multiple bias rejecting DL system based on adversarial networks for the prediction of MSI in CRC from tissue microarrays (TMAs), trained and validated in 1788 patients from EPICOLON and HGUA. The system consists of an end-to-end image preprocessing module that tile samples at multiple magnifications and a tissue classification module linked to the bias-rejecting MSI predictor. We detected three biases associated with the learned representations of a baseline model: the project of origin of samples, the patient's spot and the TMA glass where each spot was placed. The system was trained to directly avoid learning the batch effects of those variables. The learned features from the bias-ablated model achieved maximum discriminative power with respect to the task and minimal statistical mean dependence with the biases. The impact of different magnifications, types of tissues and the model performance at tile vs patient level is analyzed. The AUC at tile level, and including all three selected tissues (tumor epithelium, mucin and lymphocytic regions) and 4 magnifications, was 0.87 ± 0.03 and increased to 0.9 ± 0.03 at patient level. To the best of our knowledge, this is the first work that incorporates a multiple bias ablation technique at the DL architecture in digital pathology, and the first using TMAs for the MSI prediction task.

SDCBP Modulates Stemness and Chemoresistance in Head and Neck Squamous Cell Carcinoma through Src Activation

Simple Summary

Drug resistance is the principal limiting factor to achieving good survival rates in patients with cancer. The identification of potential biomarkers for diagnosis and prognostic prediction, as well as the design of new molecular-targeted treatments, will be essential to improving head and neck squamous cell carcinoma (HNSCC) patient outcomes. In this sense, the sensitization of resistant cells and cancer stem cells (CSCs) represents a major challenge in cancer therapy. We conducted a proteomic study involving cisplatin-resistance and CSCs with the aim to unravel the molecular and cellular mechanisms by which tumor cells acquire resistance to chemotherapy. Syntenin-1 (SDCBP) was identified as an important protein involved in the chemoresistance and stemness of HNSCC tumors.

Abstract

To characterize the mechanisms that govern chemoresistance, we performed a comparative proteomic study analyzing head and neck squamous cell carcinoma (HNSCC) cells: CCL-138 (parental), CCL-138-R (cisplatin-resistant), and cancer stem cells (CSCs). Syntenin-1 (SDCBP) was upregulated in CCL-138-R cells and CSCs over parental cells. SDCBP depletion sensitized biopsy-derived and established HNSCC cell lines to cisplatin (CDDP) and reduced CSC markers, Src activation being the main SDCBP downstream target. In mice, SDCBP-depleted cells formed tumors with decreased mitosis, Ki-67 positivity, and metastasis over controls. Moreover, the fusocellular pattern of CCL-138-R cell-derived tumors reverted to a more epithelial morphology upon SDCBP silencing. Importantly, SDCBP expression was associated with Src activation, poor differentiated tumor grade, advanced tumor stage, and shorter survival rates in a series of 382 HNSCC patients. Our results reveal that SDCBP might be a promising therapeutic target for effectively eliminating CSCs and CDDP resistance.

The early-life exposome and epigenetic age acceleration in children

The early-life exposome influences future health and accelerated biological aging has been proposed as one of the underlying biological mechanisms. We investigated the association between more than 100 exposures assessed during pregnancy and in childhood (including indoor and outdoor air pollutants, built environment, green environments, tobacco smoking, lifestyle exposures, and biomarkers of chemical pollutants), and epigenetic age acceleration in 1,173 children aged 7 years old from the Human Early-Life Exposome project. Age acceleration was calculated based on Horvath's Skin and Blood clock using child blood DNA methylation measured by Infinium HumanMethylation450 BeadChips. We performed an exposure-wide association study between prenatal and childhood exposome and age acceleration. Maternal tobacco smoking during pregnancy was nominally associated with increased age acceleration. For childhood exposures, indoor particulate matter absorbance (PM_abs) and parental smoking were nominally associated with an increase in age acceleration. Exposure to the organic pesticide dimethyl dithiophosphate and the persistent pollutant polychlorinated biphenyl-138 (inversely associated with child body mass index) were protective for age acceleration. None of the associations remained significant after multiple-testing correction. Pregnancy and childhood exposure to tobacco smoke and childhood exposure to indoor PM_abs may accelerate epigenetic aging from an early age.

Obesity-related genetic determinants of stroke

As obesity, circulating lipids and other vascular/metabolic factors influence the risk of stroke, we examined if genetic variants associated with these conditions are related to risk of stroke using a case-control study in Galicia, Spain. A selection of 200 single-nucleotide polymorphisms previously found to be related to obesity, body mass index, circulating lipids, type 2 diabetes, heart failure, obesity-related cancer and cerebral infarction were genotyped in 465 patients diagnosed with stroke and 480 population-based controls. An unsupervised Lasso regression procedure was carried out for single-nucleotide polymorphism selection based on their potential effect on stroke according to obesity. Selected genotypes were further analysed through multivariate logistic regression to study their association with risk of stroke. Using unsupervised selection procedures, nine single-nucleotide polymorphisms were found to be related to risk of stroke overall and after stratification by obesity. From these, rs10761731, rs2479409 and rs6511720 in obese subjects [odds ratio (95% confidence interval) = 0.61 (0.39-0.95) (P = 0.027); 0.54 (0.35-0.84) (P = 0.006) and 0.42 (0.22-0.80) (P = 0.0075), respectively], and rs865686 in non-obese subjects [odds ratio (95% confidence interval) = 0.67 (0.48-0.94) (P = 0.019)], were independently associated with risk of stroke after multivariate logistic regression procedures. The associations between the three single-nucleotide polymorphisms found to be associated with stroke risk in obese subjects were more pronounced among females; for rs10761731, odds ratios among obese males and females were 1.07 (0.58-1.97) (P = 0.84), and 0.31 (0.14-0.69) (P = 0.0018), respectively; for rs2479409, odd ratios were 0.66 (0.34-1.27) (P = 0.21), and 0.49 (0.24-0.99) (P = 0.04), for obese males and females, respectively; the stroke-rs6511720 association was also slightly more pronounced among obese females, odds ratios were 0.33 (0.13-0.87) (P = 0.022), and 0.28 (0.09-0.85) (P = 0.02) for obese males and females, respectively. The rs865686-stroke association was more pronounced among non-obese males [odds ratios = 0.61 (0.39-0.96) (P = 0.029) and 0.72 (0.42-1.22) (P = 0.21), for non-obese males and females, respectively]. A combined genetic score of variants rs10761731, rs2479409 and rs6511720 was highly predictive of stroke risk among obese subjects (P = 2.04 × 10^-5), particularly among females (P = 4.28 × 10^-6). In summary, single-nucleotide polymorphisms rs1076173, rs2479409 and rs6511720 were found to independently increase the risk of stroke in obese subjects after adjustment for established risk factors. A combined score with the three genomic variants was an independent predictor of risk of stroke among obese subjects in our population.

Development and Evaluation of the Ancestry Informative Marker Panel of the VISAGE Basic Tool

We detail the development of the ancestry informative single nucleotide polymorphisms (SNPs) panel forming part of the VISAGE Basic Tool (BT), which combines 41 appearance predictive SNPs and 112 ancestry predictive SNPs (three SNPs shared between sets) in one massively parallel sequencing (MPS) multiplex, whereas blood-based age analysis using methylation markers is run in a parallel MPS analysis pipeline. The selection of SNPs for the BT ancestry panel focused on established forensic markers that already have a proven track record of good sequencing performance in MPS, and the overall SNP multiplex scale closely matched that of existing forensic MPS assays. SNPs were chosen to differentiate individuals from the five main continental population groups of Africa, Europe, East Asia, America, and Oceania, extended to include differentiation of individuals from South Asia. From analysis of 1000 Genomes and HGDP-CEPH samples from these six population groups, the BT ancestry panel was shown to have no classification error using the Bayes likelihood calculators of the Snipper online analysis portal. The differentiation power of the component ancestry SNPs of BT was balanced as far as possible to avoid bias in the estimation of co-ancestry proportions in individuals with admixed backgrounds. The balancing process led to very similar cumulative population-specific divergence values for Africa, Europe, America, and Oceania, with East Asia being slightly below average, and South Asia an outlier from the other groups. Comparisons were made of the African, European, and Native American estimated co-ancestry proportions in the six admixed 1000 Genomes populations, using the BT ancestry panel SNPs and 572,000 Affymetrix Human Origins array SNPs. Very similar co-ancestry proportions were observed down to a minimum value of 10%, below which, low-level co-ancestry was not always reliably detected by BT SNPs. The Snipper analysis portal provides a comprehensive population dataset for the BT ancestry panel SNPs, comprising a 520-sample standardised reference dataset; 3445 additional samples from 1000 Genomes, HGDP-CEPH, Simons Foundation and Estonian Biocentre genome diversity projects; and 167 samples of six populations from in-house genotyping of individuals from Middle East, North and East African regions complementing those of the sampling regimes of the other diversity projects.

Variability of multi-omics profiles in a population-based child cohort

BACKGROUND

Multiple omics technologies are increasingly applied to detect early, subtle molecular responses to environmental stressors for future disease risk prevention. However, there is an urgent need for further evaluation of stability and variability of omics profiles in healthy individuals, especially during childhood.

METHODS

We aimed to estimate intra-, inter-individual and cohort variability of multi-omics profiles (blood DNA methylation, gene expression, miRNA, proteins and serum and urine metabolites) measured 6 months apart in 156 healthy children from five European countries. We further performed a multi-omics network analysis to establish clusters of co-varying omics features and assessed the contribution of key variables (including biological traits and sample collection parameters) to omics variability.

RESULTS

All omics displayed a large range of intra- and inter-individual variability depending on each omics feature, although all presented a highest median intra-individual variability. DNA methylation was the most stable profile (median 37.6% inter-individual variability) while gene expression was the least stable (6.6%). Among the least stable features, we identified 1% cross-omics co-variation between CpGs and metabolites (e.g. glucose and CpGs related to obesity and type 2 diabetes). Explanatory variables, including age and body mass index (BMI), explained up to 9% of serum metabolite variability.

CONCLUSIONS

Methylation and targeted serum metabolomics are the most reliable omics to implement in single time-point measurements in large cross-sectional studies. In the case of metabolomics, sample collection and individual traits (e.g. BMI) are important parameters to control for improved comparability, at the study design or analysis stage. This study will be valuable for the design and interpretation of epidemiological studies that aim to link omics signatures to disease, environmental exposures, or both.

Identification of Novel Regulators of Zalcitabine-Induced Neuropathic Pain

Neuropathic pain is one of the foremost adverse effects that worsens quality of life for patients undergoing an antiretroviral treatment. Currently, there are no effective analgesics for relieving it; thus, there is an urgent need to develop novel treatments for neuropathic pain. Previously, we described and validated F11 cells as a model of DRG (dorsal root ganglia) neurons. In the current work, we employed F11 cells to identify regulators of antiretroviral-induced neuropathic pain combining functional and transcriptomic analysis. The antiretroviral zalcitabine (ddC) increased the excitability of differentiated F11 cells associated with calcium signaling without morphological changes in the neuronal phenotype, mimicking the observed increase of painful signaling in patients suffering from antiretroviral-induced neuropathic pain. Employing RNA sequencing, we observed that zalcitabine treatment upregulated genes related with oxidative stress and calcium homeostasis. The functional impact of the transcriptomic changes was explored, finding that the exposure to zalcitabine significantly increased intracellular oxidative stress and reduced store-operated calcium entry (SOCE). Because the functional and transcriptomic evidence points toward fundamental changes in calcium signaling and oxidative stress upon zalcitabine exposure, we identified that NAD(P)H quinone dehydrogenase and the sarcoplasmic/endoplasmic reticulum calcium ATPase 3 were involved in zalcitabine-induced hyperexcitability of F11 cells. Overexpression of those genes increases the calcium-elicited hyperexcitability response and reduces SOCE, as well as increases intracellular ROS levels. These data do not only mimic the effects of zalcitabine but also highlight the relevance of oxidative stress and of calcium-mediated signaling in antiretroviral-induced hyperexcitability of sensory neurons, shedding light on new therapeutic targets for antiviral-induced neuropathic pain.

Genetic Susceptibility to Periodontal Disease in Down Syndrome: A Case-Control Study

Severe periodontitis is prevalent in Down syndrome (DS). This study aimed to identify genetic variations associated with periodontitis in individuals with DS. The study group was distributed into DS patients with periodontitis (n = 50) and DS patients with healthy periodontium (n = 36). All samples were genotyped with the “Axiom Spanish Biobank” array, which contains 757,836 markers. An association analysis at the individual marker level using logistic regression, as well as at the gene level applying the sequence kernel association test (SKAT) was performed. The most significant genes were included in a pathway analysis using the free DAVID software. C12orf74 (rs4315121, p = 9.85 × 10⁻⁵, OR = 8.84), LOC101930064 (rs4814890, p = 9.61 × 10⁻⁵, OR = 0.13), KBTBD12 (rs1549874, p = 8.27 × 10⁻⁵, OR = 0.08), PIWIL1 (rs11060842, p = 7.82 × 10⁻⁵, OR = 9.05) and C16orf82 (rs62030877, p = 8.92 × 10⁻⁵, OR = 0.14) showed a higher probability in the individual analysis. The analysis at the gene level highlighted PIWIL, MIR9-2, LHCGR, TPR and BCR. At the signaling pathway level, PI3K-Akt, long-term depression and FoxO achieved nominal significance (p = 1.3 × 10⁻², p = 5.1 × 10⁻³, p = 1.2 × 10⁻², respectively). In summary, various metabolic pathways are involved in the pathogenesis of periodontitis in DS, including PI3K-Akt, which regulates cell proliferation and inflammatory response.

Common variants in Alzheimer's disease and risk stratification by polygenic risk scores

Genetic discoveries of Alzheimer's disease are the drivers of our understanding, and together with polygenetic risk stratification can contribute towards planning of feasible and efficient preventive and curative clinical trials. We first perform a large genetic association study by merging all available case-control datasets and by-proxy study results (discovery n = 409,435 and validation size n = 58,190). Here, we add six variants associated with Alzheimer's disease risk (near APP, CHRNE, PRKD3/NDUFAF7, PLCG2 and two exonic variants in the SHARPIN gene). Assessment of the polygenic risk score and stratifying by APOE reveal a 4 to 5.5 years difference in median age at onset of Alzheimer's disease patients in APOE ɛ4 carriers. Because of this study, the underlying mechanisms of APP can be studied to refine the amyloid cascade and the polygenic risk score provides a tool to select individuals at high risk of Alzheimer's disease.

LIPG endothelial lipase and breast cancer risk by subtypes

Experimental data showed that endothelial lipase (LIPG) is a crucial player in breast cancer. However, very limited data exists on the role of LIPG on the risk of breast cancer in humans. We examined the LIPG-breast cancer association within our population-based case-control study from Galicia, Spain, BREOGAN (BREast Oncology GAlicia Network). Plasma LIPG and/or OxLDL were measured on 114 breast cancer cases and 82 controls from our case-control study, and were included in the present study. The risk of breast cancer increased with increasing levels of LIPG (multivariable OR for the highest category (95% CI) 2.52 (1.11-5.81), P-trend = 0.037). The LIPG-breast cancer association was restricted to Pre-menopausal breast cancer (Multivariable OR for the highest LIPG category (95% CI) 4.76 (0.94-28.77), P-trend = 0.06, and 1.79 (0.61-5.29), P-trend = 0.372, for Pre-menopausal and Post-menopausal breast cancer, respectively). The LIPG-breast cancer association was restricted to Luminal A breast cancers (Multivariable OR for the highest LIPG category (95% CI) 3.70 (1.42-10.16), P-trend = 0.015, and 2.05 (0.63-7.22), P-trend = 0.311, for Luminal A and non-Luminal A breast cancers, respectively). Subset analysis only based on HER2 receptor indicated that the LIPG-breast cancer relationship was restricted to HER2-negative breast cancers (Multivariable OR for the highest LIPG category (95% CI) 4.39 (1.70-12.03), P-trend = 0.012, and 1.10 (0.28-4.32), P-trend = 0.745, for HER2-negative and HER2-positive tumors, respectively). The LIPG-breast cancer association was restricted to women with high total cholesterol levels (Multivariable OR for the highest LIPG category (95% CI) 6.30 (2.13-20.05), P-trend = 0.018, and 0.65 (0.11-3.28), P-trend = 0.786, among women with high and low cholesterol levels, respectively). The LIPG-breast cancer association was also restricted to non-postpartum breast cancer (Multivariable OR for the highest LIPG category (95% CI) 3.83 (1.37-11.39), P-trend = 0.003, and 2.35 (0.16-63.65), P-trend = 0.396, for non-postpartum and postpartum breast cancer, respectively), although we lacked precision. The LIPG-breast cancer association was more pronounced among grades II and III than grade I breast cancers (Multivariable ORs for the highest category of LIPG (95% CI) 2.73 (1.02-7.69), P-trend = 0.057, and 1.90 (0.61-6.21), P-trend = 0.170, for grades II and III, and grade I breast cancers, respectively). No association was detected for OxLDL levels and breast cancer (Multivariable OR for the highest versus the lowest category (95% CI) 1.56 (0.56-4.32), P-trend = 0.457).

Schuurs-Hoeijmakers Syndrome (PACS1 Neurodevelopmental Disorder): Seven Novel Patients and a Review

Schuurs–Hoeijmakers syndrome (SHMS) or PACS1 Neurodevelopmental disorder is a rare disorder characterized by intellectual disability, abnormal craniofacial features and congenital malformations. SHMS is an autosomal dominant hereditary disease caused by pathogenic variants in the PACS1 gene. PACS1 is a trans-Golgi-membrane traffic regulator that directs protein cargo and several viral envelope proteins. It is upregulated during human embryonic brain development and has low expression after birth. So far, only 54 patients with SHMS have been reported. In this work, we report on seven new identified SHMS individuals with the classical c.607C > T: p.Arg206Trp PACS1 pathogenic variant and review clinical and molecular aspects of all the patients reported in the literature, providing a summary of clinical findings grouped as very frequent (≥75% of patients), frequent (50–74%), infrequent (26–49%) and rare (less than ≤25%).

Tumor Profiling at the Service of Cancer Therapy

Cancer treatment options have evolved significantly in the past few years. From the initial surgical procedures, to the latest next-generation technologies, we are now in the position to analyze and understand tumors in a one-by-one basis and use that to our advantage to provide with individualized treatment options that may increase patient survival. In this review, we will focus on how tumor profiling has evolved over the past decades to deliver more efficient and personalized treatment options, and how novel technologies can help us envisage the future of precision oncology toward a better management and, ultimately, increased survival.

CSVS, a crowdsourcing database of the Spanish population genetic variability

The knowledge of the genetic variability of the local population is of utmost importance in personalized medicine and has been revealed as a critical factor for the discovery of new disease variants. Here, we present the Collaborative Spanish Variability Server (CSVS), which currently contains more than 2000 genomes and exomes of unrelated Spanish individuals. This database has been generated in a collaborative crowdsourcing effort collecting sequencing data produced by local genomic projects and for other purposes. Sequences have been grouped by ICD10 upper categories. A web interface allows querying the database removing one or more ICD10 categories. In this way, aggregated counts of allele frequencies of the pseudo-control Spanish population can be obtained for diseases belonging to the category removed. Interestingly, in addition to pseudo-control studies, some population studies can be made, as, for example, prevalence of pharmacogenomic variants, etc. In addition, this genomic data has been used to define the first Spanish Genome Reference Panel (SGRP1.0) for imputation. This is the first local repository of variability entirely produced by a crowdsourcing effort and constitutes an example for future initiatives to characterize local variability worldwide. CSVS is also part of the GA4GH Beacon network. CSVS can be accessed at: http://csvs.babelomics.org/.

The Road so Far in Colorectal Cancer Pharmacogenomics: Are We Closer to Individualised Treatment?

In recent decades, survival rates in colorectal cancer have improved greatly due to pharmacological treatment. However, many patients end up developing adverse drug reactions that can be severe or even life threatening, and that affect their quality of life. These remain a limitation, as they may force dose reduction or treatment discontinuation, diminishing treatment efficacy. From candidate gene approaches to genome-wide analysis, pharmacogenomic knowledge has advanced greatly, yet there is still huge and unexploited potential in the use of novel technologies such as next-generation sequencing strategies. This review summarises the road of colorectal cancer pharmacogenomics so far, presents considerations and directions to be taken for further works and discusses the path towards implementation into clinical practice.

Association between DNA methylation and ADHD symptoms from birth to school age: a prospective meta-analysis

Attention-deficit and hyperactivity disorder (ADHD) is a common childhood disorder with a substantial genetic component. However, the extent to which epigenetic mechanisms play a role in the etiology of the disorder is unknown. We performed epigenome-wide association studies (EWAS) within the Pregnancy And Childhood Epigenetics (PACE) Consortium to identify DNA methylation sites associated with ADHD symptoms at two methylation assessment periods: birth and school age. We examined associations of both DNA methylation in cord blood with repeatedly assessed ADHD symptoms (age 4-15 years) in 2477 children from 5 cohorts and of DNA methylation at school age with concurrent ADHD symptoms (age 7-11 years) in 2374 children from 9 cohorts, with 3 cohorts participating at both timepoints. CpGs identified with nominal significance (p < 0.05) in either of the EWAS were correlated between timepoints (ρ = 0.30), suggesting overlap in associations; however, top signals were very different. At birth, we identified nine CpGs that predicted later ADHD symptoms (p < 1 × 10^-7), including ERC2 and CREB5. Peripheral blood DNA methylation at one of these CpGs (cg01271805 in the promoter region of ERC2, which regulates neurotransmitter release) was previously associated with brain methylation. Another (cg25520701) lies within the gene body of CREB5, which previously was associated with neurite outgrowth and an ADHD diagnosis. In contrast, at school age, no CpGs were associated with ADHD with p < 1 × 10^-7. In conclusion, we found evidence in this study that DNA methylation at birth is associated with ADHD. Future studies are needed to confirm the utility of methylation variation as biomarker and its involvement in causal pathways.

Broadening the Applicability of a Custom Multi-Platform Panel of Microhaplotypes: Bio-Geographical Ancestry Inference and Expanded Reference Data

The development of microhaplotype (MH) panels for massively parallel sequencing (MPS) platforms is gaining increasing relevance for forensic analysis. Here, we expand the applicability of a 102 autosomal and 11 X-chromosome panel of MHs, previously validated with both MiSeq and Ion S5 MPS platforms and designed for identification purposes. We have broadened reference population data for identification purposes, including data from 240 HGDP-CEPH individuals of native populations from North Africa, the Middle East, Oceania and America. Using the enhanced population data, the panel was evaluated as a marker set for bio-geographical ancestry (BGA) inference, providing a clear differentiation of the five main continental groups of Africa, Europe, East Asia, Native America, and Oceania. An informative degree of differentiation was also achieved for the population variation encompassing North Africa, Middle East, Europe, South Asia, and East Asia. In addition, we explored the potential for individual BGA inference from simple mixed DNA, by simulation of mixed profiles followed by deconvolution of mixture components.

Impact of CYP2C19 Genotype and Drug Interactions on Voriconazole Plasma Concentrations: A Spain Pharmacogenetic-Pharmacokinetic Prospective Multicenter Study

BACKGROUND

Voriconazole, a first-line agent for the treatment of invasive fungal infections, is mainly metabolized by cytochrome P450 (CYP) 2C19. A significant portion of patients fail to achieve therapeutic voriconazole trough concentrations, with a consequently increased risk of therapeutic failure.

OBJECTIVE

To show the association between subtherapeutic voriconazole concentrations and factors affecting voriconazole pharmacokinetics: CYP2C19 genotype and drug-drug interactions.

METHODS

Adults receiving voriconazole for antifungal treatment or prophylaxis were included in a multicenter prospective study conducted in Spain. The prevalence of subtherapeutic voriconazole troughs was analyzed in the rapid metabolizer and ultra-rapid metabolizer patients (RMs and UMs, respectively), and compared with the rest of the patients. The relationship between voriconazole concentration, CYP2C19 phenotype, adverse events (AEs), and drug-drug interactions was also assessed.

RESULTS

In this study 78 patients were included with a wide variability in voriconazole plasma levels with only 44.8% of patients attaining trough concentrations within the therapeutic range of 1 and 5.5 µg/ml. The allele frequency of *17 variant was found to be 29.5%. Compared with patients with other phenotypes, RMs and UMs had a lower voriconazole plasma concentration (RM/UM: 1.85 ± 0.24 µg/ml vs other phenotypes: 2.36 ± 0.26 µg/ml). Adverse events were more common in patients with higher voriconazole concentrations (p<0.05). No association between voriconazole trough concentration and other factors (age, weight, route of administration, and concomitant administration of enzyme inducer, enzyme inhibitor, glucocorticoids, or proton pump inhibitors) was found.

CONCLUSION

These results suggest the potential clinical utility of using CYP2C19 genotype-guided voriconazole dosing to achieve concentrations in the therapeutic range in the early course of therapy. Larger studies are needed to confirm the impact of pharmacogenetics on voriconazole pharmacokinetics.

In utero and childhood exposure to tobacco smoke and multi-layer molecular signatures in children

BACKGROUND

The adverse health effects of early life exposure to tobacco smoking have been widely reported. In spite of this, the underlying molecular mechanisms of in utero and postnatal exposure to tobacco smoke are only partially understood. Here, we aimed to identify multi-layer molecular signatures associated with exposure to tobacco smoke in these two exposure windows.

METHODS

We investigated the associations of maternal smoking during pregnancy and childhood secondhand smoke (SHS) exposure with molecular features measured in 1203 European children (mean age 8.1 years) from the Human Early Life Exposome (HELIX) project. Molecular features, covering 4 layers, included blood DNA methylation and gene and miRNA transcription, plasma proteins, and sera and urinary metabolites.

RESULTS

Maternal smoking during pregnancy was associated with DNA methylation changes at 18 loci in child blood. DNA methylation at 5 of these loci was related to expression of the nearby genes. However, the expression of these genes themselves was only weakly associated with maternal smoking. Conversely, childhood SHS was not associated with blood DNA methylation or transcription patterns, but with reduced levels of several serum metabolites and with increased plasma PAI1 (plasminogen activator inhibitor-1), a protein that inhibits fibrinolysis. Some of the in utero and childhood smoking-related molecular marks showed dose-response trends, with stronger effects with higher dose or longer duration of the exposure.

CONCLUSION

In this first study covering multi-layer molecular features, pregnancy and childhood exposure to tobacco smoke were associated with distinct molecular phenotypes in children. The persistent and dose-dependent changes in the methylome make CpGs good candidates to develop biomarkers of past exposure. Moreover, compared to methylation, the weak association of maternal smoking in pregnancy with gene expression suggests different reversal rates and a methylation-based memory to past exposures. Finally, certain metabolites and protein markers evidenced potential early biological effects of postnatal SHS, such as fibrinolysis.

A Comparison of Forensic Age Prediction Models Using Data From Four DNA Methylation Technologies

Individual age estimation can be applied to criminal, legal, and anthropological investigations. DNA methylation has been established as the biomarker of choice for age prediction, since it was observed that specific CpG positions in the genome show systematic changes during an individual’s lifetime, with progressive increases or decreases in methylation levels. Subsequently, several forensic age prediction models have been reported, providing average age prediction error ranges of ±3–4 years, using a broad spectrum of technologies and underlying statistical analyses. DNA methylation assessment is not categorical but quantitative. Therefore, the detection platform used plays a pivotal role, since quantitative and semi-quantitative technologies could potentially result in differences in detected DNA methylation levels. In the present study, we analyzed as a shared sample pool, 84 blood-based DNA controls ranging from 18 to 99 years old using four different technologies: EpiTYPER^®, pyrosequencing, MiSeq, and SNaPshot^TM. The DNA methylation levels detected for CpG sites from ELOVL2, FHL2, and MIR29B2 with each system were compared. A restricted three CpG-site age prediction model was rebuilt for each system, as well as for a combination of technologies, based on previous training datasets, and age predictions were calculated accordingly for all the samples detected with the previous technologies. While the DNA methylation patterns and subsequent age predictions from EpiTYPER^®, pyrosequencing, and MiSeq systems are largely comparable for the CpG sites studied, SNaPshot^TM gives bigger differences reflected in higher predictive errors. However, these differences can be reduced by applying a z-score data transformation.