Critical illness in COVID-19 is an extreme and clinically homogeneous disease phenotype that we have previously shown1 to be highly efficient for discovery of genetic associations2. Despite the advanced stage of illness at presentation, we have shown that host genetics in patients who are critically ill with COVID-19 can identify immunomodulatory therapies with strong beneficial effects in this group3. Here we analyse 24,202 cases of COVID-19 with critical illness comprising a combination of microarray genotype and whole-genome sequencing data from cases of critical illness in the international GenOMICC (11,440 cases) study, combined with other studies recruiting hospitalized patients with a strong focus on severe and critical disease: ISARIC4C (676 cases) and the SCOURGE consortium (5,934 cases). To put these results in the context of existing work, we conduct a meta-analysis of the new GenOMICC genome-wide association study (GWAS) results with previously published data. We find 49 genome-wide significant associations, of which 16 have not been reported previously. To investigate the therapeutic implications of these findings, we infer the structural consequences of protein-coding variants, and combine our GWAS results with gene expression data using a monocyte transcriptome-wide association study (TWAS) model, as well as gene and protein expression using Mendelian randomization. We identify potentially druggable targets in multiple systems, including inflammatory signalling (JAK1), monocyte–macrophage activation and endothelial permeability (PDE4A), immunometabolism (SLC2A5 and AK5), and host factors required for viral entry and replication (TMPRSS2 and RAB2A).

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

IntroductionHypertrophic cardiomyopathy (HCM) is the most prevalent inherited heart disease. Next-generation sequencing (NGS) is the preferred genetic test, but the diagnostic value of screening for minor and candidate genes, and the role of copy number variants (CNVs) deserves further evaluation.MethodsThree hundred and eighty-seven consecutive unrelated patients with HCM were screened for genetic variants in the 5 most frequent genes (MYBPC3, MYH7, TNNT2, TNNI3 and TPM1) using Sanger sequencing (N = 84) or NGS (N = 303). In the NGS cohort we analyzed 20 additional minor or candidate genes, and applied a proprietary bioinformatics algorithm for detecting CNVs. Additionally, the rate and classification of TTN variants in HCM were compared with 427 patients without structural heart disease.ResultsThe percentage of patients with pathogenic/likely pathogenic (P/LP) variants in the main genes was 33.3%, without significant differences between the Sanger sequencing and NGS cohorts. The screening for 20 additional genes revealed LP variants in ACTC1, MYL2, MYL3, TNNC1, GLA and PRKAG2 in 12 patients. This approach resulted in more inconclusive tests (36.0% vs. 9.6%, p<0.001), mostly due to variants of unknown significance (VUS) in TTN. The detection rate of rare variants in TTN was not significantly different to that found in the group of patients without structural heart disease. In the NGS cohort, 4 patients (1.3%) had pathogenic CNVs: 2 deletions in MYBPC3 and 2 deletions involving the complete coding region of PLN.ConclusionsA small percentage of HCM cases without point mutations in the 5 main genes are explained by P/LP variants in minor or candidate genes and CNVs. Screening for variants in TTN in HCM patients drastically increases the number of inconclusive tests, and shows a rate of VUS that is similar to patients without structural heart disease, suggesting that this gene should not be analyzed for clinical purposes in HCM.

The serrated polyposis syndrome (SPS) is the most common and yet underdiagnosed colorectal polyposis syndrome. It is characterized by multiple and/or large colonic serrated polyps and a higher associated risk for colorectal cancer (CRC). The main objective of this study was to identify new candidate genes involved in the germline predisposition to SPS/CRC. Thirty-nine SPS patients from 16 families (≥2 patients per family) were recruited without alterations in well-known hereditary CRC genes, and germline and somatic whole-exome sequencing were performed. Germline rare variants with plausible pathogenicity, located in genes involved in cancer development, senescence and epigenetic regulation were selected. Somatic mutational profiling and signature analysis was pursued in one sample per family, when possible. After data filtering, ANXA10, ASXL1, CFTR, DOT1L, HIC1, INO80, KLF3, MCM3AP, MCM8, PDLIM2, POLD1, TP53BP1, WNK2 and WRN were highlighted as the more promising candidate genes for SPS germline predisposition with potentially pathogenic variants shared within families. Somatic analysis characterized mutational profiles in advanced serrated polyps/tumors, revealing a high proportion of hypermutated samples, with a prevalence of clock-like mutational signatures in most samples and the presence of DNA mismatch repair-defective signatures in some cases. In conclusion, we identified new candidate genes to be involved in familial SPS. Further functional studies and replication in additional cohorts are required to confirm the selected candidates.

This is the second reported case of a child with holoprosencephaly and trisomy 21. The first case was born to a diabetic woman; in our case, there was no evidence of diabetes in the mother. Most of the distinctive facial features of Down syndrome were obscured by the presence of cyclopia and a supraorbital proboscis in this infant. The relevance of chromosome analysis in cases with holoprosencephaly is discussed.