Novel genes and sex differences in COVID-19 severity.

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.

Rearrangement of cell types in the rat carotid body neurogenic niche induced by chronic intermittent hypoxia.

The carotid body (CB) is a prototypical acute oxygen (O2 )-sensing organ that mediates reflex hyperventilation and increased cardiac output in response to hypoxaemia. CB overactivation, secondary to the repeated stimulation produced by the recurrent episodes of intermittent hypoxia, is believed to contribute to the pathogenesis of sympathetic hyperactivity present in sleep apnoea patients. Although CB functional plasticity induced by chronic intermittent hypoxia (CIH) has been demonstrated, the underlying mechanisms are not fully elucidated. Here, we show that CIH induces a small increase in CB volume and rearrangement of cell types in the CB, characterized by a mobilization of immature quiescent neuroblasts, which enter a process of differentiation into mature, O2 -sensing and neuron-like, chemoreceptor glomus cells. Prospective isolation of individual cell classes has allowed us to show that maturation of CB neuroblasts is paralleled by an upregulation in the expression of specific glomus cell genes involved in acute O2 -sensing. CIH enhances mitochondrial responsiveness to hypoxia in maturing neuroblasts as well as in glomus cells. These data provide novel perspectives on the pathogenesis of CB-mediated sympathetic overflow that may lead to the development of new pharmacological strategies of potential applicability in sleep apnoea patients. KEY POINTS: Obstructive sleep apnoea is a frequent condition in the human population that predisposes to severe cardiovascular and metabolic alterations. Activation of the carotid body, the main arterial oxygen-sensing chemoreceptor, by repeated episodes of hypoxaemia induces exacerbation of the carotid body-mediated chemoreflex and contributes to sympathetic overflow characteristic of sleep apnoea patients. In rats, chronic intermittent hypoxaemia induces fast neurogenesis in the carotid body with rapid activation of neuroblasts, which enter a process of proliferation and maturation into O2 -sensing chemoreceptor glomus cells. Maturing carotid body neuroblasts and glomus cells exposed to chronic intermittent hypoxia upregulate genes involved in acute O2 sensing and enhance mitochondrial responsiveness to hypoxia. These findings provide novel perspectives on the pathogenesis of carotid body-mediated sympathetic hyperactivation. Pharmacological modulation of carotid body fast neurogenesis could help to ameliorate the deleterious effects of chronic intermittent hypoxaemia in sleep apnoea patients.

Evaluation of a whole-exome sequencing pipeline and benchmarking of causal germline variant prioritizers.

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.

Wilson Disease Prevalence: Discrepancy Between Clinical Records, Registries and Mutation Carrier Frequency.

OBJECTIVES: Diagnosis of Wilson disease (WD) is difficult and, as early detection may prevent all symptoms, it is essential to know the exact prevalence to evaluate the cost-efficacy of a screening program. As the number of WD patients was high in our population, we wished to estimate prevalence by determining the carrier frequency for clinically relevant ATP7B mutations. METHODS: To estimate prevalence, screening for the most prevalent mutation was performed in 1661 individuals with ancestry in Gran Canaria, and the frequency of other mutations was estimated from patient records. Alternatively, ATP7B mutations were detected from exomes and genomes from 851 individuals with Canarian ancestry, 236 from Gran Canaria, and a public Spanish exome database. RESULTS: Estimated carrier frequencies in Gran Canaria ranged from 1 in 20 to 28, depending on the method used, resulting in prevalences of 1 case per 1547 to 3140 inhabitants. Alternatively, the estimated affected frequencies were 1 in 5985 to 7980 and 1 in 6278 to 16,510 in the archipelago or mainland Spain respectively. CONCLUSIONS: The number of carriers predicts much higher prevalences than reported, suggesting that WD is underdiagnosed; specific mutations may remain unnoticed due to low penetrance or no signs of disease at all; regional prevalence rather than national prevalence should be considered in cost-efficacy models to approach preventive screening in the asymptomatic population and genetic screening strategies will have to deal with the genetic heterogeneity of ATP7B in the general population and in patients.

Transactive Response DNA-Binding Protein (TARDBP/TDP-43) Regulates Cell Permissivity to HIV-1 Infection by Acting on HDAC6.

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.

Admixture mapping of asthma in southwestern Europeans with North African ancestry influences.

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.

Interactive Web-Based Resource for Annotation of Genetic Variants Causing Hereditary Angioedema (HADA): Database Development, Implementation, and Validation.

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.

Genomic Analyses of Human European Diversity at the Southwestern Edge: Isolation, African Influence and Disease Associations in the Canary Islands.

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.

PKCα-Mediated Downregulation of RhoA Activity in Depolarized Vascular Smooth Muscle: Synergistic Vasorelaxant Effect of PKCα and ROCK Inhibition.

BACKGROUND/AIMS: Protein kinase C (PKC)- and RhoA/Rho-associated kinase (ROCK) play important roles in arterial sustained contraction. Although depolarization-elicited RhoA/ROCK activation is accepted, the role of PKC in depolarized vascular smooth muscle cells (VSMCs) is a subject of controversy. Our aim was to study the role of PKC in arterial contraction and its interaction with RhoA/ROCK. METHODS: Mass spectrometry was used to identify the PKC isoenzymes. PKCα levels and RhoA activity were analyzed by western blot and G-LISA, respectively, and isometric force was measured in arterial rings. RESULTS: In depolarized VSMCs RhoA and PKCα were translocated to the plasma membrane, where they colocalize and coimmunoprecipitate. Interestingly, depolarization-induced RhoA activation was downregulated by PKCα, effect reverted by PKCα inhibition. Phorbol 12,13-dibutyrate (PDBu) induced the translocation of PKCα to the plasma membrane, increased the level of RhoA in the cytosol and reduced RhoA/ROCK activity. These effects were reverted when PKC was inhibited. Pharmacological or siRNA inhibition of PKCα synergistically potentiated the vasorelaxant effect of RhoA/ROCK inhibition. CONCLUSION: The present study provides the first evidence that RhoA activity is downregulated by PKCα in depolarized and PDBu treated freshly isolated VSMCs and arteries, with an important physiological role on arterial contractility.

Relation of RhoA in Peripheral Blood Mononuclear Cells With Severity of Aneurysmal Subarachnoid Hemorrhage and Vasospasm.

BACKGROUND AND PURPOSE: Rho-kinase, an effector of RhoA, is associated with various cardiovascular diseases in circulating blood cells. However, the role of RhoA/Rho-kinase in peripheral blood mononuclear cells from patients with spontaneous aneurysmal subarachnoid hemorrhage (aSAH) has not yet been studied in relation to the severity of this disease. Therefore, we analyzed the expression and activity of RhoA as a possible biomarker in aSAH. METHODS: Twenty-four patients with aSAH and 15 healthy subjects were examined. Peripheral blood mononuclear cells were collected, and RhoA activity and expression were determined by RhoA activation assay kit (G-LISA) and enzyme-linked immunosorbent assay tests, respectively. The severity of aSAH was determined from the World Federation of Neurological Surgeon scale, and vasospasm was evaluated using clinical symptoms, arteriography, and sonography. RESULTS: RhoA expression was significantly increased in peripheral blood mononuclear cells from patients on days 0, 2, and 4 after aSAH versus healthy subjects (P=0.036, 0.010, and 0.018, respectively, by U Mann-Whitney analysis). There was a significant correlation between RhoA expression and injury severity on days 2 and 4 (Spearman test, day 2: r=0.682, n=14, P=0.007; day 4: r=0.721, n=14, P=0.004). No significant correlation was observed on day 0 (day 0: r=0.131, n=6, P=0.805). Active RhoA was not significantly different in patients and healthy subjects on days 0, 2, and 4 (P=0.243, 0.222, and 0.600, respectively) nor did it increase significantly on days 0 and 2 in patients with vasospasm versus patients without vasospasm (P=0.064 and 0.519, respectively). In contrast, active RhoA was significantly higher on day 4 in patients who developed vasospasm versus patients without vasospasm (P=0.028). CONCLUSIONS: Our preliminary results indicate that RhoA expression and activity in peripheral blood mononuclear cells might be related with aSAH severity and cerebral vasospasm. RhoA is a potential biomarker of the risks associated with aSAH.

Plasma membrane insertion of epithelial sodium channels occurs with dual kinetics.

The epithelial sodium channel (ENaC) constitutes the rate-limiting step for Na(+) transport across electrically tight epithelia. Regulation of ENaC activity is critical for electrolyte and extracellular volume homeostasis, as well as for lung liquid clearance and colon Na(+) handling. ENaC activity is tightly controlled by a combination of mechanisms involving changes in open probability and plasma membrane abundance. The latter reflects a combination in channel biosynthesis and trafficking to and from the membrane. Studying ENaC trafficking with different techniques in a variety of expression systems has yielded inconsistent results, indicating either fast or slow rates of insertion and retrieval, which range from the order of minutes to several hours. Here, we use Xenopus oocytes as ENaC expression system to study channel insertion rate in the membrane using two different techniques under comparable conditions: (1) confocal microscopy coupled to fluorescence recovery after photobleaching (FRAP) measurements; and (2) fluorescent bungarotoxin (BTX) binding to ENaC subunits modified to include BTX binding sites (BBSs) in their extracellular domain, a technique that has not been previously used to study ENaC trafficking. Our confocal-FRAP data indicate a fast rate of ENaC incorporation to the membrane in a process conditioned by channel subunit composition. On the other hand, BTX binding experiments indicate much slower channel insertion rates, with matching slow ENaC retrieval rates. The data support a model that includes fast recycling of endocytosed ENaC with parallel incorporation of newly synthesized channels at a slower rate.

Molecular components of nitrate and nitrite efflux in yeast.

Some eukaryotes, such as plant and fungi, are capable of utilizing nitrate as the sole nitrogen source. Once transported into the cell, nitrate is reduced to ammonium by the consecutive action of nitrate and nitrite reductase. How nitrate assimilation is balanced with nitrate and nitrite efflux is unknown, as are the proteins involved. The nitrate assimilatory yeast Hansenula polymorpha was used as a model to dissect these efflux systems. We identified the sulfite transporters Ssu1 and Ssu2 as effective nitrate exporters, Ssu2 being quantitatively more important, and we characterize the Nar1 protein as a nitrate/nitrite exporter. The use of strains lacking either SSU2 or NAR1 along with the nitrate reductase gene YNR1 showed that nitrate reductase activity is not required for net nitrate uptake. Growth test experiments indicated that Ssu2 and Nar1 exporters allow yeast to cope with nitrite toxicity. We also have shown that the well-known Saccharomyces cerevisiae sulfite efflux permease Ssu1 is also able to excrete nitrite and nitrate. These results characterize for the first time essential components of the nitrate/nitrite efflux system and their impact on net nitrate uptake and its regulation.

The neuronal serum- and glucocorticoid-regulated kinase 1.1 reduces neuronal excitability and protects against seizures through upregulation of the M-current.

The M-current formed by tetramerization of Kv7.2 and Kv7.3 subunits is a neuronal voltage-gated K(+) conductance that controls resting membrane potential and cell excitability. In Xenopus laevis oocytes, an increase in Kv7.2/3 function by the serum- and glucocorticoid-regulated kinase 1 (SGK1) has been reported previously (Schuetz et al., 2008). We now show that the neuronal isoform of this kinase (SGK1.1), with distinct subcellular localization and modulation, upregulates the Kv7.2/3 current in Xenopus oocytes and mammalian human embryonic kidney HEK293 cells. In contrast to the ubiquitously expressed SGK1, the neuronal isoform SGK1.1 interacts with phosphoinositide-phosphatidylinositol 4,5-bisphosphate (PIP(2)) and is distinctly localized to the plasma membrane (Arteaga et al., 2008). An SGK1.1 mutant with disrupted PIP(2) binding sites produced no effect on Kv7.2/3 current amplitude. SGK1.1 failed to modify the voltage dependence of activation and did not change activation or deactivation kinetics of Kv7.2/3 channels. These results suggest that the kinase increases channel membrane abundance, which was confirmed with flow cytometry assays. To evaluate the effect of the kinase in neuronal excitability, we generated a transgenic mouse (Tg.sgk) expressing a constitutively active form of SGK1.1 (S515D). Superior cervical ganglion (SCG) neurons isolated from Tg.sgk mice showed a significant increase in M-current levels, paralleled by reduced excitability and more negative resting potentials. SGK1.1 effect on M-current in Tg.sgk-SCG neurons was counteracted by muscarinic receptor activation. Transgenic mice with increased SGK1.1 activity also showed diminished sensitivity to kainic acid-induced seizures. Altogether, our results unveil a novel role of SGK1.1 as a physiological regulator of the M-current and neuronal excitability.

Genomic epidemiology of the primary methicillin-resistant Staphylococcus aureus clones causing invasive infections in Paraguayan children.

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.

Novel risk loci for COVID-19 hospitalization among admixed American populations.

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.

Benchmarking of human Y-chromosomal haplogroup classifiers with whole-genome and whole-exome sequence data.

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.

Reinfection rate and disease severity of the BA.5 Omicron SARS-CoV-2 lineage compared to previously circulating variants of concern in the Canary Islands (Spain).

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

Digging into the admixture strata of current-day Canary Islanders based on mitogenomes.

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.

Bioinformatic approaches to draft the viral genome sequence of Canary Islands cases related to the multicountry mpox virus 2022-outbreak.

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.