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.

Human adaptation, demography and cattle domestication: an overview of the complexity of lactase persistence in Africa.

Lactase persistence (LP) is a genetically-determined trait that is prevalent in African, European and Arab populations with a tradition of animal herding and milk consumption. To date, genetic analyses have identified several common variants that are associated with LP. Furthermore, data have indicated that these functional alleles likely have been maintained in pastoralist populations due to the action of recent selection, exemplifying the ongoing evolution of anatomically modern humans. Additionally, demographic history has also played a role in the geographic distribution of LP and associated alleles in Africa. In particular, the migration of ancestral herders and their subsequent admixture with local populations were integral to the spread of LP alleles and the culture of pastoralism across the continent. The timing of these demographic events was often correlated with known major environmental changes and/or the ability of domesticated cattle to resist/avoid infectious diseases. This review summarizes recent advances in our understanding of the genetic basis and evolutionary history of LP, as well as the factors that influenced the origin and spread of pastoralism in Africa.

Deregulated miRNA Expression in Triple-Negative Breast Cancer of Ancestral Genomic-Characterized Latina Patients.

Among patients with triple-negative breast cancer (TNBC), several studies have suggested that deregulated microRNA (miRNA) expression may be associated with a more aggressive phenotype. Although tumor molecular signatures may be race- and/or ethnicity-specific, there is limited information on the molecular profiles in women with TNBC of Hispanic and Latin American ancestry. We simultaneously profiled TNBC biopsies for the genome-wide copy number and miRNA global expression from 28 Latina women and identified a panel of 28 miRNAs associated with copy number alterations (CNAs). Four selected miRNAs (miR-141-3p, miR-150-5p, miR-182-5p, and miR-661) were validated in a subset of tumor and adjacent non-tumor tissue samples, with miR-182-5p being the most discriminatory among tissue groups (AUC value > 0.8). MiR-141-3p up-regulation was associated with increased cancer recurrence; miR-661 down-regulation with larger tumor size; and down-regulation of miR-150-5p with larger tumor size, high p53 expression, increased cancer recurrence, presence of distant metastasis, and deceased status. This study reinforces the importance of integration analysis of CNAs and miRNAs in TNBC, allowing for the identification of interactions among molecular mechanisms. Additionally, this study emphasizes the significance of considering the patients ancestral background when examining TNBC, as it can influence the relationship between intrinsic tumor molecular characteristics and clinical manifestations of the disease.

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.

A panel of miRNAs as prognostic markers for African-American patients with triple negative breast cancer.

BACKGROUND: To investigate the global expression profile of miRNAs, their impact on cellular signaling pathways, and their association with poor prognostic parameters in African-American (AA) patients with triple negative breast cancer (TNBC). METHODS: Twenty-five samples of AA TNBC patients were profiled for global miRNA expression and stratified considering three clinical-pathological parameters: tumor size, lymph node (LN), and recurrence (REC) status. Differential miRNA expression analysis was performed for each parameter, and their discriminatory power was determined by Receiver Operating Characteristic (ROC) curve analysis. KMplotter was assessed to determine the association of the miRNAs with survival, and functional enrichment analysis to determine the main affected pathways and miRNA/mRNA target interactions. RESULTS: A panel of eight, 23 and 27 miRNAs were associated with tumor size, LN, and REC status, respectively. Combined ROC analysis of two (miR-2117, and miR-378c), seven (let-7f-5p, miR-1255b-5p, miR-1268b, miR-200c-3p, miR-520d, miR-527, and miR-518a-5p), and three (miR-1200, miR-1249-3p, and miR-1271-3p) miRNAs showed a robust discriminatory power based on tumor size (AUC = 0.917), LN (AUC = 0.945) and REC (AUC = 0.981) status, respectively. Enrichment pathway analysis revealed their involvement in proteoglycans and glycan and cancer-associated pathways. Eight miRNAs with deregulated expressions in patients with large tumor size, positive LN metastasis, and recurrence were significantly associated with lower survival rates. Finally, the construction of miRNA/mRNA networks based in experimentally validated mRNA targets, revealed nodes of critical cancer genes, such as AKT1, BCL2, CDKN1A, EZR and PTEN. CONCLUSIONS: Altogether, our data indicate that miRNA deregulated expression is a relevant biological factor that can be associated with the poor prognosis in TNBC of AA patients, by conferring to their TNBC cells aggressive phenotypes that are reflected in the clinical characteristics evaluated in this study.

Population genetic evidence for positive and purifying selection acting at the human IFN-γ locus in Africa.

Despite its critical role in the defense against microbial infection and tumor development, little is known about the range of nucleotide and haplotype variation at IFN-γ, or the evolutionary forces that have shaped patterns of diversity at this locus. To address this gap in knowledge, we examined sequence data from the IFN-γ gene in 1461 individuals from 15 worldwide populations. Our analyses uncovered novel patterns of variation in distinct African populations, including an excess of high frequency-derived alleles, unusually long haplotype structure surrounding the IFN-γ gene, and a "star-like" genealogy of African-specific haplotypes carrying variants previously associated with infectious disease. We also inferred a deep time to coalescence of variation at IFN-γ (~ 0.8 million years ago) and ancient ages for common polymorphisms predating the evolution of modern humans. Taken together, these results are congruent with a model of positive selection on standing variation in African populations. Furthermore, we inferred that common variants in intron 3 of IFN-γ are the likely targets of selection. In addition, we observed a paucity of non-synonymous substitutions relative to synonymous changes in the exons of IFN-γ in African and non-African populations, suggestive of strong purifying selection. Therefore, we contend that positive and purifying selection have influenced levels of diversity in different regions of IFN-γ, implying that these distinct genic regions are, or have been, functionally important. Overall, this study provides additional insights into the evolutionary events that have contributed to the frequency and distribution of alleles having a role in human health and disease.

Correction to: Multiple selective sweeps of ancient polymorphisms in and around LTα located in the MHC class III region on chromosome 6.

After publication of our article [1] we were notified that a few duplicate sentences were included on Figure 3 and Figure 4 legends.

Multiple selective sweeps of ancient polymorphisms in and around LTα located in the MHC class III region on chromosome 6.

BACKGROUND: Lymphotoxin-α (LTα), located in the Major Histocompatibility Complex (MHC) class III region on chromosome 6, encodes a cytotoxic protein that mediates a variety of antiviral responses among other biological functions. Furthermore, several genotypes at this gene have been implicated in the onset of a number of complex diseases, including myocardial infarction, autoimmunity, and various types of cancer. However, little is known about levels of nucleotide variation and linkage disequilibrium (LD) in and near LTα, which could also influence phenotypic variance. To address this gap in knowledge, we examined sequence variation across ~ 10 kilobases (kbs), encompassing LTα and the upstream region, in 2039 individuals from the 1000 Genomes Project originating from 21 global populations. RESULTS: Here, we observed striking patterns of diversity, including an excess of intermediate-frequency alleles, the maintenance of multiple common haplotypes and a deep coalescence time for variation (dating > 1.0 million years ago), in global populations. While these results are generally consistent with a model of balancing selection, we also uncovered a signature of positive selection in the form of long-range LD on chromosomes with derived alleles primarily in Eurasian populations. To reconcile these findings, which appear to support different models of selection, we argue that selective sweeps (particularly, soft sweeps) of multiple derived alleles in and/or near LTα occurred in non-Africans after their ancestors left Africa. Furthermore, these targets of selection were predicted to alter transcription factor binding site affinity and protein stability, suggesting they play a role in gene function. Additionally, our data also showed that a subset of these functional adaptive variants are present in archaic hominin genomes. CONCLUSIONS: Overall, this study identified candidate functional alleles in a biologically-relevant genomic region, and offers new insights into the evolutionary origins of these loci in modern human populations.

Detection of Regional Variation in Selection Intensity within Protein-Coding Genes Using DNA Sequence Polymorphism and Divergence.

Numerous approaches have been developed to infer natural selection based on the comparison of polymorphism within species and divergence between species. These methods are especially powerful for the detection of uniform selection operating across a gene. However, empirical analyses have demonstrated that regions of protein-coding genes exhibiting clusters of amino acid substitutions are subject to different levels of selection relative to other regions of the same gene. To quantify this heterogeneity of selection within coding sequences, we developed Model Averaged Site Selection via Poisson Random Field (MASS-PRF). MASS-PRF identifies an ensemble of intragenic clustering models for polymorphic and divergent sites. This ensemble of models is used within the Poisson Random Field framework to estimate selection intensity on a site-by-site basis. Using simulations, we demonstrate that MASS-PRF has high power to detect clusters of amino acid variants in small genic regions, can reliably estimate the probability of a variant occurring at each nucleotide site in sequence data and is robust to historical demographic trends and recombination. We applied MASS-PRF to human gene polymorphism derived from the 1,000 Genomes Project and divergence data from the common chimpanzee. On the basis of this analysis, we discovered striking regional variation in selection intensity, indicative of positive or negative selection, in well-defined domains of genes that have previously been associated with neurological processing, immunity, and reproduction. We suggest that amino acid-altering substitutions within these regions likely are or have been selectively advantageous in the human lineage, playing important roles in protein function.

Genetic origins of lactase persistence and the spread of pastoralism in Africa.

In humans, the ability to digest lactose, the sugar in milk, declines after weaning because of decreasing levels of the enzyme lactase-phlorizin hydrolase, encoded by LCT. However, some individuals maintain high enzyme amounts and are able to digest lactose into adulthood (i.e., they have the lactase-persistence [LP] trait). It is thought that selection has played a major role in maintaining this genetically determined phenotypic trait in different human populations that practice pastoralism. To identify variants associated with the LP trait and to study its evolutionary history in Africa, we sequenced MCM6 introns 9 and 13 and ~2 kb of the LCT promoter region in 819 individuals from 63 African populations and in 154 non-Africans from nine populations. We also genotyped four microsatellites in an ~198 kb region in a subset of 252 individuals to reconstruct the origin and spread of LP-associated variants in Africa. Additionally, we examined the association between LP and genetic variability at candidate regulatory regions in 513 individuals from eastern Africa. Our analyses confirmed the association between the LP trait and three common variants in intron 13 (C-14010, G-13907, and G-13915). Furthermore, we identified two additional LP-associated SNPs in intron 13 and the promoter region (G-12962 and T-956, respectively). Using neutrality tests based on the allele frequency spectrum and long-range linkage disequilibrium, we detected strong signatures of recent positive selection in eastern African populations and the Fulani from central Africa. In addition, haplotype analysis supported an eastern African origin of the C-14010 LP-associated mutation in southern Africa.

Origin and differential selection of allelic variation at TAS2R16 associated with salicin bitter taste sensitivity in Africa.

Bitter taste perception influences human nutrition and health, and the genetic variation underlying this trait may play a role in disease susceptibility. To better understand the genetic architecture and patterns of phenotypic variability of bitter taste perception, we sequenced a 996 bp region, encompassing the coding exon of TAS2R16, a bitter taste receptor gene, in 595 individuals from 74 African populations and in 94 non-Africans from 11 populations. We also performed genotype-phenotype association analyses of threshold levels of sensitivity to salicin, a bitter anti-inflammatory compound, in 296 individuals from Central and East Africa. In addition, we characterized TAS2R16 mutants in vitro to investigate the effects of polymorphic loci identified at this locus on receptor function. Here, we report striking signatures of positive selection, including significant Fay and Wu's H statistics predominantly in East Africa, indicating strong local adaptation and greater genetic structure among African populations than expected under neutrality. Furthermore, we observed a "star-like" phylogeny for haplotypes with the derived allele at polymorphic site 516 associated with increased bitter taste perception that is consistent with a model of selection for "high-sensitivity" variation. In contrast, haplotypes carrying the "low-sensitivity" ancestral allele at site 516 showed evidence of strong purifying selection. We also demonstrated, for the first time, the functional effect of nonsynonymous variation at site 516 on salicin phenotypic variance in vivo in diverse Africans and showed that most other nonsynonymous substitutions have weak or no effect on cell surface expression in vitro, suggesting that one main polymorphism at TAS2R16 influences salicin recognition. Additionally, we detected geographic differences in levels of bitter taste perception in Africa not previously reported and infer an East African origin for high salicin sensitivity in human populations.

Limited evidence for adaptive evolution and functional effect of allelic variation at rs702424 in the promoter of the TAS2R16 bitter taste receptor gene in Africa.

Bitter taste perception, mediated by receptors encoded by the TAS2R loci, has important roles in human health and nutrition. Prior studies have demonstrated that nonsynonymous variation at site 516 in the coding exon of TAS2R16, a bitter taste receptor gene on chromosome 7, has been subject to positive selection and is strongly correlated with differences in sensitivity to salicin, a bitter anti-inflammatory compound, in human populations. However, a recent study suggested that the derived G-allele at rs702424 in the TAS2R16 promoter has also been the target of recent selection and may have an additional effect on the levels of salicin bitter taste perception. Here, we examined alleles at rs702424 for signatures of selection using Extended Haplotype Homozygosity (EHH) and FST statistics in diverse populations from West Central, Central and East Africa. We also performed a genotype-phenotype analysis of salicin sensitivity in a subset of 135 individuals from East Africa. Based on our data, we did not find evidence for positive selection at rs702424 in African populations, suggesting that nucleotide position 516 is likely the site under selection at TAS2R16. Moreover, we did not detect a significant association between rs702424 alleles and salicin bitter taste recognition, implying that this site does not contribute to salicin phenotypic variance. Overall, this study of African diversity provides further information regarding the genetic architecture and evolutionary history of a biologically-relevant trait in humans.

Evolution of functionally diverse alleles associated with PTC bitter taste sensitivity in Africa.

Although human bitter taste perception is hypothesized to be a dietary adaptation, little is known about genetic signatures of selection and patterns of bitter taste perception variability in ethnically diverse populations with different diets, particularly from Africa. To better understand the genetic basis and evolutionary history of bitter taste sensitivity, we sequenced a 2,975 bp region encompassing TAS2R38, a bitter taste receptor gene, in 611 Africans from 57 populations in West Central and East Africa with diverse subsistence patterns, as well as in a comparative sample of 132 non-Africans. We also examined the association between genetic variability at this locus and threshold levels of phenylthiocarbamide (PTC) bitterness in 463 Africans from the above populations to determine how variation influences bitter taste perception. Here, we report striking patterns of variation at TAS2R38, including a significant excess of novel rare nonsynonymous polymorphisms that recently arose only in Africa, high frequencies of haplotypes in Africa associated with intermediate bitter taste sensitivity, a remarkably similar frequency of common haplotypes across genetically and culturally distinct Africans, and an ancient coalescence time of common variation in global populations. Additionally, several of the rare nonsynonymous substitutions significantly modified levels of PTC bitter taste sensitivity in diverse Africans. While ancient balancing selection likely maintained common haplotype variation across global populations, we suggest that recent selection pressures may have also resulted in the unusually high level of rare nonsynonymous variants in Africa, implying a complex model of selection at the TAS2R38 locus in African populations. Furthermore, the distribution of common haplotypes in Africa is not correlated with diet, raising the possibility that common variation may be under selection due to their role in nondietary biological processes. In addition, our data indicate that novel rare mutations contribute to the phenotypic variance of PTC sensitivity, illustrating the influence of rare variation on a common trait, as well as the relatively recent evolution of functionally diverse alleles at this locus.

Addressing Africa's pandemic puzzle: Perspectives on COVID-19 transmission and mortality in sub-Saharan Africa.

To date, SARS-CoV-2 (the virus that causes COVID-19) has spread to almost every region of the world, infecting millions and resulting in the deaths of hundreds of thousands of people. Although it was predicted that Africa would suffer a massive loss of life due to this pandemic, the number of COVID-19 cases has been relatively low across the continent. Researchers have speculated that several factors may be responsible for this outcome in Africa, including the extensive experience that countries have with infectious diseases and the young median age of their populations. However, it is still important for African countries to adopt aggressive and bold approaches against COVID-19, in case the nature of the pandemic changes. This short review will summarize the status of the outbreak in Africa and propose possible reasons for current trends, as well as discuss interventions aimed at preventing a rapid increase in the number of COVID-19 cases in the future.

The epidemiology, diagnosis and management of scrub typhus disease in China.

Thirty-nine years ago, scrub typhus (ST), a disease, was not among the China's notifiable diseases. However, ST has reemerged to become a growing public health issue in the southwest part of China. The major factors contributing to an increased incidence and prevalence of this disease include rapid globalization, urbanization, expansion of humans into previously uninhabited areas, and climate change. The clinical manifestation of ST also consists of high fever, headache, weakness, myalgia, rash, and an eschar. In severe cases, complications (e.g. multi-organ failure, jaundice, acute renal failure, pneumonitis, myocarditis, and even death) can occur. The diagnosis of ST is mainly based on serological identification by indirect immunofluorescence assay and other molecular methods. Furthermore, several groups of antibiotics (e.g. tetracycline, chloramphenicol, macrolides, and rifampicin) are currently effective in treating this disease. This fact suggests the need for robust early diagnostic techniques, increased surveillance, and prompt treatment, and develop future vaccine.

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

We investigated global patterns of genetic variation and signatures of natural selection at host genes relevant to SARS-CoV-2 infection (ACE2, TMPRSS2, DPP4, and LY6E). We analyzed novel data from 2,012 ethnically diverse Africans and 15,997 individuals of European and African ancestry with electronic health records, and integrated with global data from the 1000GP. At ACE2, we identified 41 non-synonymous variants that were rare in most populations, several of which impact protein function. However, three non-synonymous variants were common among Central African hunter-gatherers from Cameroon and are on haplotypes that exhibit signatures of positive selection. We identify strong 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. Genetic variants that are targets of natural selection are associated with clinical phenotypes common in patients with COVID-19.

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.

Loci associated with skin pigmentation identified in African populations.

Despite the wide range of skin pigmentation in humans, little is known about its genetic basis in global populations. Examining ethnically diverse African genomes, we identify variants in or near SLC24A5, MFSD12, DDB1, TMEM138, OCA2, and HERC2 that are significantly associated with skin pigmentation. Genetic evidence indicates that the light pigmentation variant at SLC24A5 was introduced into East Africa by gene flow from non-Africans. At all other loci, variants associated with dark pigmentation in Africans are identical by descent in South Asian and Australo-Melanesian populations. Functional analyses indicate that MFSD12 encodes a lysosomal protein that affects melanogenesis in zebrafish and mice, and that mutations in melanocyte-specific regulatory regions near DDB1/TMEM138 correlate with expression of ultraviolet response genes under selection in Eurasians.

Global population-specific variation in miRNA associated with cancer risk and clinical biomarkers.

BACKGROUND: MiRNA expression profiling is being actively investigated as a clinical biomarker and diagnostic tool to detect multiple cancer types and stages as well as other complex diseases. Initial investigations, however, have not comprehensively taken into account genetic variability affecting miRNA expression and/or function in populations of different ethnic backgrounds. Therefore, more complete surveys of miRNA genetic variability are needed to assess global patterns of miRNA variation within and between diverse human populations and their effect on clinically relevant miRNA genes. METHODS: Genetic variation in 1524 miRNA genes was examined using whole genome sequencing (60x coverage) in a panel of 69 unrelated individuals from 14 global populations, including European, Asian and African populations. RESULTS: We identified 33 previously undescribed miRNA variants, and 31 miRNA containing variants that are globally population-differentiated in frequency between African and non-African populations (PD-miRNA). The top 1% of PD-miRNA were significantly enriched for regulation of genes involved in glucose/insulin metabolism and cell division (p < 10(-7)), most significantly the mitosis pathway, which is strongly linked to cancer onset. Overall, we identify 7 PD-miRNAs that are currently implicated as cancer biomarkers or diagnostics: hsa-mir-202, hsa-mir-423, hsa-mir-196a-2, hsa-mir-520h, hsa-mir-647, hsa-mir-943, and hsa-mir-1908. Notably, hsa-mir-202, a potential breast cancer biomarker, was found to show significantly high allele frequency differentiation at SNP rs12355840, which is known to affect miRNA expression levels in vivo and subsequently breast cancer mortality. CONCLUSION: MiRNA expression profiles represent a promising new category of disease biomarkers. However, population specific genetic variation can affect the prevalence and baseline expression of these miRNAs in diverse populations. Consequently, miRNA genetic and expression level variation among ethnic groups may be contributing in part to health disparities observed in multiple forms of cancer, specifically breast cancer, and will be an essential consideration when assessing the utility of miRNA biomarkers for the clinic.