Mutational analysis of 23 autosomal short tandem repeats based on trio paternity testing in the Korean population

This study aimed to estimate A-STR mutation rates in 2,317 Korean parent-child trios by examining 20 Combined DNA Index System (CODIS) core loci (D3S1358, D5S818, D7S820, D8S1179, D13S317, D16S539, D18S51, D21S11, CSF1PO, FGA, TH01, TPOX, vWA, D1S1656, D2S441, D2S1338, D10S1248, D12S391, D19S433, and D22S1045) and three non-CODIS loci (Penta E, Penta D, and SE33). Locus-specific mutation rate estimates varied from 0.00 to 8.63 × 10-3 per generation, with an average mutation rate of 1.62 × 10-3 (95 % CI, 1.39-1.88 × 10-3). We also combined data from previous studies to obtain comprehensive genetic values for the Korean population, and the average mutation rate was 1.59 × 10-3 (95 % CI, 1.38-1.82 × 10-3). Single-step mutations (95.69 %) and double-step mutations (3.35 %) were observed in the mutation pattern analysis, and cases expected to have multi-step mutations (0.96 %) were also observed. Large-sized alleles exhibited more loss mutations than gain mutations, and paternal mutations (62.68 %) were more frequently observed than maternal mutations (19.62 %). The calculated values and features of the 23 A-STRs explored in this study are expected to play a crucial role in establishing criteria for forensic genetic interpretation.

Considerations on the application of a mutation model for Y-STR interpretation

If Y-STR profiling is to be more effective in criminal casework, the methods used to evaluate evidential weight require improvement. Many forensic scientists assign an evidential weight by estimating the number of times a Y-STR profile obtained from a questioned sample has been observed in YHRD datasets. More sophisticated models have been suggested but not yet implemented into routine casework, e.g. Andersen & Balding [1]. Mutation is inherent to STR meiosis (or inheritance) and is encountered in practice. We evaluated a mutation model that can be incorporated into a method for assigning evidential weight to Y-STR profiles, an essential part of bringing any method into practice. Since an important part of implementation to casework is communication, the article is written in an accessible format for practitioners as well as statisticians. The mutation component within the MUTEA model by Willems et al. [2] incorporates the potential for multistep mutations and a tendency for alleles to revert towards a central length, reflecting observed mutation data, e.g. [3]. We have estimated the parameters in this model and in a simplified symmetric version of this model, using sequence data from father/son pairs [4] and deep-rooted pedigrees [5]. Both datasets contain multistep mutations, which may have an effect on models based on simulations [1]. We introduce Beta-Binomial and Beta-Geometric conjugate analyses for estimating rate and step parameters for the mutation models presented here, which require only summations and multiplications. We proved mathematically that the parameters can be estimated independently. We show the importance of reporting the variability of the parameters and not only a point estimate. The parameters can be easily incorporated into statistical models, and updated sequentially as more data becomes available. We recommend fuller publication of data to enable the development and evaluation of a wider range of mutation models.

The paradoxes of Mycobacterium tuberculosis molecular evolution and consequences for the inference of tuberculosis emergence date

The date of Mycobacterium tuberculosis complex emergence has been the subject of long debates. New studies joining archaeological efforts with sequencing methods raise high hopes for solving whether this emergence is closer to 70,000 or to 6000 years before present. Inferring the date of emergence of this pathogen based on sequence data requires a molecular clock. Several clocks inferred from different types of loci and/or different samples, using both sound reasoning and reliable data, are actually very different, which we refer to as the paradoxes of M. tuberculosis molecular evolution. After having presented these paradoxes, we will remind the limits of the molecular clocks used in the different studies such as the assumption of homogeneous substitution rate. We will then review recent results that shed new light on the characteristics of M. tuberculosis molecular evolution: traces of diverse selection pressures, the impact of host dynamics, etc. We provide some ideas on what to do next to get nearer to a reliable dating of Mycobacterium tuberculosis complex emergence. Among them, the collection of additional remains from ancient tuberculosis seems still essential.

Methods for two nonstandard problems arising from the Luria-Delbrück experiment

The fluctuation experiment, devised by Luria and Delbrück in 1943, remains the method of choice for measuring microbial mutation rates in the laboratory. While most inference problems commonly encountered in a fluctuation experiment can be tackled by existing standard algorithms, investigators from time to time run into nonstandard problems not amenable to any existing algorithms. A major obstacle to solving these nonstandard problems is the construction of confidence intervals for mutation rates. This note describes methods for two important categories of nonstandard problems, namely, pooling data from separate experiments and analyzing grouped mutant count data, focusing on the construction of likelihood ratio confidence intervals. In addition to illustrative examples using real-world data, simulation results are presented to help assess the proposed methods.

Signification and Application of Mutator and Antimutator Phenotype-Induced Genetic Variations in Evolutionary Adaptation and Cancer Therapeutics

Mutations present a dichotomy in their implications for cellular processes. They primarily arise from DNA replication errors or damage repair processes induced by environmental challenges. Cumulative mutations underlie genetic variations and drive evolution, yet also contribute to degenerative diseases such as cancer and aging. The mutator phenotype elucidates the heightened mutation rates observed in malignant tumors. Evolutionary adaptation, analogous to bacterial and eukaryotic systems, manifests through mutator phenotypes during changing environmental conditions, highlighting the delicate balance between advantageous mutations and their potentially detrimental consequences. Leveraging the genetic tractability of Saccharomyces cerevisiae offers unique insights into mutator phenotypes and genome instability akin to human cancers. Innovative reporter assays in yeast model organisms enable the detection of diverse genome alterations, aiding a comprehensive analysis of mutator phenotypes. Despite significant advancements, our understanding of the intricate mechanisms governing spontaneous mutation rates and preserving genetic integrity remains incomplete. This review outlines various cellular pathways affecting mutation rates and explores the role of mutator genes and mutation-derived phenotypes, particularly prevalent in malignant tumor cells. An in-depth comprehension of mutator and antimutator activities in yeast and higher eukaryotes holds promise for effective cancer control strategies.

Gene birth in a model of non-genic adaptation

BACKGROUND

Over evolutionary timescales, genomic loci can switch between functional and non-functional states through processes such as pseudogenization and de novo gene birth. Particularly, de novo gene birth is a widespread process, and many examples continue to be discovered across diverse evolutionary lineages. However, the general mechanisms that lead to functionalization are poorly understood, and estimated rates of de novo gene birth remain contentious. Here, we address this problem within a model that takes into account mutations and structural variation, allowing us to estimate the likelihood of emergence of new functions at non-functional loci.

RESULTS

Assuming biologically reasonable mutation rates and mutational effects, we find that functionalization of non-genic loci requires the realization of strict conditions. This is in line with the observation that most de novo genes are localized to the vicinity of established genes. Our model also provides an explanation for the empirical observation that emerging proto-genes are often lost despite showing signs of adaptation.

CONCLUSIONS

Our work elucidates the properties of non-genic loci that make them fertile for adaptation, and our results offer mechanistic insights into the process of de novo gene birth.

Intra-host mutation rate of acute SARS-CoV-2 infection during the initial pandemic wave

SARS-CoV-2 mutation is minimized through a proofreading function encoded by NSP-14. Most estimates of the SARS-CoV-2 mutation rate are derived from population based sequence data. Our understanding of SARS-CoV-2 evolution might be enhanced through analysis of intra-host viral mutation rates in specific populations. Viral genome analysis was performed between paired samples and mutations quantified at allele frequencies (AF) ≥ 0.25, ≥ 0.5 and ≥ 0.75. Mutation rate was determined employing F81 and JC69 evolution models and compared between isolates with (ΔNSP-14) and without (wtNSP-14) non-synonymous mutations in NSP-14 and by patient comorbidity. Forty paired samples with median interval of 13 days [IQR 8.5-20] were analyzed. The estimated mutation rate by F81 modeling was 93.6 (95%CI 90.8-96.4], 40.7 (95%CI 38.9-42.6) and 34.7 (95%CI 33.0-36.4) substitutions/genome/year at AF ≥ 0.25, ≥ 0.5, ≥ 0.75 respectively. Mutation rate in ΔNSP-14 were significantly elevated at AF ≥ 0.25 vs wtNSP-14. Patients with immune comorbidities had higher mutation rate at all allele frequencies. Intra-host SARS-CoV-2 mutation rates are substantially higher than those reported through population analysis. Virus strains with altered NSP-14 have accelerated mutation rate at low AF. Immunosuppressed patients have elevated mutation rate at all AF. Understanding intra-host virus evolution will aid in current and future pandemic modeling.

Sequence dependencies and mutation rates of localized mutational processes in cancer

BACKGROUND

Cancer mutations accumulate through replication errors and DNA damage coupled with incomplete repair. Individual mutational processes often show nucleotide sequence and functional region preferences. As a result, some sequence contexts mutate at much higher rates than others, with additional variation found between functional regions. Mutational hotspots, with recurrent mutations across cancer samples, represent genomic positions with elevated mutation rates, often caused by highly localized mutational processes.

METHODS

We count the 11-mer genomic sequences across the genome, and using the PCAWG set of 2583 pan-cancer whole genomes, we associate 11-mers with mutational signatures, hotspots of single nucleotide variants, and specific genomic regions. We evaluate the mutation rates of individual and combined sets of 11-mers and derive mutational sequence motifs.

RESULTS

We show that hotspots generally identify highly mutable sequence contexts. Using these, we show that some mutational signatures are enriched in hotspot sequence contexts, corresponding to well-defined sequence preferences for the underlying localized mutational processes. This includes signature 17b (of unknown etiology) and signatures 62 (POLE deficiency), 7a (UV), and 72 (linked to lymphomas). In some cases, the mutation rate and sequence preference increase further when focusing on certain genomic regions, such as signature 62 in transcribed regions, where the mutation rate is increased up to 9-folds over cancer type and mutational signature average.

CONCLUSIONS

We summarize our findings in a catalog of localized mutational processes, their sequence preferences, and their estimated mutation rates.

The high mutation rate at the D614G hotspot-furin cleavage site region increases the priming efficiency of the Spike protein by furin protease: analysis of Indonesian SARS-CoV-2 G614 variants obtained during the early COVID-19 pandemic

D614G mutation plays a significant role in the transmissibility of SARS-CoV-2. Identification of other mutations related to D614G mutation within the Spike protein is pivotal as they might contribute to the pathogenicity of SARS-CoV-2. This study aims to analyze the mutation rate of furin cleavage site (FCS) region of Indonesian origin SARS-CoV-2 and to predict the effect of mutation against Spike priming efficiency by furin. A total of 375 sequences of Indonesian isolates obtained during the early pandemic were used for mutation analysis. Mutation analysis includes mutation pattern, variability, frequency of mutation, amino acid conservation, and mutation rate. The effect of mutation against Spike priming efficiency by furin protease from eight sequences with mutation in the FCS region was analyzed by protein-protein docking. We showed that mutations related to the G614 variant were increasing through time, in contrast to the D614 variant. The FCS region at the position 675-692 contained the most variable (66.67%) as well as the highest mutation frequency (85.92%) and has been observed to be the hotspot mutations linked to the D614G mutation. The D614G hotspot-FCS region (residue 600-700) had the highest amino acid change per site (20.8%) as well as the highest mutation rate as 1.34 × 10^-2 substitution per site per year (95% CI 1.79 × 10^-3-2.74 × 10^-2), compared with other Spike protein regions. Mutations in the FCS region were the most common mutation found after the D614G mutation. These mutations were predicted to increase the Spike priming efficiency by furin. Thus, this study elucidates the importance of D614G mutation to other mutations located in the FCS region and their significance to Spike priming efficiency by furin.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13337-023-00827-w.

Analysis of mutation rates and haplotypes of 23 Y-chromosomal STRs in Korean father-son pairs

Y-chromosomal short tandem repeats (Y-STRs) have been widely used in forensic genetics, and accurate knowledge of mutation rates at Y-STR loci is essential in kinship analysis. The main aim of this study was to estimate Y-STR mutation rates in Korean males. To obtain locus-specific mutations and haplotypes at 23 Y-STRs, we analyzed samples from 620 Korean father-son pairs. In addition, we also analyzed 476 unrelated individuals using the PowerPlex® Y23 System, with the aim of augmenting the available data for the Korean population. The PowerPlex® Y23 system facilitates analysis of the 23 Y-STR loci (DYS576, DYS570, DYS458, DYS635, DYS389 II, DYS549, DYS385, DYS481, DYS439, DYS456, DYS389 I, DYS19, DYS393, DYS391, DYS533, DYS437, DYS390, Y GATA H4, DYS448, DYS438, DYS392, and DYS643). Locus-specific mutation rate estimates varied from 0.00 to 8.06 × 10^-3 per generation, with an average mutation rate of 2.17 × 10^-3 (95% CI, 1.5-3.1 × 10^-3). To obtain comprehensive genetic values for the Korean population, we combined data obtained in this study with previously reported values, thereby enabling us to estimate the locus-specific mutation rates regarding 22,711 allele transmissions. By combining these data, we obtained an overall average mutation rate of 2.91 × 10^-3 (95% CI, 2.3-3.7 × 10^-3). In addition, among the 476 unrelated Korean males, we detected 467 different haplotypes, with an overall haplotype diversity value of 0.9999. By extracting haplotypes of Y-STRs described in previous literature on 23 Y-STR reported in Korea, we obtained gene diversities for 1133 Korean individuals. We believe that the values and characteristics of the 23 Y-STRs analyzed in this study will contribute to establishing criteria for forensic genetic interpretation, including kinship analysis.

The how of counter-defense: viral evolution to combat host immunity

Viruses are locked in an evolutionary arms race with their hosts. What ultimately determines viral evolvability, or capacity for adaptive evolution, is their ability to efficiently explore and expand sequence space while under the selective regime imposed by their ecology, which includes innate and adaptive host defenses. Viral genomes have significantly higher evolutionary rates than their host counterparts and should have advantages relative to their slower-evolving hosts. However, functional constraints on virus evolutionary landscapes along with the modularity and mutational tolerance of host defense proteins may help offset the advantage conferred to viruses by high evolutionary rates. Additionally, cellular life forms from all domains of life possess many highly complex defense mechanisms that act as hurdles to viral replication. Consequently, viruses constantly probe sequence space through mutation and genetic exchange and are under pressure to optimize diverse counter-defense strategies.

Efficient, robust, and versatile fluctuation data analysis using MLE MUtation Rate calculator (mlemur)

The fluctuation assay remains an important tool for analyzing the levels of mutagenesis in microbial populations. The mutant counts originating from some average number of mutations are usually assumed to obey the Luria-Delbrück distribution. While several tools for estimating mutation rates are available, they sometimes lack accuracy or versatility under non-standard conditions. In this work, extensions to the Luria-Delbrück protocol to account for phenotypic lag and cellular death with either perfect or partial plating were developed. Hence, the novel MLE MUtation Rate calculator, or mlemur, is the first tool that provides a user-friendly graphical interface allowing the researchers to model their data with consideration for partial plating, differential growth of mutants and non-mutants, phenotypic lag, cellular death, variability of the final number of cells, post-exponential-phase mutations, and the size of the inoculum. Additionally, mlemur allows the users to incorporate most of these special conditions at the same time to obtain highly accurate estimates of mutation rates and P values, confidence intervals for an arbitrary function of data (such as fold), and perform power analysis and sample size determination for the likelihood ratio test. The accuracy of point and interval estimates produced by mlemur against historical and simulated fluctuation experiments are assessed. Both mlemur and the analyses in this work might be of great help when evaluating fluctuation experiments and increase the awareness of the limitations of the widely-used Lea-Coulson formulation of the Luria-Delbrück distribution in the more realistic biological contexts.

[Polymorphism and Mutation Rate of 20 Rapidly Mutating Y-Chromosomal Short Tandem Repeats in Chinese Han Population of Sichuan Province]

Objective

To explore the applicability of 20 rapidly mutating Y-chromosomal short tandem repeats (RM Y-STRs) in Chinese Han population of Sichuan province.

Methods

Two RM Y-STR multiple amplification systems (RM1, including DYF404S1, DYF399S1, DYS547, DYS526a/b, DYS626, DYF403S1a/b, and DYS612, and RM2, including DYS1003, DYS1007, DYR88, DYS712, DYS711, DYS724, and DYF1002, with 14 RM Y-STR loci in total) and Y41SE-V1.2 (including 6 RM Y-STR loci of DYS627, DYS576, DYF387S1, DYS518, DYS570, and DYS449, 30 ordinary Y-chromosomal short tandem repeats [Y-STR] loci, and 1 Indel locus) were used for the amplification and typing of 200 unrelated males and 260 father-son pairs. The polymorphisms and mutation rates of 20 RM Y-STRs and 30 ordinary Y-STRs in Chinese Han population of Sichuan province were investigated and compared.

Results

In the 200 unrelated males, the gene diversity (GD) of 20 RM Y-STR loci ranged from 0.7910 to 0.9975, and there were 200 haplotypes. Haplotype diversity (HD) was 1 and the discriminative capacity (DC) was 1. A total of 198 haplotypes were found in Y41se-v1.2 (the 30 Y-STRs), with 4 cases sharing two haplotypes, the haplotype diversity being 0.9999, and the discriminative capacity being 0.99. A total of 68 mutations were found at the 20 RM Y-STRs loci in the 260 father-son pairs, and there was slightly more increase than decrease of allele repeats (1.19∶1), with the mutation rate ranging from <3.85×10 -3 (95% C I: 0.00-1.41×10 -2) to 2.69×10 -2 (95% CI: 1.09×10 -2-5.47×10 -2), and the average mutation rate being 1.19×10 -2 (95% CI: 9.20×10 -3-1.51×10 -2). The 20 RM Y-STRs and the Y41SE-V1.2 (the 30 Y-STRs) could be used to distinguish 22.3% and 13.8% father-son pairs, respectively.

Conclusion

The 20 RM Y-STRs have high gene and haplotype diversity and paternal lineage differentiation rate in Chinese Han population of Sichuan province, showing great potential for application in Chinese Han population of Sichuan province.

Development of a sensitive double TaqMan Probe-based qPCR Angle-Degree method to detect mutation frequencies

We designed a method to examine the mutation frequencies of the A3243G mutation of mitochondrial DNA (mtDNA) in patients with mitochondrial encephalopathy, lactic acidosis, and stroke-like episodes (MELAS) syndrome. We performed a qPCR assay using the FAM and VIC TaqMan probes, which detect the 3243G (mutated) and 3243A (wild-type) sequences of mtDNA, respectively. The results obtained by "degree" in a series of differential mutation frequencies were used to plot a standard curve of the mutation frequency. The standard curve was then applied for qPCR assays of the desired samples. The standard deviation (%) of the samples calculated using the standard curve for the TaqMan probe was 2.4 ± 1.5%. This method could be used to examine mutation frequencies in the context of diabetes, aging, cancer, and neurodegenerative diseases.

Very Low Rates of Spontaneous Gene Deletions and Gene Duplications in Dictyostelium discoideum

The study of spontaneous mutation rates has revealed a wide range of heritable point mutation rates across species, but there are comparatively few estimates for large-scale deletion and duplication rates. The handful of studies that have directly calculated spontaneous rates of deletion and duplication using mutation accumulation lines have estimated that genes are duplicated and deleted at orders of magnitude greater rates than the spontaneous point mutation rate. In our study, we tested whether spontaneous gene deletion and gene duplication rates are also high in Dictyostelium discoideum, a eukaryote with among the lowest point mutation rates (2.5 × 10-11 per site per generation) and an AT-rich genome (GC content of 22%). We calculated mutation rates of gene deletions and duplications using whole-genome sequencing data originating from a mutation accumulation experiment and determined the association between the copy number mutations and GC content. Overall, we estimated an average of 3.93 × 10-8 gene deletions and 1.18 × 10-8 gene duplications per gene per generation. While orders of magnitude greater than their point mutation rate, these rates are much lower compared to gene deletion and duplication rates estimated from mutation accumulation lines in other organisms (that are on the order of ~ 10-6 per gene/generation). The deletions and duplications were enriched in regions that were AT-rich even compared to the genomic background, in contrast to our expectations if low GC content was contributing to low mutation rates. The low deletion and duplication mutation rates in D. discoideum compared to other eukaryotes mirror their low point mutation rates, supporting previous work suggesting that this organism has high replication fidelity and effective molecular machinery to avoid the accumulation of mutations in their genome.

Estimation of Rates of Non-neutral Mutations When Bacteria are Exposed to Subinhibitory Levels of Antibiotics

The increasingly widespread use of antibiotics has raised an urgent public health question: how does nonlethal exposure of antibiotics increase rates of bacterial mutations that confer antibiotic resistance? Quantitative studies of this sort are rare, and interested investigators tend to neglect scenarios where the mutation inducer and the selective agent are the same. One possible reason for this puzzling neglect is that such an investigation may lead to results that are prone to a Lamarckian interpretation, but a more tangible reason is that successful estimation of mutation rates in such scenarios requires new computational methods. This study presents computational methods tailored for a mutation model in which some wide-type cells may be killed by nonlethal exposure to an antibiotic, but in which mutants proliferate unimpeded. Methods for computing maximum likelihood estimates of mutation rates and their corresponding confidence intervals are devised and evaluated by simulations.

Biocides used as material preservatives modify rates of de novo mutation and horizontal gene transfer in bacteria

Antimicrobial resistance (AMR) is a global health problem with the environment being an important compartment for the evolution and transmission of AMR. Previous studies showed that de-novo mutagenesis and horizontal gene transfer (HGT) by conjugation or transformation - important processes underlying resistance evolution and spread - are affected by antibiotics, metals and pesticides. However, natural microbial communities are also frequently exposed to biocides used as material preservatives, but it is unknown if these substances induce mutagenesis and HGT. Here, we show that active substances used in material preservatives can increase rates of mutation and conjugation in a species- and substance-dependent manner, while rates of transformation are not increased. The bisbiguanide chlorhexidine digluconate, the quaternary ammonium compound didecyldimethylammonium chloride, the metal copper, the pyrethroid-insecticide permethrin, and the azole-fungicide propiconazole increase mutation rates in Escherichia coli, whereas no increases were identified for Bacillus subtilis and Acinetobacter baylyi. Benzalkonium chloride, chlorhexidine and permethrin increased conjugation in E. coli. Moreover, our results show a connection between the RpoS-mediated general stress and the RecA-linked SOS response with increased rates of mutation and conjugation, but not for all biocides. Taken together, our data show the importance of assessing the contribution of material preservatives on AMR evolution and spread.

RMplex reveals population differences in RM Y-STR mutation rates and provides improved father-son differentiation in Japanese

Rapidly mutating Y chromosomal short tandem repeat markers (RM Y-STRs) -characterized by at least one mutation per 100 generations- are suitable for differentiating both related and unrelated males. The recently introduced multiplex method RMplex allows for the efficient analysis of 30 Y-STRs with increased mutation rates, including all 26 currently known RM Y-STRs. While currently available RM Y-STR mutation rates were established mostly from European individuals, here we applied RMplex to DNA samples of 178 genetically confirmed father-son pairs from East Asia. For several Y-STRs, we found significantly higher mutation rates in Japanese compared to previous estimates. The consequent father-son differentiation rate based on RMplex was significantly higher (52%) in Japanese than previously reported for Europeans (42%), and much higher than with Yfiler Plus in both sample sets (14% and 13%, respectively). Further analysis suggests that the higher mutation and relative differentiation rates in Japanese can in part be explained by on average longer Y-STR alleles relative to Europeans. Moreover, we show that the most striking difference, which was found in DYS712, could be linked to a Y-SNP haplogroup (O1b2-P49) that is common in Japanese and rare in other populations. We encourage the forensic Y-STR community to generate more RMplex data from more population samples of sufficiently large sample size in combination with Y-SNP data to further investigate population effects on mutation and relative differentiation rates. Until more RMplex data from more populations become available, caution shall be placed when applying RM Y-STR mutation rate estimates established in one population, such as Europeans, to forensic casework involving male suspects of paternal origin from other populations, such as non-Europeans.

Frequency of TERT Promoter Mutations in Real-World Analysis of 2,092 Thyroid Carcinoma Patients

BACKGRUOUND

Telomerase reverse transcriptase (TERT) promoter mutations are associated with increased recurrence and mortality in patients with thyroid carcinoma. Previous studies on TERT promoter mutations were retrospectively conducted on a limited number of patients.

METHODS

We prospectively collected data on all consecutive patients who underwent thyroid carcinoma surgery between January 2019 and December 2020 at the Samsung Medical Center, Seoul, Korea. We included 2,092 patients with thyroid carcinoma.

RESULTS

Of 2,092 patients, 72 patients (3.4%) had TERT promoter mutations. However, the frequency of TERT promoter mutations was 0.5% in papillary thyroid microcarcinoma (PTMC) ≤1 cm and it was 5.8% in papillary thyroid carcinoma (PTC) >1 cm. The frequency of TERT promoter mutations was significantly associated with older age at diagnosis (odds ratio [OR], 1.12; P<0.001), larger primary tumor size (OR, 2.02; P<0.001), and aggressive histological type (OR, 7.78 in follicular thyroid carcinoma; OR, 10.33 in poorly differentiated thyroid carcinoma; OR, 45.92 in anaplastic thyroid carcinoma; P<0.001). Advanced T stage, advanced N stage, and distant metastasis at diagnosis were highly prevalent in mutated thyroid cancers. However, initial distant metastasis was not present in patients with TERT promoter mutations in PTMC. Although the C228T mutation was more highly detected than the C250T mutation (64 cases vs. 7 cases), there were no significant clinicopathological differences.

CONCLUSION

This study is the first attempt to investigate the frequency of TERT promoter mutations in a real-world setting. The frequency of TERT promoter mutations in PTC was lower than expected, and in PTMC, young patients, and female patients, the frequency was very low.

Emerging SARS-CoV-2 variants: Why, how, and what's next?

The emergence of the SARS-CoV-2 Omicron variant poses a striking threat to human society. More than 30 mutations in the Spike protein of the Omicron variant severely compromised the protective immunity elicited by either vaccination or prior infection. The persistent viral evolutionary trajectory generates Omicron-associated lineages, such as BA.1 and BA.2. Moreover, the virus recombination upon Delta and Omicron co-infections has been reported lately, although the impact remains to be assessed. This minireview summarizes the characteristics, evolution and mutation control, and immune evasion mechanisms of SARS-CoV-2 variants, which will be helpful for the in-depth understanding of the SARS-CoV-2 variants and policy-making related to COVID-19 pandemic control.

Large-scale discovery of novel neurodevelopmental disorder-related genes through a unified analysis of single-nucleotide and copy number variants

Background

Previous large-scale studies of de novo variants identified a number of genes associated with neurodevelopmental disorders (NDDs); however, it was also predicted that many NDD-associated genes await discovery. Such genes can be discovered by integrating copy number variants (CNVs), which have not been fully considered in previous studies, and increasing the sample size.

Methods

We first constructed a model estimating the rates of de novo CNVs per gene from several factors such as gene length and number of exons. Second, we compiled a comprehensive list of de novo single-nucleotide variants (SNVs) in 41,165 individuals and de novo CNVs in 3675 individuals with NDDs by aggregating our own and publicly available datasets, including denovo-db and the Deciphering Developmental Disorders study data. Third, summing up the de novo CNV rates that we estimated and SNV rates previously established, gene-based enrichment of de novo deleterious SNVs and CNVs were assessed in the 41,165 cases. Significantly enriched genes were further prioritized according to their similarity to known NDD genes using a deep learning model that considers functional characteristics (e.g., gene ontology and expression patterns).

Results

We identified a total of 380 genes achieving statistical significance (5% false discovery rate), including 31 genes affected by de novo CNVs. Of the 380 genes, 52 have not previously been reported as NDD genes, and the data of de novo CNVs contributed to the significance of three genes (GLTSCR1, MARK2, and UBR3). Among the 52 genes, we reasonably excluded 18 genes [a number almost identical to the theoretically expected false positives (i.e., 380 × 0.05 = 19)] given their constraints against deleterious variants and extracted 34 “plausible” candidate genes. Their validity as NDD genes was consistently supported by their similarity in function and gene expression patterns to known NDD genes. Quantifying the overall similarity using deep learning, we identified 11 high-confidence (> 90% true-positive probabilities) candidate genes: HDAC2, SUPT16H, HECTD4, CHD5, XPO1, GSK3B, NLGN2, ADGRB1, CTR9, BRD3, and MARK2.

Conclusions

We identified dozens of new candidates for NDD genes. Both the methods and the resources developed here will contribute to the further identification of novel NDD-associated genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-022-01042-w.

Genome-wide identification and prediction of SARS-CoV-2 mutations show an abundance of variants: Integrated study of bioinformatics and deep neural learning

Genomic data analysis is a fundamental system for monitoring pathogen evolution and the outbreak of infectious diseases. Based on bioinformatics and deep learning, this study was designed to identify the genomic variability of SARS-CoV-2 worldwide and predict the impending mutation rate. Analysis of 259044 SARS-CoV-2 isolates identified 3334545 mutations with an average of 14.01 mutations per isolate. Globally, single nucleotide polymorphism (SNP) is the most prevalent mutational event. The prevalence of C > T (52.67%) was noticed as a major alteration across the world followed by the G > T (14.59%) and A > G (11.13%). Strains from India showed the highest number of mutations (48) followed by Scotland, USA, Netherlands, Norway, and France having up to 36 mutations. D416G, F106F, P314L, UTR:C241T, L93L, A222V, A199A, V30L, and A220V mutations were found as the most frequent mutations. D1118H, S194L, R262H, M809L, P314L, A8D, S220G, A890D, G1433C, T1456I, R233C, F263S, L111K, A54T, A74V, L183A, A316T, V212F, L46C, V48G, Q57H, W131R, G172V, Q185H, and Y206S missense mutations were found to largely decrease the structural stability of the corresponding proteins. Conversely, D3L, L5F, and S97I were found to largely increase the structural stability of the corresponding proteins. Multi-nucleotide mutations GGG > AAC, CC > TT, TG > CA, and AT > TA have come up in our analysis which are in the top 20 mutational cohort. Future mutation rate analysis predicts a 17%, 7%, and 3% increment of C > T, A > G, and A > T, respectively in the future. Conversely, 7%, 7%, and 6% decrement is estimated for T > C, G > A, and G > T mutations, respectively. T > G\A, C > G\A, and A > T\C are not anticipated in the future. Since SARS-CoV-2 is mutating continuously, our findings will facilitate the tracking of mutations and help to map the progression of the COVID-19 intensity worldwide.

Multivariate phase-type theory for the site frequency spectrum

Linear functions of the site frequency spectrum (SFS) play a major role for understanding and investigating genetic diversity. Estimators of the mutation rate (e.g. based on the total number of segregating sites or average of the pairwise differences) and tests for neutrality (e.g. Tajima's D) are perhaps the most well-known examples. The distribution of linear functions of the SFS is important for constructing confidence intervals for the estimators, and to determine significance thresholds for neutrality tests. These distributions are often approximated using simulation procedures. In this paper we use multivariate phase-type theory to specify, characterize and calculate the distribution of linear functions of the site frequency spectrum. In particular, we show that many of the classical estimators of the mutation rate are distributed according to a discrete phase-type distribution. Neutrality tests, however, are generally not discrete phase-type distributed. For neutrality tests we derive the probability generating function using continuous multivariate phase-type theory, and numerically invert the function to obtain the distribution. A main result is an analytically tractable formula for the probability generating function of the SFS. Software implementation of the phase-type methodology is available in the R package PhaseTypeR, and R code for the reproduction of our results is available as an accompanying vignette.

Sex without crossing over in the yeast Saccharomycodes ludwigii

BACKGROUND

Intermixing of genomes through meiotic reassortment and recombination of homologous chromosomes is a unifying theme of sexual reproduction in eukaryotic organisms and is considered crucial for their adaptive evolution. Previous studies of the budding yeast species Saccharomycodes ludwigii suggested that meiotic crossing over might be absent from its sexual life cycle, which is predominated by fertilization within the meiotic tetrad.

RESULTS

We demonstrate that recombination is extremely suppressed during meiosis in Sd. ludwigii. DNA double-strand break formation by the conserved transesterase Spo11, processing and repair involving interhomolog interactions are required for normal meiosis but do not lead to crossing over. Although the species has retained an intact meiotic gene repertoire, genetic and population analyses suggest the exceptionally rare occurrence of meiotic crossovers in its genome. A strong AT bias of spontaneous mutations and the absence of recombination are likely responsible for its unusually low genomic GC level.

CONCLUSIONS

Sd. ludwigii has followed a unique evolutionary trajectory that possibly derives fitness benefits from the combination of frequent mating between products of the same meiotic event with the extreme suppression of meiotic recombination. This life style ensures preservation of heterozygosity throughout its genome and may enable the species to adapt to its environment and survive with only minimal levels of rare meiotic recombination. We propose Sd. ludwigii as an excellent natural forum for the study of genome evolution and recombination rates.

Rates and patterns of molecular evolution in bryophyte genomes, with focus on complex thalloid liverworts, Marchantiopsida

Plants commonly referred to as "bryophytes" belong to three major lineages of non-vascular plants: the liverworts, the hornworts and the mosses. They are unique among land plants in having a dominant haploid generation and a short-lived diploid sporophytic generation. The dynamics of selection acting on a haploid genome differs from those acting on a diploid genome: new mutations are directly exposed to selection. The general aim of this paper is to investigate the diversification rateof bryophytes - measured as silent site substitution rate representing neutral evolution (mutation rate) and the nonsynonymous to synonymous substitution rate ratio (dN/dS) representing selective evolution - and compare it with earlier studies on vascular plants. Results show that the silent site substitution rate is lower for liverworts as compared to angiosperms, but not as low as for gymnosperms. The selection pressure, measured as dN/dS, isnot remarkably lower for bryophytes as compared to other diploid dominant plants as would be expected by the masking hypothesis, indicating that other factors are more important than ploidy.

Mutation analysis of 28 autosomal short tandem repeats in the Chinese Han population

Short tandem repeats (STRs) have been extensively used in forensic genetics. However, according to previous studies, the mutation rates of STRs are relatively high and are affected by many factors. Therefore, it is important to analyze STR mutations and determine the influence of underlying factors on STR mutation rates. Mutation rates of 28 autosomal STRs were determined from 8708 paternity testing cases in the Chinese Han population, and the relationships between STR mutation rates and population, sex, age, allele length and heterozygosity were investigated. A total of 279 mutations were observed at 27 loci in a total of 233,530 meiosis cases, including 273 (97.8%) one-step, 5 (1.8%) two-step and 1 (0.4%) three-step mutations. The overall average mutation rate was 1.19 × 10^-3 (95% CI 1.06 × 10^-3 - 1.34 × 10^-3) ranging from 0 (TPOX) to 2.79 × 10^-3 (D13S325). Mutation rate comparisons revealed statistically significant differences at several STRs among populations. Paternal mutations occurred more frequently than maternal mutations, at a ratio of 6.04:1, and the mutation rate tended to increase with paternal age. Moreover, our study revealed a bias towards contraction mutations for long alleles and expansion mutations for short alleles. No obvious bias was observed in the overall mutation direction. In addition, STR loci with higher expected heterozygosity (H_exp) tended to have higher mutation rates. This work revealed the relationships between STR mutation rates and several influencing factors, providing useful data and information for further research on STR mutations in forensic genetics.

The Y-STR landscape of coastal southeastern Han: Forensic characteristics, haplotype analyses, mutation rates, and population genetics

The Y-STR landscape of Coastal Southeastern Han (CSEH) living in Chinese southeast areas (including Guangdong, Fujian, and Zhejiang provinces) is still unclear. We investigated 62 Y-STR markers in a reasonably large number of 1021 unrelated males and 1027 DNA-confirmed father-son pairs to broaden the genetic backgrounds of CSEH. In total, 85 null alleles, 121 off-ladder alleles, and 95 copy number variants were observed, and 1012 distinct haplotypes were determined with the overall HD and DC values of 0.999974 and 0.9912. We observed 369 mutations in 76 099 meiotic transfers, and the average estimated Y-STR mutation rate was 4.85 × 10^-3 (95% CI, 4.4 × 10^-3 -5.4 × 10^-3 ). The Spearman correlation analyses indicated that GD values (R² = 0.6548) and average allele sizes (R² = 0.5989) have positive correlations with Y-STR mutation rates. Our RM Y-STR set including 8 candidate RM Y-STRs, of which DYS534, DYS630, and DYS713 are new candidates in CSEH, distinguished 18.52% of father-son pairs. This study also clarified the population structures of CSEH which isolated in population-mixed South China relatively. The strategy, SM Y-STRs for familial searching and RM Y-STRs for individual identification regionally, could be applicable based on enough knowledge of the Y-STR mutability of different populations.

Base-substitution mutation rate across the nuclear genome of Alpheus snapping shrimp and the timing of isolation by the Isthmus of Panama

Background

The formation of the Isthmus of Panama and final closure of the Central American Seaway (CAS) provides an independent calibration point for examining the rate of DNA substitutions. This vicariant event has been widely used to estimate the substitution rate across mitochondrial genomes and to date evolutionary events in other taxonomic groups. Nuclear sequence data is increasingly being used to complement mitochondrial datasets for phylogenetic and evolutionary investigations; these studies would benefit from information regarding the rate and pattern of DNA substitutions derived from the nuclear genome.

Results

To estimate the genome-wide neutral mutation rate (µ), genotype-by-sequencing (GBS) datasets were generated for three transisthmian species pairs in Alpheus snapping shrimp. A range of bioinformatic filtering parameters were evaluated in order to minimize potential bias in mutation rate estimates that may result from SNP filtering. Using a Bayesian coalescent approach (G-PhoCS) applied to 44,960 GBS loci, we estimated µ to be 2.64E−9 substitutions/site/year, when calibrated with the closure of the CAS at 3 Ma. Post-divergence gene flow was detected in one species pair. Failure to account for this post-split migration inflates our substitution rate estimates, emphasizing the importance of demographic methods that can accommodate gene flow.

Conclusions

Results from our study, both parameter estimates and bioinformatic explorations, have broad-ranging implications for phylogeographic studies in other non-model taxa using reduced representation datasets. Our best estimate of µ that accounts for coalescent and demographic processes is remarkably similar to experimentally derived mutation rates in model arthropod systems. These results contradicted recent suggestions that the closure of the Isthmus was completed much earlier (around 10 Ma), as mutation rates based on an early calibration resulted in uncharacteristically low genomic mutation rates. Also, stricter filtering parameters resulted in biased datasets that generated lower mutation rate estimates and influenced demographic parameters, serving as a cautionary tale for the adherence to conservative bioinformatic strategies when generating reduced-representation datasets at the species level. To our knowledge this is the first use of transisthmian species pairs to calibrate the rate of molecular evolution from GBS data.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01836-3.

Genome-wide DNA mutations in Arabidopsis plants after multigenerational exposure to high temperatures

BACKGROUND

Elevated temperatures can cause physiological, biochemical, and molecular responses in plants that can greatly affect their growth and development. Mutations are the most fundamental force driving biological evolution. However, how long-term elevations in temperature influence the accumulation of mutations in plants remains unknown.

RESULTS

Multigenerational exposure of Arabidopsis MA (mutation accumulation) lines and MA populations to extreme heat and moderate warming results in significantly increased mutation rates in single-nucleotide variants (SNVs) and small indels. We observe distinctive mutational spectra under extreme and moderately elevated temperatures, with significant increases in transition and transversion frequencies. Mutation occurs more frequently in intergenic regions, coding regions, and transposable elements in plants grown under elevated temperatures. At elevated temperatures, more mutations accumulate in genes associated with defense responses, DNA repair, and signaling. Notably, the distribution patterns of mutations among all progeny differ between MA populations and MA lines, suggesting that stronger selection effects occurred in populations. Methylation is observed more frequently at mutation sites, indicating its contribution to the mutation process at elevated temperatures. Mutations occurring within the same genome under elevated temperatures are significantly biased toward low gene density regions, special trinucleotides, tandem repeats, and adjacent simple repeats. Additionally, mutations found in all progeny overlap significantly with genetic variations reported in 1001 Genomes, suggesting non-uniform distribution of de novo mutations through the genome.

CONCLUSION

Collectively, our results suggest that elevated temperatures can accelerate the accumulation, and alter the molecular profiles, of DNA mutations in plants, thus providing significant insight into how environmental temperatures fuel plant evolution.

Response of Escherichia coli minimal ter operon to UVC and auto-aggregation: pilot study

Aim

The study of minimal ter operon as a determinant of tellurium resistance (Te^R) is important for the purpose of confirming the relationship of these genes to the pathogenicity of microorganisms. The ter operon is widespread among bacterial species and pathogens, implicated also in phage inhibition, oxidative stress and colicin resistance. So far, there is no experimental evidence for the role of the Escherichia coli (E. coli) minimal ter operon in ultraviolet C (UVC) resistance, biofilm formation and auto-aggregation. To identify connection with UVC resistance of the minimal ter operon, matched pairs of Ter-positive and -negative E. coli cells were stressed and differences in survival and whole genome sequence analysis were performed. This study was aimed also to identify differences in phenotype of cells induced by environmental stress.

Methods

In the current study, a minimal ter operon(terBCDEΔF) originating from the uropathogenic strain E. coli KL53 was used. Clonogenic assay was the method of choice to determine cell reproductive death after treatment with UVC irradiation at certain time intervals. Bacterial suspensions were irradiated with 254 nm UVC-light (germicidal lamp in biological safety cabinet) in vitro. UVC irradiance output was 2.5 mW/cm² (calculated at the UVC device aperture) and plate-lamp distance of 60 cm. DNA damage analysis was performed using shotgun sequencing on Illumina MiSeq platform. Biofilm formation was measured by a crystal violet retention assay. Auto-aggregation assay was performed according to the Ghane, Babaeekhou & Ketabi (2020).

Results

A large fraction of Ter-positive E. coli cells survived treatment with 120-s UVC light (300 mJ/cm²) compared to matched Ter-negative cells; ∼5-fold higher resistance of Ter-positive cells to UVC dose (p = 0.0007). Moreover, UVC surviving Ter-positive cells showed smaller mutation rate as Ter-negative cells. The study demonstrated that a 1200-s exposure to UVC (3,000 mJ/cm²) was sufficient for 100% inhibition of growth for all the Ter-positive and -negative E. coli cells. The Ter-positive strain exhibited of 26% higher auto-aggregation activities and was able to inhibit biofilm formation over than Ter- negative strain (**** P < 0.0001).

Conclusion

Our study shows that Ter-positive cells display lower sensitivity to UVC radiation, corresponding to a presence in minimal ter operon. In addition, our study suggests that also auto-aggregation ability is related to minimal ter operon. The role of the minimal ter operon (terBCDEΔF) in resistance behavior of E. coli under environmental stress is evident.

A game-theoretical dynamic imitation model on networks

A game-theoretical model is constructed to capture the effect of imitation on the evolution of cooperation. This imitation describes the case where successful individuals are more likely to be imitated by newcomers who will employ their strategies and social networks. Two classical repeated strategies 'always defect (ALLD)' and 'tit-for-tat (TFT)' are adopted. Mathematical analyses are mainly conducted by the method of coalescence theory. Under the assumption of a large population size and weak selection, the results show that the evolution of cooperation is promoted in this dynamic network. As we observed that the critical benefit-to-cost ratio is smaller compared to that in well-mixed populations. The critical benefit-to-cost ratio approaches a specific value which depends on three parameters, the repeated rounds of the game, the effective strategy mutation rate, and the effective link mutation rate. Specifically, for a very high value of the effective link mutation rate, the critical benefit-to-cost ratio approaches 1. Remarkably, for a low value of the effective link mutation rate, by letting the effective strategy mutation is nearly equal to zero, the critical benefit-to-cost ratio approaches [Formula: see text] for the resulting highly connected networks, which allows TFT to be evolutionary stable. It illustrates that dominance of TFTs is associated with more connected networks. This research can enrich the theory of the coevolution of game strategy and network structure with dynamic imitation.

Mutation analysis for 25 Y-STR markers in Japanese population

In this study, we analyzed DNA samples from 213 Japanese father son pairs with 25 Y-chromosome short tandem repeat (Y-STR) (DYS576, DYS389I, DYS635, DYS389II, DYS627, DYS460, DYS458, DYS19, YGATAH4, DYS448, DYS391, DYS456, DYS390, DYS438, DYS392, DYS518, DYS570, DYS437, DYS385, DYS449, DYS393, DYS439, DYS481, DYF387S1, and DYS533) markers using the Yfiler™ Plus PCR amplification kit. We calculated Y-STR mutation rates for each locus to evaluate the efficacy of the 25 Y-STR markers for paternity testing and forensic identification using samples from male relatives. Six rapidly mutating Y-STR markers (DYS576, DYS627, DYS518, DYS570, DYS449 and DYF387S1), previously reported to have high mutation rates (>1.0 × 10^-2), are included in the 25 Y-STR markers, but our findings revealed that the mutation rates for all Y-STR markers except for DYS576 and DYS458 were lower than 1.0 × 10^-2. Therefore, the use of these 25 Y-STR markers may be useful for forensic identification in the Japanese population.

Sum of High-Risk Gene Mutation (SHGM): A Novel Attempt to Assist Differential Diagnosis for Adrenocortical Carcinoma with Benign Adenoma, Based on Detection of Mutations of Nine Target Genes

There has been no research on applying gene detection to differential diagnosis of adrenocortical carcinoma (ACC). We attempted to explore a novel auxiliary method for differential diagnosis between ACC with benign adrenocortical adenoma (ACA), based on mutations of target genes in tissues. Nine genes were chosen as target genes, including TP53, CTNNB1, ARMC5, PRKAR1A, ZNRF3, RB1, APC, MEN1, and RPL22. Exons sequencing of target genes were performed in 98 cases of tissue samples by FastTarget technology, including 41 ACC tissues, 32 ACA tissues, and 25 normal adrenal gland tissues. Significant mutations were detected and identified, and the clinical information was collected, for further comparative analysis and application to assist differential diagnosis of ACC. We identified 132 significant gene mutations and 227 significant mutation sites in 37 ACC tissues, much more than ACA and normal adrenal gland tissues. Mutation rates of 6 genes in ACC tissues were obviously higher than ACA tissues, including ZNRF3, ARMC5, TP53, APC, RB1, and PRKAR1A, regarded as high-risk genes. The sum of mutated high-risk genes detected in each sample was denominated sum of high-risk gene mutation (SHGM), and the rates of SHGM > 0 and SHGM > 1 in ACC tissues were 73.0% and 62.2%, respectively, both obviously higher than those in ACA tissues, with significant statistic differences. Especially for 8 cases of ACC with diameter < 5 cm, SHGM > 0 and SHGM > 1 were found in 6 samples (75%) and 4 samples (50%), respectively. However, no relevance was found between SHGM and clinical characteristics of ACC. We identified 6 high-risk genes in ACC tissues, with significantly higher mutation rates than ACA or normal adrenal gland tissues. The sum of mutated high-risk genes detected in ACC tissues was denominated SHGM, which was potential to assist the differential diagnosis of ACC with ACA, especially for the small-size ACC.

Genetic Variation and Evolution of the 2019 Novel Coronavirus

INTRODUCTION

SARS-CoV-2 is a new type of coronavirus causing a pandemic severe acute respiratory syndrome (SARS-2). Coronaviruses are very diverting genetically and mutate so often periodically. The natural selection of viral mutations may cause host infection selectivity and infectivity.

METHODS

This study was aimed to indicate the diversity between human and animal coronaviruses through finding the rate of mutation in each of the spike, nucleocapsid, envelope, and membrane proteins.

RESULTS

The mutation rate is abundant in all 4 structural proteins. The most number of statistically significant amino acid mutations were found in spike receptor-binding domain (RBD) which may be because it is responsible for a corresponding receptor binding in a broad range of hosts and host selectivity to infect. Among 17 previously known amino acids which are important for binding of spike to angiotensin-converting enzyme 2 (ACE2) receptor, all of them are conservative among human coronaviruses, but only 3 of them significantly are mutated in animal coronaviruses. A single amino acid aspartate-454, that causes dissociation of the RBD of the spike and ACE2, and F486 which gives the strength of binding with ACE2 remain intact in all coronaviruses.

DISCUSSION/CONCLUSION

Observations of this study provided evidence of the genetic diversity and rapid evolution of SARS-CoV-2 as well as other human and animal coronaviruses.

Competition for nutritional resources masks the true frequency of bacterial mutants

Background

It is widely assumed that all mutant microorganisms present in a culture are able to grow and form colonies, provided that they express the features required for selection. Unlike wild-type Escherichia coli, PHO-constitutive mutants overexpress alkaline phosphatase and hence can hydrolyze glycerol-2-phosphate (G2P) to glycerol and form colonies on plates having G2P as the sole carbon source. These mutations mostly occur in the pst operon. However, the frequency of PHO-constitutive colonies on the G2P selective plate is exceptionally low.

Results

We show that the rate in which spontaneous PHO-constitutive mutations emerge is about 8.0 × 10⁻⁶/generation, a relatively high rate, but the growth of most existing mutants is inhibited by their neighboring wild-type cells. This inhibition is elicited only by non-mutant viable bacteria that can take up and metabolize glycerol formed by the mutants. Evidence indicates that the few mutants that do form colonies derive from microclusters of mutants on the selective plate. A mathematical model that describes the fate of the wild-type and mutant populations under these circumstances supports these results.

Conclusion

This scenario in which neither the wild-type nor the majority of the mutants are able to grow resembles an unavoidable “tragedy of the commons” case which results in the collapse of the majority of the population. Cooperation between rare adjacent mutants enables them to overcome the competition and eventually form mutant colonies. The inhibition of PHO-constitutive mutants provides an example of mutant frequency masked by orders of magnitude due to a competition between mutants and their ancestral wild-type cells. Similar “tragedy of the commons-like” cases may occur in other settings and should be taken into consideration while estimating true mutant frequencies and mutation rates.

Supplementary Information

The online version contains supplementary material available at (doi:10.1186/s12915-020-00913-1).

Estimation of the SNP Mutation Rate in Two Vegetatively Propagating Species of Duckweed

Mutation rate estimates for vegetatively reproducing organisms are rare, despite their frequent occurrence across the tree of life. Here we report mutation rate estimates in two vegetatively reproducing duckweed species, Lemna minor and Spirodela polyrhiza We use a modified approach to estimating mutation rates by taking into account the reduction in mutation detection power that occurs when new individuals are produced from multiple cell lineages. We estimate an extremely low per generation mutation rate in both species of duckweed and note that allelic coverage at de novo mutation sites is very skewed. We also find no substantial difference in mutation rate between mutation accumulation lines propagated under benign conditions and those grown under salt stress. Finally, we discuss the implications of interpreting mutation rate estimates in vegetatively propagating organisms.

Mutations of SARS-CoV-2 nsp14 exhibit strong association with increased genome-wide mutation load

SARS-CoV-2 is a betacoronavirus responsible for COVID-19, a pandemic with global impact that first emerged in late 2019. Since then, the viral genome has shown considerable variance as the disease spread across the world, in part due to the zoonotic origins of the virus and the human host adaptation process. As a virus with an RNA genome that codes for its own genomic replication proteins, mutations in these proteins can significantly impact the variance rate of the genome, affecting both the survival and infection rate of the virus, and attempts at combating the disease. In this study, we analyzed the mutation densities of viral isolates carrying frequently observed mutations for four proteins in the RNA synthesis complex over time in comparison to wildtype isolates. Our observations suggest mutations in nsp14, an error-correcting exonuclease protein, have the strongest association with increased mutation load without selective pressure and across the genome, compared to nsp7, nsp8 and nsp12, which form the core polymerase complex. We propose nsp14 as a priority research target for understanding genomic variance rate in SARS-CoV-2 isolates and nsp14 mutations as potential predictors for high mutability strains.

A genome assembly and the somatic genetic and epigenetic mutation rate in a wild long-lived perennial Populus trichocarpa

BACKGROUND

Plants can transmit somatic mutations and epimutations to offspring, which in turn can affect fitness. Knowledge of the rate at which these variations arise is necessary to understand how plant development contributes to local adaption in an ecoevolutionary context, particularly in long-lived perennials.

RESULTS

Here, we generate a new high-quality reference genome from the oldest branch of a wild Populus trichocarpa tree with two dominant stems which have been evolving independently for 330 years. By sampling multiple, age-estimated branches of this tree, we use a multi-omics approach to quantify age-related somatic changes at the genetic, epigenetic, and transcriptional level. We show that the per-year somatic mutation and epimutation rates are lower than in annuals and that transcriptional variation is mainly independent of age divergence and cytosine methylation. Furthermore, a detailed analysis of the somatic epimutation spectrum indicates that transgenerationally heritable epimutations originate mainly from DNA methylation maintenance errors during mitotic rather than during meiotic cell divisions.

CONCLUSION

Taken together, our study provides unprecedented insights into the origin of nucleotide and functional variation in a long-lived perennial plant.

Time series prediction of COVID-19 by mutation rate analysis using recurrent neural network-based LSTM model

SARS-CoV-2, a novel coronavirus mostly known as COVID-19 has created a global pandemic. The world is now immobilized by this infectious RNA virus. As of June 15, already more than 7.9 million people have been infected and 432k people died. This RNA virus has the ability to do the mutation in the human body. Accurate determination of mutation rates is essential to comprehend the evolution of this virus and to determine the risk of emergent infectious disease. This study explores the mutation rate of the whole genomic sequence gathered from the patient's dataset of different countries. The collected dataset is processed to determine the nucleotide mutation and codon mutation separately. Furthermore, based on the size of the dataset, the determined mutation rate is categorized for four different regions: China, Australia, the United States, and the rest of the World. It has been found that a huge amount of Thymine (T) and Adenine (A) are mutated to other nucleotides for all regions, but codons are not frequently mutating like nucleotides. A recurrent neural network-based Long Short Term Memory (LSTM) model has been applied to predict the future mutation rate of this virus. The LSTM model gives Root Mean Square Error (RMSE) of 0.06 in testing and 0.04 in training, which is an optimized value. Using this train and testing process, the nucleotide mutation rate of 400^th patient in future time has been predicted. About 0.1% increment in mutation rate is found for mutating of nucleotides from T to C and G, C to G and G to T. While a decrement of 0.1% is seen for mutating of T to A, and A to C. It is found that this model can be used to predict day basis mutation rates if more patient data is available in updated time.

RdRp mutations are associated with SARS-CoV-2 genome evolution

COVID-19, caused by the novel SARS-CoV-2 virus, started in China in late 2019, and soon became a global pandemic. With the help of thousands of viral genome sequences that have been accumulating, it has become possible to track the evolution of the viral genome over time as it spread across the world. An important question that still needs to be answered is whether any of the common mutations affect the viral properties, and therefore the disease characteristics. Therefore, we sought to understand the effects of mutations in RNA-dependent RNA polymerase (RdRp), particularly the common 14408C>T mutation, on mutation rate and viral spread. By focusing on mutations in the slowly evolving M or E genes, we aimed to minimize the effects of selective pressure. Our results indicate that 14408C>T mutation increases the mutation rate, while the third-most common RdRp mutation, 15324C>T, has the opposite effect. It is possible that 14408C>T mutation may have contributed to the dominance of its co-mutations in Europe and elsewhere.

Counterbalancing the time-dependent effect on the human mitochondrial DNA molecular clock

BACKGROUND

The molecular clock is an important genetic tool for estimating evolutionary timescales. However, the detection of a time-dependent effect on substitution rate estimates complicates its application. It has been suggested that demographic processes could be the main cause of this confounding effect. In the present study, I propose a new algorithm for estimating the coalescent age of phylogenetically related sequences, taking into account the observed time-dependent effect on the molecular rate detected by others.

RESULTS

By applying this method to real human mitochondrial DNA trees with shallow and deep topologies, I obtained significantly older molecular ages for the main events of human evolution than were previously estimated. These ages are in close agreement with the most recent archaeological and paleontological records favoring the emergence of early anatomically modern humans in Africa 315 ± 34 thousand years ago (kya) and the presence of recent modern humans outside of Africa as early as 174 ± 48 thousand years ago. Furthermore, during the implementation process, I demonstrated that in a population with fluctuating sizes, the probability of fixation of a new neutral mutant depends on the effective population size, which is in better accordance with the fact that under the neutral theory of molecular evolution, the fate of a molecular mutation is mainly determined by random drift.

CONCLUSIONS

I suggest that the demographic history of populations has a more decisive effect than purifying selection and/or mutational saturation on the time-dependent effect observed for the substitution rate, and I propose a new method that corrects for this effect.

Highly prevalent BRAF V600E and low-frequency TERT promoter mutations underlie papillary thyroid carcinoma in Koreans

BACKGROUND

The presence of telomerase reverse transcriptase (TERT) promoter mutations have been associated with a poor prognosis in patients with papillary thyroid carcinomas (PTC). The frequency of TERT promoter mutations varies widely depending on the population and the nature of the study.

METHODS

Data were prospectively collected in 724 consecutive patients who underwent thyroidectomy for PTC from 2018 to 2019. Molecular testing for BRAF V600E and TERT promoter mutations was performed in all cases.

RESULTS

TERT promoter alterations in two hotspots (C228T and C250T) and C216T were found in 16 (2.2%) and 4 (0.6%) of all PTCs, respectively. The hotspot mutations were significantly associated with older age at diagnosis, larger tumor size, extrathyroidal extension, higher pathologic T category, lateral lymph node metastasis, and higher American Thyroid Association recurrence risk. The patients with C216T variant were younger and had a lower American Thyroid Association recurrence risk than those with hotspot mutations. Concurrent BRAF V600E was found in 19 of 20 cases with TERT promoter mutations. Of 518 microcarcinomas measuring ≤1.0 cm in size, hotspot mutations and C216T variants were detected in five (1.0%) and three (0.6%) cases, respectively.

CONCLUSIONS

Our study indicates a low frequency of TERT promoter mutations in Korean patients with PTC and supports previous findings that TERT promoter mutations are more common in older patients with unfavorable clinicopathologic features and BRAF V600E. TERT promoter mutations in patients with microcarcinoma are uncommon and may have a limited role in risk stratification. The C216T variant seems to have no clinicopathologic effect on PTC.

Analysis of HBsAg mutations in the 25 years after the implementation of the hepatitis B vaccination plan in China

Since 1992, China has promoted hepatitis B vaccination. Concurrently, during this period, increasing use of immunoglobulins and nucleoside analogues might have exerted selective pressure on the hepatitis B virus (HBV) S gene, driving mutations in the HBsAg and changed the subtype. Using the National Center for Biotechnology Information database, we obtained gene sequence information for HBV strains from China and analysed changes in HBsAg subtypes and substitution mutations in HBsAg in 5-year intervals over 25 years to identify potential challenges to the prevention and treatment of hepatitis B. Most HBV sequences from China were genotype C (1996/2833, 70.46%) or B (706/2833, 24.92%). During the implementation of hepatitis B vaccination (recombinant hepatitis B vaccine was subgenotype A2 and HBsAg subtype adw2), the proportion of subtypes ayw1 and adw3 in genotype B and ayw2 in genotype C increased over the programme period. The overall mutation rate in HBsAg tended to decrease for genotype B, whereas, for genotype C, the rate increased gradually and then decreased slightly. Moreover, the mutation rate at some HBsAg amino acid sites (such as sG145 of genotype B and sG130 and sK141 of genotype C) is gradually increasing. HBV strains with internal stop codons of HBsAg (e.g., sC69*) and additional N-glycosylation (e.g., sG130N) mutations should be studied extensively to prevent them from becoming dominant circulating strains. The development of HBV vaccines and antiviral immunoglobulins and use of antiviral drugs may require making corresponding changes.

The construction and application of a new 17-plex Y-STR system using universal fluorescent PCR

Y-chromosomal short tandem repeat (Y-STR) polymorphisms are useful in forensic identification, population genetics, and human structures. However, the current Y-STR systems are limited in discriminating distant relatives in a family with a low discrimination power. Increasing the capacity of detecting Y chromosomal polymorphisms will drastically narrow down the matching number of genealogy populations or pedigrees. In this study, we developed a system containing 17 Y-STRs that are complementary to the current commercially available Y-STR kits. This system was constructed by multiplex PCR with expected sizes of 126-400 bp labeled by different fluorescence molecules (DYS715, DYS709, DYS716, DYS713, and DYS607 labeled by FAM; DYS718, DYS723, DYS708, and DYS714 labeled by JOE; DYS712, DYS717, DYS721, and DYS605 labeled by TAMRA; and DYS719, DYS726, DYS598, and DYS722 labeled by ROX). The system was extensively tested for sensitivity, male specificity, species specificity, mixture, population genetics, and mutation rates following the Scientific Working Group on DNA Analysis Methods (SWGDAM) guidelines. The genetic data were obtained from eight populations with a total of 1260 individuals. Our results showed that all the 17 Y-STRs are human- and male-specific and include only one copy of the Y-chromosome. The 17 Y-STR system detects 143 alleles and has a high discrimination power (0.996031746). Mutation rates were different among the 17 Y-STRs, ranging from 0.30 to 3.03%. In conclusion, our study provides a robust, sensitive, and cost-effective genotyping method for human identification, which will be beneficial for narrowing the search scope when applied to genealogy searching with the Y-STR DNA databank.

Evaluating the contributions of purifying selection and progeny-skew in dictating within-host Mycobacterium tuberculosis evolution

The within-host evolutionary dynamics of tuberculosis (TB) remain unclear, and underlying biological characteristics render standard population genetic approaches based upon the Wright-Fisher model largely inappropriate. In addition, the compact genome combined with an absence of recombination is expected to result in strong purifying selection effects. Thus, it is imperative to establish a biologically relevant evolutionary framework incorporating these factors in order to enable an accurate study of this important human pathogen. Further, such a model is critical for inferring fundamental evolutionary parameters related to patient treatment, including mutation rates and the severity of infection bottlenecks. We here implement such a model and infer the underlying evolutionary parameters governing within-patient evolutionary dynamics. Results demonstrate that the progeny skew associated with the clonal nature of TB severely reduces genetic diversity and that the neglect of this parameter in previous studies has led to significant mis-inference of mutation rates. As such, our results suggest an underlying de novo mutation rate that is considerably faster than previously inferred, and a progeny distribution differing significantly from Wright-Fisher assumptions. This inference represents a more appropriate evolutionary null model, against which the periodic effects of positive selection, associated with drug-resistance for example, may be better assessed.

Diverse roles of nucleoside diphosphate kinase in genome stability and growth fitness

Nucleoside diphosphate kinase (NDK), a ubiquitous enzyme, catalyses reversible transfer of the γ phosphate from nucleoside triphosphates to nucleoside diphosphates and functions to maintain the pools of ribonucleotides and deoxyribonucleotides in the cell. As even a minor imbalance in the nucleotide pools can be mutagenic, NDK plays an antimutator role in maintaining genome integrity. However, the mechanism of the antimutator roles of NDK is not completely understood. In addition, NDKs play important roles in the host-pathogen interactions, metastasis, gene regulation, and various cellular metabolic processes. To add to these diverse roles of NDK in cells, a recent study now reveals that NDK may even confer mutator phenotypes to the cell by acting on the damaged deoxyribonucleoside diphosphates that may be formed during the oxidative stress. In this review, we discuss the roles of NDK in homeostasis of the nucleotide pools and genome integrity, and its possible implications in conferring growth/survival fitness to the organisms in the changing environmental niches.

Paternal and maternal mutations in X-STRs: A GHEP-ISFG collaborative study

The GHEP-ISFG organized a collaborative study to estimate mutation rates for the markers included in the Investigator Argus X-12 QS kit Qiagen. A total of 16 laboratories gathered data from 1,612 father/mother/daughter trios, which were used to estimate both maternal and paternal mutation rates, when pooled together with other already published data. Data on fathers and mothers' age at the time of birth of the daughter were also available for ∼93 % of the cases. Population analyses were computed considering the genetic information of a subset of 1,327 unrelated daughters, corresponding to 2,654 haplotypes from residents in several regions of five countries: Argentina, Brazil, Ecuador, Portugal and Spain. Genetic differentiation analyses between the population samples from the same country did not reveal signs of significant stratification, although results from Hardy-Weinberg and linkage disequilibrium tests indicated the need of larger studies for Ecuador and Brazilian populations. The high genetic diversity of the markers resulted in a large number of haplotype combinations, showing the need of huge databases for reliable estimates of their frequencies. It should also be noted the high number of new alleles found, many of them not included in the allelic ladders provided with the kit, as very diverse populations were analyzed. The overall estimates for locus specific average mutation rates varied between 7.5E-04 (for DXS7423) and 1.1E-02 (for DXS10135), the latter being a troublesome figure for kinship analyses. Most of the found mutations (∼92 %) are compatible with the gain or loss of a single repeat. Paternal mutation rates showed to be 5.2 times higher than maternal ones. We also found that older fathers were more prone to transmit mutated alleles, having this trend not been observed in the case of the mothers.

Focusing on DNA Repair and Damage Tolerance Mechanisms in Mycobacterium tuberculosis: An Emerging Therapeutic Theme

Tuberculosis (TB) is one such disease that has become a nuisance in the world scenario and one of the most deadly diseases of the current times. The etiological agent of tuberculosis, Mycobacterium tuberculosis (M. tb) kills millions of people each year. Not only 1.7 million people worldwide are estimated to harbor M. tb in the latent form but also 5 to 15 percent of which are expected to acquire an infection during a lifetime. Though curable, a long duration of drug regimen and expense leads to low patient adherence. The emergence of multi-, extensive- and total- drug-resistant strains of M. tb further complicates the situation. Owing to high TB burden, scientists worldwide are trying to design novel therapeutics to combat this disease. Therefore, to identify new drug targets, there is a growing interest in targeting DNA repair pathways to fight this infection. Thus, this review aims to explore DNA repair and damage tolerance as an efficient target for drug development by understanding M. tb DNA repair and tolerance machinery and its regulation, its role in pathogenesis and survival, mutagenesis, and consequently, in the development of drug resistance.

Epigenomic landscapes and their relationship to variation, fitness, and evolution

Doctor Tan's paper insightfully describes how epigenomic investigations of twins can help clarifying the mediation between environmental pressures, the genome, and behaviour. Here, I will expand on two emerging implications and their connections to the twin method. The first proposition deals with the transgenerational inheritance of epigenetic marks. The second examines the relationship between DNA methylation (DNAmeth) and mutation rate (MR), and their evolutionary implications.

Mutation rate and efflux response of bacteria exposed to a novel antimicrobial iodo-thiocyanate complex

OBJECTIVES

Antimicrobials, at sub-lethal concentrations, can act as selectors and promoters of resistance by increasing mutation rates. We measured the rate of Escherichia coli mutation from levofloxacin (LVX) sensitivity to resistance when it was grown under the near-lethal challenge of the novel biocidal iodo-thiocyanate complex (ITC). Another relevant factor affecting the emergence of antimicrobial resistance is the role of efflux pumps. Consequently, we evaluated whether ITC could potentially be a substrate for efflux pumps, and thus that efflux-mediated resistance could arise towards ITC.

METHODS

The mutation rate was measured by fluctuation analysis, when multiple parallel E. coli cultures were grown in the absence and presence of ITC. Then the mutational events, which occurred independently in each culture, were scored by plating the fraction of the culture in LVX-selective solid media and compared with the total cell number. To detect if ITC is a substrate for efflux pumps, minimum inhibitory concentrations (MICs) were determined against Pseudomonas aeruginosa in the absence and presence of the efflux pump inhibitor (EPI).

RESULTS

We have found that the E. coli exposed to the near-lethal level of ITC had a slight, but not significant, increase in mutation rate compared with unexposed cultures. Furthermore, the presence of EPI decreased the MIC of ITC by a modest 2-fold, showing that ITC was not a target for efflux pumps.

CONCLUSIONS

ITC usage most likely will not promote resistance development via increased mutation rates, and efflux-mediated resistance emergence to it is less likely than for some other antimicrobials.