Dps-dependent in vivo mutation enhances long-term host adaptation in Vibrio cholerae.

As one of the most successful pathogenic organisms, Vibrio cholerae (V. cholerae) has evolved sophisticated regulatory mechanisms to overcome host stress. During long-term colonization by V. cholerae in adult mice, many spontaneous nonmotile mutants (approximately 10% at the fifth day post-infection) were identified. These mutations occurred primarily in conserved regions of the flagellar regulator genes flrA, flrC, and rpoN, as shown by Sanger and next-generation sequencing, and significantly increased fitness during colonization in adult mice. Intriguingly, instead of key genes in DNA repair systems (mutS, nfo, xthA, uvrA) or ROS and RNS scavenging systems (katG, prxA, hmpA), which were generally thought to be associated with bacterial mutagenesis, we found that deletion of the cyclin gene dps significantly increased the mutation rate (up to 53% at the fifth day post-infection) in V. cholerae. We further determined that the dpsD65A and dpsF46E point mutants showed a similar mutagenesis profile as the Δdps mutant during long-term colonization in mice, which strongly indicated that the antioxidative function of Dps directly contributes to the development of V. cholerae nonmotile mutants. Methionine metabolism pathway may be one of the mechanism for ΔflrA, ΔflrC and ΔrpoN mutant increased colonization in adult mice. Our results revealed a new phenotype in which increased fitness of V. cholerae in the host gut via spontaneous production nonmotile mutants regulated by cyclin Dps, which may represent a novel adaptation strategy for directed evolution of pathogens in the host.

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.

Classification of MLH1 Missense VUS Using Protein Structure-Based Deep Learning-Ramachandran Plot-Molecular Dynamics Simulations Method.

Pathogenic variation in DNA mismatch repair (MMR) gene MLH1 is associated with Lynch syndrome (LS), an autosomal dominant hereditary cancer. Of the 3798 MLH1 germline variants collected in the ClinVar database, 38.7% (1469) were missense variants, of which 81.6% (1199) were classified as Variants of Uncertain Significance (VUS) due to the lack of functional evidence. Further determination of the impact of VUS on MLH1 function is important for the VUS carriers to take preventive action. We recently developed a protein structure-based method named "Deep Learning-Ramachandran Plot-Molecular Dynamics Simulation (DL-RP-MDS)" to evaluate the deleteriousness of MLH1 missense VUS. The method extracts protein structural information by using the Ramachandran plot-molecular dynamics simulation (RP-MDS) method, then combines the variation data with an unsupervised learning model composed of auto-encoder and neural network classifier to identify the variants causing significant change in protein structure. In this report, we applied the method to classify 447 MLH1 missense VUS. We predicted 126/447 (28.2%) MLH1 missense VUS were deleterious. Our study demonstrates that DL-RP-MDS is able to classify the missense VUS based solely on their impact on protein structure.

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.

CBS-derived H2S facilitates host colonization of Vibrio cholerae by promoting the iron-dependent catalase activity of KatB.

Sensing and resisting oxidative stress is critical for Vibrio cholerae to survive in either the aquatic environment or the gastrointestinal tract. Previous studies mainly focused on the mechanisms of oxidative stress response regulation that rely on enzymatic antioxidant systems, while functions of non-enzymatic antioxidants are rarely discussed in V. cholerae. For the first time, we investigated the role of hydrogen sulfide (H2S), the simplest thiol compound, in protecting V. cholerae against oxidative stress. We found that degradation of L-cysteine by putative cystathionine ß-synthase (CBS) is the major source of endogenous H2S in V. cholerae. Our results indicate that intracellular H2S level has a positive correlation with cbs expression, while the enhanced H2S production can render V. cholerae cells less susceptible to H2O2 in vitro. Using proteome analysis and real-time qPCR assay, we found that cbs expression could stimulate the expression of several enzymatic antioxidants, including reactive oxygen species (ROS) detoxifying enzymes SodB, KatG and AhpC, the DNA protective protein DPS and the protein redox regulator Trx1. Assays of ROS detoxification capacities revealed that CBS-derived H2S could promote catalase activity at the post-translational level, especially for KatB, which serves as an important way that endogenous H2S participates in H2O2 detoxification. The enhancement of catalase activity by H2S is achieved through facilitating the uptake of iron. Adult mice experiments showed that cbs mutant has colonization defect, while either complementation of cbs or exogenous supplement of N-Acetyl-L-Cysteine restores its fitness in the host environment. Herein, we proposed that V. cholerae regulates CBS-dependent H2S production for better survival and proliferation under ROS stress.

Prevalence and spectrum of DNA mismatch repair gene variation in the general Chinese population.

BACKGROUND: Identifying genetic disease-susceptible individuals through population screening is considered as a promising approach for disease prevention. DNA mismatch repair (MMR) genes including MLH1, MSH2, MSH6 and PMS2 play essential roles in maintaining microsatellite stability through DNA mismatch repair, and pathogenic variation in MMR genes causes microsatellite instability and is the genetic predisposition for cancer as represented by the Lynch syndrome. While the prevalence and spectrum of MMR variation has been extensively studied in cancer, it remains largely elusive in the general population. Lack of the knowledge prevents effective prevention for MMR variation-caused cancer. In the current study, we addressed the issue by using the Chinese population as a model. METHODS: We performed extensive data mining to collect MMR variant data from 18 844 ethnic Chinese individuals and comprehensive analyses for the collected MMR variants to determine its prevalence, spectrum and features of the MMR data in the Chinese population. RESULTS: We identified 17 687 distinct MMR variants. We observed substantial differences of MMR variation between the general Chinese population and Chinese patients with cancer, identified highly Chinese-specific MMR variation through comparing MMR data between Chinese and non-Chinese populations, predicted the enrichment of deleterious variants in the unclassified Chinese-specific MMR variants, determined MMR pathogenic prevalence of 0.18% in the general Chinese population and determined that MMR variation in the general Chinese population is evolutionarily neutral. CONCLUSION: Our study provides a comprehensive view of MMR variation in the general Chinese population, a resource for biological study of human MMR variation, and a reference for MMR-related cancer applications.

Ethnic-specific BRCA1/2 variation within Asia population: evidence from over 78 000 cancer and 40 000 non-cancer cases of Indian, Chinese, Korean and Japanese populations.

BACKGROUND: Germline mutation in BRCA1 and BRCA2 (BRCA) is genetic predisposition for breast and ovarian cancer. Identification of mutation carriers is a critical step to prevent and treat the cancer in the mutation carriers. Human BRCA variation has been well determined as ethnic-specific by studies in Ashkenazi Jewish, Polish and Icelandic populations in the 1990s. However, sufficient evidence is lacking to determine if ethnic-specific BRCA variation is also present in Asia population, which is the largest and the most diversified in modern humans. Our current study aims to investigate ethnic-specific BRCA variation in Asian population. METHODS: We performed a comprehensive data mining to collect BRCA variation data in Indian, Chinese, Korean and Japanese populations derived from over 78 000 cancer and 40 000 non-cancer cases. We standardised all BRCA variation data following the international standard. We made a systematic comparison between the datasets including variant composition, variation spectrum, variant type, clinical class, founder mutation and high-frequent variants. RESULTS: Our analysis showed that over half of the Asian BRCA variants were Asian-specific, and significant differences were present between the four Asia populations in each category analysed. CONCLUSION: Data from our study reveal that ethnic-specific BRCA variation is commonly present in Asia population as existing in non-Asian populations. Our study indicates that ethnicity should be an important factor to consider in prevention and treatment of BRCA mutation-related cancer in the Asia population. We recommend that the current BRCA variation databases should include ethnic variation information in order to function as true global BRCA references.

Can population BRCA screening be applied in non-Ashkenazi Jewish populations? Experience in Macau population.

BACKGROUND: Pathogenic mutation in BRCA genes causes high cancer risk. Identifying the mutation carriers plays key roles in preventing BRCA mutation-related cancer. Population screening has demonstrated its power for comprehensive identification of the mutation carriers. However, it is only recommended for the Ashkenazi Jewish population with high prevalence of three founder mutations, but not for non-Ashkenazi Jewish populations as the cost-effectiveness could be too low due to their lower mutation prevalence and lack of founder mutation. Population screening would not benefit the majority of the human population for BRCA mutation-related cancer prevention. METHODS: We used population BRCA screening in 6000 residents, 1% of the Macau population, an ethnic Chinese population with unique genetic, linguistic and cultural features, and its BRCA mutation has not been analysed before. RESULTS: We called BRCA variants, identified 18 carriers with 14 pathogenic mutations and determined the prevalence of 0.29% in the population (95% CI 0.15% to 0.42%). We compared the testing cost between the Ashkenazi Jewish population, the Sephardi Jewish population and the Macau population, and observed only a few fold differences. CONCLUSION: Our study shows that testing cost is not the most important factor in considering population BRCA screening, at least for the populations in the developed countries/regions, regardless of the status of mutation prevalence and founder mutation.

Comprehensive Identification of Deleterious TP53 Missense VUS Variants Based on Their Impact on TP53 Structural Stability.

TP53 plays critical roles in maintaining genome stability. Deleterious genetic variants damage the function of TP53, causing genome instability and increased cancer risk. Of the large quantity of genetic variants identified in TP53, however, many remain functionally unclassified as variants of unknown significance (VUS) due to the lack of evidence. This is reflected by the presence of 749 (42%) VUS of the 1785 germline variants collected in the ClinVar database. In this study, we addressed the deleteriousness of TP53 missense VUS. Utilizing the protein structure-based Ramachandran Plot-Molecular Dynamics Simulation (RPMDS) method that we developed, we measured the effects of missense VUS on TP53 structural stability. Of the 340 missense VUS tested, we observed deleterious evidence for 193 VUS, as reflected by the TP53 structural changes caused by the VUS-substituted residues. We compared the results from RPMDS with those from other in silico methods and observed higher specificity of RPMDS in classification of TP53 missense VUS than these methods. Data from our current study address a long-standing challenge in classifying the missense VUS in TP53, one of the most important tumor suppressor genes.

Germline variation in BRCA1/2 is highly ethnic-specific: Evidence from over 30,000 Chinese hereditary breast and ovarian cancer patients.

BRCA1 and BRCA2 play essential roles in maintaining the genome stability. Pathogenic germline mutations in these two genes disrupt their function, lead to genome instability and increase the risk of developing breast and ovarian cancers. BRCA mutations have been extensively screened in Caucasian populations, and the resulting information are used globally as the standard reference in clinical diagnosis, treatment and prevention of BRCA-related cancers. Recent studies suggest that BRCA mutations can be ethnic-specific, raising the question whether a Caucasian-based BRCA mutation information can be used as a universal standard worldwide, or whether an ethnicity-based BRCA mutation information system need to be developed for the corresponding ethnic populations. In this study, we used Chinese population as a model to test ethnicity-specific BRCA mutations considering that China has one of the latest numbers of breast cancer patients therefore BRCA mutation carriers. Through comprehensive data mining, standardization and annotation, we collected 1,088 distinct BRCA variants derived from over 30,000 Chinese individuals, one of the largest BRCA data set from a non-Caucasian population covering nearly all known BRCA variants in the Chinese population (https://dbBRCA-Chinese.fhs.umac.mo). Using this data, we performed multi-layered analyses to determine the similarities and differences of BRCA variation between Chinese and non-Chinese ethnic populations. The results show the substantial differences of BRCA data between Chinese and non-Chinese ethnicities. Our study indicates that the current Caucasian population-based BRCA data is not adequate to represent the BRCA status in non-Caucasian populations. Therefore, ethnic-based BRCA standards need to be established to serve for the non-Caucasian populations.

[Kuntai capsule combined with gonadotropin releasing hormone agonist in treatment of moderate-severe endometriosis: a clinical observation].

OBJECTIVE: To observe the effect of Kuntai Capsule (KC), a Chinese patent medicine, in add-back therapy for gonadotropin-releasing hormone agonist (GnRH-a) treatment for moderate-severe endometriosis (EM). METHODS: Totally 100 patients suffering from stage III/IV EM, who were confirmed by laparoscopic surgery were randomly assigned to the GnRH-a group (A) and the KC combined GnRH-a group (B), 50 in each group. Patients in Group A were hypodermically injected with goserelin (3.6 mg), once per 4 weeks. Those in Group B additionally took KC, 4 pills each time, three times per day. The therapeutic course for all was 12 weeks. Serum levels of estradiol (E2), follicle stimulating hormone (FSH), bone gamma-carboxyglutamic-acid-containing proteins (BGP) were measured respectively. Kupperman Menopausal Index (KMI) and bone mineral density (BMD) of the lumbar vertebra were also compared between the two groups. RESULTS: Serum levels of E2 and FSH both significantly decreased in the two groups at week 12 of the treatment (P < 0.05), when compared with pre-treatment. Compared with before treatment in the same group, KMI increased in the two groups (P < 0.05). Compared with before treatment in the same group, BMI decreased in the two groups with no statistical difference (P > 0.05). Serum BGP increased after 12-week treatment (P < 0.05). Compared with Group A after treatment, serum levels of E2 and FSH both significantly increased in Group B (P < 0.05). There was no statistical difference in KMI between the two groups (P > 0.05). As for the incidence of menopausal symptoms, better effects in improving symptoms such as hot flashes, sleep disorders, and vaginal dryness were obtained in Group B than in Group A (P < 0.05). There was no significant difference in the post-pre-treatment difference of BMI between the two groups, but with statistical post-pre-treatment difference in the BGP level (P < 0.05). CONCLUSIONS: HKC combined GnRH-a could effectively reduce GnRH-a treatment induced partial low estrogen symptoms, improve increased serum BGP levels after GnRH-a therapy.

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/.

Thiol-Based Modification of MarR Protein VnrR Regulates Resistance Toward Nitrofuran in Vibrio cholerae By Promoting the Expression of a Novel Nitroreductase VnrA and of NO-Detoxifying Enzyme HmpA.

Aims: Epidemiological investigations have indicated low resistance toward nitrofuran in clinical isolates, suggesting its potential application in the treatment of multidrug-resistant bacteria. Therefore, it is valuable to explore the mechanism of bacterial resistance to nitrofuran. Results: Through phenotypic screening of ten multiple antibiotic resistance regulator (MarR) proteins in Vibrio cholerae, we discovered that the regulator VnrR (VCA1058) plays a crucial role in defending against nitrofuran, specifically furazolidone (FZ). Our findings demonstrate that VnrR responds to FZ metabolites, such as hydroxylamine, methylglyoxal, hydrogen peroxide (H2O2), ß-hydroxyethylhydrazine. Notably, VnrR exhibits reversible responses to the addition of H2O2 through three cysteine residues (Cys180, Cys223, Cys247), leading to the derepression of its upstream gene, vnrA (vca1057). Gene vnrA encodes a novel nitroreductase, which directly contributes to the degradation of FZ. Our study reveals that V. cholerae metabolizes FZ via the vnrR-vnrA system and achieves resistance to FZ with the assistance of the classical reactive oxygen/nitrogen species scavenging pathway. Innovation and Conclusion: This study represents a significant advancement in understanding the antibiotic resistance mechanisms of V. cholerae and other pathogens. Our findings demonstrate that the MarR family regulator, VnrR, responds to the FZ metabolite H2O2, facilitating the degradation and detoxification of this antibiotic in a thiol-dependent manner. These insights not only enrich our knowledge of antibiotic resistance but also provide new perspectives for the control and prevention of multidrug-resistant bacteria.

Magnetic chitin beads (MCB) coated with Vibrio cholerae reveals transcriptome dynamics in adult mice with a complex gut microbiota.

Vibrio cholerae adapts to the host environment by altering gene expression. Because of the complexity of the gut microbiome, current in vivo V. cholerae transcriptome studies have focused on microbiota-undeveloped conditions, neglecting the interaction between the host's commensal gut microbiota and V. cholerae. In this study, we analyzed the transcriptome of fully colonized adult mice in vivo using V. cholerae coated-magnetic chitin beads (vcMCB). This provides a simple yet powerful method for obtaining high-quality RNA from V. cholerae during colonization in mice. The transcriptome of V. cholerae recovered from adult mice infected with vcMCB shows differential expression of several genes when compared to V. cholerae recovered from the infant mouse and infant rabbit model. Some of these genes were also observed to be differentially expressed in previous studies of V. cholera recovered from human infection when compared to V. cholerae grown in vitro. In particular, we confirmed that V. cholerae resists the inhibitory effects of low pH and formic acid from gut microbiota, such as Anaerostipes caccae and Dorea formicigenerans, by downregulating vc1080. We propose that the vc1080 product may protect V. cholerae from formic acid stress through a novel acid tolerance response mechanism. Transcriptomic data obtained using the vcMCB system provide new perspectives on the interaction between V. cholerae and the gut microbiota, and this approach can also be applied to studies of other pathogenic bacteria.

Cross-Platform Transcriptomic Data Integration, Profiling, and Mining in Vibrio cholerae.

A large number of transcriptome studies generate important data and information for the study of pathogenic mechanisms of pathogens, including Vibrio cholerae. V. cholerae transcriptome data include RNA-seq and microarray: microarray data mainly include clinical human and environmental samples, and RNA-seq data mainly focus on laboratory processing conditions, including different stresses and experimental animals in vivo. In this study, we integrated the data sets of both platforms using Rank-in and the Limma R package normalized Between Arrays function, achieving the first cross-platform transcriptome data integration of V. cholerae. By integrating the entire transcriptome data, we obtained the profiles of the most active or silent genes. By transferring the integrated expression profiles into the weighted correlation network analysis (WGCNA) pipeline, we identified the important functional modules of V. cholerae in vitro stress treatment, gene manipulation, and in vitro culture as DNA transposon, chemotaxis and signaling, signal transduction, and secondary metabolic pathways, respectively. The analysis of functional module hub genes revealed the uniqueness of clinical human samples; however, under specific expression patterning, the Δhns, ΔoxyR1 strains, and tobramycin treatment group showed high expression profile similarity with human samples. By constructing a protein-protein interaction (PPI) interaction network, we discovered several unreported novel protein interactions within transposon functional modules. IMPORTANCE We used two techniques to integrate RNA-seq data for laboratory studies with clinical microarray data for the first time. The interactions between V. cholerae genes were obtained from a global perspective, as well as comparing the similarity between clinical human samples and the current experimental conditions, and uncovering the functional modules that play a major role under different conditions. We believe that this data integration can provide us with some insight and basis for elucidating the pathogenesis and clinical control of V. cholerae.

Distinct landscapes of deleterious variants in DNA damage repair system 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.

Identification of deleterious variants of uncertain significance in BRCA2 BRC4 repeat through molecular dynamics simulations.

Large quantity of variants of uncertain significance (VUS) has been identified in cancer predisposition genes, but classification of VUS remains a big challenge. We proposed that the impact of VUS on protein structure stability can be used to identify these with deleterious effects by using molecular dynamics simulation (MDS)-based approach and developed a MDS-based method for missense VUS classification. In the current study, we applied the system to classify the missense VUS in BRCA2. BRCA2 plays an important role in maintaining genome stability by repairing double-strand DNA damage through homologous recombination. BRCA2 BRC repeats bring RAD51 from cytoplasm to the break sites in nucleus to initiate the repairing process. Missense variants in BRCA2 BRC repeats can interfere the interaction between BRCA2 and RAD51, impair double-strand break repair, cause genome instability and increase cancer risk. We characterized the missense VUS in BRCA2 BRC4 repeat, the primary site of BRCA2 interacting with RAD51. Based on the well-determined BRC4 structure, we applied MDS to measure the impact of BRC4 missense VUS on the stability of BRC4 structure by testing the equilibrium state, flexibility, compactness, hydrogen bonds and surface accessibility. Of the 46 missense VUS analyzed, we were able to differentiate them into 24 Deleterious and 22 Tolerated variants. Comparison between the MDS-based and other 24 existing computational methods for variant classification showed that the MDS-based approach is highly sensitive and specific for classifying missense VUS in cancer predisposition genes.

Human BRCA pathogenic variants were originated during recent human history.

BRCA1 and BRCA2 (BRCA) play essential roles in maintaining genome stability. BRCA germline pathogenic variants increase cancer risk. However, the evolutionary origin of human BRCA pathogenic variants remains largely elusive. We tested the 2,972 human BRCA1 and 3,652 human BRCA2 pathogenic variants from ClinVar database in 100 vertebrates across eight clades, but failed to find evidence to show cross-species evolution conservation as the origin; we searched the variants in 2,792 ancient human genome data, and identified 28 BRCA1 and 22 BRCA2 pathogenic variants in 44 cases dated from 45,000 to 300 yr ago; we analyzed the haplotype-dated human BRCA pathogenic founder variants, and observed that they were mostly arisen within the past 3,000 yr; we traced ethnic distribution of human BRCA pathogenic variants, and found that the majority were present in single or a few ethnic populations. Based on the data, we propose that human BRCA pathogenic variants were highly likely arisen in recent human history after the latest out-of-Africa migration, and the expansion of modern human population could largely increase the variation spectrum.

Mr.Vc v2: An updated version of database with increased data of transcriptome and experimental validated interactions.

Mr.Vc is a database of curated Vibrio cholerae transcriptome data and annotated information. The main objective is to facilitate the accessibility and reusability of the rapidly growing Vibrio cholerae omics data and relevant annotation. To achieve these goals, we performed manual curation on the transcriptome data and organized the datasets in an experiment-centric manner. We collected unknown operons annotated through text-mining analysis that would provide more clues about how Vibrio cholerae modulates gene regulation. Meanwhile, to understand the relationship between genes or experiments, we performed gene co-expression analysis and experiment-experiment correlation analysis. In additional, functional module named "Interactions" which dedicates to collecting experimentally validated interactions about Vibrio cholerae from public databases, MEDLINE documents and literature in life science journals. To date, Mr.Vc v2, which is significantly increased from the previous version, contains 107 microarray experiments, 106 RNA-seq experiments, and 3 Tn-seq projects, covering 56,839 entries of DEGs (Differentially Expressed Genes) from transcriptomes and 7,463 related genes from Tn-seq, respectively. and a total of 270,129 gene co-expression entries and 11,990 entries of experiment-experiment correlation was obtained, in total 1,316 entries of interactions were collected, including 496 protein-chemical signaling molecule interactions, 472 protein-protein interactions, 306 TF (Transcription Factor)-gene interactions and 42 Vibrio cholerae-virus interactions, most of which obtained from 402 literature through text-mining analysis. To make the information easier to access, Mr.Vc v2 is equipped with a search widget, enabling users to query what they are interested in. Mr.Vc v2 is freely available at http://mrvcv2.biownmc.info.

Antibiotic-Induced Primary Biles Inhibit SARS-CoV-2 Endoribonuclease Nsp15 Activity in Mouse Gut.

The gut microbiome profile of COVID-19 patients was found to correlate with a viral load of SARS-CoV-2, COVID-19 severity, and dysfunctional immune responses, suggesting that gut microbiota may be involved in anti-infection. In order to investigate the role of gut microbiota in anti-infection against SARS-CoV-2, we established a high-throughput in vitro screening system for COVID-19 therapeutics by targeting the endoribonuclease (Nsp15). We also evaluated the activity inhibition of the target by substances of intestinal origin, using a mouse model in an attempt to explore the interactions between gut microbiota and SARS-CoV-2. The results unexpectedly revealed that antibiotic treatment induced the appearance of substances with Nsp15 activity inhibition in the intestine of mice. Comprehensive analysis based on functional profiling of the fecal metagenomes and endoribonuclease assay of antibiotic-enriched bacteria and metabolites demonstrated that the Nsp15 inhibitors were the primary bile acids that accumulated in the gut as a result of antibiotic-induced deficiency of bile acid metabolizing microbes. This study provides a new perspective on the development of COVID-19 therapeutics using primary bile acids.