The major genetic risk factor for severe COVID-19 is inherited from Neanderthals

Data-driven mathematical modeling approaches for COVID-19: A survey

In this review, we successively present the methods for phenomenological modeling of the evolution of reported and unreported cases of COVID-19, both in the exponential phase of growth and then in a complete epidemic wave. After the case of an isolated wave, we present the modeling of several successive waves separated by endemic stationary periods. Then, we treat the case of multi-compartmental models without or with age structure. Eventually, we review the literature, based on 260 articles selected in 11 sections, ranging from the medical survey of hospital cases to forecasting the dynamics of new cases in the general population. This review favors the phenomenological approach over the mechanistic approach in the choice of references and provides simulations of the evolution of the number of observed cases of COVID-19 for 10 states (California, China, France, India, Israel, Japan, New York, Peru, Spain and United Kingdom).

Recurrent gene flow between Neanderthals and modern humans over the past 200,000 years

Although it is well known that the ancestors of modern humans and Neanderthals admixed, the effects of gene flow on the Neanderthal genome are not well understood. We develop methods to estimate the amount of human-introgressed sequences in Neanderthals and apply it to whole-genome sequence data from 2000 modern humans and three Neanderthals. We estimate that Neanderthals have 2.5 to 3.7% human ancestry, and we leverage human-introgressed sequences in Neanderthals to revise estimates of Neanderthal ancestry in modern humans, show that Neanderthal population sizes were significantly smaller than previously estimated, and identify two distinct waves of modern human gene flow into Neanderthals. Our data provide insights into the genetic legacy of recurrent gene flow between modern humans and Neanderthals.

Genomic Landscape of Susceptibility to Severe COVID-19 in the Slovenian Population

Determining the genetic contribution of susceptibility to severe SARS-CoV-2 infection outcomes is important for public health measures and individualized treatment. Through intense research on this topic, several hundred genes have been implicated as possibly contributing to the severe infection phenotype(s); however, the findings are complex and appear to be population-dependent. We aimed to determine the contribution of human rare genetic variants associated with a severe outcome of SARS-CoV-2 infections and their burden in the Slovenian population. A panel of 517 genes associated with severe SARS-CoV-2 infection were obtained by combining an extensive review of the literature, target genes identified by the COVID-19 Host Genetic Initiative, and the curated Research COVID-19 associated genes from PanelApp, England Genomics. Whole genome sequencing was performed using PCR-free WGS on DNA from 60 patients hospitalized due to severe COVID-19 disease, and the identified rare genomic variants were analyzed and classified according to the ACMG criteria. Background prevalence in the general Slovenian population was determined by comparison with sequencing data from 8025 individuals included in the Slovenian genomic database (SGDB). Results show that several rare pathogenic/likely pathogenic genomic variants in genes CFTR, MASP2, MEFV, TNFRSF13B, and RNASEL likely contribute to the severe infection outcomes in our patient cohort. These results represent an insight into the Slovenian genomic diversity associated with a severe COVID-19 outcome.

Basic implications on three pathways associated with SARS-CoV-2

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) interacts between the host and virus and govern induction, resulting in multiorgan impacts. Its pathophysiology involves the followings: 1) the angiotensin-converting enzyme (ACE2) and Toll-like receptor (TLR) pathways: 2) the neuropilin (NRP) pathway: 3) the spike protein pathway. Therefore, it is necessary to block the pathological course with modulating innate lymphoid cells against diverse corona variants in the future.

Critically-ill COVID-19 susceptibility gene CCR3 shows natural selection in sub-Saharan Africans

The prevalence of COVID-19 critical illness varies across ethnicities, with recent studies suggesting that genetic factors may contribute to this variation. The aim of this study was to investigate natural selection signals of genes associated with critically-ill COVID-19 in sub-Saharan Africans. Severe COVID-19 SNPs were obtained from the HGI website. Selection signals were assessed in 661 sub-Sahara Africans from 1000 Genomes Project using integrated haplotype score (iHS), cross-population extended haplotype homozygosity (XP-EHH), and fixation index (Fst). Allele frequency trajectory analysis of ancient DNA samples were used to validate the existing of selection in sub-Sahara Africans. We also used Mendelian randomization to decipher the correlation between natural selection and critically-ill COVID-19. We identified that CCR3 exhibited significant natural selection signals in sub-Sahara Africans. Within the CCR3 gene, rs17217831-A showed both high iHS (Standardized iHS = 2) and high XP-EHH (Standardized XP-EHH = 2.5) in sub-Sahara Africans. Allele frequency trajectory of CCR3 rs17217831-A revealed natural selection occurring in the recent 1,500 years. Natural selection resulted in increased CCR3 expression in sub-Sahara Africans. Mendelian Randomization provided evidence that increased blood CCR3 expression and eosinophil counts lowered the risk of critically ill COVID-19. Our findings suggest that sub-Saharan Africans are resistant to critically ill COVID-19 due to natural selection and identify CCR3 as a potential novel therapeutic target.

The association of rs17713054 with Neanderthal origin at 3p21.31 locus with the severity of COVID-19 in Iranian patients

Since the COVID-19 pandemic, the diversity of clinical manifestations in patients has been a tremendous challenge. It seems that genetic variations, as one of the players, contribute to the variety of symptoms. Genome-wide association studies have demonstrated the influence of certain genomic regions on the disease prognosis. Particularly, a haplotype at 3p21.31 locus, inherited from Neanderthals, showed an association with COVID-19 severity. Despite several studies regarding this haplotype, some key variants are not sufficiently addressed. In the present study, we investigated the association of rs17713054 at 3p21.31 with COVID-19 severity. We analyzed the genotype of 251 Iranian COVID-19 patients (151 patients with asymptomatic to mild form as control and 100 patients with severe to critical symptoms without any comorbidities as case group) using the ARMS-PCR method. Results demonstrated that the A allele confers an almost twofold increased risk for COVID-19 severity (P value = 0.008). The AA genotype also raises the risk by more than 11 times following the recessive model (P value = 0.013). In conclusion, the A allele in rs17713054 was a risk allele in Iranian patients and was independently associated with COVID-19 severity. More studies are beneficial to confirm these findings in other populations and to develop strategies for risk assessment, prevention, and personalized medicine.

Difference in TMPRSS2 usage by Delta and Omicron variants of SARS-CoV-2: Implication for a sudden increase among children

It has been postulated from a combination of evidence that a sudden increase in COVID-19 cases among pediatric patients after onset of the Omicron wave was attributed to a reduced requirement for TMPRSS2-mediated entry in pediatric airways with lower expression levels of TMPRSS2. Epidemic strains were isolated from the indigenous population in an area, and the levels of TMPRSS2 required for Delta and Omicron variants were assessed. As a result, Delta variants proliferated fully in cultures of TMPRSS2-positive Vero cells but not in TMPRSS2-negative Vero cell culture (350-fold, Delta vs 9.6-fold, Omicron). There was no obvious age-dependent selection of Omicron strains affected by the TMPRSS2 (9.6-fold, Adults vs. 12-fold, Children). A phylogenetic tree was generated and Blast searches (up to 100 references) for the spread of strains in the study area showed that each strain had almost identical homology (>99.5%) with foreign isolates, although indigenous strains had obvious differences from each other. This suggested that the differences had been present abroad for a long period. Therefore, the lower requirement for TMPRSS2 by Omicron strains might be applicable to epidemic strains globally. In conclusion, the property of TMPRSS2-independent cleavage makes Omicron proliferate with ease and allows epidemics among children with fewer TMPRSS2 on epithelial surfaces of the respiratory organs.

Genomic analysis of severe COVID-19 considering or not asthma comorbidity: GWAS insights from the BQC19 cohort

BACKGROUND

The severity of COVID-19 is influenced by various factors including the presence of respiratory diseases. Studies have indicated a potential relationship between asthma and COVID-19 severity.

OBJECTIVE

This study aimed to conduct a genome-wide association study (GWAS) to identify genetic and clinical variants associated with the severity of COVID-19, both among patients with and without asthma.

METHODS

We analyzed data from 2131 samples sourced from the Biobanque québécoise de la COVID-19 (BQC19), with 1499 samples from patients who tested positive for COVID-19. Among these, 1110 exhibited mild-to-moderate symptoms, 389 had severe symptoms, and 58 had asthma. We conducted a comparative analysis of clinical data from individuals in these three groups and GWAS using a logistic regression model. Phenotypic data analysis resulted in the refined covariates integrated into logistic models for genetic studies.

RESULTS

Considering a significance threshold of 1 × 10-6, seven genetic variants were associated with severe COVID-19. These variants were located proximal to five genes: sodium voltage-gated channel alpha subunit 1 (SCN10A), desmoplakin (DSP), RP1 axonemal microtubule associated (RP1), IGF like family member 1 (IGFL1), and docking protein 5 (DOK5). The GWAS comparing individuals with severe COVID-19 with asthma to those without asthma revealed four genetic variants in transmembrane protein with EGF like and two follistatin like domains 2 (TMEFF2) and huntingtin interacting protein-1 (HIP1) genes.

CONCLUSION

This study provides significant insights into the genetic profiles of patients with severe forms of the disease, whether accompanied by asthma or not. These findings enhance our comprehension of the genetic factors that affect COVID-19 severity.

KEY MESSAGES

Seven genetic variants were associated with the severe form of COVID-19; Four genetic variants were associated with the severe form of COVID-19 in individuals with comorbid asthma; These findings help define the genetic component of the severe form of COVID-19 in relation to asthma as a comorbidity.

Evaluation of IFNAR2 and TYK2 transcripts' prognostic role in COVID-19 patients: a retrospective study

Background and objectives

This study aimed to investigate the possible prognostic significance of interferon alpha-beta receptor subunit 2 (IFNAR2) and tyrosine kinase 2 (TYK2) expressions.

Methods

We conducted a retrospective study including COVID-19 adult patients. All blood samples were collected before any interventions. The expressions of IFNAR2 and TYK2 were assessed using real-time PCR in venous blood samples of 54 cases and 56 controls. The transcript quantities of IFNAR2 and TYK2 genes were assessed using a Delta-Ct method.

Results

Our findings show no significant differences in gene expression levels for IFNAR2 and TYK2 between patients who required oxygen (O2) therapy and those who did not (p-value = 0.732 and p-value = 0.629, respectively). Likewise, there were no significant differences in IFNAR2 and TYK2 expressions between patients hospitalized for less than 7 days and those hospitalized for 7 days or more (p-value = 0.455 and p-value = 0.626, respectively). We also observed a weak correlation between IFNAR2 expression and CRP (p-value = 0.045, r = 0.192). There was a negative correlation between the expression levels of IFNAR2 and TYK2 transcripts in COVID-19 patients (p-value = 0.044; partial correlation coefficient = -0.283). Additionally, IFNAR2 and TYK2 were significantly downregulated in the COVID-19 group compared to healthy subjects (p-value = 0.002 and p-value = 0.028, respectively). However, neither IFNAR2 nor TYK2 expression was significantly different between the case subgroups based on COVID-19 severity. The IFNAR2 ΔΔCt (B = -0.184, 95% CI: -0.524-0.157, p-value = 0.275) and the TYK2 ΔΔCt (B = 0.114, 95% CI: -0.268-0.496, p-value = 0.543) were not found to be significant predictors of hospitalization duration. The area under the curve (AUC) for IFNAR2 expression is 0.655 (p-value = 0.005, 95% CI: 0.554-0.757), suggesting its poor discriminative value.

Conclusion

We were unable to comment definitively on the prognostic power of IFNAR2 and TYK2 expressions in COVID-19 patients, and larger-scale studies are needed. The principal limitations of this study included the lack of longitudinal analysis and limited sample size.

Using Portuguese BRCA pathogenic variation as a model to study the impact of human admixture on human health

BACKGROUND

Admixture occurs between different ethnic human populations. The global colonization in recent centuries by Europeans led to the most significant admixture in human history. While admixture may enhance genetic diversity for better fitness, it may also impact on human health by transmitting genetic variants for disease susceptibility in the admixture population. The admixture by Portuguese global exploration initiated in the 15th century has reached over 20 million of Portuguese-heritage population worldwide. It provides a valuable model to study the impact of admixture on human health. BRCA1 and BRCA2 (BRCA) are two of the important tumor suppressor genes. The pathogenic variation (PV) in BRCA is well determined to cause high risk of hereditary breast and ovarian cancer. Tracing the distribution of Portuguese BRCA PV in Portuguese-heritage population will help to understand the impact of admixture on cancer susceptibility in modern humans. In this study, we analyzed the distribution of the Portuguese-originated BRCA variation in Brazilian population, which has high degree Portuguese-heritage.

METHODS

By comprehensive data mining, standardization and annotation, we generated a Portuguese-derived BRCA variation dataset and a Brazilian-derived BRCA variation dataset. We compared the two BRCA variation datasets to identify the BRCA variants shared between the two populations.

RESULTS

The Portuguese-derived BRCA variation dataset consists of 220 BRCA variants including 78 PVs from 11,482 Portuguese cancer patients, 93 (42.2%) in BRCA1 and 127 (57.7%) in BRCA2. Of the 556 Portuguese BRCA PV carriers carrying the 78 PVs, 331 (59.5%) carried the three Portuguese-BRCA founder PVs of BRCA1 c.2037delinsCC, BRCA1 c.3331_3334del and BRCA2 c.156_157insAlu. The Brazilian-derived BRCA variation dataset consists of 255 BRCA PVs from 7,711 cancer patients, 136 (53.3%) in BRCA1 and 119 (46.6%) in BRCA2. We developed an open database named dbBRCA-Portuguese ( https://genemutation.fhs.um.edu.mo/dbbrca-portuguese/ ) and an open database named dbBRCA-Brazilian ( https://genemutation.fhs.um.edu.mo/dbbrca-brazilian ) to host the BRCA variation data from Portuguese and Brazilian populations. We compared the BRCA PV datasets between Portuguese and Brazilian populations, and identified 29 Portuguese-specific BRCA PVs shared between Portuguese and Brazilian populations, 14 in BRCA1 including the Portuguese founder BRCA1 c.3331_3334del and BRCA1 c.2037delinsCC, and 15 in BRCA2 including the Portuguese founder BRCA2 c.156_157insAlu. Searching the 78 Portuguese BRCA PVs in over 5,000 ancient human genomes identified evolution origin for only 8 PVs in Europeans dated between 37,470 and 3,818 years before present, confirming the Portuguese-specificity of Portuguese BRCA PVs; comparing the 78 Portuguese BRCA PVs Portuguese, 255 Brazilian BRCA PVs, and 134 African BRCA PVs showed little overlapping, ruling out the possibility that the BRCA PVs shared between Portuguese and Brazilian may also be contributed by African.

CONCLUSION

Our study provides evidence that the admixture in recent human history contributed to cancer susceptibility in modern humans.

Prevalence of the protective OAS1 rs10774671-G allele against severe COVID-19 in Moroccans: implications for a North African Neanderthal connection

The clinical presentation of COVID-19 shows high variability among individuals, which is partly due to genetic factors. The OAS1/2/3 cluster has been found to be strongly associated with COVID-19 severity. We examined this locus in the Moroccan population for the occurrence of the critical variant rs10774671 and its respective haplotype blocks. The frequency of single-nucleotide polymorphisms (SNPs) in the cluster of OAS immunity genes in 157 unrelated individuals of Moroccan origin was determined using an in-house exome database. OAS1 exon 6 of 71 SARS-CoV-2-positive individuals with asymptomatic/mild disease and 74 with moderate/severe disease was sequenced by the Sanger method. The genotypic, allelic, and haplotype frequencies of three SNPs were compared between these two groups. Finally, males in our COVID-19 series were genotyped for the Berber-specific marker E-M81. The prevalence of the OAS1 rs10774671-G allele in present-day Moroccans was found to be 40.4%, which is similar to that found in Europeans. However, it was found equally in both the Neanderthal GGG haplotype and the African GAC haplotype, with a frequency of 20% each. These two haplotypes, and hence the rs10774671-G allele, were significantly associated with protection against severe COVID-19 (p = 0.034, p = 0.041, and p = 0.008, respectively). Surprisingly, in men with the Berber-specific uniparental markers, the African haplotype was absent, while the prevalence of the Neanderthal haplotype was similar to that in Europeans. The protective rs10774671-G allele of OAS1 was found only in the Neanderthal haplotype in Berbers, the indigenous people of North Africa, suggesting that this region may have served as a stepping-stone for the passage of hominids to other continents.

Decoding triancestral origins, archaic introgression, and natural selection in the Japanese population by whole-genome sequencing

We generated Japanese Encyclopedia of Whole-Genome/Exome Sequencing Library (JEWEL), a high-depth whole-genome sequencing dataset comprising 3256 individuals from across Japan. Analysis of JEWEL revealed genetic characteristics of the Japanese population that were not discernible using microarray data. First, rare variant-based analysis revealed an unprecedented fine-scale genetic structure. Together with population genetics analysis, the present-day Japanese can be decomposed into three ancestral components. Second, we identified unreported loss-of-function (LoF) variants and observed that for specific genes, LoF variants appeared to be restricted to a more limited set of transcripts than would be expected by chance, with PTPRD as a notable example. Third, we identified 44 archaic segments linked to complex traits, including a Denisovan-derived segment at NKX6-1 associated with type 2 diabetes. Most of these segments are specific to East Asians. Fourth, we identified candidate genetic loci under recent natural selection. Overall, our work provided insights into genetic characteristics of the Japanese population.

A Neanderthal haplotype introgressed into the human genome confers protection against membranous nephropathy

Class 2 HLA and PLA2R1 alleles are exceptionally strong genetic risk factors for membranous nephropathy (MN), leading, through an unknown mechanism, to a targeted autoimmune response. Introgressed archaic haplotypes (introduced from an archaic human genome into the modern human genome) might influence phenotypes through gene dysregulation. Here, we investigated the genomic region surrounding the PLA2R1 gene. We reconstructed the phylogeny of Neanderthal and modern haplotypes in this region and calculated the probability of the observed clustering being the result of introgression or common descent. We imputed variants for the participants in our previous genome-wide association study and we compared the distribution of Neanderthal variants between MN cases and controls. The region associated with the lead MN risk locus in the PLA2R1 gene was confirmed and showed that, within a 507 kb region enriched in introgressed sequence, a stringently defined 105 kb haplotype, intersecting the coding regions for PLA2R1 and ITGB6, is inherited from Neanderthals. Thus, introgressed Neanderthal haplotypes overlapping PLA2R1 are differentially represented in MN cases and controls, with enrichment In controls suggesting a protective effect.

The inflammatory microenvironment of the lung at the time of infection governs innate control of SARS-CoV-2 replication

SARS-CoV-2 infection leads to vastly divergent clinical outcomes ranging from asymptomatic infection to fatal disease. Co-morbidities, sex, age, host genetics and vaccine status are known to affect disease severity. Yet, how the inflammatory milieu of the lung at the time of SARS-CoV-2 exposure impacts the control of viral replication remains poorly understood. We demonstrate here that immune events in the mouse lung closely preceding SARS-CoV-2 infection significantly impact viral control and we identify key innate immune pathways required to limit viral replication. A diverse set of pulmonary inflammatory stimuli, including resolved antecedent respiratory infections with S. aureus or influenza, ongoing pulmonary M. tuberculosis infection, ovalbumin/alum-induced asthma or airway administration of defined TLR ligands and recombinant cytokines, all establish an antiviral state in the lung that restricts SARS-CoV-2 replication upon infection. In addition to antiviral type I interferons, the broadly inducible inflammatory cytokines TNFα and IL-1 precondition the lung for enhanced viral control. Collectively, our work shows that SARS-CoV-2 may benefit from an immunologically quiescent lung microenvironment and suggests that heterogeneity in pulmonary inflammation that precedes or accompanies SARS-CoV-2 exposure may be a significant factor contributing to the population-wide variability in COVID-19 disease outcomes.

Legacy of a magic gene-CCR5-∆32: From discovery to clinical benefit in a generation

The discovery of the 32-bp deletion allele of the chemokine receptor gene CCR5 showed that homozygous carriers display near-complete resistance to HIV infection, irrespective of exposure. Algorithms of molecular evolutionary theory suggested that the CCR5-∆32 mutation occurred but once in the last millennium and rose by strong selective pressure relatively recently to a ~10% allele frequency in Europeans. Several lines of evidence support the hypothesis that CCR5-∆32 was selected due to its protective influence to resist Yersinia pestis, the agent of the Black Death/bubonic plague of the 14th century. Powerful anti-AIDS entry inhibitors targeting CCR5 were developed as a treatment for HIV patients, particularly those whose systems had developed resistance to powerful anti-retroviral therapies. Homozygous CCR5-∆32/∆32 stem cell transplant donors were used to produce HIV-cleared AIDS patients in at least five "cures" of HIV infection. CCR5 has also been implicated in regulating infection with Staphylococcus aureus, in recovery from stroke, and in ablation of the fatal graft versus host disease (GVHD) in cancer transplant patients. While homozygous CCR5-∆32/32 carriers block HIV infection, alternatively they display an increased risk for encephalomyelitis and death when infected with the West Nile virus.

[Genetic factors associated with long COVID]

INTRODUCTION

The variability in expression and evolution of COVID is not completely explained by clinical factors. In fact, genetic factors play an important role. Moreover, it is unknown whether the genetic factor that contribute to susceptibility and severity are also involved in the onset and evolution of long-COVID. The objective of this review is to gather information from literature to understand which genetic factors are involved in the onset of persistent COVID.

MATERIAL AND METHODS

Systematic review in PubMed and bioRxiv and medRxiv repositories based on MeSH-descriptors and MeSH-terms related to COVID and genetic factors. Using these terms 2715 articles were pooled. An initial screening performed by authors independently, selected 205 articles of interest. A final deeper screening a total of 85 articles were chosen for complete reading and summarized in this review.

RESULTS

Although ACE2 and TMPSS6 are involved in COVID susceptibility, their involvement in long-COVID has not been found. On the other hand, the severity of the disease and the onset of long-COVID has been associated with different genes involved in the inflammatory and immune response. Particularly interesting has been the association found with the FOXP4 locus.

CONCLUSIONS

Although studies on long-COVID are insufficient to fully comprehend the cause, it is clear that the current identified genetic factors do not fully explain the progression and onset of long-COVID. Other factors such as polygenic action, pleiotropic genes, the microbiota and epigenetic changes must be considered and studied.

Host genetic variants associated with COVID-19 reconsidered in a Slovak cohort

We present the results of an association study involving hospitalized coronavirus disease 2019 (COVID-19) patients with a clinical background during the 3rd pandemic wave of COVID-19 in Slovakia. Seventeen single nucleotide variants (SNVs) in the eleven most relevant genes, according to the COVID-19 Host Genetics Initiative, were investigated. Our study confirms the validity of the influence of LZTFL1 and 2'-5'-oligoadenylate synthetase (OAS)1/OAS3 genetic variants on the severity of COVID-19. For two LZTFL1 SNVs in complete linkage disequilibrium, rs17713054 and rs73064425, the odds ratios of baseline allelic associations and logistic regressions (LR) adjusted for age and sex ranged in the four tested designs from 2.04 to 2.41 and from 2.05 to 3.98, respectively. The OAS1/OAS3 haplotype 'gttg' carrying a functional allele G of splice-acceptor variant rs10774671 manifested its protective function in the Delta pandemic wave. Significant baseline allelic associations of two DPP9 variants in all tested designs and two IFNAR2 variants in the Omicron pandemic wave were not confirmed by adjusted LR. Nevertheless, adjusted LR showed significant associations of NOTCH4 rs3131294 and TYK2 rs2304256 variants with severity of COVID-19. Hospitalized patients' reported comorbidities were not correlated with genetic variants, except for obesity, smoking (IFNAR2), and hypertension (NOTCH4). The results of our study suggest that host genetic variations have an impact on the severity and duration of acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. Considering the differences in allelic associations between pandemic waves, they support the hypothesis that every new SARS-CoV-2 variant may modify the host immune response by reconfiguring involved pathways.

Biogeographic Perspectives on Human Genetic Diversification

Modern humans originated in Africa 300,000 yr ago, and before leaving their continent of origin, they underwent a process of intense diversification involving complex demographic dynamics. Upon exiting Africa, different populations emerged on the four other inhabited continents, shaped by the interplay of various evolutionary processes, such as migrations, founder effects, and natural selection. Within each region, continental populations, in turn, diversified and evolved almost independently for millennia. As a backdrop to this diversification, introgressions from archaic species contributed to establishing different patterns of genetic diversity in different geographic regions, reshaping our understanding of our species' variability. With the increasing availability of genomic data, it has become possible to delineate the subcontinental human population structure precisely. However, the bias toward the genomic research focused on populations from the global North has limited our understanding of the real diversity of our species and the processes and events that guided different human groups throughout their evolutionary history. This perspective is part of a series of articles celebrating 40 yr since our journal, Molecular Biology and Evolution, was founded (Russo et al. 2024). The perspective is accompanied by virtual issues, a selection of papers on human diversification published by Genome Biology and Evolution and Molecular Biology and Evolution.

The genetic changes that shaped Neandertals, Denisovans, and modern humans

Modern human ancestors diverged from the ancestors of Neandertals and Denisovans about 600,000 years ago. Until about 40,000 years ago, these three groups existed in parallel, occasionally met, and exchanged genes. A critical question is why modern humans, and not the other two groups, survived, became numerous, and developed complex cultures. Here, we discuss genetic differences among the groups and some of their functional consequences. As more present-day genome sequences become available from diverse groups, we predict that very few, if any, differences will distinguish all modern humans from all Neandertals and Denisovans. We propose that the genetic basis of what constitutes a modern human is best thought of as a combination of genetic features, where perhaps none of them is present in each and every present-day individual.

Altered expression of serum lncRNA CASC2 and miRNA-21-5p in COVID-19 patients

Infection by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) that causes coronavirus disease 2019 (COVID-19) has a high incidence of spread. On January 30, 2020, the World Health Organization proclaimed a public health emergency of worldwide concern. More than 6.9 million deaths and more than 768 million confirmed cases had been reported worldwide as of June 18, 2023. This study included 51 patients and 50 age- and sex-matched healthy subjects. The present study aimed to identify the expression levels of lncRNA CASC2 and miRNA-21-5p (also known as miRNA-21) in COVID-19 patients and their relation to the clinicopathological characteristics of the disease. The expression levels of noncoding RNAs were measured by RT-PCR technique. Results detected that CASC2 was significantly downregulated while miRNA-21-5p was significantly upregulated in COVID-19 patients compared to healthy subjects. A significant negative correlation was found between CASC2 and miRNA-21-5p. ROC curve analysis used to distinguish COVID-19 patients from controls. MiRNA-21-p serum expression level had a significant positive association with temperature and PO2 (p = 0.04 for each). These findings indicate that CASC2 and miRNA-21-p might be used as potential diagnostic and therapeutic biomarkers in COVID-19.

Increased prevalence of the COVID-19 associated Neanderthal mutations in the Central European Roma population

BACKGROUND

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection and subsequent COVID-19 has spread world-wide and become pandemic with about 7 million deaths reported so far. Interethnic variability of the disease has been described, but a significant part of the differences remain unexplained and may be attributable to genetic factors.

AIM

To analyse genetic factors potentially influencing COVID-19 susceptibility and severity in European Roma minority.

SUBJECTS AND METHODS

Two genetic determinants, within OAS-1 (2-prime,5-prime-oligoadenylate synthetase 1, a key protein in the defence against viral infection; it activates RNases that degrade viral RNAs; rs4767027 has been analysed) and LZTFL1 (leucine zipper transcription factor-like 1, expressed in the lung respiratory epithelium; rs35044562 has been analysed) genes were screened in a population-sample of Czech Roma (N = 302) and majority population (N = 2,559).

RESULTS

For both polymorphisms, Roma subjects were more likely carriers of at least one risky allele for both rs4767027-C (p < 0.001) and rs35044562-G (p < 0.00001) polymorphism. There were only 5.3% Roma subjects without at least one risky allele in comparison with 10.1% in the majority population (p < 0.01).

CONCLUSIONS

It is possible that different genetic background plays an important role in increased prevalence of COVID-19 in the Roma minority.

Worldwide distribution of genetic factors related to severity of COVID-19 infection

BACKGROUND

Genome-wide association studies of COVID-19 severity have been carried out mostly on European or East Asian populations with small representation of other world regions. Here we explore the worldwide distribution and linkage disequilibrium (LD) patterns of genetic variants previously associated with COVID-19 severity.

METHODS

We followed up the results of a large Spanish genome-wide meta-analysis on 26 populations from the 1000 Genomes Project by calculating allele frequencies and LD scores of the nine most significant SNPs. We also used the entire set of summary statistics to compute polygenic risk scores (PRSs) and carried out comparisons at the population and continental level.

RESULTS

We observed the strongest differences among continental regions for the five top SNPs in chromosome 3. European, American, and South Asian populations showed similar LD patterns. Average PRSs in South Asian and American populations were consistently higher than those observed in Europeans. While PRS distributions were similar among South Asians, the American populations showed striking differences among them.

CONCLUSIONS

Considering the caveats of PRS transferability across ethnicities, our analysis showed that American populations present the highest genetic risk score, hence potentially higher propensity, for COVID-19 severity. Independent validation is warranted with additional summary statistics and phenotype data.

Evolutionary origin of germline pathogenic variants in human DNA mismatch repair genes

BACKGROUND

Mismatch repair (MMR) system is evolutionarily conserved for genome stability maintenance. Germline pathogenic variants (PVs) in MMR genes that lead to MMR functional deficiency are associated with high cancer risk. Knowing the evolutionary origin of germline PVs in human MMR genes will facilitate understanding the biological base of MMR deficiency in cancer. However, systematic knowledge is lacking to address the issue. In this study, we performed a comprehensive analysis to know the evolutionary origin of human MMR PVs.

METHODS

We retrieved MMR gene variants from the ClinVar database. The genomes of 100 vertebrates were collected from the UCSC genome browser and ancient human sequencing data were obtained through comprehensive data mining. Cross-species conservation analysis was performed based on the phylogenetic relationship among 100 vertebrates. Rescaled ancient sequencing data were used to perform variant calling for archeological analysis.

RESULTS

Using the phylogenetic approach, we traced the 3369 MMR PVs identified in modern humans in 99 non-human vertebrate genomes but found no evidence for cross-species conservation as the source for human MMR PVs. Using the archeological approach, we searched the human MMR PVs in over 5000 ancient human genomes dated from 45,045 to 100 years before present and identified a group of MMR PVs shared between modern and ancient humans mostly within 10,000 years with similar quantitative patterns.

CONCLUSION

Our study reveals that MMR PVs in modern humans were arisen within the recent human evolutionary history.

The Inheritance of Hearing Loss and Deafness: A Historical Perspective

If the term "genetics" is a relatively recent proposition, introduced in 1905 by English biologist William Bateson, who rediscovered and spread in the scientific community Mendel's principles of inheritance, since the dawn of human civilization the influence of heredity has been recognized, especially in agricultural crops and animal breeding. And, later, in familial dynasties. In this concise review, we outline the evolution of the idea of hereditary hearing loss, up to the current knowledge of molecular genetics and epigenetics.

The Human Genome Organisation (HUGO) and a vision for Ecogenomics: the Ecological Genome Project

BACKGROUND

The following outlines ethical reasons for widening the Human Genome Organisation's (HUGO) mandate to include ecological genomics. MAIN: The environment influences an organism's genome through ambient factors in the biosphere (e.g. climate and UV radiation), as well as the agents it comes into contact with, i.e. the epigenetic and mutagenic effects of inanimate chemicals and pollution, and pathogenic organisms. Emerging scientific consensus is that social determinants of health, environmental conditions and genetic factors work together to influence the risk of many complex illnesses. That paradigm can also explain the environmental and ecological determinants of health as factors that underlie the (un)healthy ecosystems on which communities rely. We suggest that The Ecological Genome Project is an aspirational opportunity to explore connections between the human genome and nature. We propose consolidating a view of Ecogenomics to provide a blueprint to respond to the environmental challenges that societies face. This can only be achieved by interdisciplinary engagement between genomics and the broad field of ecology and related practice of conservation. In this respect, the One Health approach is a model for environmental orientated work. The idea of Ecogenomics-a term that has been used to relate to a scientific field of ecological genomics-becomes the conceptual study of genomes within the social and natural environment.

CONCLUSION

The HUGO Committee on Ethics, Law and Society (CELS) recommends that an interdisciplinary One Health approach should be adopted in genomic sciences to promote ethical environmentalism. This perspective has been reviewed and endorsed by the HUGO CELS and the HUGO Executive Board.

ACE2 knockout hinders SARS-CoV-2 propagation in iPS cell-derived airway and alveolar epithelial cells

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), the causative agent of COVID-19, continues to spread around the world with serious cases and deaths. It has also been suggested that different genetic variants in the human genome affect both the susceptibility to infection and severity of disease in COVID-19 patients. Angiotensin-converting enzyme 2 (ACE2) has been identified as a cell surface receptor for SARS-CoV and SARS-CoV-2 entry into cells. The construction of an experimental model system using human iPS cells would enable further studies of the association between viral characteristics and genetic variants. Airway and alveolar epithelial cells are cell types of the lung that express high levels of ACE2 and are suitable for in vitro infection experiments. Here, we show that human iPS cell-derived airway and alveolar epithelial cells are highly susceptible to viral infection of SARS-CoV-2. Using gene knockout with CRISPR-Cas9 in human iPS cells we demonstrate that ACE2 plays an essential role in the airway and alveolar epithelial cell entry of SARS-CoV-2 in vitro. Replication of SARS-CoV-2 was strongly suppressed in ACE2 knockout (KO) lung cells. Our model system based on human iPS cell-derived lung cells may be applied to understand the molecular biology regulating viral respiratory infection leading to potential therapeutic developments for COVID-19 and the prevention of future pandemics.

Possible effect of OAS1 and TMPRSS6 but not DPP4 and ZNF335 polymorphisms on COVID-19 severity in the Czech population

OBJECTIVES

The acute respiratory syndrome, known as COVID-19, is characterised by high morbidity and increased mortality. Genetic factors may partially explain the differences in susceptibility to and severity of COVID-19.

METHODS

We have analysed common functional polymorphisms within the OAS1 (rs4767027), TMPRSS6 (rs855791), DPP4 (rs3788979), and ZNF335 (rs3848719) genes in SARS-CoV-2 positive subjects (n = 521, different disease severity) and in population controls (n = 2,559 subjects, COVID-19 status unknown).

RESULTS

Neither DPP4 nor ZNF335 were associated with disease susceptibility or severity in the Czech population in any of the models used for calculation. T allele carriers of the OAS1 polymorphism seem to be protective against symptomatic COVID-19 (p = 0.002 calculated for trend; asymptomatic, symptomatic, hospitalised). Similarly, within the TMPRSS6, minor TT homozygotes associated with lower plasma Fe concentrations were underrepresented in the overall patient group (p = 0.044; OR = 0.77, 95% CI: 0.59-0.99), and the difference was mainly driven by the severe COVID-19 subjects. In general, risky homozygotes of these two polymorphisms were less frequent than expected in the group of hospitalised COVID-19 survivors.

CONCLUSIONS

Common variants within OAS1 (rs4767027) and TMPRSS6 (rs855791) play some role in COVID-19 pathology in the Czech Caucasian population. Whether the depletion of minor allele carriers of these two variants is associated with increased COVID-19 mortality, needs to be analysed in an external confirmatory study.

Envisioning a new era: Complete genetic information from routine, telomere-to-telomere genomes

Advances in long-read sequencing and assembly now mean that individual labs can generate phased genomes that are more accurate and more contiguous than the original human reference genome. With declining costs and increasing democratization of technology, we suggest that complete genome assemblies, where both parental haplotypes are phased telomere to telomere, will become standard in human genetics. Soon, even in clinical settings where rigorous sample-handling standards must be met, affected individuals could have reference-grade genomes fully sequenced and assembled in just a few hours given advances in technology, computational processing, and annotation. Complete genetic variant discovery will transform how we map, catalog, and associate variation with human disease and fundamentally change our understanding of the genetic diversity of all humans.

Association between ABO blood group/genotype and COVID-19 in a Japanese population

An association between coronavirus disease 2019 (COVID-19) and the ABO blood group has been reported. However, such an association has not been studied in the Japanese population on a large scale. Little is known about the association between COVID-19 and ABO genotype. This study investigated the association between COVID-19 and ABO blood group/genotype in a large Japanese population. All Japanese patients diagnosed with COVID-19 were recruited through the Japan COVID-19 Task Force between February 2020 and October 2021. We conducted a retrospective cohort study involving 1790 Japanese COVID-19 patients whose DNA was used for a genome-wide association study. We compared the ABO blood group/genotype in a healthy population (n = 611, control) and COVID-19 patients and then analyzed their associations and clinical outcomes. Blood group A was significantly more prevalent (41.6% vs. 36.8%; P = 0.038), and group O was significantly less prevalent (26.2% vs. 30.8%; P = 0.028) in the COVID-19 group than in the control group. Moreover, genotype OO was significantly less common in the COVID-19 group. Furthermore, blood group AB was identified as an independent risk factor for most severe diseases compared with blood group O [aOR (95% CI) = 1.84 (1.00-3.37)]. In ABO genotype analysis, only genotype AB was an independent risk factor for most severe diseases compared with genotype OO. Blood group O is protective, whereas group A is associated with the risk of infection. Moreover, blood group AB is associated with the risk of the "most" severe disease.

Intronic Variants of the Angiotensin-Converting Enzyme 2 Gene Modulate Plasma ACE2 Levels and Possibly Confer Protection against Severe COVID-19

Membrane-bound angiotensin-converting enzyme 2 (ACE2) receptor acts as the entry point for the novel coronavirus, SARS-CoV-2. Polymorphisms in the ACE2 gene may alter viral binding, regulate the expression of ACE2, and thus, affect disease severity. In this study, 68 COVID-19 patients with varying degrees of severity and 40 healthy controls were enrolled. The genetic landscape of the ACE2 gene was explored by whole exome sequencing of 29 individuals, and specific regions of ACE2 were analyzed for the rest of the participants via PCR, followed by barcode-tagged sequencing. The mean soluble ACE2 level in the plasma of healthy controls and patients did not vary significantly but was higher in the patient group (3.77 ± 1.55 ng/mL vs. 3.94 ± 1.42 ng/mL). Analysis of exon 1, exon 2, and exon 8 of the ACE2 gene revealed that these regions are highly conserved in our population. Investigation of exon 11 and its flanking intronic region revealed that deletions in a stretch of 18T nucleotides in the noncoding region significantly decrease ACE2 levels in plasma, as individuals harboring wild-type variants had higher plasma ACE2 levels compared to those harboring T1del, T2del, and T3del variants. However, the intronic variants were not found to be significantly associated with disease severity.

Past human expansions shaped the spatial pattern of Neanderthal ancestry

The worldwide expansion of modern humans (Homo sapiens) started before the extinction of Neanderthals (Homo neanderthalensis). Both species coexisted and interbred, leading to slightly higher introgression in East Asians than in Europeans. This distinct ancestry level has been argued to result from selection, but range expansions of modern humans could provide an alternative explanation. This hypothesis would lead to spatial introgression gradients, increasing with distance from the expansion source. We investigate the presence of Neanderthal introgression gradients after past human expansions by analyzing Eurasian paleogenomes. We show that the out-of-Africa expansion resulted in spatial gradients of Neanderthal ancestry that persisted through time. While keeping the same gradient orientation, the expansion of early Neolithic farmers contributed decisively to reducing the Neanderthal introgression in European populations compared to Asian populations. This is because Neolithic farmers carried less Neanderthal DNA than preceding Paleolithic hunter-gatherers. This study shows that inferences about past human population dynamics can be made from the spatiotemporal variation in archaic introgression.

Genetic signature detected in T cell receptors from patients with severe COVID-19

Characterization of host genetic factors contributing to COVID-19 severity promises advances on drug discovery to fight the disease. Most genetic analyses to date have identified genome-wide significant associations involving loss-of-function variants for immune response pathways. Despite accumulating evidence supporting a role for T cells in COVID-19 severity, no definitive genetic markers have been found to support an involvement of T cell responses. We analyzed 205 whole exomes from both a well-characterized cohort of hospitalized severe COVID-19 patients and controls. Significantly enriched high impact alleles were found for 25 variants within the T cell receptor beta (TRB) locus on chromosome 7. Although most of these alleles were found in heterozygosis, at least three or more in TRBV6-5, TRBV7-3, TRBV7-6, TRBV7-7, and TRBV10-1 suggested a possible TRB loss of function via compound heterozygosis. This loss-of-function in TRB genes supports suboptimal or dysfunctional T cell responses as a major contributor to severe COVID-19 pathogenesis.

A GWAS in the pandemic epicenter highlights the severe COVID-19 risk locus introgressed by Neanderthals

Large GWAS indicated that genetic factors influence the response to SARS-CoV-2. However, sex, age, concomitant diseases, differences in ancestry, and uneven exposure to the virus impacted the interpretation of data. We aimed to perform a GWAS of COVID-19 outcome in a homogeneous population who experienced a high exposure to the virus and with a known infection status. We recruited inhabitants of Bergamo province-that in spring 2020 was the epicenter of the SARS-Cov-2 pandemic in Europe-via an online questionnaire followed by personal interviews. Cases and controls were matched by age, sex and risk factors. We genotyped 1195 individuals and replicated the association at the 3p21.31 locus with severity, but with a stronger effect size that further increased in gravely ill patients. Transcriptome-wide association study highlighted eQTLs for LZTFL1 and CCR9. We also identified 17 loci not previously reported, suggestive for an association with either COVID-19 severity or susceptibility.

COVID-19-Related Age Profiles for SARS-CoV-2 Variants in England and Wales and States of the USA (2020 to 2022): Impact on All-Cause Mortality

Since 2020, COVID-19 has caused serious mortality around the world. Given the ambiguity in establishing COVID-19 as the direct cause of death, we first investigate the effects of age and sex on all-cause mortality during 2020 and 2021 in England and Wales. Since infectious agents have their own unique age profile for death, we use a 9-year time series and several different methods to adjust single-year-of-age deaths in England and Wales during 2019 (the pre-COVID-19 base year) to a pathogen-neutral single-year-of-age baseline. This adjusted base year is then used to confirm the widely reported higher deaths in males for most ages above 43 in both 2020 and 2021. During 2020 (+COVID-19 but no vaccination), both male and female population-adjusted deaths significantly increased above age 35. A significant reduction in all-cause mortality among both males and females aged 75+ could be demonstrated in 2021 during the widespread COVID-19 vaccination period; however, deaths below age 75 progressively increased. This finding arises from a mix of vaccination coverage and year-of-age profiles of deaths for the different SARS-CoV-2 variants. In addition, specific effects of age around puberty were demonstrated, where females had higher deaths than males. There is evidence that year-of-birth cohorts may also be involved, indicating that immune priming to specific pathogen outbreaks in the past may have led to lower deaths for some birth cohorts. To specifically identify the age profile for the COVID-19 variants from 2020 to 2023, we employ the proportion of total deaths at each age that are potentially due to or 'with' COVID-19. The original Wuhan strain and the Alpha variant show somewhat limited divergence in the age profile, with the Alpha variant shifting to a moderately higher proportion of deaths below age 84. The Delta variant specifically targeted individuals below age 65. The Omicron variants showed a significantly lower proportion of overall mortality, with a markedly higher relative proportion of deaths above age 65, steeply increasing with age to a maximum around 100 years of age. A similar age profile for the variants can be seen in the age-banded deaths in US states, although they are slightly obscured by using age bands rather than single years of age. However, the US data shows that higher male deaths are greatly dependent on age and the COVID variant. Deaths assessed to be 'due to' COVID-19 (as opposed to 'involving' COVID-19) in England and Wales were especially overestimated in 2021 relative to the change in all-cause mortality. This arose as a by-product of an increase in COVID-19 testing capacity in late 2020. Potential structure-function mechanisms for the age-specificity of SARS-CoV-2 variants are discussed, along with potential roles for small noncoding RNAs (miRNAs). Using data from England, it is possible to show that the unvaccinated do indeed have a unique age profile for death from each variant and that vaccination alters the shape of the age profile in a manner dependent on age, sex, and the variant. The question is posed as to whether vaccines based on different variants carry a specific age profile.

COVID-19 host genetic risk study conducted at community pharmacies: Implications for public health, research and pharmacists' scope of practice

Community pharmacists serve a large, diverse population of patients, resulting in the potential to utilize community pharmacies as recruitment sites for clinical research. Beyond traditional roles as one of the most accessible health care professionals in the US healthcare system, pharmacists have played a major role in the response to the COVID-19 pandemic, administering hundreds of thousands of vaccines and tests. However, less emphasis is placed on the ability to leverage community pharmacies as research-focused partners for clinical studies. In this study, we demonstrate the feasibility and workflow of recruiting study participants from community pharmacies and confirm genetic markers of COVID-19 susceptibility. Specific genetic markers include those associated with COVID-19 infection risk (ACE2, TMEM27, and RAVER1), difficulty breathing (NOTCH4), and hospitalization (OAS3). In addition, collaboration with a clinical laboratory allowed for a more seamless consenting process without substantial training needs or workflow disruption at the community pharmacy site. The COVID-19 pandemic has demonstrated that the expansion of pharmacists' scope of practice is a key factor in managing the population health crisis; this study demonstrates that pharmacies can also advance clinical research studies by serving as sites for patient recruitment from a large, diverse, and ambulatory study population.

Human-specific genetics: new tools to explore the molecular and cellular basis of human evolution

Our ancestors acquired morphological, cognitive and metabolic modifications that enabled humans to colonize diverse habitats, develop extraordinary technologies and reshape the biosphere. Understanding the genetic, developmental and molecular bases for these changes will provide insights into how we became human. Connecting human-specific genetic changes to species differences has been challenging owing to an abundance of low-effect size genetic changes, limited descriptions of phenotypic differences across development at the level of cell types and lack of experimental models. Emerging approaches for single-cell sequencing, genetic manipulation and stem cell culture now support descriptive and functional studies in defined cell types with a human or ape genetic background. In this Review, we describe how the sequencing of genomes from modern and archaic hominins, great apes and other primates is revealing human-specific genetic changes and how new molecular and cellular approaches - including cell atlases and organoids - are enabling exploration of the candidate causal factors that underlie human-specific traits.

Gomariz R, Muñoz-Calleja C, Fernández-Ruiz E, González-Álvaro I, Cardeñoso L. Occurrence of SARS-CoV-2 viremia is associated with genetic variants of genes related to COVID-19 pathogenesis

Introduction

SARS-CoV-2 viral load has been related to COVID-19 severity. The main aim of this study was to evaluate the relationship between SARS-CoV-2 viremia and SNPs in genes previously studied by our group as predictors of COVID-19 severity.

Materials and methods

Retrospective observational study including 340 patients hospitalized for COVID-19 in the University Hospital La Princesa between March 2020 and December 2021, with at least one viremia determination. Positive viremia was considered when viral load was above the quantifiable threshold (20 copies/ml). A total of 38 SNPs were genotyped. To study their association with viremia a multivariate logistic regression was performed.

Results

The mean age of the studied population was 64.5 years (SD 16.6), 60.9% patients were male and 79.4% white non-Hispanic. Only 126 patients (37.1%) had at least one positive viremia. After adjustment by confounders, the presence of the minor alleles of rs2071746 (HMOX1; T/T genotype OR 9.9 p < 0.0001), rs78958998 (probably associated with SERPING1 expression; A/T genotype OR 2.3, p = 0.04 and T/T genotype OR 12.9, p < 0.0001), and rs713400 (eQTL for TMPRSS2; C/T + T/T genotype OR 1.86, p = 0.10) were associated with higher risk of viremia, whereas the minor alleles of rs11052877 (CD69; A/G genotype OR 0.5, p = 0.04 and G/G genotype OR 0.3, p = 0.01), rs2660 (OAS1; A/G genotype OR 0.6, p = 0.08), rs896 (VIPR1; T/T genotype OR 0.4, p = 0.02) and rs33980500 (TRAF3IP2; C/T + T/T genotype OR 0.3, p = 0.01) were associated with lower risk of viremia.

Conclusion

Genetic variants in HMOX1 (rs2071746), SERPING1 (rs78958998), TMPRSS2 (rs713400), CD69 (rs11052877), TRAF3IP2 (rs33980500), OAS1 (rs2660) and VIPR1 (rs896) could explain heterogeneity in SARS-CoV-2 viremia in our population.

Dissecting human population variation in single-cell responses to SARS-CoV-2

Humans display substantial interindividual clinical variability after SARS-CoV-2 infection1-3, the genetic and immunological basis of which has begun to be deciphered4. However, the extent and drivers of population differences in immune responses to SARS-CoV-2 remain unclear. Here we report single-cell RNA-sequencing data for peripheral blood mononuclear cells-from 222 healthy donors of diverse ancestries-that were stimulated with SARS-CoV-2 or influenza A virus. We show that SARS-CoV-2 induces weaker, but more heterogeneous, interferon-stimulated gene activity compared with influenza A virus, and a unique pro-inflammatory signature in myeloid cells. Transcriptional responses to viruses display marked population differences, primarily driven by changes in cell abundance including increased lymphoid differentiation associated with latent cytomegalovirus infection. Expression quantitative trait loci and mediation analyses reveal a broad effect of cell composition on population disparities in immune responses, with genetic variants exerting a strong effect on specific loci. Furthermore, we show that natural selection has increased population differences in immune responses, particularly for variants associated with SARS-CoV-2 response in East Asians, and document the cellular and molecular mechanisms by which Neanderthal introgression has altered immune functions, such as the response of myeloid cells to viruses. Finally, colocalization and transcriptome-wide association analyses reveal an overlap between the genetic basis of immune responses to SARS-CoV-2 and COVID-19 severity, providing insights into the factors contributing to current disparities in COVID-19 risk.

Smoke and Spike: Benzo[a]pyrene Enhances SARS-CoV-2 Infection by Boosting NR4A2-Induced ACE2 and TMPRSS2 Expression

Cigarette smoke aggravates severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. However, the underlying mechanisms remain unclear. Here, they show that benzo[a]pyrene in cigarette smoke extract facilitates SARS-CoV-2 infection via upregulating angiotensin-converting enzyme 2 (ACE2) and transmembrane protease serine 2 (TMPRSS2). Benzo[a]pyrene trans-activates the promoters of ACE2 and TMPRSS2 by upregulating nuclear receptor subfamily 4 A number 2 (NR4A2) and promoting its binding of NR4A2 to their promoters, which is independent of functional genetic polymorphisms in ACE2 and TMPRSS2. Benzo[a]pyrene increases the susceptibility of lung epithelial cells to SARS-CoV-2 pseudoviruses and facilitates the infection of authentic Omicron BA.5 in primary human alveolar type II cells, lung organoids, and lung and testis of hamsters. Increased expression of Nr4a2, Ace2, and Tmprss2, as well as decreased methylation of CpG islands at the Nr4a2 promoter are observed in aged mice compared to their younger counterparts. NR4A2 knockdown or interferon-λ2/λ3 stimulation downregulates the expression of NR4A2, ACE2, and TMPRSS2, thereby inhibiting the infection. In conclusion, benzo[a]pyrene enhances SARS-CoV-2 infection by boosting NR4A2-induced ACE2 and TMPRSS2 expression. This study elucidates the mechanisms underlying the detrimental effects of cigarette smoking on SARS-CoV-2 infection and provides prophylactic options for coronavirus disease 2019, particularly for the elderly population.

Overview of immunological & virological factors driving the evolution & global spread of SARS-CoV-2 variants

The SARS-CoV-2, a highly infectious positive strand RNA virus first identified in December 2019, has produced multiple genetic variants that have rapidly and sequentially spread worldwide during the coronavirus disease 2019 (COVID-19) pandemic. Genetic changes in SARS-CoV-2 for greater infectivity, replication and transmission were selected during the early stages of the pandemic. More recently, after widespread infection and vaccination, SARS-CoV-2 variants that evade antigen-specific adaptive immunity, have begun to be selected. This article provides an overview of the molecular immunological and virological factors underlying the origin and global spread of important SARS-CoV-2 variant lineages.

Association of 3p21.31 Locus (CXCR6 and LZTFL1) with COVID-19 Outcomes in Brazilian Hospitalyzed Subjects

The 3p21.31 locus has been associated with severe COVID-19 prognosis in GWAS studies. Here, we evaluated whether three polymorphisms (LZTFL1 rs10490770, CXCR6 rs2234355 and rs2234358) in the reported locus were associated with the need for mechanical ventilation, hospitalization length and death in 102 COVID-19 hospitalized Brazilian subjects. No genetic association was found with the need for mechanical ventilation and hospitalization length. CXCR6 rs2234355 was associated with mortality under the codominance model, with carriers of the A/A genotype having a greater chance of death than A/G (OR: 10.5; 95% CI: 1.55-70.76). Our results further suggest that the CXCR6 genetic variant contributes to COVID-19 outcomes.

Human Genomics of COVID-19 Pneumonia: Contributions of Rare and Common Variants

SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) infection is silent or benign in most infected individuals, but causes hypoxemic COVID-19 pneumonia in about 10% of cases. We review studies of the human genetics of life-threatening COVID-19 pneumonia, focusing on both rare and common variants. Large-scale genome-wide association studies have identified more than 20 common loci robustly associated with COVID-19 pneumonia with modest effect sizes, some implicating genes expressed in the lungs or leukocytes. The most robust association, on chromosome 3, concerns a haplotype inherited from Neanderthals. Sequencing studies focusing on rare variants with a strong effect have been particularly successful, identifying inborn errors of type I interferon (IFN) immunity in 1-5% of unvaccinated patients with critical pneumonia, and their autoimmune phenocopy, autoantibodies against type I IFN, in another 15-20% of cases. Our growing understanding of the impact of human genetic variation on immunity to SARS-CoV-2 is enabling health systems to improve protection for individuals and populations.

Seeking old wisdoms for new AMP discovery

In this issue of Cell Host & Microbe, Maasch et al. built a deep-learning model for proteolytic site prediction and applied this approach to modern and extinct human proteomes, revealing antimicrobial peptides. The de-extinction of archaic encrypted peptides provides a revolutionary approach for discovering functional peptides with therapeutic potential.

Molecular de-extinction of ancient antimicrobial peptides enabled by machine learning

Molecular de-extinction could offer avenues for drug discovery by reintroducing bioactive molecules that are no longer encoded by extant organisms. To prospect for antimicrobial peptides encrypted within extinct and extant human proteins, we introduce the panCleave random forest model for proteome-wide cleavage site prediction. Our model outperformed multiple protease-specific cleavage site classifiers for three modern human caspases, despite its pan-protease design. Antimicrobial activity was observed in vitro for modern and archaic protein fragments identified with panCleave. Lead peptides showed resistance to proteolysis and exhibited variable membrane permeabilization. Additionally, representative modern and archaic protein fragments showed anti-infective efficacy against A. baumannii in both a skin abscess infection model and a preclinical murine thigh infection model. These results suggest that machine-learning-based encrypted peptide prospection can identify stable, nontoxic peptide antibiotics. Moreover, we establish molecular de-extinction through paleoproteome mining as a framework for antibacterial drug discovery.

Human genetic and immunological determinants of SARS-CoV-2 and Epstein-Barr virus diseases in childhood: Insightful contrasts

There is growing evidence to suggest that severe disease in children infected with common viruses that are typically benign in other children can result from inborn errors of immunity or their phenocopies. Infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a cytolytic respiratory RNA virus, can lead to acute hypoxemic COVID-19 pneumonia in children with inborn errors of type I interferon (IFN) immunity or autoantibodies against IFNs. These patients do not appear to be prone to severe disease during infection with Epstein-Barr virus (EBV), a leukocyte-tropic DNA virus that can establish latency. By contrast, various forms of severe EBV disease, ranging from acute hemophagocytosis to chronic or long-term illnesses, such as agammaglobulinemia and lymphoma, can manifest in children with inborn errors disrupting specific molecular bridges involved in the control of EBV-infected B cells by cytotoxic T cells. The patients with these disorders do not seem to be prone to severe COVID-19 pneumonia. These experiments of nature reveal surprising levels of redundancy of two different arms of immunity, with type I IFN being essential for host defense against SARS-CoV-2 in respiratory epithelial cells, and certain surface molecules on cytotoxic T cells essential for host defense against EBV in B lymphocytes.

Learning from cancer to address COVID-19

Patients with cancer have been disproportionately affected by the novel coronavirus disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Knowledge collected during the last three decades of cancer research has helped the medical research community worldwide to respond to many of the challenges raised by COVID-19, during the pandemic. The review, briefly summarizes the underlying biology and risk factors of COVID-19 and cancer, and aims to present recent evidence on cellular and molecular relationship between the two diseases, with a focus on those that are related to the hallmarks of cancer and uncovered in the first less than three years of the pandemic (2020-2022). This may not only help answer the question "Why cancer patients are considered to be at a particularly high risk of developing severe COVID-19 illness?", but also helped treatments of patients during the COVID-19 pandemic. The last session highlights the pioneering mRNA studies and the breakthrough discovery on nucleoside-modifications of mRNA by Katalin Karikó, which led to the innovation and development of the mRNA-based SARSCoV-2 vaccines saving lives of millions and also opened the door for a new era of vaccines and a new class of therapeutics.

Comparative physiological anthropogeny: exploring molecular underpinnings of distinctly human phenotypes

Anthropogeny is a classic term encompassing transdisciplinary investigations of the origins of the human species. Comparative anthropogeny is a systematic comparison of humans and other living nonhuman hominids (so-called "great apes"), aiming to identify distinctly human features in health and disease, with the overall goal of explaining human origins. We begin with a historical perspective, briefly describing how the field progressed from the earliest evolutionary insights to the current emphasis on in-depth molecular and genomic investigations of "human-specific" biology and an increased appreciation for cultural impacts on human biology. While many such genetic differences between humans and other hominids have been revealed over the last two decades, this information remains insufficient to explain the most distinctive phenotypic traits distinguishing humans from other living hominids. Here we undertake a complementary approach of "comparative physiological anthropogeny," along the lines of the preclinical medical curriculum, i.e., beginning with anatomy and considering each physiological system and in each case considering genetic and molecular components that are relevant. What is ultimately needed is a systematic comparative approach at all levels from molecular to physiological to sociocultural, building networks of related information, drawing inferences, and generating testable hypotheses. The concluding section will touch on distinctive considerations in the study of human evolution, including the importance of gene-culture interactions.