Impact of natural selection on global patterns of genetic variation and association with clinical phenotypes at genes involved in SARS-CoV-2 infection.

Human genomic diversity has been shaped by both ancient and ongoing challenges from viruses. The current coronavirus disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has had a devastating impact on population health. However, genetic diversity and evolutionary forces impacting host genes related to SARS-CoV-2 infection are not well understood. We investigated global patterns of genetic variation and signatures of natural selection at host genes relevant to SARS-CoV-2 infection (angiotensin converting enzyme 2 [ACE2], transmembrane protease serine 2 [TMPRSS2], dipeptidyl peptidase 4 [DPP4], and lymphocyte antigen 6 complex locus E [LY6E]). We analyzed data from 2,012 ethnically diverse Africans and 15,977 individuals of European and African ancestry with electronic health records and integrated with global data from the 1000 Genomes Project. At ACE2, we identified 41 nonsynonymous variants that were rare in most populations, several of which impact protein function. However, three nonsynonymous variants (rs138390800, rs147311723, and rs145437639) were common among central African hunter-gatherers from Cameroon (minor allele frequency 0.083 to 0.164) and are on haplotypes that exhibit signatures of positive selection. We identify signatures of selection impacting variation at regulatory regions influencing ACE2 expression in multiple African populations. At TMPRSS2, we identified 13 amino acid changes that are adaptive and specific to the human lineage compared with the chimpanzee genome. Genetic variants that are targets of natural selection are associated with clinical phenotypes common in patients with COVID-19. Our study provides insights into global variation at host genes related to SARS-CoV-2 infection, which have been shaped by natural selection in some populations, possibly due to prior viral infections.

Redondovirus Diversity and Evolution on Global, Individual, and Molecular Scales.

Redondoviridae is a newly established family of circular Rep-encoding single-stranded (CRESS) DNA viruses found in the human ororespiratory tract. Redondoviruses were previously found in ∼15% of respiratory specimens from U.S. urban subjects; levels were elevated in individuals with periodontitis or critical illness. Here, we report higher redondovirus prevalence in saliva samples: four rural African populations showed 61 to 82% prevalence, and an urban U.S. population showed 32% prevalence. Longitudinal, limiting-dilution single-genome sequencing revealed diverse strains of both redondovirus species (Brisavirus and Vientovirus) in single individuals, persistence over time, and evidence of intergenomic recombination. Computational analysis of viral genomes identified a recombination hot spot associated with a conserved potential DNA stem-loop structure. To assess the possible role of this site in recombination, we carried out in vitro studies which showed that this potential stem-loop was cleaved by the virus-encoded Rep protein. In addition, in reconstructed reactions, a Rep-DNA covalent intermediate was shown to mediate DNA strand transfer at this site. Thus, redondoviruses are highly prevalent in humans, found in individuals on multiple continents, heterogeneous even within individuals and encode a Rep protein implicated in facilitating recombination. IMPORTANCERedondoviridae is a recently established family of DNA viruses predominantly found in the human respiratory tract and associated with multiple clinical conditions. In this study, we found high redondovirus prevalence in saliva from urban North American individuals and nonindustrialized African populations in Botswana, Cameroon, Ethiopia, and Tanzania. Individuals on both continents harbored both known redondovirus species. Global prevalence of both species suggests that redondoviruses have long been associated with humans but have remained undetected until recently due to their divergent genomes. By sequencing single redondovirus genomes in longitudinally sampled humans, we found that redondoviruses persisted over time within subjects and likely evolve by recombination. The Rep protein encoded by redondoviruses catalyzes multiple reactions in vitro, consistent with a role in mediating DNA replication and recombination. In summary, we identify high redondovirus prevalence in humans across multiple continents, longitudinal heterogeneity and persistence, and potential mechanisms of redondovirus evolution by recombination.

Mitochondrial genetic variation in human bioenergetics, adaptation, and adult disease.

OBJECTIVES: Mitochondria are critical for the survival of eukaryotic organisms due to their ability to produce cellular energy, which drives virtually all aspects of host biology. However, the effects of mitochondrial DNA (mtDNA) variation in relation to disease etiology and adaptation within contemporary global human populations remains incompletely understood. METHODS: To develop a more holistic understanding of the role of mtDNA diversity in human adaptation, health, and disease, we investigated mitochondrial biology and bioenergetics. More specifically, we synthesized details from studies of mitochondrial function and variation in the context of haplogroup background, climatic adaptation, and oxidative disease. RESULTS: The majority of studies show that mtDNA variation arose during modern human dispersal around the world. Some of these variants appear to have been positively selected for their adaptiveness in colder climates, with these sequence changes having implications for tissue-specific function and thermogenic capacity. In addition, many variants modulating energy production are also associated with damaging metabolic byproducts and mitochondrial dysfunction, which, in turn, are implicated in the onset and severity of several different adult mitochondrial diseases. Thus, mtDNA variation that governs bioenergetics, metabolism, and thermoregulation may potentially have adverse consequences for human health, depending on the genetic background and context in which it occurs. CONCLUSIONS: Our review suggests that the mitochondrial research field would benefit from independently replicating mtDNA haplogroup-phenotype associations across global populations, incorporating potentially confounding environmental, demographic, and disease covariates into studies of mtDNA variation, and extending association-based studies to include analyses of complete mitogenomes and assays of mitochondrial function.

Facing the facts: The Runx2 gene is associated with variation in facial morphology in primates.

The phylogenetic and adaptive factors that cause variation in primate facial form-including differences among the major primate clades and variation related to feeding and/or social behavior-are relatively well understood. However, comparatively little is known about the genetic mechanisms that underlie diversity in facial form in primates. Because it is essential for osteoblastic differentiation and skeletal development, the runt-related transcription factor 2 (Runx2) is one gene that may play a role in these genetic mechanisms. Specifically, polymorphisms in the QA ratio (determined by the ratio of the number of polyglutamines to polyalanines in one functional domain of Runx2) have been shown to be correlated with variation in facial length and orientation in other mammal groups. However, to date, the relationship between variation in this gene and variation in facial form in primates has not been explicitly tested. To test the hypothesis that the QA ratio is correlated with facial form in primates, the current study quantified the QA ratio, facial length, and facial angle in a sample of 33 primate species and tested for correlation using phylogenetic generalized least squares. The results indicate that the QA ratio of the Runx2 gene is positively correlated with variation in relative facial length in anthropoid primates. However, no correlation was found in strepsirrhines, and there was no correlation between facial angle and the QA ratio in any groups. These results suggest that, in primates, the QA ratio of the Runx2 gene may play a role in modulating facial size, but not facial orientation. This study therefore provides important clues about the genetic and developmental mechanisms that may underlie variation in facial form in primates.

Regulatory element copy number differences shape primate expression profiles.

Gene expression differences are shaped by selective pressures and contribute to phenotypic differences between species. We identified 964 copy number differences (CNDs) of conserved sequences across three primate species and examined their potential effects on gene expression profiles. Samples with copy number different genes had significantly different expression than samples with neutral copy number. Genes encoding regulatory molecules differed in copy number and were associated with significant expression differences. Additionally, we identified 127 CNDs that were processed pseudogenes and some of which were expressed. Furthermore, there were copy number-different regulatory regions such as ultraconserved elements and long intergenic noncoding RNAs with the potential to affect expression. We postulate that CNDs of these conserved sequences fine-tune developmental pathways by altering the levels of RNA.

Radiologist-supervised Transfer Learning: Improving Radiographic Localization of Pneumonia and Prognostication of Patients With COVID-19.

PURPOSE: To assess the potential of a transfer learning strategy leveraging radiologist supervision to enhance convolutional neural network-based (CNN) localization of pneumonia on radiographs and to further assess the prognostic value of CNN severity quantification on patients evaluated for COVID-19 pneumonia, for whom severity on the presenting radiograph is a known predictor of mortality and intubation. MATERIALS AND METHODS: We obtained an initial CNN previously trained to localize pneumonia along with 25,684 radiographs used for its training. We additionally curated 1466 radiographs from patients who had a computed tomography (CT) performed on the same day. Regional likelihoods of pneumonia were then annotated by cardiothoracic radiologists, referencing these CTs. Combining data, a preexisting CNN was fine-tuned using transfer learning. Whole-image and regional performance of the updated CNN was assessed using receiver-operating characteristic area under the curve and Dice. Finally, the value of CNN measurements was assessed with survival analysis on 203 patients with COVID-19 and compared against modified radiographic assessment of lung edema (mRALE) score. RESULTS: Pneumonia detection area under the curve improved on both internal (0.756 to 0.841) and external (0.864 to 0.876) validation data. Dice overlap also improved, particularly in the lung bases (R: 0.121 to 0.433, L: 0.111 to 0.486). There was strong correlation between radiologist mRALE score and CNN fractional area of involvement (ρ=0.85). Survival analysis showed similar, strong prognostic ability of the CNN and mRALE for mortality, likelihood of intubation, and duration of hospitalization among patients with COVID-19. CONCLUSIONS: Radiologist-supervised transfer learning can enhance the ability of CNNs to localize and quantify the severity of disease. Closed-loop systems incorporating radiologists may be beneficial for continued improvement of artificial intelligence algorithms.

Lifestyle and the presence of helminths is associated with gut microbiome composition in Cameroonians.

BACKGROUND: African populations provide a unique opportunity to interrogate host-microbe co-evolution and its impact on adaptive phenotypes due to their genomic, phenotypic, and cultural diversity. We integrate gut microbiome 16S rRNA amplicon and shotgun metagenomic sequence data with quantification of pathogen burden and measures of immune parameters for 575 ethnically diverse Africans from Cameroon. Subjects followed pastoralist, agropastoralist, and hunter-gatherer lifestyles and were compared to an urban US population from Philadelphia. RESULTS: We observe significant differences in gut microbiome composition across populations that correlate with subsistence strategy and country. After these, the variable most strongly associated with gut microbiome structure in Cameroonians is the presence of gut parasites. Hunter-gatherers have high frequencies of parasites relative to agropastoralists and pastoralists. Ascaris lumbricoides, Necator americanus, Trichuris trichiura, and Strongyloides stercoralis soil-transmitted helminths ("ANTS" parasites) significantly co-occur, and increased frequency of gut parasites correlates with increased gut microbial diversity. Gut microbiome composition predicts ANTS positivity with 80% accuracy. Colonization with ANTS, in turn, is associated with elevated levels of TH1, TH2, and proinflammatory cytokines, indicating an association with multiple immune mechanisms. The unprecedented size of this dataset allowed interrogation of additional questions-for example, we find that Fulani pastoralists, who consume high levels of milk, possess an enrichment of gut bacteria that catabolize galactose, an end product of lactose metabolism, and of bacteria that metabolize lipids. CONCLUSIONS: These data document associations of bacterial microbiota and eukaryotic parasites with each other and with host immune responses; each of these is further correlated with subsistence practices.

Population structure of human gut bacteria in a diverse cohort from rural Tanzania and Botswana.

BACKGROUND: Gut microbiota from individuals in rural, non-industrialized societies differ from those in individuals from industrialized societies. Here, we use 16S rRNA sequencing to survey the gut bacteria of seven non-industrialized populations from Tanzania and Botswana. These include populations practicing traditional hunter-gatherer, pastoralist, and agropastoralist subsistence lifestyles and a comparative urban cohort from the greater Philadelphia region. RESULTS: We find that bacterial diversity per individual and within-population phylogenetic dissimilarity differs between Botswanan and Tanzanian populations, with Tanzania generally having higher diversity per individual and lower dissimilarity between individuals. Among subsistence groups, the gut bacteria of hunter-gatherers are phylogenetically distinct from both agropastoralists and pastoralists, but that of agropastoralists and pastoralists were not significantly different from each other. Nearly half of the Bantu-speaking agropastoralists from Botswana have gut bacteria that are very similar to the Philadelphian cohort. Based on imputed metagenomic content, US samples have a relative enrichment of genes found in pathways for degradation of several common industrial pollutants. Within two African populations, we find evidence that bacterial composition correlates with the genetic relatedness between individuals. CONCLUSIONS: Across the cohort, similarity in bacterial presence/absence compositions between people increases with both geographic proximity and genetic relatedness, while abundance weighted bacterial composition varies more significantly with geographic proximity than with genetic relatedness.

Investigating zoonotic infection barriers to ape Plasmodium parasites using faecal DNA analysis.

African apes are endemically infected with numerous Plasmodium spp. including close relatives of human Plasmodium falciparum, Plasmodium vivax, Plasmodium ovale, and Plasmodium malariae. Although these ape parasites are not believed to pose a zoonotic threat, their ability to colonise humans has not been fully explored. In particular, it remains unknown whether ape parasites are able to initiate exo-erythrocytic replication in human hepatocytes following the bite of an infective mosquito. Since animal studies have shown that liver stage infection can result in the excretion of parasite nucleic acids into the bile, we screened faecal samples from 504 rural Cameroonians for Plasmodium DNA. Using pan-Laverania as well as P. malariae- and P. vivax-specific primer sets, we amplified human P. falciparum (n = 14), P. malariae (n = 1), and P. ovale wallikeri (n = 1) mitochondrial sequences from faecal DNA of 15 individuals. However, despite using an intensified PCR screening approach we failed to detect ape Laverania, ape P. vivax or ape P. malariae parasites in these same subjects. One faecal sample from a hunter-gatherer contained a sequence closely related to the porcupine parasite Plasmodium atheruri. Since this same faecal sample also contained porcupine mitochondrial DNA, but a matching blood sample was Plasmodium-negative, it is likely that this hunter-gatherer consumed Plasmodium-infected bushmeat. Faecal Plasmodium detection was not secondary to intestinal bleeding and/or infection with gastrointestinal parasites, but indicative of blood parasitaemia. Quantitative PCR identified 26-fold more parasite DNA in the blood of faecal Plasmodium-positive than faecal Plasmodium-negative individuals (P = 0.01). However, among blood-positive individuals only 10% - 20% had detectable Plasmodium sequences in their stool. Thus, faecal screening of rural Cameroonians failed to uncover abortive ape Plasmodium infections, but detected infection with human parasites, albeit with reduced sensitivity compared with blood analysis.

Inferences of African evolutionary history from genomic data.

Africa is the origin of anatomically modern humans and a continent of linguistic, cultural, environmental, phenotypic, and genetic diversity. However, African populations remain underrepresented in genetic studies, which have largely focused on individuals with European and Asian ancestry. The expansion of high-throughput 'omic' technologies to interrogate multiple tissue types across many biomolecules-DNA, proteins, epigenetic modifications, metabolites, and others-has heralded a new era of investigation into African history. In this review, we summarize how some of these recent advances have been applied to contemporary sub-Saharan African populations to inform studies on human origins and adaptation.