Studies showed that the gastrointestinal (GI) tract is one of the most important pathways for SARS-CoV-2 infection and coronavirus disease 2019 (COVID-19). As SARS-CoV-2 cellular entry depends on the ACE2 receptor and TMPRSS2 priming of the spike protein, it is important to understand the molecular mechanisms through which these two proteins and their cognate transcripts interact and influence the pathogenesis of COVID-19. In this study, we quantified the expression, associations, genetic modulators, and molecular pathways for Tmprss2 and Ace2 mRNA expressions in GI tissues using a systems genetics approach and the expanded family of highly diverse BXD mouse strains. The results showed that both Tmprss2 and Ace2 are highly expressed in GI tissues with significant covariation. We identified a significant expression quantitative trait locus on chromosome 7 that controls the expression of both Tmprss2 and Ace2. Dhx32 was found to be the strongest candidate in this interval. Co-expression network analysis demonstrated that both Tmprss2 and Ace2 were located at the same module that is significantly associated with other GI-related traits. Protein-protein interaction analysis indicated that hub genes in this module are linked to circadian rhythms. Collectively, our data suggested that genes with circadian rhythms of expression may have an impact on COVID-19 disease, with implications related to the timing and treatment of COVID-19.
Keywords: BXD mice; COVID-19; SARS-CoV-2; co-expression; gastrointestinal tract; microbiota; systems genetics; transcriptome.