Human BRCA pathogenic variants were originated during recent human history

Pathogenic variants in human DNA damage repair genes mostly arose after the latest human out-of-Africa migration

Introduction

The DNA damage repair (DDR) system in human genome is pivotal in maintaining genomic integrity. Pathogenic variation (PV) in DDR genes impairs their function, leading to genome instability and increased susceptibility to diseases, especially cancer. Understanding the evolution origin and arising time of DDR PV is crucial for comprehending disease susceptibility in modern humans.

Methods

We used big data approach to identify the PVs in DDR genes in modern humans. We mined multiple genomic databases derived from 251,214 modern humans of African and non-Africans. We compared the DDR PVs between African and non-African. We also mined the DDR PVs in the genomic data derived from 5,031 ancient humans. We used the DDR PVs from ancient humans as the intermediate to further the DDR PVs between African and non-African.

Results and discussion

We identified 1,060 single-base DDR PVs across 77 DDR genes in modern humans of African and non-African. Direct comparison of the DDR PVs between African and non-African showed that 82.1% of the non-African PVs were not present in African. We further identified 397 single-base DDR PVs in 56 DDR genes in the 5,031 ancient humans dated between 45,045 and 100 years before present (BP) lived in Eurasian continent therefore the descendants of the latest out-of-Africa human migrants occurred 50,000-60,000 years ago. By referring to the ancient DDR PVs, we observed that 276 of the 397 (70.3%) ancient DDR PVs were exclusive in non-African, 106 (26.7%) were shared between non-African and African, and only 15 (3.8%) were exclusive in African. We further validated the distribution pattern by testing the PVs in BRCA and TP53, two of the important genes in genome stability maintenance, in African, non-African, and Ancient humans. Our study revealed that DDR PVs in modern humans mostly emerged after the latest out-of-Africa migration. The data provides a foundation to understand the evolutionary basis of disease susceptibility, in particular cancer, in modern humans.

Using Portuguese BRCA pathogenic variation as a model to study the impact of human admixture on human health

BACKGROUND

Admixture occurs between different ethnic human populations. The global colonization in recent centuries by Europeans led to the most significant admixture in human history. While admixture may enhance genetic diversity for better fitness, it may also impact on human health by transmitting genetic variants for disease susceptibility in the admixture population. The admixture by Portuguese global exploration initiated in the 15th century has reached over 20 million of Portuguese-heritage population worldwide. It provides a valuable model to study the impact of admixture on human health. BRCA1 and BRCA2 (BRCA) are two of the important tumor suppressor genes. The pathogenic variation (PV) in BRCA is well determined to cause high risk of hereditary breast and ovarian cancer. Tracing the distribution of Portuguese BRCA PV in Portuguese-heritage population will help to understand the impact of admixture on cancer susceptibility in modern humans. In this study, we analyzed the distribution of the Portuguese-originated BRCA variation in Brazilian population, which has high degree Portuguese-heritage.

METHODS

By comprehensive data mining, standardization and annotation, we generated a Portuguese-derived BRCA variation dataset and a Brazilian-derived BRCA variation dataset. We compared the two BRCA variation datasets to identify the BRCA variants shared between the two populations.

RESULTS

The Portuguese-derived BRCA variation dataset consists of 220 BRCA variants including 78 PVs from 11,482 Portuguese cancer patients, 93 (42.2%) in BRCA1 and 127 (57.7%) in BRCA2. Of the 556 Portuguese BRCA PV carriers carrying the 78 PVs, 331 (59.5%) carried the three Portuguese-BRCA founder PVs of BRCA1 c.2037delinsCC, BRCA1 c.3331_3334del and BRCA2 c.156_157insAlu. The Brazilian-derived BRCA variation dataset consists of 255 BRCA PVs from 7,711 cancer patients, 136 (53.3%) in BRCA1 and 119 (46.6%) in BRCA2. We developed an open database named dbBRCA-Portuguese ( https://genemutation.fhs.um.edu.mo/dbbrca-portuguese/ ) and an open database named dbBRCA-Brazilian ( https://genemutation.fhs.um.edu.mo/dbbrca-brazilian ) to host the BRCA variation data from Portuguese and Brazilian populations. We compared the BRCA PV datasets between Portuguese and Brazilian populations, and identified 29 Portuguese-specific BRCA PVs shared between Portuguese and Brazilian populations, 14 in BRCA1 including the Portuguese founder BRCA1 c.3331_3334del and BRCA1 c.2037delinsCC, and 15 in BRCA2 including the Portuguese founder BRCA2 c.156_157insAlu. Searching the 78 Portuguese BRCA PVs in over 5,000 ancient human genomes identified evolution origin for only 8 PVs in Europeans dated between 37,470 and 3,818 years before present, confirming the Portuguese-specificity of Portuguese BRCA PVs; comparing the 78 Portuguese BRCA PVs Portuguese, 255 Brazilian BRCA PVs, and 134 African BRCA PVs showed little overlapping, ruling out the possibility that the BRCA PVs shared between Portuguese and Brazilian may also be contributed by African.

CONCLUSION

Our study provides evidence that the admixture in recent human history contributed to cancer susceptibility in modern humans.

Pathogenic variants in human DNA damage repair genes mostly arose in recent human history

BACKGROUND

Genome stability is maintained by the DNA damage repair (DDR) system composed of multiple DNA repair pathways of hundreds of genes. Germline pathogenic variation (PV) in DDR genes damages function of the affected DDR genes, leading to genome instability and high risk of diseases, in particular, cancer. Knowing evolutionary origin of the PVs in human DDR genes is essential to understand the etiology of human diseases. However, answer to the issue remains largely elusive. In this study, we analyzed evolutionary origin for the PVs in human DDR genes.

METHODS

We identified 169 DDR genes by referring to various databases and identified PVs in the DDR genes of modern humans from ClinVar database. We performed a phylogenetic analysis to analyze the conservation of human DDR PVs in 100 vertebrates through cross-species genomic data comparison using the phyloFit program of the PHAST package and visualized the results using the GraphPad Prism software and the ggplot module. We identified DDR PVs from over 5000 ancient humans developed a database to host the DDR PVs ( https://genemutation.fhs.um.edu.mo/dbDDR-AncientHumans ). Using the PV data, we performed a molecular archeological analysis to compare the DDR PVs between modern humans and ancient humans. We analyzed evolution selection of DDR genes across 20 vertebrates using the CodeML in PAML for phylogenetic analysis.

RESULTS

Our phylogenic analysis ruled out cross-species conservation as the origin of human DDR PVs. Our archeological approach identified rich DDR PVs shared between modern and ancient humans, which were mostly dated within the last 5000 years. We also observed similar pattern of quantitative PV distribution between modern and ancient humans. We further detected a set of ATM, BRCA2 and CHEK2 PVs shared between human and Neanderthals.

CONCLUSIONS

Our study reveals that human DDR PVs mostly arose in recent human history. We propose that human high cancer risk caused by DDR PVs can be a by-product of human evolution.

Population-based BRCA germline mutation screening in the Han Chinese identifies individuals at risk of BRCA mutation-related cancer: experience from a clinical diagnostic center from greater Shanghai area

BACKGROUND

Deleterious BRCA1/2 (BRCA) mutation raises the risk for BRCA mutation-related malignancies, including breast, ovarian, prostate, and pancreatic cancer. Germline variation of BRCA exhibits substantial ethnical diversity. However, there is limited research on the Chinese Han population, constraining the development of strategies for BRCA mutation screening in this large ethnic group.

METHODS

We profile the BRCA mutational spectrum, including single nucleotide variation, insertion/deletion, and large genomic rearrangements in 2,080 apparently healthy Chinese Han individuals and 522 patients with BRCA mutation-related cancer, to determine the BRCA genetic background of the Chinese Han population, especially of the East Han. Incident cancer events were monitored in 1,005 participants from the healthy group, comprising 11 BRCA pathogenic/likely pathogenic (PLP) variant carriers and 994 PLP-free individuals, including 3 LGR carriers.

RESULTS

Healthy Chinese Han individuals demonstrated a distinct BRCA mutational spectrum compared to cancer patients, with a 0.53% (1 in 189) prevalence of pathogenic/likely pathogenic (PLP) variant, alongside a 3 in 2,080 occurrence of LGR. BRCA1 c. 5470_5477del demonstrated high prevalence (0.44%) in the North Han Chinese and penetrance for breast cancer. None of the 3 LGR carriers developed cancer during the follow-up. We calculated a relative risk of 135.55 (95% CI 25.07 to 732.88) for the development of BRCA mutation-related cancers in the BRCA PLP variant carriers (mean age 42.91 years, median follow-up 10 months) compared to PLP-free individuals (mean age 48.47 years, median follow-up 16 months).

CONCLUSION

The unique BRCA mutational profile in the Chinese Han highlights the potential for standardized population-based BRCA variant screening to enhance BRCA mutation-related cancer prevention and treatment.

Onco-Breastomics: An Eco-Evo-Devo Holistic Approach

Known as a diverse collection of neoplastic diseases, breast cancer (BC) can be hyperbolically characterized as a dynamic pseudo-organ, a living organism able to build a complex, open, hierarchically organized, self-sustainable, and self-renewable tumor system, a population, a species, a local community, a biocenosis, or an evolving dynamical ecosystem (i.e., immune or metabolic ecosystem) that emphasizes both developmental continuity and spatio-temporal change. Moreover, a cancer cell community, also known as an oncobiota, has been described as non-sexually reproducing species, as well as a migratory or invasive species that expresses intelligent behavior, or an endangered or parasite species that fights to survive, to optimize its features inside the host's ecosystem, or that is able to exploit or to disrupt its host circadian cycle for improving the own proliferation and spreading. BC tumorigenesis has also been compared with the early embryo and placenta development that may suggest new strategies for research and therapy. Furthermore, BC has also been characterized as an environmental disease or as an ecological disorder. Many mechanisms of cancer progression have been explained by principles of ecology, developmental biology, and evolutionary paradigms. Many authors have discussed ecological, developmental, and evolutionary strategies for more successful anti-cancer therapies, or for understanding the ecological, developmental, and evolutionary bases of BC exploitable vulnerabilities. Herein, we used the integrated framework of three well known ecological theories: the Bronfenbrenner's theory of human development, the Vannote's River Continuum Concept (RCC), and the Ecological Evolutionary Developmental Biology (Eco-Evo-Devo) theory, to explain and understand several eco-evo-devo-based principles that govern BC progression. Multi-omics fields, taken together as onco-breastomics, offer better opportunities to integrate, analyze, and interpret large amounts of complex heterogeneous data, such as various and big-omics data obtained by multiple investigative modalities, for understanding the eco-evo-devo-based principles that drive BC progression and treatment. These integrative eco-evo-devo theories can help clinicians better diagnose and treat BC, for example, by using non-invasive biomarkers in liquid-biopsies that have emerged from integrated omics-based data that accurately reflect the biomolecular landscape of the primary tumor in order to avoid mutilating preventive surgery, like bilateral mastectomy. From the perspective of preventive, personalized, and participatory medicine, these hypotheses may help patients to think about this disease as a process governed by natural rules, to understand the possible causes of the disease, and to gain control on their own health.

Oncogenic BRCA1,2 Mutations in the Human Lineage-A By-Product of Sexual Selection?

In this review, we discuss the long-known problem of tissue-specific carcinogenesis in BRCA1 and BRCA2 mutation carriers: while the genes are expressed ubiquitously, increased cancer risk is observed mostly in the breast and ovaries, and to a much lesser extent, in some other tissues such as the prostate or pancreas. We reevaluate hypotheses on the evolutionary origin of these mutations in humans. Also, we align together the reports that at least some great apes have much lower risks of epithelial cancers in general and breast cancer in particular with the fact that humans have more voluminous breast tissue as compared to their closest extant relatives, particularly chimpanzees and bonobos. We conjecture that this disparity may be a consequence of sexual selection, augmented via selection for enhanced lactation. Further, we argue that there is an organ-specific enigma similar to the Peto paradox: breast cancer risk in humans is only minimally correlated with breast size. These considerations lead to the hypothesis that, along with the evolutionary development of larger breasts in humans, additional changes have played a balancing role in suppressing breast cancer. These yet-to-be-discovered mechanisms, while purely speculative, may be valuable to understanding human breast cancer, though they may not be exclusive to the mammary gland epithelial cells. Combining these themes, we review some anti-carcinogenesis preventive strategies and prospects of new interventions against breast cancer.

Classification of PTEN missense VUS through exascale simulations

Phosphatase and tensin homolog (PTEN), a tumor suppressor with dual phosphatase properties, is a key factor in PI3K/AKT signaling pathway. Pathogenic germline variation in PTEN can abrogate its ability to dephosphorylate, causing high cancer risk. Lack of functional evidence lets numerous PTEN variants be classified as variants of uncertain significance (VUS). Utilizing Molecular Dynamics (MD) simulations, we performed a thorough evaluation for 147 PTEN missense VUS, sorting them into 66 deleterious and 81 tolerated variants. Utilizing replica exchange molecular dynamic (REMD) simulations, we further assessed the variants situated in the catalytic core of PTEN's phosphatase domain and uncovered conformational alterations influencing the structural stability of the phosphatase domain. There was a high degree of agreement between our results and the variants classified by Variant Abundance by Massively Parallel Sequencing, saturation mutagenesis, multiplexed functional data and experimental assays. Our extensive analysis of PTEN missense VUS should benefit their clinical applications in PTEN-related cancer.

SIGNIFICANCE STATEMENT

Classification of PTEN variants affecting its lipid phosphatase activity is important for understanding the roles of PTEN variation in the pathogenesis of hereditary and sporadic malignancies. Of the 3000 variants identified in PTEN, 1296 (43%) were assigned as VUS. Here, we applied MD and REMD simulations to investigate the effects of PTEN missense VUS on the structural integrity of the PTEN phosphatase domain consisting the WPD, P and TI active sites. We classified a total of 147 missense VUS into 66 deleterious and 81 tolerated variants by referring to the control group comprising 54 pathogenic and 12 benign variants. The classification was largely in concordance with these classified by experimental approaches.

TP53 germline pathogenic variants in modern humans were likely originated during recent human history

TP53 is crucial for maintaining genome stability and preventing oncogenesis. Germline pathogenic variation in TP53 damages its function, causing genome instability and increased cancer risk. Despite extensive study in TP53, the evolutionary origin of the human TP53 germline pathogenic variants remains largely unclear. In this study, we applied phylogenetic and archaeological approaches to identify the evolutionary origin of TP53 germline pathogenic variants in modern humans. In the phylogenic analysis, we searched 406 human TP53 germline pathogenic variants in 99 vertebrates distributed in eight clades of Primate, Euarchontoglires, Laurasiatheria, Afrotheria, Mammal, Aves, Sarcopterygii and Fish, but we observed no direct evidence for the cross-species conservation as the origin; in the archaeological analysis, we searched the variants in 5031 ancient human genomes dated between 45045 and 100 years before present, and identified 45 pathogenic variants in 62 ancient humans dated mostly within the last 8000 years; we also identified 6 pathogenic variants in 3 Neanderthals dated 44000 to 38515 years before present and 1 Denisovan dated 158 550 years before present. Our study reveals that TP53 germline pathogenic variants in modern humans were likely originated in recent human history and partially inherited from the extinct Neanderthals and Denisovans.

Evolutionary Origin of Human PALB2 Germline Pathogenic Variants

PALB2 (Partner and localizer of BRCA2) is crucial for repairing DNA double-stranded breaks (DSBs) through homologous recombination (HR). Germline pathogenic variation in PALB2 disrupts DNA damage repair and increases the risk of Fanconi Anemia, breast cancer, and ovarian cancer. Determination of the evolutionary origin of human PALB2 variants will promote a deeper understanding of the biological basis of PALB2 germline variation and its roles in human diseases. We tested the evolution origin for 1444 human PALB2 germline variants, including 484 pathogenic and 960 benign variants. We performed a phylogenic analysis by tracing the variants in 100 vertebrates. However, we found no evidence to show that cross-species conservation was the origin of PALB2 germline pathogenic variants, but it is indeed a rich source for PALB2 germline benign variants. We performed a paleoanthropological analysis by tracing the variants in over 5000 ancient humans. We identified 50 pathogenic in 71 ancient humans dated from 32,895 to 689 before the present, of which 90.1% were dated within the recent 10,000 years. PALB2 benign variants were also highly shared with ancient humans. Data from our study reveal that human PALB2 pathogenic variants mostly arose in recent human history.

Integration of deep learning with Ramachandran plot molecular dynamics simulation for genetic variant classification

Functional classification of genetic variants is a key for their clinical applications in patient care. However, abundant variant data generated by the next-generation DNA sequencing technologies limit the use of experimental methods for their classification. Here, we developed a protein structure and deep learning (DL)-based system for genetic variant classification, DL-RP-MDS, which comprises two principles: 1) Extracting protein structural and thermodynamics information using the Ramachandran plot-molecular dynamics simulation (RP-MDS) method, 2) combining those data with an unsupervised learning model of auto-encoder and a neural network classifier to identify the statistical significance patterns of the structural changes. We observed that DL-RP-MDS provided higher specificity than over 20 widely used in silico methods in classifying the variants of three DNA damage repair genes: TP53, MLH1, and MSH2. DL-RP-MDS offers a powerful platform for high-throughput genetic variant classification. The software and online application are available at https://genemutation.fhs.um.edu.mo/DL-RP-MDS/.

Ethnic-specificity, evolution origin and deleteriousness of Asian BRCA variation revealed by over 7500 BRCA variants derived from Asian population

Pathogenic variation in BRCA1 and BRCA2 (BRCA) causes high risk of breast and ovarian cancer, and BRCA variation data are important markers for BRCA-related clinical cancer applications. However, comprehensive BRCA variation data are lacking from the Asian population despite its large population size, heterogenous genetic background and diversified living environment across the Asia continent. We performed a systematic study on BRCA variation in Asian population including extensive data mining, standardization, annotation and characterization. We identified 7587 BRCA variants from 685 592 Asian individuals in 40 Asia countries and regions, including 1762 clinically actionable pathogenic variants and 4915 functionally unknown variants (https://genemutation.fhs.um.edu.mo/Asian-BRCA/). We observed the highly ethnic-specific nature of Asian BRCA variants between Asian and non-Asian populations and within Asian populations, highlighting that the current European descendant population-based BRCA data is inadequate to reflect BRCA variation in the Asian population. We also provided archeological evidence for the evolutionary origin and arising time of Asian BRCA variation. We further provided structural-based evidence for the deleterious variants enriched within the functionally unknown Asian BRCA variants. The data from our study provide a current view of BRCA variation in the Asian population and a rich resource to guide clinical applications of BRCA-related cancer for the Asian population.

A global perspective on the ethnic-specific BRCA variation and its implication in clinical application

Pathogenic BRCA1 and BRCA2 (BRCA) variation is the genetic predisposition for high cancer risk affecting mostly breast and ovarian. BRCA variation information is widely used in clinical diagnosis, treatment, and prevention of BRCA-related cancer. The positive selection imposed on human BRCA leads to highly ethnic-specific BRCA variation to adapt different living environment on earth. Most of the human BRCA variants identified so far were from the European descendant populations and used as the standard reference for global human populations, whereas BRCA variation in other ethnic populations remains poorly characterized. This review addresses the origin of ethnic-specific BRCA variation, the importance of ethnic-specific BRCA variation in clinical application, the limitation of current BRCA variation data, and potential solutions to fill the gap.

Evolutionary Origin of Germline Pathogenic MUTYH Variations in Modern Humans

MUTYH plays an essential role in preventing oxidation-caused DNA damage. Pathogenic germline variations in MUTYH damage its function, causing intestinal polyposis and colorectal cancer. Determination of the evolutionary origin of the variation is essential to understanding the etiological relationship between MUTYH variation and cancer development. In this study, we analyzed the origins of pathogenic germline variants in human MUTYH. Using a phylogenic approach, we searched pathogenic MUTYH variants in modern humans in the MUTYH of 99 vertebrates across eight clades. We did not find pathogenic variants shared between modern humans and the non-human vertebrates following the evolutionary tree, ruling out the possibility of cross-species conservation as the origin of human pathogenic variants in MUTYH. We then searched the variants in the MUTYH of 5031 ancient humans and extinct Neanderthals and Denisovans. We identified 24 pathogenic variants in 42 ancient humans dated between 30,570 and 480 years before present (BP), and three pathogenic variants in Neanderthals dated between 65,000 and 38,310 years BP. Data from our study revealed that human pathogenic MUTYH variants mostly arose in recent human history and were partially inherited from Neanderthals.

The Etiology of Parkinson's Disease: New Perspectives from Gene-Environment Interactions

Parkinson's disease is now the most rapidly growing neurodegenerative disease worldwide. It is therefore critical to identify which factors, and to what extent, contribute to the multifactorial etiology of Parkinson's disease. Here, we address two interesting elements from the perspective of genetics, namely (a) the estimated age of several genetic risk factors related to Parkinson's disease; and (b) the relative contribution of genetics to the etiology of Parkinson's disease, as derived from twin studies. Based on these two perspectives, we argue that most genetic risk factors are by themselves insufficient to explain the majority of Parkinson's disease, and that environmental factors are required for these genetic factors to become pathophysiologically relevant.

Heterozygotic Brca1 mutation initiates mouse genome instability at embryonic stage

BRCA1 mutation is the genetic predisposition in causing genome instability towards cancer. BRCA1 mutation is predominantly germline inherited at the fertilization. However, when the inherited mutation initiates genome instability in the mutation carriers remains largely elusive. We used a heterozygotic Brca1-knockout mouse as a model to investigate the issue. Through whole-genome sequencing and bioinformatics analysis, we monitored genome status across the developmental stages from embryo to adulthood in the mouse model. We observed that genome instability as reflected by structural variation, indel and copy number variation already appeared at 10.5-day embryo and progressively towards adulthood. We also observed that the genome instability was not linearly accumulated but dynamically changed along the developmental process, affecting many oncogenic genes and pathways including DNA damage repair, estrogen signaling, and oncogenesis. We further observed that many genome abnormalities in the cancer caused by Brca1 mutation were originated at embryonic stage, and Trp53 (TP53) mutation was not essential for the Brca1 mutation-caused genome instability in the non-cancer cells. Our study revealed that heterozygotic Brca1 mutation alone can cause genome instability at embryonic stage, highlighting that prevention of BRCA1 mutation-related cancer in humans may need to start earlier than currently considered.

Distinct landscapes of deleterious variants in DNA damage repair system in ethnic human populations

Deleterious variants in the DNA damage repair system can cause genome instability and increase cancer risk. The highly ethnic-specific DDR deleterious variation from this study suggests its potential relationship with different disease susceptibility in ethnic human populations.

Deleterious variants in DNA damage repair (DDR) system can cause genome instability and increase cancer risk. In this study, we analyzed the deleterious variants in DDR system in 16 ethnic human populations. From the genetic variants in 169 DDR genes involved in nine DDR pathways collected from 158,612 individuals of different ethnic background, we identified 1,781 deleterious variants in 81 DDR genes in eight DDR pathways (https://genemutation.fhs.um.edu.mo/dbddr-global/). Our analysis showed although the quantity of deleterious variants was loaded at a similar level, the landscape of the variants differed substantially among different populations that two-third of the variants were present in single ethnic populations, and the rest was mostly shared between the populations with closer geographic and genetic relationship. The highly ethnic-specific DDR deleterious variation suggests its potential relationship with different disease susceptibility in ethnic human populations.