Genes for covid-19 resilience: identification of DNA markers corresponding to coronavirus resistance and susceptibility
 
Coronaviruses (CoVs) (order Nidovirales, family Coronaviridae,
                                    subfamily
                                    Coronavirinae) are responsible for respiratory disease outbreaks in many
                                    vertebrate species.
                                    They are a large family of single-stranded RNA enveloped viruses (+ssRNA) that can
                                    be isolated in
                                    different animal species. They have genome sizes ranging between 26 to 32 kilobases
                                    (kb) in length,
                                    being the largest genomes for RNA viruses (consequently increasing the effectiveness
                                    of facemasks).
                                    COVID-19 also known as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2),
                                    or "novel
                                    coronavirus 2019" is a new virus
                                    and we are just beginning to understand resistance and susceptibility in humans. 
COVID-19 is similar to severe acute respiratory syndrome (SARS) in the respect
                                    that both viruses
                                    infect their human hosts via the same receptor, the angiotensin-converting enzyme 2
                                    (ACE2 receptor),
                                    and cause similar clinical and pathological features.
                                    Interestingly the spike protein which is responsible for receptor binding is highly
                                    similar between
                                    2019-nCoV and SARS-CoV, this is a result of significant selection for the same
                                    receptor (Wu.,
                                        2020). Research into how
                                    our bodies defend against SARS might reveal how our bodies could defend against
                                    covid-19.
                                
Several recent Genome Wide Association Studies (GWAS) have provided a much deeper
                                    insight into
                                    the genetic variations that can help to explain why some individuals are virtually
                                    unaffected by
                                    COVID-19, and for others the virus is life threatening or even fatal. 
In this post, we provide a review of the peer reviewed literature and present information about candidate genes for SARS-CoV resistance. If you have taken an at home DNA test like those available from 23andMe, Ancestry DNA, Dante labs, you can evaluate your raw DNA data and see how your DNA sequence compares with the research findings.
How to analyze your DNA for coronavirus resistance or susceptibility?
                                
Step 1) Download your raw autosomal DNA file and save it to a safe and secure location
To analyze your DNA data, start by downloading your raw autosomal DNA and save it to a safe location. Here are instructions to download your raw DNA file from: 23andMe, Ancestry DNA, Family Tree DNA, Dante Labs, My Heritage, Genes For Good, Vitagene, and Living DNA.
Step 2) Analyze your raw DNA file
Search your raw DNA data using a text editor such as "text wrangler" or "notepad" using the "find" function, or using command line.
Open your raw DNA file and you will notice the headers of unique SNP ID (rs# or i#), chromosome, position and genotype. The formats differ slightly between each direct to consumer DNA testing company.
To evaluate your risk of poor recovery from the COVID-19 have a look for these DNA markers described below:
Several Genome Wide Association Study (GWAS) were recently published that describe loci associated with respiratory failure in patients infected with SARS-CoV-2 and three studies identified SNP markers in the same ~50 kb genomic segment that is inherited from Neandertals (Ellinghaus D et al., 2020, Zeberg & Pääbo, and Blokland et al., 2020). In addition these GWAS studies also identified a number of other DNA markers that are associated with COVID-19, and each of these are presented in the tabel below.
In addition, other DNA markers covered in this post include rs4804803 which was associated with SARS, and those positioned in the angiotensin-converting enzyme-2 (ACE2) receptor which was proved to be the same receptor for the human respiratory coronavirus NL63, , SARS-coronavirus (SARS-CoV), and the novel coronavirus 2019-nCoV/SARS-CoV (Li et al., 2017; Lu et al., 2019). Since the spike protein of the coronavirus has evolved to match the ACE2 receptor, it's likely that individuals with variations that alter the protein sequence would result in a degree of resistance to covid-19. Below are non-synonymous SNPs from the ACE2 transcript NM_021804.2 and of particular interest are SNPs that cause major changes like rs199951323 which results in a premature stop codon.
| Gene | dbsnp | Chromosome (GRCh37) | POS | REF | ALT | Risk Alleles | Marker Effect | Reference | 
|---|---|---|---|---|---|---|---|---|
| IVNS1ABP | rs6668622 | 1 | 185414582 | T | C | Susceptibile Variant T:T and T:C in males | odds ratio 1.44 | Roberts., 2020; | 
| SRRM1 | rs111972040 | 1 | 24999361 | A | G | risk genotypes G:G and A:G, 3_prime_UTR_variant | odds ratio for hospitalization = 8.29 | |
| LZTFL1 | rs35044562 | 3 | 45909024 | A | G | risk genotypes G:G and A:G, intron_variant,genic_upstream_transcript_variant | odds ratio 1.60 | Blokland et al., 2020; Zeberg & Pääbo | 
| LZTFL1 | rs11385942 | 3 | 45876460 | A | - or A or AAA | InDel, A:A and A:- have higher respiratory failure susceptibility, intron_variant | odds ratio 1.77 | Ellinghaus et al., 2020; Roberts., 2020 | 
| LZTFL1 | rs10490770 *LD with rs11385942 | 3 | 45864732 | T | C | risk genotypes T:C and C:C, intron_variant | odds ratio for heterozygous carriers 1.7 | Zeberg and Pääbo., 2020; | 
| LZTFL1 | rs67959919 *LD with rs11385942 | 3 | 45871908 | G | A | risk genotypes G:A and A:A, intron_variant | odds ratio for heterozygous carriers 1.7 | |
| LZTFL1 | rs35624553 *LD with rs11385942 | 3 | 45867440 | A | G | risk genotypes G:A and G:G, intron_variant | odds ratio for heterozygous carriers 1.7 | |
| LZTFL1 | rs71325088 *LD with rs11385942 | 3 | 45862952 | T | C | risk genotypes C:T and C:C, intron_variant | odds ratio for heterozygous carriers 1.7 | |
| ABO | rs657152 | 9 | 136139265 | A | C or T | A risk allele, intron_variant | odds ratio 1.77 | Ellinghaus et al., 2020; Roberts., 2020 | 
| Intergenic | rs5798227 | 12 | 53120100 | C | - | Risk allele is the deletion | p = 2.2x10-7 | Blokland et al., 2020; | 
| IGHV3-7 | rs11844522 | 14 | 106522576 | C | T | Susceptibile variations T:T, C:T | p=1.9x10-7 | |
| Immunoglobulin Lambda Locus (IGL) | rs73166864 | 22 | 23340580 | T | C or G | Susceptibile variations T:T and T:C | odds ratio 1.7 | Roberts., 2020; | 
| TLR7 | rs200553089 | ChrX | 12906010 | G | T | risk genotypes T:G and T:T, missense_variant | Made et al., 2020; | |
| Synonymous SNPs positioned in ACE2 | ||||||||
| ACE2 | rs373153165 | chrX | 15580093 | C | T or A | missense_variant | p.Asp785Asn/c.2353G>A | Cao et al., 2020 | 
| ACE2 | rs140016715 | chrX | 15582154 | G | A | missense_variant | p.Arg768Trp/c.2302C>T | |
| ACE2 | rs147311723 | chrX | 15582265 | G | A | missense_variant | p.Leu731Phe/c.2191C>T | |
| ACE2 | rs41303171 | chrX | 15582298 | T | C | missense_variant | p.Asn720Asp/c.2158A>G | |
| ACE2 | rs370187012 | chrX | 15582327 | C | T | missense_variant | p.Arg710His/c.2129G>A | |
| ACE2 | rs776995986 | chrX | 15582334 | G | A | missense_variant | p.Arg708Trp/c.2122C>T | |
| ACE2 | rs149039346 | chrX | 15584416 | A | G | missense_variant | p.Ser692Pro/c.2074T>C | |
| ACE2 | rs200180615 | chrX | 15584488 | C | T | missense_variant | p.Glu668Lys/c.2002G>A | |
| ACE2 | * | chrX | 15585879 | A | C | stop_gained | p.Leu656*/c.1967T>G | |
| ACE2 | rs183135788 | chrX | 15585933 | T | C | missense_variant | p.Asn638Ser/c.1913A>G | |
| ACE2 | rs748163894 | chrX | 15588434 | G | A | missense_variant | ||
| ACE2 | rs202137736 | chrX | 15591485 | T | C | splice_region_variant+intron_variant | c.1541+5A>G | |
| ACE2 | rs140473595 | chrX | 15591530 | C | T | missense_variant | p.Ala501Thr/c.1501G>A | |
| ACE2 | rs191860450 | chrX | 15593829 | T | C | missense_variant | p.Ile468Val/c.1402A>G | |
| ACE2 | rs758142853 | chrX | 15609868 | A | G | missense_variant | p.Val184Ala/c.551T>C | |
| ACE2 | rs754511501 | chrX | 15609902 | C | T | missense_variant | p.Gly173Ser/c.517G>A | |
| ACE2 | rs746034076 | chrX | 15609943 | T | C | missense_variant | p.Asn159Ser/c.476A>G | |
| ACE2 | rs373252182 | chrX | 15609973 | T | C | missense_variant | p.Asn149Ser/c.446A>G | |
| ACE2 | rs2285666 | chrX | 15610348 | C | T | splice_region_variant+intron_variant | c.439+4G>A | |
| ACE2 | rs768736934 | chrX | 15612963 | C | T | splice_region_variant+intron_variant | c.345+5G>A | |
| ACE2 | rs4646116 | chrX | 15618958 | T | C | missense_variant | p.Lys26Arg/c.77A>G | |
| ACE2 | rs73635825 | chrX | 15618980 | A | G | missense_variant | p.Ser19Pro/c.55T>C | |
| SNPs associated with SARS | ||||||||
| CD209 | rs4804803 | 19 | 7812733 | A | G | Susceptibile genotype A:A, upstream_transcript_variant | NC_000019.10:7747846 | Sakuntabhai et al., 2005; Chan et al., 2010 | 
