Characterization of the Common Genetic Variation in the Spanish Population of Navarre.

Large-scale genomic studies have significantly increased our knowledge of genetic variability across populations. Regional genetic profiling is essential for distinguishing common benign variants from disease-causing ones. To this end, we conducted a comprehensive characterization of exonic variants in the population of Navarre (Spain), utilizing whole genome sequencing data from 358 unrelated individuals of Spanish origin. Our analysis revealed 61,410 biallelic single nucleotide variants (SNV) within the Navarrese cohort, with 35% classified as common (MAF > 1%). By comparing allele frequency data from 1000 Genome Project (excluding the Iberian cohort of Spain, IBS), Genome Aggregation Database, and a Spanish cohort (including IBS individuals and data from Medical Genome Project), we identified 1069 SNVs common in Navarre but rare (MAF ≤ 1%) in all other populations. We further corroborated this observation with a second regional cohort of 239 unrelated exomes, which confirmed 676 of the 1069 SNVs as common in Navarre. In conclusion, this study highlights the importance of population-specific characterization of genetic variation to improve allele frequency filtering in sequencing data analysis to identify disease-causing variants.

Mutational landscape of risk variants in comorbid depression and obesity: a next-generation sequencing approach.

Major depression (MD) and obesity are complex genetic disorders that are frequently comorbid. However, the study of both diseases concurrently remains poorly addressed and therefore the underlying genetic mechanisms involved in this comorbidity remain largely unknown. Here we examine the contribution of common and rare variants to this comorbidity through a next-generation sequencing (NGS) approach. Specific genomic regions of interest in MD and obesity were sequenced in a group of 654 individuals from the PISMA-ep epidemiological study. We obtained variants across the entire frequency spectrum and assessed their association with comorbid MD and obesity, both at variant and gene levels. We identified 55 independent common variants and a burden of rare variants in 4 genes (PARK2, FGF21, HIST1H3D and RSRC1) associated with the comorbid phenotype. Follow-up analyses revealed significantly enriched gene-sets associated with biological processes and pathways involved in metabolic dysregulation, hormone signaling and cell cycle regulation. Our results suggest that, while risk variants specific to the comorbid phenotype have been identified, the genes functionally impacted by the risk variants share cell biological processes and signaling pathways with MD and obesity phenotypes separately. To the best of our knowledge, this is the first study involving a targeted sequencing approach toward the study of the comorbid MD and obesity. The framework presented here allowed a deep characterization of the genetics of the co-occurring MD and obesity, revealing insights into the mutational and functional profile that underlies this comorbidity and contributing to a better understanding of the relationship between these two disabling disorders.

Primary Ciliary Dyskinesia and Retinitis Pigmentosa: Novel RPGR Variant and Possible Modifier Gene.

We report a novel RPGR missense variant co-segregated with a familial X-linked retinitis pigmentosa (XLRP) case. The brothers were hemizygous for this variant, but only the proband presented with primary ciliary dyskinesia (PCD). Thus, we aimed to elucidate the role of the RPGR variant and other modifier genes in the phenotypic variability observed in the family and its impact on motile cilia. The pathogenicity of the variant on the RPGR protein was evaluated by in vitro studies transiently transfecting the mutated RPGR gene, and immunofluorescence analysis on nasal brushing samples. Whole-exome sequencing was conducted to identify potential modifier variants. In vitro studies showed that the mutated RPGR protein could not localise to the cilium and impaired cilium formation. Accordingly, RPGR was abnormally distributed in the siblings' nasal brushing samples. In addition, a missense variant in CEP290 was identified. The concurrent RPGR variant influenced ciliary mislocalisation of the protein. We provide a comprehensive characterisation of motile cilia in this XLRP family, with only the proband presenting PCD symptoms. The variant's pathogenicity was confirmed, although it alone does not explain the respiratory symptoms. Finally, the CEP290 gene may be a potential modifier for respiratory symptoms in patients with RPGR mutations.

[Effect of frailty in health resource use in aged cancer patients]. / Efecto de la fragilidad en el consumo de recursos sanitarios en ancianos con cáncer.

OBJECTIVE: To evaluate the effect of frailty on health resource use in aged population with cancer. METHOD: Population-based cohort study with retrospective data collection and follow-up from January 2018 to December 2019 in people ≥65 years with cancer. RESULTS: Overall, 996 individuals were included, with a prevalence of frailty of 22.1%. Mortality at 2 years was 14.1% in the frail and 9.0% in the non-frail (p=0.028). Independently of age and sex, frailty increased the number of urgent hospitalizations (168%) and planned hospitalizations (64%), visits to the emergency room (111%), outpatient consultations (59%), day hospital sessions (30%) and visits to primary care (114%). CONCLUSIONS: Frailty is more prevalent, determines a poorer prognostic and is associated with higher health resource use in aged population with cancer.

High life satisfaction reported among small-scale societies with low incomes.

Global polls have shown that people in high-income countries generally report being more satisfied with their lives than people in low-income countries. The persistence of this correlation, and its similarity to correlations between income and life satisfaction within countries, could lead to the impression that high levels of life satisfaction can only be achieved in wealthy societies. However, global polls have typically overlooked small-scale, nonindustrialized societies, which can provide an alternative test of the consistency of this relationship. Here, we present results from a survey of 2,966 members of Indigenous Peoples and local communities among 19 globally distributed sites. We find that high average levels of life satisfaction, comparable to those of wealthy countries, are reported for numerous populations that have very low monetary incomes. Our results are consistent with the notion that human societies can support very satisfying lives for their members without necessarily requiring high degrees of monetary wealth.

Evidence of the association between increased use of direct oral anticoagulants and a reduction in the rate of atrial fibrillation-related stroke and major bleeding at the population level (2012-2019).

BACKGROUND: The introduction of direct-acting oral anticoagulants (DOACs) has shown to decrease atrial fibrillation (AF)-related stroke and bleeding rates in clinical studies, but there is no certain evidence about their effects at the population level. Our aim was to assess changes in AF-related stroke and major bleeding rates between 2012 and 2019 in Andalusia (Spain), and the association between DOACs use and events rates at the population level. METHODS: All patients with an AF diagnosis from 2012 to 2019 were identified using the Andalusian Health Population Base, that provides clinical information on all Andalusian people. Annual ischemic and hemorrhagic stroke, major bleeding rates, and used antithrombotic treatments were determined. Marginal hazard ratios (HR) were calculated for each treatment. RESULTS: A total of 95,085 patients with an AF diagnosis were identified. Mean age was 76.1±10.2 years (49.7% women). An increase in the use of DOACs was observed throughout the study period in both males and females (p<0.001). The annual rate of ischemic stroke decreased by one third, while that of hemorrhagic stroke and major bleeding decreased 2-3-fold from 2012 to 2019. Marginal HR was lower than 0.50 for DOACs compared to VKA for all ischemic or hemorrhagic events. CONCLUSIONS: In this contemporary population-based study using clinical and administrative databases in Andalusia, a significant reduction in the incidence of AF-related ischemic and hemorrhagic stroke and major bleeding was observed between 2012 and 2019. The increased use of DOACs seems to be associated with this reduction.

Indigenous Peoples' rights in national climate governance: An analysis of Nationally Determined Contributions (NDCs).

Although the recognition of Indigenous Peoples' contributions to climate governance by the international community has gradually increased, a rights-based approach in national climate action is still largely absent. This article analyses the recognition of Indigenous Peoples' rights in Nationally Determined Contributions (NDCs) under the Paris Agreement. We conducted a content analysis of all NDCs submitted between 2016 and May 2022. Through a five-pronged framework of sustainable self-determination, we assessed how the NDCs recognise: i. Indigenous Peoples as rights-holders; ii. Indigenous jurisdiction over land; iii. Indigenous knowledge systems; iv. Indigenous Peoples' right to full and effective participation in climate governance; and v. the legacy of colonialism. NDCs with references related to Indigenous Peoples are increasing. However, questions remain regarding their sincerity and commitment to implementation. States must therefore make more significant efforts to ensure that the NDCs take a rights-based approach and contribute to strengthening Indigenous Peoples' role and say in climate governance.

Real-world evidence with a retrospective cohort of 15,968 COVID-19 hospitalized patients suggests 21 new effective treatments.

PURPOSE: Despite the extensive vaccination campaigns in many countries, COVID-19 is still a major worldwide health problem because of its associated morbidity and mortality. Therefore, finding efficient treatments as fast as possible is a pressing need. Drug repurposing constitutes a convenient alternative when the need for new drugs in an unexpected medical scenario is urgent, as is the case with COVID-19. METHODS: Using data from a central registry of electronic health records (the Andalusian Population Health Database), the effect of prior consumption of drugs for other indications previous to the hospitalization with respect to patient outcomes, including survival and lymphocyte progression, was studied on a retrospective cohort of 15,968 individuals, comprising all COVID-19 patients hospitalized in Andalusia between January and November 2020. RESULTS: Covariate-adjusted hazard ratios and analysis of lymphocyte progression curves support a significant association between consumption of 21 different drugs and better patient survival. Contrarily, one drug, furosemide, displayed a significant increase in patient mortality. CONCLUSIONS: In this study we have taken advantage of the availability of a regional clinical database to study the effect of drugs, which patients were taking for other indications, on their survival. The large size of the database allowed us to control covariates effectively.

Fractal approaches to scaling transformations to sustainability.

Responses to sustainability challenges are not delivering results at the scale and speed called for by science, international agreements, and concerned citizens. Yet there is a tendency to underestimate the large-scale impacts of small-scale, local, and contextualized actions, and particularly the role of individuals in scaling transformations. Here, we explore a fractal approach to scaling sustainability transformations based on "universal values." Universal values are proposed as intrinsic characteristics that connect humans and nature in a coherent, acausal way. Drawing on the Three Spheres of Transformation framework, we consider how enacting universal values can generate fractal-like patterns of sustainability that repeat recursively across scales. Fractal approaches shift the focus from scaling through "things" (e.g., technologies, behaviors, projects) to scaling through a quality of agency based on values that apply to all. We discuss practical steps involved in fractal approaches to scaling transformations to sustainability, provide examples, and conclude with questions for future research.

Placing diverse knowledge systems at the core of transformative climate research.

We argue that solutions-based research must avoid treating climate change as a merely technical problem, recognizing instead that it is symptomatic of the history of European and North American colonialism. It must therefore be addressed by decolonizing the research process and transforming relations between scientific expertise and the knowledge systems of Indigenous Peoples and of local communities. Partnership across diverse knowledge systems can be a path to transformative change only if those systems are respected in their entirety, as indivisible cultural wholes of knowledge, practices, values, and worldviews. This argument grounds our specific recommendations for governance at the local, national, and international scales. As concrete mechanisms to guide collaboration across knowledge systems, we propose a set of instruments based on the principles of consent, intellectual and cultural autonomy, and justice. We recommend these instruments as tools to ensure that collaborations across knowledge systems embody just partnerships in support of a decolonial transformation of relations between human communities and between humanity and the more-than-human world.

A crowdsourcing database for the copy-number variation of the Spanish population.

BACKGROUND: Despite being a very common type of genetic variation, the distribution of copy-number variations (CNVs) in the population is still poorly understood. The knowledge of the genetic variability, especially at the level of the local population, is a critical factor for distinguishing pathogenic from non-pathogenic variation in the discovery of new disease variants. RESULTS: Here, we present the SPAnish Copy Number Alterations Collaborative Server (SPACNACS), which currently contains copy number variation profiles obtained from more than 400 genomes and exomes of unrelated Spanish individuals. By means of a collaborative crowdsourcing effort whole genome and whole exome sequencing data, produced by local genomic projects and for other purposes, is continuously collected. Once checked both, the Spanish ancestry and the lack of kinship with other individuals in the SPACNACS, the CNVs are inferred for these sequences and they are used to populate the database. A web interface allows querying the database with different filters that include ICD10 upper categories. This allows discarding samples from the disease under study and obtaining pseudo-control CNV profiles from the local population. We also show here additional studies on the local impact of CNVs in some phenotypes and on pharmacogenomic variants. SPACNACS can be accessed at: http://csvs.clinbioinfosspa.es/spacnacs/ . CONCLUSION: SPACNACS facilitates disease gene discovery by providing detailed information of the local variability of the population and exemplifies how to reuse genomic data produced for other purposes to build a local reference database.

Assessing the Impact of SARS-CoV-2 Lineages and Mutations on Patient Survival.

OBJECTIVES: More than two years into the COVID-19 pandemic, SARS-CoV-2 still remains a global public health problem. Successive waves of infection have produced new SARS-CoV-2 variants with new mutations for which the impact on COVID-19 severity and patient survival is uncertain. METHODS: A total of 764 SARS-CoV-2 genomes, sequenced from COVID-19 patients, hospitalized from 19th February 2020 to 30 April 2021, along with their clinical data, were used for survival analysis. RESULTS: A significant association of B.1.1.7, the alpha lineage, with patient mortality (log hazard ratio (LHR) = 0.51, C.I. = [0.14,0.88]) was found upon adjustment by all the covariates known to affect COVID-19 prognosis. Moreover, survival analysis of mutations in the SARS-CoV-2 genome revealed 27 of them were significantly associated with higher mortality of patients. Most of these mutations were located in the genes coding for the S, ORF8, and N proteins. CONCLUSIONS: This study illustrates how a combination of genomic and clinical data can provide solid evidence for the impact of viral lineage on patient survival.

Protein and functional isoform levels and genetic variants of the BAFF and APRIL pathway components in systemic lupus erythematosus.

Systemic lupus erythematosus (SLE) is the prototype of an autoimmune disease. Belimumab, a monoclonal antibody targets BAFF, is the only biologic approved for SLE and active lupus nephritis. BAFF is a cytokine with a key-regulatory role in the B cell homeostasis, which acts by binding to three receptors: BAFF-R, TACI and BCMA. TACI and BCMA also bind APRIL. Many studies reported elevated soluble BAFF and APRIL levels in the sera of SLE patients, but other questions about the role of this system in the disease remain open. The study aimed to investigate the utility of the cytokine levels in serum and urine as biomarkers, the role of non-functional isoforms, and the association of gene variants with the disease. This case-control study includes a cohort (women, 18-60 years old) of 100 patients (48% with nephritis) and 100 healthy controls. We used ELISA assays to measure the cytokine concentrations in serum (sBAFF and sAPRIL) and urine (uBAFF and uAPRIL); TaqMan Gene Expression Assays to quantify the relative mRNA expression of ΔBAFF, ßAPRIL, and εAPRIL, and next-generation sequencing to genotype the cytokine (TNFSF13 and TNFSF13B) and receptor (TNFRSF13B, TNFRSF17 and TNFRSF13C) genes. The statistical tests used were: Kruskal-Wallis (qualitative variables), the Spearman Rho coefficient (correlations), the Chi-square and SKAT (association of common and rare genetic variants, respectively). As expected, sBAFF and sAPRIL levels were higher in patients than in controls (p ≤ 0.001) but found differences between patient subgroups. sBAFF and sAPRIL significantly correlated only in patients with nephritis (rs = 0.67, p ≤ 0.001) and ßAPRIL levels were lower in patients with nephritis (p = 0.04), and ΔBAFF levels were lower in patients with dsDNA antibodies (p = 0.04). Rare variants of TNFSF13 and TNFRSF13B and TNFSF13 p.Gly67Arg and TNFRSF13B p.Val220Ala were associated with SLE. Our study supports differences among SLE patient subgroups with diverse clinical features in the BAFF/APRIL pathway. In addition, it suggests the involvement of genetic variants in the susceptibility to the disease.

Incidence and Prevalence of Children's Diffuse Lung Disease in Spain.

BACKGROUND: Children's diffuse lung disease, also known as children's Interstitial Lung Diseases (chILD), are a heterogeneous group of rare diseases with relevant morbidity and mortality, which diagnosis and classification are very complex. Epidemiological data are scarce. The aim of this study was to analyse incidence and prevalence of chILD in Spain. METHODS: Multicentre observational prospective study in patients from 0 to 18 years of age with chILD to analyse its incidence and prevalence in Spain, based on data reported in 2018 and 2019. RESULTS: A total of 381 cases with chILD were notified from 51 paediatric pulmonology units all over Spain, covering the 91.7% of the paediatric population. The average incidence of chILD was 8.18 (CI 95% 6.28-10.48) new cases/million of children per year. The average prevalence of chILD was 46.53 (CI 95% 41.81-51.62) cases/million of children. The age group with the highest prevalence were children under 1 year of age. Different types of disorders were seen in children 2-18 years of age compared with children 0-2 years of age. Most frequent cases were: primary pulmonary interstitial glycogenosis in neonates (17/65), neuroendocrine cell hyperplasia of infancy in infants from 1 to 12 months (44/144), idiopathic pulmonary haemosiderosis in children from 1 to 5 years old (13/74), hypersensitivity pneumonitis in children from 5 to 10 years old (9/51), and scleroderma in older than 10 years old (8/47). CONCLUSIONS: We found a higher incidence and prevalence of chILD than previously described probably due to greater understanding and increased clinician awareness of these rare diseases.

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/.

SMN1 copy-number and sequence variant analysis from next-generation sequencing data.

Spinal muscular atrophy (SMA) is a severe neuromuscular autosomal recessive disorder affecting 1/10,000 live births. Most SMA patients present homozygous deletion of SMN1, while the vast majority of SMA carriers present only a single SMN1 copy. The sequence similarity between SMN1 and SMN2, and the complexity of the SMN locus makes the estimation of the SMN1 copy-number by next-generation sequencing (NGS) very difficult. Here, we present SMAca, the first python tool to detect SMA carriers and estimate the absolute SMN1 copy-number using NGS data. Moreover, SMAca takes advantage of the knowledge of certain variants specific to SMN1 duplication to also identify silent carriers. This tool has been validated with a cohort of 326 samples from the Navarra 1000 Genomes Project (NAGEN1000). SMAca was developed with a focus on execution speed and easy installation. This combination makes it especially suitable to be integrated into production NGS pipelines. Source code and documentation are available at https://www.github.com/babelomics/SMAca.

Hypothalamic miR-30 regulates puberty onset via repression of the puberty-suppressing factor, Mkrn3.

Mkrn3, the maternally imprinted gene encoding the makorin RING-finger protein-3, has recently emerged as putative pubertal repressor, as evidenced by central precocity caused by MKRN3 mutations in humans; yet, the molecular underpinnings of this key regulatory action remain largely unexplored. We report herein that the microRNA, miR-30, with three binding sites in a highly conserved region of its 3' UTR, operates as repressor of Mkrn3 to control pubertal onset. Hypothalamic miR-30b expression increased, while Mkrn3 mRNA and protein content decreased, during rat postnatal maturation. Neonatal estrogen exposure, causing pubertal alterations, enhanced hypothalamic Mkrn3 and suppressed miR-30b expression in female rats. Functional in vitro analyses demonstrated a strong repressive action of miR-30b on Mkrn3 3' UTR. Moreover, central infusion during the juvenile period of target site blockers, tailored to prevent miR-30 binding to Mkrn3 3' UTR, reversed the prepubertal down-regulation of hypothalamic Mkrn3 protein and delayed female puberty. Collectively, our data unveil a novel hypothalamic miRNA pathway, involving miR-30, with a prominent role in the control of puberty via Mkrn3 repression. These findings expand our current understanding of the molecular basis of puberty and its disease states.

TransFlow: a modular framework for assembling and assessing accurate de novo transcriptomes in non-model organisms.

BACKGROUND: The advances in high-throughput sequencing technologies are allowing more and more de novo assembling of transcriptomes from many new organisms. Some degree of automation and evaluation is required to warrant reproducibility, repetitivity and the selection of the best possible transcriptome. Workflows and pipelines are becoming an absolute requirement for such a purpose, but the issue of assembling evaluation for de novo transcriptomes in organisms lacking a sequenced genome remains unsolved. An automated, reproducible and flexible framework called TransFlow to accomplish this task is described. RESULTS: TransFlow with its five independent modules was designed to build different workflows depending on the nature of the original reads. This architecture enables different combinations of Illumina and Roche/454 sequencing data, and can be extended to other sequencing platforms. Its capabilities are illustrated with the selection of reliable plant reference transcriptomes and the assembling six transcriptomes (three case studies for grapevine leaves, olive tree pollen, and chestnut stem, and other three for haustorium, epiphytic structures and their combination for the phytopathogenic fungus Podosphaera xanthii). Arabidopsis and poplar transcriptomes revealed to be the best references. A common result regarding de novo assemblies is that Illumina paired-end reads of 100 nt in length assembled with OASES can provide reliable transcriptomes, while the contribution of longer reads is noticeable only when they complement a set of short, single-reads. CONCLUSIONS: TransFlow can handle up to 181 different assembling strategies. Evaluation based on principal component analyses allows its self-adaptation to different sets of reads to provide a suitable transcriptome for each combination of reads and assemblers. As a result, each case study has its own behaviour, prioritises evaluation parameters, and gives an objective and automated way for detecting the best transcriptome within a pool of them. Sequencing data type and quantity (preferably several hundred millions of 2×100 nt or longer), assemblers (OASES for Illumina, MIRA4 and EULER-SR reconciled with CAP3 for Roche/454) and strategy (preferably scaffolding with OASES, and probably merging with Roche/454 when available) arise as the most impacting factors.

Bioinformatic Prediction of S-Nitrosylation Sites in Large Protein Datasets.

S-nitrosylation is an essential and reversible posttranslational modification of proteins involved in numerous biological processes. The experimental determination of S-nitrosylation sites is laborious and time-consuming. Therefore, the use of computational prediction tools of this modification represents a convenient first approach to generate useful information for subsequent experimental verification. Here we describe an in silico analysis pipeline to integrate the use of several bionformatic tools while dealing with big query protein sets.