key: cord-0892939-t7o5g5mx
authors: Koch, Benjamin Florian
title: SARS-CoV-2 and human retroelements: a case for molecular mimicry?
date: 2022-04-08
journal: BMC Genom Data
DOI: 10.1186/s12863-022-01040-2
sha: 65e0cf8212aafa49c269efee682cd53902a2d3d1
doc_id: 892939
cord_uid: t7o5g5mx

BACKGROUND: The factors driving the late phase of COVID-19 are still poorly understood. However, autoimmunity is an evolving theme in COVID-19’s pathogenesis. Additionally, deregulation of human retroelements (RE) is found in many viral infections, and has also been reported in COVID-19. RESULTS: Unexpectedly, coronaviruses (CoV) – including SARS-CoV-2 – harbour many RE-identical sequences (up to 35 base pairs), and some of these sequences are part of SARS-CoV-2 epitopes associated to COVID-19 severity. Furthermore, RE are expressed in healthy controls and human cells and become deregulated after SARS-CoV-2 infection, showing mainly changes in long interspersed nuclear element (LINE1) expression, but also in endogenous retroviruses. CONCLUSION: CoV and human RE share coding sequences, which are targeted by antibodies in COVID-19 and thus could induce an autoimmune loop by molecular mimicry. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1186/s12863-022-01040-2.

At the end of 2019, a severe acute respiratory syndrome (SARS)-like disease was noted in eastern China and a novel coronavirus (later designated SARS-CoV-2) recognized as the factor for the disease, COVID-19 [1] . By the spring of 2022, 447 million people have been infected globally, with 6 million casualties [2] . COVID-19 can be divided into an early viral replication phase and a late stage of organ failure [3, 4] . While the inhibition of SARS-CoV-2 replication has already been achieved [5] [6] [7] [8] [9] [10] , the factors driving the late phase of the disease are poorly understood [11, 12] . However, it has been reported that autoimmunity [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] and deregulation of human retroelements (RE) might contribute to the outcome of COVID-19 patients [28] [29] [30] [31] .

The RE share a reverse transcriptase as a common denominator. Together with an endonuclease, they can move by "copy and paste. " Based on the presence of an envelope gene, they can be divided into long terminal repeat (LTR) positive and LTR negative retrotransposons. The former and endogenous retroviruses (ERV) belong to LTR positive elements. Long interspersed nuclear elements (LINE), short interspersed nuclear elements (SINE) and SVA elements (SINE-R, VNTR and Alu) belong to LTR negative elements [32] [33] [34] [35] . The LINE contain at least two open reading frames (ORFs), ORF1, coding for a nucleic acid binding protein with chaperone activity (ORF1p) and ORF2, which codes for a reverse transcriptase/endonuclease (ORF2p) [35, 36] . Importantly, RE make up 50 -70% of the human genome [37, 38] . About 20% of the genome is made up from LINE sequences (c. 500,000 copies), of which more than 100 LINE1 family members are still intact and about 68 active in humans. The LINE1 show strong interpersonal differences [39, 40] and an age-dependent expression pattern [41] [42] [43] . By comparison, ERV make up about 8% of the human genome. Despite -similar to LINE -predominant inactivation, there are still hundreds of intact viral promoters and open reading frames from which the expression of ERV transcripts and proteins is possible [44] [45] [46] . The RE activation is known from many viral infections, such as HIV [47] , dengue [48] , influenza A [48] , Zika virus [48] , West Nile virus [48] , measles [48] , Epstein-Barr virus [49] and cytomegalovirus [50] . Therefore, I looked for the relationship of coronaviruses (CoV) to human RE based on genome, transcriptome, epitope and peptide array data. Here, transcriptome analysis coincidentally revealed many RE-identical sequences and shared epitopes in the CoV family members investigated, such as SARS-CoV-2, MERS-CoV and HKU1. To the best of my knowledge, these findings have never been reported. Importantly, epitopes are shared between human LINE1-and SARS-CoV-2 proteins and antibodies against some of these epitopes have been found to be correlated to COVID-19's severity. In addition, RE are expressed in healthy controls and deregulated in COVID-19 patients, as well as in SARS-CoV-2-infected human cells.

The CoV genomes harbour a large number of RE-identical sequences. Several of these sequences represent shared RE-SARS-CoV-2 epitopes. Importantly, antibodies against some of these epitopes are correlated to the severity of COVID-19. In addition, RE are widely expressed in healthy controls and deregulated in COVID-19 patients, as well as in SARS-CoV-2-infected human cells.

A sequence identity (≥12 bp, range 12 -35 bp, Fig. 1A ) of human RE sequences to CoV genomes from SARS-CoV-2, SARS-CoV-1, MERS-CoV, NL63, 229E, OC43, HKU1, bat CoV RA13591, bat CoV RATG13 and bat CoV RSSHC014 was found by sequence alignment of human RE sequences and different CoV genomes (Figs. 1 and 2, Table 1 ). Very high counts of RE-identical sequences in CoV were seen at ≥12, ≥ 15 and ≥ 18 bp (Table 1) .

A cut-off ≥18 bp (correlating to potential epitopes of at least 6 aa) was chosen for downstream analysis for sensitivity and epitope size reasons. A 6 aa cut-off corresponds well to a known immuno-relevant linear epitope length of 4 -12 aa, as about 50% of them have a length ≤ 8 aa (about 25% ≤ 6 aa, and only a few of 4 aa) [51] . At this cut-off point, the majority of RE-identical sequences are seen in HKU1 (332), followed by NL63 (206) and SARS-CoV-2 (191) ( Fig. 2A and B, Table 1 ). SARS-CoV-2 and RE sequence data were further explored by "LAST" in order to allow single nucleotide polymorphisms to be included, thereby alignments to RE sequences up to 35 bp were seen (Supplementary Table 2 ). In the RE-CoV data, LINE1 represent the majority of all shared sequences, while alignment to ERV sequences is a relevant minority and includes the 35 bp hits (Fig. 1B, Supplementary   Fig. 2 Sequence alignments of CoV genomes to retroelements by nucmer (cut-off ≥18 bp). A. Proportion of LINE1 (L1) and endogenous retrovirus sequences, showing a dominance of L1 sequences in all virus genomes (nucmer) analysed. B. Dot plot of shared RE sequences in CoV genomes, showing the highest RE-identical sequences in HKU1, followed by NL63 and SARS-CoV-2 (nucmer). Each dot represents an ≥18 bp retroelement sequence also found in the respective CoV genome Tables 1 and 2 ). In conclusion, genome analysis revealed the presence of many short RE-identical sequences in CoV genomes, including SARS-CoV-2.

Subsequently, all RE-identical sequences ≥18 bp were compared to the coding regions of the genome of SARS-CoV-2. Accordingly, 70 sequences showing identical aa sequences in CoV and RE were identified (Supplementary Table 1 ). These sequences were then compared to results from a peptide array, which investigated epitope signatures in COVID-19 patients (severe vs. mild) [52] . An overlap of human LINE1 proteins to SARS-CoV-2 epitopes from the RNA-dependent RNA polymerase (RdRp), helicase and 2′-O-ribose methyltransferase was detected for epitopes targeted with > 2-fold elevated antibody levels in severe cases (Fig. 3) . Importantly, antibodies targeting an epitope of the SARS-CoV-2 RdRp polymerase, which is identical to an epitope of the LINE1 ORF2p endonuclease domain, were 39-fold elevated in severely compared to only mildly affected COVID-19 patients (Fig. 3A) . The same is seen with antibodies targeting the shared CoV-RE epitopes from the 2′-O-ribose methyltransferase (Fig. 3C ) and helicase (Fig. 3D) . The latter is also a known B cell epitope, aa "PARA-RVECFDKFKV" (the known B cell epitope is depicted in bold) [53] . Many other shared RE-CoV peptides (similar to those displayed in Fig. 3B) were not targeted by antibodies in severe vs. mild COVID-19 (Supplementary  Table 2 ), but some are known as T cell epitopes, such as the one present in all three chains of the spike protein shown in Fig. 3B (aa VKQIYKTPPIKDF, the known T cell epitope sequence is depicted in bold) [54] . Taken together, SARS-CoV-2 and RE share peptide sequences, of which some are epitopes correlated to COVID-19 severity.

An RE analysis of COVID-19 patient data (bronchoalveolar lavage fluid, BALF), SARS-CoV-2 infected lung epithelial cells and SARS-CoV-2 infected macrophages was performed to explore the presence of and changes in RE expression after SARS-CoV-2 infection. Infection resulted in a highly significant (adjusted p-value ≤0.05) and relevant (fold change ≥2) deregulation of human RE in all samples. Transcriptome data from COVID-19 patients' BALF compared to healthy controls shows an upregulation of 2035 and downregulation of 3144 RE (Fig. 4A) . Among the top deregulated RE are mainly LINE1 (Fig. 4D) . SARS-CoV-2-infected epithelial lung cells (Calu-3) show 34 up-and 29 downregulated RE (Fig. 4E) , while infected human macrophages have 8 upand 24 downregulated RE. Among the top de-regulated RE for both are also mainly LINE1 (Fig. 4E, F) .

In conclusion, RE are expressed in COVID-19 patients and human cells and become deregulated after SARS-CoV-2 infection, showing mainly changes in LINE1 expression.

The factors driving the late phase of COVID-19 are still not fully understood [11, 12] . However, there is evidence that autoantibodies and autoreactive lymphocytes could contribute to the disease's final outcome [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] . Therefore, the question of autoantibody formation in COVID-19 has to be asked. The employment of a comprehensive RE database revealed many RE-identical sequences in ten CoV family members investigated, such as in SARS-CoV-2, MERS-CoV and HKU1 ( Figs. 1 and 2) . Crucially, it was found that the LINE1 proteins ORF1p and ORF2p have peptides identical to SARS-CoV-2 epitopes (Fig. 3) , and that some of these epitopes are associated with COVID-19's severity, as shown by correlation to COVID-19 patients' antibody titres (Fig. 3 ). In addition, RE are deregulated in COVID-19 patients (Fig. 4A) , as well as SARS-CoV-2-infected human epithelial lung cells and macrophages ( Fig. 4B and C) , which has occasionally been reported in the last few months for cell lines and patients [28] [29] [30] [31] . Among the analysed REs, LINE1 are strongly represented in all results (Figs. 2, 3 and 4, Supplementary Table 1 and 2) . The LINE1 code for at least a nucleic acid binding protein with chaperone activity (ORF1p) and a reverse transcriptase/endonuclease (ORF2p). Importantly, autoantibodies targeting the LINE1 ORF2p endonuclease domain have been reported in 41% of SARS-CoV-1 patients [55] . The RE are also targeted by autoantibodies in several connective tissue diseases, for example, antibodies against LINE1's ORF1p or ERV HERV-K's envelope protein have been described in patients with systemic lupus erythematosus, lupus nephritis, rheumatoid arthritis, Sjogren's syndrome and mixed connective tissue disease [56] [57] [58] [59] [60] [61] [62] [63] [64] [65] . Relating to SARS, the autoantibodies' target, LINE1 ORF2p, was prominently stained post-mortem in lung macrophages (residing in blood vessels), leading the authors to suspect a build-up of autoreactive CD4+ Th cells and, thus, an autoimmune loop in SARS [55] . Importantly, there is also increasing evidence for an autoimmune pathogenesis in severe COVID-19 [13-27, 66, 67] . One explanation for autoantibody formation is by molecular mimicry, i.e. shared epitopes between pathogens and hosts [68] [69] [70] [71] [72] . The evolution of mimicry epitopes in pathogens could be based on chance. However, although the RE-identical sequences in CoV observed are short (12 -35 bp) , the sequence lengths observed make formation by chance highly unlikely. Exemplarily, taking the genetic code (A, T, C, G) raised to a sequence of 18 bp (4 18 ) results in 68,719,476,736 possible bp combinations, thus, the chance of getting one identical sequence is 1:69 billion. Additionally, a myriad of 12 bp events (Table 1) occurring by chance is stochastically very unlikely (4 12 = 16,777,216) at more than 18,000 events. Moreover, an observed 35 bp hit such as ERVL_Xq21.31b (4 35 ) corresponds to 1.18 E21 possible bp combinations, thus, the chance of getting an identical sequence is 1:1.1 trilliardwithout accounting for all the other matching sequences. Therefore, recombination activities more probably account for the phenomena observed. The exchange of genetic material by recombination in RNA viruses is generally associated with virulence, host range and host response [73] . It is known that recombination in CoV can take place during co-infections at a high frequency by homologous and non-homologous recombination [74] [75] [76] . Mechanistically, an explanation could be the switching of the RdRp between multiple available RNA strands during replication [77] . This could have happened in a CoV host/ancestor with relevant LINE1 expression, as this is possible in some bat species. The black-bearded tomb bat (Taphozous melanopogon), for example, harbours two active LINE families [78] and shows relevant SARS-CoV-2 infection efficiency [79] . Moreover, lots of ERV families also reside in bats [80] . Therefore, serial acquisition of RE sequences, possibly taken from CoV in host animals (starting many million years ago) is a feasible scenario. Relating to the rather short sequence lengths observed, there might be an evolutionary functional constraint working against the uptake of longer RE sequences, but a benefit for the virus by coating itself with host self-antigens ("self-peptide coat"). This would dampen the innate and adaptive immune response by the presentation of "viral but self-like" peptides. The consequence of this hypothesis is in line with the view of autoimmune disease as a breakdown of self-tolerance [81, 82] . Based on the findings, autoantibodies targeting human RE could be a factor in CoV-induced disease, like COVID-19. However, this report has limitations, as the data basis for a more extensive analysis of anti-RE autoantibodies in COVID-19 still does not exist.

In conclusion, it was found that CoV -including SARS-CoV-2 -harbour many RE-identical sequences, and that some of these sequences are part of SARS-CoV-2 epitopes associated with COVID-19 severity. MG916904.1 = Ra1359) were downloaded from GenBank (https:// www. ncbi. nlm. nih. gov/ genba nk/). Retro.hg38.v1 (https:// github. com/ mlben dall/ teles cope_ annot ation_ db/ tree/ master/ builds) was employed as an RE database. The database contains 28.513 RE and is made of "RepeatMasker" hits for 60 HERV families (RepeatMasker Open-4.0, http:// www. repea tmask er. org/) and all LINE elements from "L1base v2" (https:// l1base. chari te. de/) [83] . Alignment of the retro. hg38.v1 database to CoV genomes was done by the genome sequence aligner "nucmer" [84] (4.0.0beta2) on galaxy. org [85] and a local installation of "LAST" (v1250), a programme for genome scale sequence comparison [86] . The minimum sequence length cut-off (with 100% sequence identity) was stepwise chosen at 12, 15, 18, 21, 24 , and ≥ 27, based on an immuno-relevant epitope size of about 4 -12 amino acids (aa) (many epitopes are less than 8 aa, about 25% ≤ 6 aa, but only a few at 4 aa [51] ). The nucmer "-b" and "-L" variables were used accordingly, and "Show-Coords" as well as "Mummerplot" from the "MUMmer 4" package [84] were employed to extract and plot data. Regarding to "LAST, " firstly, an RE database was built ("lastdb -uNEAR -c RE_ db retro.hg38.v1.fa") and then CoV genomes were compared to the RE database ("lastal -D100 RE_db CoV_genome.fa > RE_db_CoV.maf").

The SARS-CoV-2 epitope-specific antibody data (IgG) in severely vs. mildly affected COVID-19 patients are from Schwarz et al. [52] "Peptide microarray data -severe vs. mild -IgG, " with the peptides: 1060 (NSP12, QTVK-PGNFNKDFYDF, LogFC 5.3, p-value 2.4E-04, FDR-adj. p-value 2.8E-02), 1243 (NSP16, ENDSKEGFFTYICGF, LogFC 2.2, p-value 4.0E-02, FDR-adj. p-value 5.2E-01), 1227 (NSP13, IPARARVECFDKFKV, LogFC − 0.9, p-value 3.2E-01, FDR-adj. p-value 5.3E-01) and 1690 (Spike, AQVKQIYKTPPIKDF, LogFC 0.2, p-value 8.3E-01, FDR-adj. p-value 8.5E-01). "L1base v2" was used for comparison with coding LINE1 sequences (https:// l1base. chari te. de/) [83] . Known SARS-CoV-2 B-and T-cell epitopes are from Phan et al. [53] and Griffoni et al. [54] . The PDB data for the SARS-CoV-2 RdRp (PDB ID: 7BW4), helicase (PDB ID: 7NNG), 2′-O-ribose methyltransferase (PDB ID: 7JYY) and -spike protein (PDB ID: 7LSS) were downloaded from https:// www. rcsb. org and epitopes displayed by "UCSF Chimera v1.15" (for Mac OS) [87] .

Total RNA sequencing data from SARS-CoV-2-infected macrophages (BioProject ID PRJNA637580, Sequence Read Archive (SRA) ID mock: SRR11934391, SRR11934392, SRR11934393, infected: SRR11934394, SRR11934395, SRR11934396) [88] , Calu-3 adrenocarcinomic lung epithelial cells (PRJNA615032, mock: SRR11517744, SRR11517745, SRR11517746, infected: SRR11517747, SRR11517748, SRR11517749) [89] and bronchoalveolar lavage (BALF) samples from intensive care COVID-19 patients (PRJNA605983SRA, SRA: SRR11092056, SRR11092057, SRR11092058, SRR11 092059, SRR11092060, SRR11092061, SRR11092062, SRR11092063, SRR11092064) [90] compared to healthy controls (PRJNA316136, SRA: SRR3286988, SRR3286989, SRR3286990, SRR3286991, SRR5515942, SRR5515943, SRR5515944) [91] were downloaded from SRA (https:// www. ncbi. nlm. nih. gov/ sra), quality controlled by FastQC (Babraham Institute, Cambridge, UK, http:// www. bioin forma tics. babra ham. ac. uk/ proje cts/ fastqc/) and Illumina adapters trimmed by Trimmomatic [92] . Salmon [93] and DESeq2 [94] were employed for differential RE analysis, with standard parameters after indexing the retro.hg38.v1 database ("salmon index -t retro.hg38.v1.fa -i retro.hg38.v1_ index -k 31"). Heatmaps were done by iDEP v0.92 [95] and graphs by GraphPad Prism software version 8.0 for OS X (GraphPad Software Inc., USA). 

A novel coronavirus from patients with pneumonia in China

An interactive web-based dashboard to track COVID-19 in real time

Virology, transmission, and pathogenesis of SARS-CoV-2

A biological profile for diagnosis and outcome of COVID-19 patients

A SARS-CoV-2 cytopathicity dataset generated by high-content screening of a large drug repurposing collection

Aprotinin inhibits SARS-CoV-2 replication

Growth factor receptor signaling inhibition prevents SARS-CoV-2 replication

Proteomics of SARS-CoV-2-infected host cells reveals therapy targets

Identification of antiviral drug candidates against SARS-CoV-2 from FDA-approved drugs

FDA-approved drugs with potent in vitro antiviral activity against severe acute respiratory syndrome coronavirus 2

Postmortem findings associated with SARS-CoV-2

Postmortem examination of COVID-19 patients reveals diffuse alveolar damage with severe capillary congestion and variegated findings in lungs and other organs suggesting vascular dysfunction

Rogue antibodies could be driving severe COVID-19

Autoantibodies against type I IFNs in patients with life-threatening COVID-19

Diverse functional autoantibodies in patients with COVID-19

COVID-19: infection or autoimmunity

Autoimmunity to the lung protective phospholipid-binding protein Annexin A2 predicts mortality among hospitalized COVID-19 patients

Covid-19 and autoimmunity

Prothrombotic autoantibodies in serum from patients hospitalized with COVID-19

Can molecular mimicry explain the cytokine storm of SARS-CoV-2?: an in silico approach

Guillain-Barré syndrome, SARS-CoV-2 and molecular mimicry

High frequency of cerebrospinal fluid autoantibodies in COVID-19 patients with neurological symptoms

Reaction of human monoclonal antibodies to SARS-CoV-2 proteins with tissue antigens: implications for autoimmune diseases

The JANUS of chronic inflammatory and autoimmune diseases onset during COVID-19 -a systematic review of the literature

SARS-CoV-2 genome-wide T cell epitope mapping reveals immunodominance and substantial CD8+ T cell activation in COVID-19 patients

SARS-Cov-2 fulminant myocarditis: an autopsy and histopathological case study

Alterations in lymphocyte subsets and monocytes in patients diagnosed with SARS-CoV-2 pneumonia: a mini review of the literature

Evidence of the pathogenic HERV-W envelope expression in T lymphocytes in association with the respiratory outcome of COVID-19 patients

Genomic study of COVID-19

Corona virus excludes its origin from recombination or characterized biological sources and suggests a role for HERVS in its wide range symptoms

Human endogenous retrovirus K activation in the lower respiratory tract of severe COVID-19 patients associates with early mortality

Alu retrotransposons and COVID-19 susceptibility and morbidity

Cas9 targeted enrichment of mobile elements using nanopore sequencing

Variable number tandem repeats -their emerging role in sickness and health

New insights into the functional role of retrotransposon dynamics in mammalian somatic cells

The long (LINEs) and the short (SINEs) of it: altered methylation as a precursor to toxicity

Initial sequencing and analysis of the human genome

Repetitive elements may comprise over two-thirds of the human genome

Sequencing, identification and mapping of primed L1 elements (SIMPLE) reveals significant variation in full length L1 elements between individuals

LINE1 Derepression in aged wild-type and SIRT6-deficient mice drives inflammation

Aging-associated distinctive DNA methylation changes of LINE-1 retrotransposons in pure cell-free DNA from human blood

Skeletal muscle LINE-1 ORF1 mRNA is higher in older humans but decreases with endurance exercise and is negatively associated with higher physical activity

Identification of endogenous retroviral reading frames in the human genome

Endogenous viruses: insights into viral evolution and impact on host biology

HERV-K(HML-2), the best preserved family of HERVs: Endogenization, expression, and implications in health and disease

Modulation of human endogenous retrovirus (HERV) transcription during persistent and de novo HIV-1 infection

Transcriptome analyses implicate endogenous retroviruses involved in the host antiviral immune system through the interferon pathway

Epstein-Barr virus Transactivates the human endogenous retrovirus HERV-K18 that encodes a Superantigen

Human cytomegalovirus (HCMV) induces human endogenous retrovirus (HERV) transcription. -PubMed -NCBI

High-resolution mapping of linear antibody epitopes using ultrahighdensity peptide microarrays*

SARS-CoV-2 proteome-wide analysis revealed significant epitope signatures in COVID-19 patients

In silico detection of SARS-CoV-2 specific B-cell epitopes and validation in ELISA for serological diagnosis of COVID-19

A sequence homology and Bioinformatic approach can predict candidate targets for immune responses to SARS-CoV-2

Human LINE1 endonuclease domain as a putative target of SARS-associated autoantibodies involved in the pathogenesis of severe acute respiratory syndrome

A role for human endogenous retrovirus-K (HML-2) in rheumatoid arthritis: investigating mechanisms of pathogenesis

Detection of serum antibodies to retroviral proteins in patients with primary sjögren's syndrome (autoimmune exocrinopathy)

Serum antibody to retroviral gag proteins in systemic sclerosis

A conserved idiotype and antibodies to retroviral proteins in systemic lupus erythematosus

Selective antibody reactivity with peptides from human endogenous retroviruses and nonviral poly(amino acids) in patients with systemic lupus erythematosus

Detection of antibodies to a recombinant gag protein derived from human endogenous retrovirus clone 4-1 in autoimmune diseases

How retroviruses and Retrotransposons in our genome may contribute to autoimmunity in Rheumatological conditions

High prevalence and disease correlation of autoantibodies against p40 encoded by long interspersed nuclear elements in systemic lupus Erythematosus

The Ro60 autoantigen binds endogenous retroelements and regulates inflammatory gene expression

Expression of long interspersed nuclear element 1 Retroelements and induction of type I interferon in patients with systemic autoimmune disease

COVID-19 Immunobiology: lessons learned

Adaptive immunity and the risk of autoreactivity in COVID-19

Mimicking the way to autoimmunity: an evolving theory of sequence and structural homology

Molecular mimicry and immune-mediated diseases

Molecular mimicry: a mechanism for autoimmune injury1

Molecular mimicry in biological adaptation

Molecular mimicry revisited

RNA recombination enhances adaptability and is required for virus spread and virulence

The molecular biology of coronaviruses

Recombination, reservoirs, and the modular spike: mechanisms of coronavirus cross-species transmission▿

High-frequency RNA recombination of murine coronaviruses

Tracing the origins of SARS-COV-2 in coronavirus phylogenies: a review

Flying around in the genome: characterization of LINE-1 in Chiroptera

ACE2 receptor usage reveals variation in susceptibility to SARS-CoV and SARS-CoV-2 infection among bat species

Retroviruses of bats: a threat waiting in the wings?

Breakdown of self-tolerance and the pathogenesis of autoimmunity

L1Base 2: more retrotransposition-active LINE-1s, more mammalian genomes

MUMmer4: A fast and versatile genome alignment system

Using galaxy to perform large-scale interactive data analyses-an update

Adaptive seeds tame genomic sequence comparison

UCSF Chimera--a visualization system for exploratory research and analysis

Cardiomyocytes recruit monocytes upon SARS-CoV-2 infection by secreting CCL2. Stem Cell Rep

Imbalanced host response to SARS-CoV-2 drives development of COVID-19

A pneumonia outbreak associated with a new coronavirus of probable bat origin

Bronchoalveolar lavage (BAL) cells in idiopathic pulmonary fibrosis express a complex pro-inflammatory, pro-repair, angiogenic activation pattern, likely associated with macrophage iron accumulation

Trimmomatic: a flexible trimmer for Illumina sequence data

Salmon provides fast and bias-aware quantification of transcript expression

Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2

iDEP: an integrated web application for differential expression and pathway analysis of RNA-Seq data

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations

The author declares no competing interests.