Search

Background: Community-acquired pneumonia (CAP) is associated with high morbidity and hospitalization rate. In infectious diseases, host genetics plays a critical role in susceptibility and immune response, and the immune pathways involved are highly dependent on the microorganism and its route of infection. Here we aimed to identify genetic risk loci for CAP using a case-control genome-wide association study (GWAS)., Methods: We performed a GWAS on 3,765 Spanish individuals, including 257 adult patients hospitalized with CAP and 3,508 population controls. Pneumococcal CAP was documented in 30% of patients; the remaining 70% were selected among patients with unidentified microbiological etiology. We tested 7,6 million imputed genotypes using logistic regressions. UK Biobank GWAS of bacterial pneumonia were used for results validation. Subsequently, we prioritized genes and likely causal variants based on Bayesian fine mapping and functional evidence. Imputation and association of classical HLA alleles and amino acids were also conducted., Results: Six independent sentinel variants reached the genome-wide significance (p < 5 × 10 -8 ), three on chromosome 6p21.32, and one for each of the chromosomes 4q28.2, 11p12, and 20q11.22. Only one variant at 6p21.32 was validated in independent GWAS of bacterial and pneumococcal pneumonia. Our analyses prioritized C4orf33 on 4q28.2, TAPBP on 6p21.32, and ZNF341 on 20q11.22. Interestingly, genetic defects of TAPBP and ZNF341 are previously known inborn errors of immunity predisposing to bacterial pneumonia, including pneumococcus and Haemophilus influenzae. Associations were all non-significant for the classical HLA alleles., Conclusions: We completed a GWAS of CAP and identified four novel risk loci involved in CAP susceptibility., (© 2024. The Author(s).)

Background: The clinical course of idiopathic pulmonary fibrosis (IPF) is highly variable and unpredictable, with multiple genetic variants influencing IPF outcomes. Notably, rare pathogenic variants in telomere-related genes are associated with poorer clinical outcomes in these patients. Here we assessed whether rare qualifying variants (QVs) in monogenic adult-onset pulmonary fibrosis (PF) genes are associated with IPF survival. Using polygenic risk scores (PRS), we also evaluated the influence of common IPF risk variants in individuals carrying these QVs., Methods: We identified QVs in telomere and non-telomere genes linked to monogenic PF forms using whole-genome sequences (WGS) from 888 Pulmonary Fibrosis Foundation Patient Registry (PFFPR) individuals. We also derived a PRS for IPF (PRS-IPF) from 19 previously published common sentinel IPF variants. Using regression models, we then examined the mutual relationships of QVs and PRS-IPF and their association with survival. Validation of results was sought in WGS from an independent IPF study (PROFILE, n=472), and results from the two cohorts were meta-analyzed., Results: Carriers of QVs in monogenic adult-onset PF genes, representing nearly 1 out of 6 IPF patients, were associated with lower PRS-IPF (Odds Ratio [OR]: 1.79; 95% Confidence Interval [CI]: 1.15-2.81; p=0.010) and shorter survival (Hazard Ratio [HR]: 1.53; 95% CI: 1.12-2.10; p=7.3×10 -3 ). Notably, carriers of pathogenic variants at telomere genes showed the strongest association with survival (HR: 1.76; 95% CI: 1.13-2.76; p=0.013). The meta-analysis of the results showed a consistent direction of effect across both cohorts., Conclusions: We revealed the opposite effects of QVs and PRS-IPF on IPF survival. Thus, a distinct IPF molecular subtype might be defined by QVs in monogenic adult-onset PF genes. Assessing the carrier status for QVs and modelling PRS-IPF promises to further contribute to predicting disease progression among IPF patients.

The genetic basis of severe COVID-19 has been thoroughly studied, and many genetic risk factors shared between populations have been identified. However, reduced sample sizes from non-European groups have limited the discovery of population-specific common risk loci. In this second study nested in the SCOURGE consortium, we conducted a genome-wide association study (GWAS) for COVID-19 hospitalization in admixed Americans, comprising a total of 4702 hospitalized cases recruited by SCOURGE and seven other participating studies in the COVID-19 Host Genetic Initiative. We identified four genome-wide significant associations, two of which constitute novel loci and were first discovered in Latin American populations ( BAZ2B and DDIAS ). A trans-ethnic meta-analysis revealed another novel cross-population risk locus in CREBBP . Finally, we assessed the performance of a cross-ancestry polygenic risk score in the SCOURGE admixed American cohort. This study constitutes the largest GWAS for COVID-19 hospitalization in admixed Latin Americans conducted to date. This allowed to reveal novel risk loci and emphasize the need of considering the diversity of populations in genomic research., Competing Interests: SD, RC, AL, JL, Md, IQ, RG, VN, MP, JT, JN, JA, CA, SA, VA, BA, NA, ÁA, EA, CA, MA, MA, RB, MB, MB, JB, MB, LB, MB, RB, MB, EC, CC, LC, JC, RC, MC, JC, MC, MD, AD, Vd, Ad, ED, RE, MF, MF, UF, AF, TF, MG, BG, JG, BÁ, FB, AG, JG, JG, MH, AH, MJ, ML, AB, RL, EM, CM, VM, OM, IM, MM, AM, JM, EM, PM, VC, SO, EO, MP, EP, NS, PP, PP, MP, TP, MP, GP, EP, GP, AP, SL, SR, MF, ER, FR, JR, FR, JR, PR, JS, JS, ET, AT, JT, AT, LT, MU, JV, ZY, RZ, Id, AR, PS, LR, EG, CA, EP, JR, AR, CF, PL, ÁC No competing interests declared, (© 2024, Diz-de Almeida, Cruz et al.)

Background: Monoallelic loss-of-function IKZF1 (IKAROS) variants cause B-cell deficiency or combined immunodeficiency, whereas monoallelic gain-of-function (GOF) IKZF1 variants have recently been reported to cause hypergammaglobulinemia, abnormal plasma cell differentiation, autoimmune and allergic manifestations, and infections., Objective: We studied 7 relatives with autoimmune/inflammatory and lymphoproliferative manifestations to identify the immunologic disturbances and the genetic cause of their disease., Methods: We analyzed biopsy results and performed whole-exome sequencing and immunologic studies., Results: Disease onset occurred at a mean age of 25.2 years (range, 10-64, years). Six patients suffered from autoimmune/inflammatory diseases, 4 had confirmed IG4-related disease (IgG4-RD), and 5 developed B-cell malignancies: lymphoma in 4 and multiple myeloma in the remaining patient. Patients without immunosuppression were not particularly prone to infectious diseases. Three patients suffered from life-threatening coronavirus disease 2019 pneumonia, of whom 1 had autoantibodies neutralizing IFN-α. The recently described IKZF1 GOF p.R183H variant was found in the 5 affected relatives tested and in a 6-year-old asymptomatic girl. Immunologic analysis revealed hypergammaglobulinemia and high frequencies of certain lymphocyte subsets (exhausted B cells, effector memory CD4 T cells, effector memory CD4 T cells that have regained surface expression of CD45RA and CD28 - CD57 + CD4 + and CD8 + T cells, T H 2, and Tfh2 cells) attesting to immune dysregulation. Partial clinical responses to rituximab and corticosteroids were observed, and treatment with lenalidomide, which promotes IKAROS degradation, was initiated in 3 patients., Conclusions: Heterozygosity for GOF IKZF1 variants underlies autoimmunity/inflammatory diseases, IgG4-RD, and B-cell malignancies, the onset of which may occur in adulthood. Clinical and immunologic data are similar to those for patients with unexplained IgG4-RD. Patients may therefore benefit from treatments inhibiting pathways displaying IKAROS-mediated overactivity., (Copyright © 2024. Published by Elsevier Inc.)

Background: North African human populations present a complex demographic scenario due to the presence of an autochthonous genetic component and population substructure, plus extensive gene flow from the Middle East, Europe, and sub-Saharan Africa., Results: We conducted a comprehensive analysis of 364 genomes to construct detailed demographic models for the North African region, encompassing its two primary ethnic groups, the Arab and Amazigh populations. This was achieved through an Approximate Bayesian Computation with Deep Learning (ABC-DL) framework and a novel algorithm called Genetic Programming for Population Genetics (GP4PG). This innovative approach enabled us to effectively model intricate demographic scenarios, utilizing a subset of 16 whole genomes at > 30X coverage. The demographic model suggested by GP4PG exhibited a closer alignment with the observed data compared to the ABC-DL model. Both point to a back-to-Africa origin of North African individuals and a close relationship with Eurasian populations. Results support different origins for Amazigh and Arab populations, with Amazigh populations originating back in Epipaleolithic times, while GP4PG supports Arabization as the main source of Middle Eastern ancestry. The GP4PG model includes population substructure in surrounding populations (sub-Saharan Africa and Middle East) with continuous decaying gene flow after population split. Contrary to ABC-DL, the best GP4PG model does not require pulses of admixture from surrounding populations into North Africa pointing to soft splits as drivers of divergence in North Africa., Conclusions: We have built a demographic model on North Africa that points to a back-to-Africa expansion and a differential origin between Arab and Amazigh populations., (© 2024. The Author(s).)

Methicillin-resistant Staphylococcus aureus (MRSA) is one of the major human pathogens. It could carry numerous resistance genes and virulence factors in its genome, some of which are related to the severity of the infection. An observational, descriptive, cross-sectional study was designed to molecularly analyze MRSA isolates that cause invasive infections in Paraguayan children from 2009 to 2013. Ten representative MRSA isolates of the main clonal complex identified were analyzed with short-read paired-end sequencing and assessed for the virulome, resistome, and phylogenetic relationships. All the genetically linked MRSA isolates were recovered from diverse clinical sources, patients, and hospitals at broad gap periods. The pan-genomic analysis of these clones revealed three major and different clonal complexes (CC30, CC5, and CC8), each composed of clones closely related to each other. The CC30 genomes prove to be a successful clone, strongly installed and disseminated throughout our country, and closely related to other CC30 public genomes from the region and the world. The CC5 shows the highest genetic variability, and the CC8 carried the complete arginine catabolic mobile element (ACME), closely related to the USA300-NAE-ACME+, identified as the major cause of CA-MRSA infections in North America. Multiple virulence and resistance genes were identified for the first time in this study, highlighting the complex virulence profiles of MRSA circulating in the country. This study opens a wide range of new possibilities for future projects and trials to improve the existing knowledge on the epidemiology of MRSA circulating in Paraguay., Importance: The increasing prevalence of methicillin-resistant Staphylococcus aureus (MRSA) is a public health problem worldwide. The most frequent MRSA clones identified in Paraguay in previous studies (including community and hospital acquired) were the Pediatric (CC5-ST5-IV), the Cordobes-Chilean (CC5-ST5-I), the SouthWest Pacific (CC30-ST30-IV), and the Brazilian (CC8-ST239-III) clones. In this study, the pan-genomic analysis of the most representative MRSA clones circulating in invasive infection in Paraguayan children over the years 2009-2013, such as the CC30-ST30-IV, CC5-ST5-IV, and CC8-ST8-IV, was carried out to evaluate their genetic diversity, their repertoire of virulence factors, and antimicrobial resistance determinants. This revealed multiple virulence and resistance genes, highlighting the complex virulence profiles of MRSA circulating in Paraguay. Our work is the first genomic study of MRSA in Paraguay and will contribute to the development of genomic surveillance in the region and our understanding of the global epidemiology of this pathogen., Competing Interests: The authors declare no conflict of interest.

ABSTRACT The emergence of the Omicron SARS-CoV-2 variant of concern has changed the COVID-19 scenario as this variant is characterized by high transmissibility and immune evasion ability. To evaluate the impact of this variant on the Canary Islands (Spain) population, we determined the reinfection rates and disease severity associated with the Omicron sublineages and the previously circulating variants of concern. We performed a retrospective observational study on 21,745 SARS-CoV-2 viral genomes collected from December 2020 to July 2022 in the Canary Islands (Spain). We compared the reinfection rates between lineages using pairwise proportion and Fisher's exact tests. To assess disease severity, we studied the association of Alpha, Delta, BA.1, BA.2, BA.5, and other risk factors on 28-day hospital mortality using logistic regression and Cox proportional hazard models. We observed 127 bona fide reinfection cases throughout the study period. We found that BA.5 had the highest reinfection rate compared to other lineages (vs. Delta p =  2.89 × 10 -25 ; vs. BA.1 p =  5.17 × 10 -11 ; vs. BA.2 p =  0.002). Among the 1,094 hospitalized patients, multivariate logistic regression showed that Alpha (Odds Ratio [OR] =  0.45, 95% Confidence Interval [CI] =  0.23-0.87, p =  0.02), BA.2 (OR =   0.38, 95% CI = 0.22-0.63, p =  1.91 × 10 -4 ), and BA.5 (OR = 0.30, 95% CI = 0.16-0.55, p =  1.05 × 10 -4 ) had lower 28-day hospital mortality compared to Delta. These results were confirmed by using Cox proportional hazard models. Omicron lineages, and in particular BA.5, were associated with higher reinfection rates and lower disease severity (28-day hospital mortality) than previously circulating variants of concern.

In anthropological, medical, and forensic studies, the nonrecombinant region of the human Y chromosome (NRY) enables accurate reconstruction of pedigree relationships and retrieval of ancestral information. Using high-throughput sequencing (HTS) data, we present a benchmarking analysis of command-line tools for NRY haplogroup classification. The evaluation was performed using paired Illumina data from whole-genome sequencing (WGS) and whole-exome sequencing (WES) experiments from 50 unrelated donors. Additionally, as a validation, we also used paired WGS/WES datasets of 54 individuals from the 1000 Genomes Project. Finally, we evaluated the tools on data from third-generation HTS obtained from a subset of donors and one reference sample. Our results show that WES, despite typically offering less genealogical resolution than WGS, is an effective method for determining the NRY haplogroup. Y-LineageTracker and Yleaf showed the highest accuracy for WGS data, classifying precisely 98% and 96% of the samples, respectively. Yleaf outperforms all benchmarked tools in the WES data, classifying approximately 90% of the samples. Yleaf, Y-LineageTracker, and pathPhynder can correctly classify most samples (88%) sequenced with third-generation HTS. As a result, Yleaf provides the best performance for applications that use WGS and WES. Overall, our study offers researchers with a guide that allows them to select the most appropriate tool to analyze the NRY region using both second- and third-generation HTS data., Competing Interests: The authors declare no conflict of interest., (© 2023 The Authors.)

Rationale: Idiopathic pulmonary fibrosis (IPF) is a devastating disease characterized by limited treatment options and high mortality. A better understanding of the molecular drivers of IPF progression is needed. Objectives: To identify and validate molecular determinants of IPF survival. Methods: A staged genome-wide association study was performed using paired genomic and survival data. Stage I cases were drawn from centers across the United States and Europe and stage II cases from Vanderbilt University. Cox proportional hazards regression was used to identify gene variants associated with differential transplantation-free survival (TFS). Stage I variants with nominal significance ( P  < 5 × 10 -5 ) were advanced for stage II testing and meta-analyzed to identify those reaching genome-wide significance ( P  < 5 × 10 -8 ). Downstream analyses were performed for genes and proteins associated with variants reaching genome-wide significance. Measurements and Main Results: After quality controls, 1,481 stage I cases and 397 stage II cases were included in the analysis. After filtering, 9,075,629 variants were tested in stage I, with 158 meeting advancement criteria. Four variants associated with TFS with consistent effect direction were identified in stage II, including one in an intron of PCSK6 (proprotein convertase subtilisin/kexin type 6) reaching genome-wide significance (hazard ratio, 4.11 [95% confidence interval, 2.54-6.67]; P  = 9.45 × 10 -9 ). PCSK6 protein was highly expressed in IPF lung parenchyma. PCSK6 lung staining intensity, peripheral blood gene expression, and plasma concentration were associated with reduced TFS. Conclusions: We identified four novel variants associated with IPF survival, including one in PCSK6 that reached genome-wide significance. Downstream analyses suggested that PCSK6 protein plays a potentially important role in IPF progression.

X-linked recessive deficiency of TLR7, a MyD88- and IRAK-4-dependent endosomal ssRNA sensor, impairs SARS-CoV-2 recognition and type I IFN production in plasmacytoid dendritic cells (pDCs), thereby underlying hypoxemic COVID-19 pneumonia with high penetrance. We report 22 unvaccinated patients with autosomal recessive MyD88 or IRAK-4 deficiency infected with SARS-CoV-2 (mean age: 10.9 yr; 2 mo to 24 yr), originating from 17 kindreds from eight countries on three continents. 16 patients were hospitalized: six with moderate, four with severe, and six with critical pneumonia, one of whom died. The risk of hypoxemic pneumonia increased with age. The risk of invasive mechanical ventilation was also much greater than in age-matched controls from the general population (OR: 74.7, 95% CI: 26.8-207.8, P < 0.001). The patients' susceptibility to SARS-CoV-2 can be attributed to impaired TLR7-dependent type I IFN production by pDCs, which do not sense SARS-CoV-2 correctly. Patients with inherited MyD88 or IRAK-4 deficiency were long thought to be selectively vulnerable to pyogenic bacteria, but also have a high risk of hypoxemic COVID-19 pneumonia., (© 2023 García García et al.)

On July 23, 2022, monkeypox disease (mpox) was declared a Public Emergency of International Concern (PHEIC) by the World Health Organization (WHO) due to a multicountry outbreak. In Europe, several cases of mpox virus (MPXV) infection related to this outbreak were detected in the Canary Islands (Spain). Here we describe the combination of viral DNA sequencing and bioinformatic approaches, including methods for de novo genome assembly and short- and long-read technologies, used to reconstruct the first MPXV genome isolated in the Canary Islands on the 31st of May 2022 from a male adult patient with mild symptoms. The same sequencing and bioinformatic approaches were then validated with three other positive cases of MPXV infection from the same mpox outbreak. We obtained the best results using a reference-based approach with short reads, evidencing 46-79 nucleotide variants against viral sequences from the 2018-2019 mpox outbreak and placing the viral sequences in the new B.1 sublineage of clade IIb of the MPXV classification. This study of MPXV demonstrates the potential of metagenomics sequencing for rapid and precise pathogen identification., Competing Interests: The authors declare no conflict of interest., (© 2023 The Authors.)

The conquest of the Canary Islands by Europeans began at the beginning of the 15th century and culminated in 1496 with the surrender of the aborigines. The collapse of the aboriginal population during the conquest and the arrival of settlers caused a drastic change in the demographic composition of the archipelago. To shed light on this historical process, we analyzed 896 mitogenomes of current inhabitants from the seven main islands. Our findings confirm the continuity of aboriginal maternal contributions and the persistence of their genetic footprints in the current population, even at higher levels (>60% on average) than previously evidenced. Moreover, the age estimates for most autochthonous founder lineages support a first aboriginal arrival to the islands at the beginning of the first millennium. We also revealed for the first time that the main recognizable genetic influences from Europe are from Portuguese and Galicians., Competing Interests: The authors declare no potential conflicts of interest with respect to the authorship and/or publication of this article. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results., (© 2022 The Author(s).)

Most causal variants of Mendelian diseases are exonic. Whole-exome sequencing (WES) has become the diagnostic gold standard, but causative variant prioritization constitutes a bottleneck. Here we assessed an in-house sample-to-sequence pipeline and benchmarked free prioritization tools for germline causal variants from WES data. WES of 61 unselected patients with a known genetic disease cause was obtained. Variant prioritizations were performed by diverse tools and recorded to obtain a diagnostic yield when the causal variant was present in the first, fifth, and 10th top rankings. A fraction of causal variants was not captured by WES (8.2%) or did not pass the quality control criteria (13.1%). Most of the applications inspected were unavailable or had technical limitations, leaving nine tools for complete evaluation. Exomiser performed best in the top first rankings, while LIRICAL led in the top fifth rankings. Based on the more conservative top 10th rankings, Xrare had the highest diagnostic yield, followed by a three-way tie among Exomiser, LIRICAL, and PhenIX, then followed by AMELIE, TAPES, Phen-Gen,  AIVar, and VarNote-PAT. Xrare, Exomiser, LIRICAL, and PhenIX are the most efficient options for variant prioritization in real patient WES data., (© 2022 Wiley Periodicals LLC.)

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection and the associated coronavirus disease 2019 (COVID-19), which severely affect the respiratory system and several organs and tissues, and may lead to death, have shown how science can respond when challenged by a global emergency, offering as a response a myriad of rapid technological developments. Development of vaccines at lightning speed is one of them. SARS-CoV-2 outbreaks have stressed healthcare systems, questioning patients care by using standard non-adapted therapies and diagnostic tools. In this scenario, nanotechnology has offered new tools, techniques and opportunities for prevention, for rapid, accurate and sensitive diagnosis and treatment of COVID-19. In this review, we focus on the nanotechnological applications and nano-based materials (i.e., personal protective equipment) to combat SARS-CoV-2 transmission, infection, organ damage and for the development of new tools for virosurveillance, diagnose and immune protection by mRNA and other nano-based vaccines. All the nano-based developed tools have allowed a historical, unprecedented, real time epidemiological surveillance and diagnosis of SARS-CoV-2 infection, at community and international levels. The nano-based technology has help to predict and detect how this Sarbecovirus is mutating and the severity of the associated COVID-19 disease, thereby assisting the administration and public health services to make decisions and measures for preparedness against the emerging variants of SARS-CoV-2 and severe or lethal COVID-19., (Copyright © 2022 Valenzuela-Fernández, Cabrera-Rodriguez, Ciuffreda, Perez-Yanes, Estevez-Herrera, González-Montelongo, Alcoba-Florez, Trujillo-González, García-Martínez de Artola, Gil-Campesino, Díez-Gil, Lorenzo-Salazar, Flores and Garcia-Luis.)

Next-generation sequencing (NGS) applications have flourished in the last decade, permitting the identification of cancer driver genes and profoundly expanding the possibilities of genomic studies of cancer, including melanoma. Here we aimed to present a technical review across many of the methodological approaches brought by the use of NGS applications with a focus on assessing germline and somatic sequence variation. We provide cautionary notes and discuss key technical details involved in library preparation, the most common problems with the samples, and guidance to circumvent them. We also provide an overview of the sequence-based methods for cancer genomics, exposing the pros and cons of targeted sequencing vs. exome or whole-genome sequencing (WGS), the fundamentals of the most common commercial platforms, and a comparison of throughputs and key applications. Details of the steps and the main software involved in the bioinformatics processing of the sequencing results, from preprocessing to variant prioritization and filtering, are also provided in the context of the full spectrum of genetic variation (SNVs, indels, CNVs, structural variation, and gene fusions). Finally, we put the emphasis on selected bioinformatic pipelines behind (a) short-read WGS identification of small germline and somatic variants, (b) detection of gene fusions from transcriptomes, and (c) de novo assembly of genomes from long-read WGS data. Overall, we provide comprehensive guidance across the main methodological procedures involved in obtaining sequencing results for the most common short- and long-read NGS platforms, highlighting key applications in melanoma research.

Here, we describe the results of a genome-wide study conducted in 11 939 coronavirus disease 2019 (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 < 5 × 10-8) was crossed for variants in 3p21.31 and 21q22.11 loci only among males (P = 1.3 × 10-22 and P = 8.1 × 10-12, respectively), and for variants in 9q21.32 near TLE1 only among females (P = 4.4 × 10-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.7 × 10-8) and ARHGAP33 (P = 1.3 × 10-8), respectively. The meta-analysis of both phases with four European studies stratified by sex from the Host Genetics Initiative (HGI) 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.1 × 10-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 Author(s) 2022. Published by Oxford University Press.)

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., (© 2022. The Author(s).)

The current inhabitants of the Canary Islands have a unique genetic makeup in the European diversity landscape due to the existence of African footprints from recent admixture events, especially of North African components (> 20%). The underrepresentation of non-Europeans in genetic studies and the sizable North African ancestry, which is nearly absent from all existing catalogs of worldwide genetic diversity, justify the need to develop CIRdb, a population-specific reference catalog of natural genetic variation in the Canary Islanders. Based on array genotyping of the selected unrelated donors and comparisons against available datasets from European, sub-Saharan, and North African populations, we illustrate the intermediate genetic differentiation of Canary Islanders between Europeans and North Africans and the existence of within-population differences that are likely driven by genetic isolation. Here we describe the overall design and the methods that are being implemented to further develop CIRdb. This resource will help to strengthen the implementation of Precision Medicine in this population by contributing to increase the diversity in genetic studies. Among others, this will translate into improved ability to fine map disease genes and simplify the identification of causal variants and estimate the prevalence of unattended Mendelian diseases., (© 2022. The Author(s).)

Hereditary angioedema (HAE) is a rare disease where known causes involve C1 inhibitor dysfunction or dysregulation of the kinin cascade. The updated HAE management guidelines recommend performing genetic tests to reach a precise diagnosis. Unfortunately, genetic tests are still uncommon in the diagnosis routine. Here, we characterized for the first time the genetic causes of HAE in affected families from the Canary Islands (Spain). Whole-exome sequencing data was obtained from 41 affected patients and unaffected relatives from 29 unrelated families identified in the archipelago. The Hereditary Angioedema Database Annotation (HADA) tool was used for pathogenicity classification and causal variant prioritization among the genes known to cause HAE. Manual reclassification of prioritized variants was used in those families lacking known causal variants. We detected a total of eight different variants causing HAE in this patient series, affecting essentially SERPING1 and F12 genes, one of them being a novel SERPING1 variant (c.686-12A>G) with a predicted splicing effect which was reclassified as likely pathogenic in one family. Altogether, the diagnostic yield by assessing previously reported causal genes and considering variant reclassifications according to the American College of Medical Genetics guidelines reached 66.7% (95% Confidence Interval [CI]: 30.1-91.0) in families with more than one affected member and 10.0% (95% CI: 1.8-33.1) among cases without family information for the disease. Despite the genetic causes of many patients remain to be identified, our results reinforce the need of genetic tests as first-tier diagnostic tool in this disease, as recommended by the international WAO/EAACI guidelines for the management of HAE., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Mendoza-Alvarez, Tosco-Herrera, Muñoz-Barrera, Rubio-Rodríguez, Alonso-Gonzalez, Corrales, Iñigo-Campos, Almeida-Quintana, Martin-Fernandez, Martinez-Beltran, Perez-Rodriguez, Callero, Garcia-Robaina, González-Montelongo, Marcelino-Rodriguez, Lorenzo-Salazar and Flores.)

Several variants of concern (VOCs) explain most of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) epidemic waves in Europe. We aimed to dissect the spread of the SARS-CoV-2 VOCs in the Canary Islands (Spain) between December 2020 and September 2021 at a micro-geographical level. We sequenced the viral genome of 8,224 respiratory samples collected in the archipelago. We observed that Alpha (B.1.1.7) and Delta (B.1.617.2 and sublineages) were ubiquitously present in the islands, while Beta (B.1.351) and Gamma (P.1/P.1.1) had a heterogeneous distribution and were responsible for fewer and more controlled outbreaks. This work represents the largest effort for viral genomic surveillance in the Canary Islands so far, helping the public health bodies in decision-making throughout the pandemic., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Ciuffreda, González-Montelongo, Alcoba-Florez, García-Martínez de Artola, Gil-Campesino, Rodríguez-Pérez, Íñigo-Campos, De Miguel-Martínez, Tosco-Nuñez, Díez-Gil, Valenzuela-Fernández, Lorenzo-Salazar and Flores.)

The transactive response DNA-binding protein (TARDBP/TDP-43) influences the processing of diverse transcripts, including that of histone deacetylase 6 (HDAC6). Here, we assessed TDP-43 activity in terms of regulating CD4+ T-cell permissivity to HIV-1 infection. We observed that overexpression of wt-TDP-43 increased both mRNA and protein levels of HDAC6, resulting in impaired HIV-1 infection independently of the viral envelope glycoprotein complex (Env) tropism. Consistently, using an HIV-1 Env-mediated cell-to-cell fusion model, the overexpression of TDP-43 levels negatively affected viral Env fusion capacity. Silencing of endogenous TDP-43 significantly decreased HDAC6 levels and increased the fusogenic and infection activities of the HIV-1 Env. Using pseudovirus bearing primary viral Envs from HIV-1 individuals, overexpression of wt-TDP-43 strongly reduced the infection activity of Envs from viremic non-progressors (VNP) and rapid progressors (RP) patients down to the levels of the inefficient HIV-1 Envs observed in long-term non-progressor elite controllers (LTNP-EC). On the contrary, silencing endogenous TDP-43 significantly favored the infectivity of primary Envs from VNP and RP individuals, and notably increased the infection of those from LTNP-EC. Taken together, our results indicate that TDP-43 shapes cell permissivity to HIV-1 infection, affecting viral Env fusion and infection capacities by altering the HDAC6 levels and associated tubulin-deacetylase anti-HIV-1 activity.

The precisionFDA Truth Challenge V2 aimed to assess the state of the art of variant calling in challenging genomic regions. Starting with FASTQs, 20 challenge participants applied their variant-calling pipelines and submitted 64 variant call sets for one or more sequencing technologies (Illumina, PacBio HiFi, and Oxford Nanopore Technologies). Submissions were evaluated following best practices for benchmarking small variants with updated Genome in a Bottle benchmark sets and genome stratifications. Challenge submissions included numerous innovative methods, with graph-based and machine learning methods scoring best for short-read and long-read datasets, respectively. With machine learning approaches, combining multiple sequencing technologies performed particularly well. Recent developments in sequencing and variant calling have enabled benchmarking variants in challenging genomic regions, paving the way for the identification of previously unknown clinically relevant variants., Competing Interests: DECLARATION OF INTERESTS C.B. is an employee and shareholder of SAGA Diagnostics AB. A.C., P.-C.C., A.K., M.N., G.B., S.G., and H.Y. are employees of Google, and A.C. is a shareholder. S.D.-B., D.K.-Z., D.T., Ö.K., G.B., K.N., E.A., R.B., I.J.J., A.D., V.S., A.J., and H.S.T. are employees of Seven Bridges Genomics. O.S. and S. T.W. are employees of DNAnexus. G.L., C.M., L.T.F., Y.D., and S.Z. are employees of Genetalks. V.J., M.R., B.L., C.R., S.C., and R.M. are employees of Illumina. S.M.E.S. and M.M. are employees of Roche. C.H. is an employee of Wasai Technology. H.F., Z.L., and L.C. are employees of Sentieon.

The repetitive nature and complexity of some medically relevant genes poses a challenge for their accurate analysis in a clinical setting. The Genome in a Bottle Consortium has provided variant benchmark sets, but these exclude nearly 400 medically relevant genes due to their repetitiveness or polymorphic complexity. Here, we characterize 273 of these 395 challenging autosomal genes using a haplotype-resolved whole-genome assembly. This curated benchmark reports over 17,000 single-nucleotide variations, 3,600 insertions and deletions and 200 structural variations each for human genome reference GRCh37 and GRCh38 across HG002. We show that false duplications in either GRCh37 or GRCh38 result in reference-specific, missed variants for short- and long-read technologies in medically relevant genes, including CBS, CRYAA and KCNE1. When masking these false duplications, variant recall can improve from 8% to 100%. Forming benchmarks from a haplotype-resolved whole-genome assembly may become a prototype for future benchmarks covering the whole genome., (© 2022. This is a U.S. government work and not under copyright protection in the U.S.; foreign copyright protection may apply.)

Sepsis is a severe systemic inflammatory response to infections that is accompanied by organ dysfunction. Although the ancestral genetic background is a relevant factor for sepsis susceptibility, there is a lack of studies using the genetic singularities of a recently admixed population to identify loci involved in sepsis susceptibility. Here we aimed to discover new sepsis loci by completing the first admixture mapping study of sepsis in Canary Islanders, leveraging their distinctive genetic makeup as a mixture of Europeans and African ancestries. We used a case-control approach and inferred local ancestry blocks from genome-wide data from 113,414 polymorphisms genotyped in 343 patients with sepsis and 410 unrelated controls, all ascertained for grandparental origin in the Canary Islands (Spain). Deviations in local ancestries between cases and controls were tested using logistic regressions, followed by fine-mapping analyses based on imputed genotypes, in silico functional assessments, and gene expression analysis centered on the region of interest. The admixture mapping analysis detected that local European ancestry in a locus spanning 1.2 megabases of chromosome 8p23.1 was associated with sepsis (lowest p = 1.37 × 10 -4 ; Odds Ratio [OR] = 0.51; 95%CI = 0.40-0.66). Fine-mapping studies prioritized the variant rs13249564 within intron 1 of MFHAS1 gene associated with sepsis ( p = 9.94 × 10 -4 ; OR = 0.65; 95%CI = 0.50-0.84). Functional and gene expression analyses focused on 8p23.1 allowed us to identify alternative genes with possible biological plausibility such as defensins, which are well-known effector molecules of innate immunity. By completing the first admixture mapping study of sepsis, our results revealed a new genetic locus (8p23.1) harboring a number of genes with plausible implications in sepsis susceptibility., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Hernandez-Beeftink, Marcelino-Rodríguez, Guillen-Guio, Rodríguez-Pérez, Lorenzo-Salazar, Corrales, Díaz-de Usera, González-Montelongo, Domínguez, Espinosa, Villar and Flores.)

Despite asthma has a considerable genetic component, an important proportion of genetic risks remain unknown, especially for non-European populations. Canary Islanders have the largest African genetic ancestry observed among Southwestern Europeans and the highest asthma prevalence in Spain. Here we examined broad chromosomal regions previously associated with an excess of African genetic ancestry in Canary Islanders, with the aim of identifying novel risk variants associated with asthma susceptibility. In a two-stage cases-control study, we revealed a variant within HLA-DQB1 significantly associated with asthma risk (rs1049213, meta-analysis p = 1.30 × 10 -7 , OR [95% CI] = 1.74 [1.41-2.13]) previously associated with asthma and broad allergic phenotype. Subsequent fine-mapping analyses of classical HLA alleles revealed a novel allele significantly associated with asthma protection (HLA-DQA1*01:02, meta-analysis p = 3.98 × 10 -4 , OR [95% CI] = 0.64 [0.50-0.82]) that had been linked to infectious and autoimmune diseases, and peanut allergy. HLA haplotype analyses revealed a novel haplotype DQA1*01:02-DQB1*06:04 conferring asthma protection (meta-analysis p = 4.71 × 10 -4 , OR [95% CI] = 0.47 [0.29- 0.73])., (© 2021. The Author(s).)

Chagas disease is an infection caused by the parasite Trypanosoma cruzi, endemic in Latino America. Leveraging the three-way admixture between Native American (AMR), European (EUR) and African (AFR) populations in Latin Americans, we aimed to better understand the genetic basis of Chagas disease by performing an admixture mapping study in a Colombian population. A two-stage study was conducted, and subjects were classified as seropositive and seronegative for T. cruzi. In stage 1, global and local ancestries were estimated using reference data from the 1000 Genomes Project (1KGP), and local ancestry associations were performed by logistic regression models. The AMR ancestry showed a protective association with Chagas disease within the major histocompatibility complex region [Odds ratio (OR) = 0.74, 95% confidence interval (CI) = 0.66-0.83, lowest P-value = 4.53 × 10-8]. The fine mapping assessment on imputed genotypes combining data from stage 1 and 2 from an independent Colombian cohort, revealed nominally associated variants in high linkage disequilibrium with the top signal (rs2032134, OR = 0.93, 95% CI = 0.90-0.97, P-value = 3.54 × 10-4) in the previously associated locus. To assess ancestry-specific adaptive signals, a selective sweep scan in an AMR reference population from 1KGP together with an in silico functional analysis highlighted the Tripartite Motif family and the human leukocyte antigen genes, with crucial role in the immune response against pathogens. Furthermore, these analyses emphasized the macrophages, neutrophils and eosinophils, as key players in the defense against T. cruzi. This first admixture mapping study in Chagas disease provided novel insights underlying the host immune response in the pathogenesis of this neglected disease., (© The Author(s) 2021. Published by Oxford University Press.)

The mitochondrial genome (mtDNA) is of interest for a range of fields including evolutionary, forensic, and medical genetics. Human mitogenomes can be classified into evolutionary related haplogroups that provide ancestral information and pedigree relationships. Because of this and the advent of high-throughput sequencing (HTS) technology, there is a diversity of bioinformatic tools for haplogroup classification. We present a benchmarking of the 11 most salient tools for human mtDNA classification using empirical whole-genome (WGS) and whole-exome (WES) short-read sequencing data from 36 unrelated donors. We also assessed the best performing tool in third-generation long noisy read WGS data obtained with nanopore technology for a subset of the donors. We found that, for short-read WGS, most of the tools exhibit high accuracy for haplogroup classification irrespective of the input file used for the analysis. However, for short-read WES, Haplocheck and MixEmt were the most accurate tools. Based on the performance shown for WGS and WES, and the accompanying qualitative assessment, Haplocheck stands out as the most complete tool. For third-generation HTS data, we also showed that Haplocheck was able to accurately retrieve mtDNA haplogroups for all samples assessed, although only after following assembly-based approaches (either based on a referenced-based assembly or a hybrid de novo assembly). Taken together, our results provide guidance for researchers to select the most suitable tool to conduct the mtDNA analyses from HTS data., (© 2021. The Author(s).)

Acute respiratory distress syndrome (ARDS) is an inflammatory process of the lungs that develops primarily in response to pulmonary or systemic sepsis, resulting in a disproportionate death toll in intensive care units (ICUs). Given its role as a critical activator of the inflammatory and innate immune responses, previous studies have reported that an increase of circulating cell-free mitochondrial DNA (mtDNA) is a biomarker for fatal outcome in the ICU. Here we analyzed the association of whole-blood mtDNA (wb-mtDNA) copies with 28-day survival from sepsis and sepsis-associated ARDS. We analyzed mtDNA data from 687 peripheral whole-blood samples within 24 h of sepsis diagnosis from unrelated Spanish patients with sepsis (264 with ARDS) included in the GEN-SEP study. The wb-mtDNA copies were obtained from the array intensities of selected probes, with 100% identity with mtDNA and with the largest number of mismatches with the nuclear sequences, and normalized across the individual-probe intensities. We used Cox regression models for testing the association with 28-day survival. We observed that wb-mtDNA copies were significantly associated with 28-day survival in ARDS patients (hazard ratio = 3.65, 95% confidence interval = 1.39-9.59, p = 0.009) but not in non-ARDS patients. Our findings support that wb-mtDNA copies at sepsis diagnosis could be considered an early prognostic biomarker in sepsis-associated ARDS patients. Future studies will be needed to evaluate the mechanistic links of this observation with the pathogenesis of ARDS., Competing Interests: MG-G was employed by SPECTRUM, LLC. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Hernández-Beeftink, Guillen-Guio, Rodríguez-Pérez, Marcelino-Rodríguez, Lorenzo-Salazar, Corrales, Prieto-González, Rodríguez-Pérez, Carriedo, Blanco, Ambrós, González-Higueras, Casanova, González-Garay, Espinosa, Muriel, Domínguez, de Lorenzo, Añón, Soro, Belda, Garcia, Villar and Flores.)

Introduction: Fibrotic interstitial lung diseases (ILDs) are the first indication for lung transplantation (LT). Telomere dysfunction has been associated with poor post-transplant outcomes. The aim of the study was to evaluate the morbi-mortality and quality of life in fibrotic ILDs after lung transplant depending on telomere biology. Methods: Fibrotic ILD patients that underwent lung transplant were allocated to two arms; with or without telomere dysfunction at diagnosis based on the telomere length and telomerase related gene mutations revealed by whole-exome sequencing. Post-transplant evaluation included: (1) short and long-term mortality and complications and (2) quality of life. Results: Fifty-five percent of patients that underwent LT carried rare coding mutations in telomerase-related genes. Patients with telomere shortening more frequently needed extracorporeal circulation and presented a higher rate of early post-transplant hematological complications, longer stay in the intensive care unit (ICU), and a higher number of long-term hospital admissions. However, post-transplant 1-year survival was higher than 80% regardless of telomere dysfunction, with improvement in the quality of life and oxygen therapy withdrawal. Conclusions: Post-transplant morbidity is higher in patients with telomere dysfunction and differs according to elapsed time from transplantation. However, lung transplant improves survival and quality of life and the associated complications are manageable., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Planas-Cerezales, Arias-Salgado, Berastegui, Montes-Worboys, González-Montelongo, Lorenzo-Salazar, Vicens-Zygmunt, Garcia-Moyano, Dorca, Flores, Perona, Román and Molina-Molina.)

Whole-exome sequencing has become a popular technique in research and clinical settings, assisting in disease diagnosis and increasing the understanding of disease pathogenesis. In this study, we aimed to compare common enrichment capture solutions available in the market. Peripheral blood-purified DNA samples were enriched with SureSelect QXT V6 (Agilent) and various Illumina solutions: TruSeq DNA Nano, TruSeq DNA Exome, Nextera DNA Exome, and Illumina DNA Prep with Enrichment, and sequenced on a HiSeq 4000. We found that their percentage of duplicate reads was as much as 2 times higher than previously reported values for the previous HiSeq series. SureSelect QXT and Illumina DNA Prep with Enrichment showed the best average on-target coverage, which improved when off-target regions were included. At high coverage levels and in shared bases, these two solutions and TruSeq DNA Exome provided three of the best performances. With respect to the number of small variants detected, SureSelect QXT presented the lowest number of detected variants in target regions. When off-target regions were considered, its ability equalized to other solutions. Our results show SureSelect QXT and Illumina DNA Prep with Enrichment to be the best enrichment capture solutions.

Background: Hereditary angioedema is a rare genetic condition caused by C1 esterase inhibitor deficiency, dysfunction, or kinin cascade dysregulation, leading to an increased bradykinin plasma concentration. Hereditary angioedema is a poorly recognized clinical entity and is very often misdiagnosed as a histaminergic angioedema. Despite its genetic nature, first-line genetic screening is not integrated in routine diagnosis. Consequently, a delay in the diagnosis, and inaccurate or incomplete diagnosis and treatment of hereditary angioedema are common., Objective: In agreement with recent recommendations from the International Consensus on the Use of Genetics in the Management of Hereditary Angioedema, to facilitate the clinical diagnosis and adapt it to the paradigm of precision medicine and next-generation sequencing-based genetic tests, we aimed to develop a genetic annotation tool, termed Hereditary Angioedema Database Annotation (HADA)., Methods: HADA is built on top of a database of known variants affecting function, including precomputed pathogenic assessment of each variant and a ranked classification according to the current guidelines from the American College of Medical Genetics and Genomics., Results: HADA is provided as a freely accessible, user-friendly web-based interface with versatility for the entry of genetic information. The underlying database can also be incorporated into automated command-line stand-alone annotation tools., Conclusions: HADA can achieve the rapid detection of variants affecting function for different hereditary angioedema types, and further integrates useful information to reduce the diagnosis odyssey and improve its delay., (©Alejandro Mendoza-Alvarez, Adrián Muñoz-Barrera, Luis Alberto Rubio-Rodríguez, Itahisa Marcelino-Rodriguez, Almudena Corrales, Antonio Iñigo-Campos, Ariel Callero, Eva Perez-Rodriguez, Jose Carlos Garcia-Robaina, Rafaela González-Montelongo, Jose Miguel Lorenzo-Salazar, Carlos Flores. Originally published in the Journal of Medical Internet Research (http://www.jmir.org), 09.10.2020.)

The prevalence of asthma symptoms in Canary Islanders, a southwestern European population from Spain, is almost three times higher than the country average. Because the genetic risks identified so far explain <5% of asthma heritability, here we aimed to discover new asthma loci by completing the first admixture mapping study in Canary Islanders leveraging their distinctive genetic makeup, where significant northwest African influences coexist in the European genetic diversity landscape. A 2-stage study was conducted in 1,491 unrelated individuals self-declaring having a Canary Islands origin for the 4 grandparents. Local ancestry estimates were obtained for the shared positions with reference data from putative ancestral populations from Europe, North Africa, and sub-Saharan Africa. Case-control deviations in local ancestry were tested for each ancestry separately using logistic regressions adjusted for principal components, followed by fine-mapping analyses based on imputed genotypes and analyses of the likely deleterious exonic variants. The admixture mapping analysis of asthma detected that local North African ancestry in a locus spanning 365 kb of chromosome 16q23.3 was associated with asthma risk at study-wide significance [lowest P = 1.12 × 10 -4 ; odds ratio (OR) = 2.05; 95% confidence interval (CI) = 1.41-3.00]. Fine-mapping studies identified a variant associated with asthma, and results were replicated in independent samples (rs3852738, OR = 1.34; 95% CI = 1.13-1.59, P = 7.58 × 10 -4 ). Whole exome sequencing data from a subset of individuals revealed an enrichment of likely deleterious variants among asthma cases in 16q23.3, particularly in the phospholipase Cγ2 ( PLCG2 ) gene ( P = 3.67 × 10 -4 ). By completing the first mapping study of asthma in admixed populations from Europe, our results revealed a new plausible asthma locus.

Background: Acute respiratory distress syndrome (ARDS) is a lung inflammatory process caused mainly by sepsis. Most previous studies that identified genetic risks for ARDS focused on candidates with biological relevance. We aimed to identify novel genetic variants associated with ARDS susceptibility and to provide complementary functional evidence of their effect in gene regulation., Methods: We did a case-control genome-wide association study (GWAS) of 1935 European individuals, using patients with sepsis-associated ARDS as cases and patients with sepsis without ARDS as controls. The discovery stage included 672 patients admitted into a network of Spanish intensive care units between January, 2002, and January, 2017. The replication stage comprised 1345 individuals from two independent datasets from the MESSI cohort study (Sep 22, 2008-Nov 30, 2017; USA) and the VISEP (April 1, 2003-June 30, 2005) and MAXSEP (Oct 1, 2007-March 31, 2010) trials of the SepNet study (Germany). Results from discovery and replication stages were meta-analysed to identify association signals. We then used RNA sequencing data from lung biopsies, in-silico analyses, and luciferase reporter assays to assess the functionallity of associated variants., Findings: We identified a novel genome-wide significant association with sepsis-associated ARDS susceptibility (rs9508032, odds ratio [OR] 0·61, 95% CI 0·41-0·91, p=5·18 × 10 -8 ) located within the Fms-related tyrosine kinase 1 (FLT1) gene, which encodes vascular endothelial growth factor receptor 1 (VEGFR-1). The region containing the sentinel variant and its best proxies acted as a silencer for the FLT1 promoter, and alleles with protective effects in ARDS further reduced promoter activity (p=0·0047). A literature mining of all previously described ARDS genes validated the association of vascular endothelial growth factor A (VEGFA; OR 0·55, 95% CI 0·41-0·73; p=4·69 × 10 -5 )., Interpretation: A common variant within the FLT1 gene is associated with sepsis-associated ARDS. Our findings support a role for the vascular endothelial growth factor signalling pathway in ARDS pathogenesis and identify VEGFR-1 as a potential therapeutic target., Funding: Instituto de Salud Carlos III, European Regional Development Funds, Instituto Tecnológico y de Energías Renovables., (Copyright © 2020 Elsevier Ltd. All rights reserved.)

Background: Specific common and rare single nucleotide variants (SNVs) increase the likelihood of developing sporadic idiopathic pulmonary fibrosis (IPF). We performed target-enriched sequencing on three loci previously identified by a genome-wide association study to gain a deeper understanding of the full spectrum of IPF genetic risk and performed a two-stage case-control association study., Methods: A total of 1.7 Mb of DNA from 181 IPF patients was deep sequenced (>100×) across 11p15.5, 14q21.3 and 17q21.31 loci. Comparisons were performed against 501 unrelated controls and replication studies were assessed in 3968 subjects., Results: 36 SNVs were associated with IPF susceptibility in the discovery stage (p<5.0×10 -8 ). After meta-analysis, the strongest association corresponded to rs35705950 (p=9.27×10 -57 ) located upstream from the mucin 5B gene ( MUC5B ). Additionally, a novel association was found for two co-inherited low-frequency SNVs (<5%) in MUC5AC , predicting a missense amino acid change in mucin 5AC (lowest p=2.27×10 -22 ). Conditional and haplotype analyses in 11p15.5 supported the existence of an additional contribution of MUC5AC variants to IPF risk., Conclusions: This study reinforces the significant IPF associations of these loci and implicates MUC5AC as another key player in IPF susceptibility., Competing Interests: Conflict of interest: J.M. Lorenzo-Salazar has nothing to disclose. Conflict of interest: S-F. Ma has nothing to disclose. Conflict of interest: J. Jou has nothing to disclose. Conflict of interest: P-C. Hou has nothing to disclose. Conflict of interest: B. Guillen-Guio has nothing to disclose. Conflict of interest: R.J. Allen has nothing to disclose. Conflict of interest: R.G. Jenkins reports grants from GlaxoSmithKline, during the conduct of the study; grants from Biogen and Galecto, personal fees from Boehringer Ingelheim, Galapagos, Heptares, Pliant and Roche, grants and personal fees from GlaxoSmithKline and MedImmune, and has been on advisory boards for NuMedii and Redex, outside the submitted work. He is a trustee of the British Thoracic Society and Action for Pulmonary Fibrosis, and is an NIHR Research Professor (RP-2017-08-014). Conflict of interest: L.V. Wain reports grants from GlaxoSmithKline, outside the submitted work, and holds a GlaxoSmithKline/British Lung Foundation Chair in Respiratory Research. Conflict of interest: J.M. Oldham has nothing to disclose. Conflict of interest: I. Noth reports personal fees and consultancy fees from Boehringer Ingelheim, and personal fees from Genentech/Hoffman La Roche, Global Blood therapeutics, Sanofi Aventis and Zambon, outside the submitted work. In addition, I. Noth has a patent TOLLIP and NAC in IPF pending. Conflict of interest: C. Flores reports grants from Instituto de Salud Carlos III and Instituto Tecnológico y de Energías Renovables (ITER), during the conduct of the study.

Clear cell renal cell carcinoma (ccRCC) is among the most aggressive histologic subtypes of kidney cancer, representing about 3% of all human cancers. Patients at stage IV have nearly 60% of mortality in 2-3 years after diagnosis. To date, most ccRCC studies have used DNA microarrays and targeted sequencing of a small set of well-established, commonly altered genes. An exception is the large multi-omics study of The Cancer Genome Atlas Kidney Renal Clear Cell Carcinoma (TCGA-KIRC), which identified new ccRCC genes based on whole exome-sequencing (WES) data, and molecular prognostic signatures based on transcriptomics, epigenetics and proteomics data. Applying WES to simultaneously interrogate virtually all exons in the human genome for somatic variation, here we analyzed the burden of coding somatic mutations in metastatic ccRCC primary tumors, and its association with patient mortality from cancer, in patients who received VEGF receptor-targeting drugs as the first-line therapy. To this end, we sequenced the exomes of ten tumor-normal pairs of ccRCC patient tissues from primary biopsies at >100× mean depth and called somatic coding variation. Mutation burden analysis prioritized 138 genes linked to patient mortality. A gene set enrichment analysis evidenced strong statistical support for the abundance of genes involved in the development of kidney cancer ( p = 2.31 × 10 -9 ) and carcinoma ( p = 1.22 × 10 -5 ), with 49 genes having direct links with kidney cancer according to the published records. Two of these genes, SIPA1L2 and EIF3A , demonstrated independent associations with mortality in TCGA-KIRC project data. Besides, three mutational signatures were found to be operative in the tumor exomes, one of which (COSMIC signature 12) has not been previously reported in ccRCC. Selection analysis yielded no detectable evidence of overall positive or negative selection, with the exome-wide number of nonsynonymous substitutions per synonymous site reflecting largely neutral tumor evolution. Despite the limited sample size, our results provide evidence for candidate genes where somatic mutation burden is tentatively associated with patient mortality in metastatic ccRCC, offering new potential pharmacological targets and a basis for further validation studies.

Background: The Oxford Nanopore Technologies (ONT) MinION portable sequencer makes it possible to use cutting-edge genomic technologies in the field and the academic classroom., Results: We present NanoDJ, a Jupyter notebook integration of tools for simplified manipulation and assembly of DNA sequences produced by ONT devices. It integrates basecalling, read trimming and quality control, simulation and plotting routines with a variety of widely used aligners and assemblers, including procedures for hybrid assembly., Conclusions: With the use of Jupyter-facilitated access to self-explanatory contents of applications and the interactive visualization of results, as well as by its distribution into a Docker software container, NanoDJ is aimed to simplify and make more reproducible ONT DNA sequence analysis. The NanoDJ package code, documentation and installation instructions are freely available at https://github.com/genomicsITER/NanoDJ .

The Canary Islands' indigenous people have been the subject of substantial archaeological, anthropological, linguistic and genetic research pointing to a most probable North African Berber source. However, neither agreement about the exact point of origin nor a model for the indigenous colonization of the islands has been established. To shed light on these questions, we analyzed 48 ancient mitogenomes from 25 archaeological sites from the seven main islands. Most lineages observed in the ancient samples have a Mediterranean distribution, and belong to lineages associated with the Neolithic expansion in the Near East and Europe (T2c, J2a, X3a…). This phylogeographic analysis of Canarian ancient mitogenomes, the first of its kind, shows that some lineages are restricted to Central North Africa (H1cf, J2a2d and T2c1d3), while others have a wider distribution, including both West and Central North Africa, and, in some cases, Europe and the Near East (U6a1a1, U6a7a1, U6b, X3a, U6c1). In addition, we identify four new Canarian-specific lineages (H1e1a9, H4a1e, J2a2d1a and L3b1a12) whose coalescence dates correlate with the estimated time for the colonization of the islands (1st millennia CE). Additionally, we observe an asymmetrical distribution of mtDNA haplogroups in the ancient population, with certain haplogroups appearing more frequently in the islands closer to the continent. This reinforces results based on modern mtDNA and Y-chromosome data, and archaeological evidence suggesting the existence of two distinct migrations. Comparisons between insular populations show that some populations had high genetic diversity, while others were probably affected by genetic drift and/or bottlenecks. In spite of observing interinsular differences in the survival of indigenous lineages, modern populations, with the sole exception of La Gomera, are homogenous across the islands, supporting the theory of extensive human mobility after the European conquest., Competing Interests: Tibicena Arqueología y Patrimonio provided support in the form of salaries for authors V.A. and M.A.M.B and radiocarbon dating of four archaeological samples. This does not alter our adherence to PLOS ONE policies on sharing data and materials.

Despite the genetic resemblance of Canary Islanders to other southern European populations, their geographical isolation and the historical admixture of aborigines (from North Africa) with sub-Saharan Africans and Europeans have shaped a distinctive genetic makeup that likely affects disease susceptibility and health disparities. Based on single nucleotide polymorphism array data and whole genome sequencing (30×), we inferred that the last African admixture took place ∼14 generations ago and estimated that up to 34% of the Canary Islander genome is of recent African descent. The length of regions in homozygosis and the ancestry-related mosaic organization of the Canary Islander genome support the view that isolation has been strongest on the two smallest islands. Furthermore, several genomic regions showed significant and large deviations in African or European ancestry and were significantly enriched in genes involved in prevalent diseases in this community, such as diabetes, asthma, and allergy. The most prominent of these regions were located near LCT and the HLA, two well-known targets of selection, at which 40‒50% of the Canarian genome is of recent African descent according to our estimates. Putative selective signals were also identified in these regions near the SLC6A11-SLC6A1, KCNMB2, and PCDH20-PCDH9 genes. Taken together, our findings provide solid evidence of a significant recent African admixture, population isolation, and adaptation in this part of Europe, with the favoring of African alleles in some chromosome regions. These findings may have medical implications for populations of recent African ancestry.