Hearing silence: non-neutral evolution at synonymous sites in mammals

Expanding the CRISPR base editing toolbox in Drosophila melanogaster

CRISPR base editors can introduce point mutations into DNA precisely, and cytosine base editors (CBEs) catalyze C to T transitions. While CBEs have been thoroughly explored in cell culture and organisms such as mice, little is known about DNA base editing in insects. In this study, we evaluated germline editing rates of three different CBEs expressed under actin (ubiquitous) or nanos (germline) promoters utilizing Drosophila melanogaster. The original Rattus norvegicus-derived cytosine deaminase APOBEC1 (rAPO-1) displayed high base editing rates (~99%) with undetectable indel formation. Additionally, we show that base editors can be used for generating male sterility and female lethality. Overall, this study highlights the importance of promoter choice and sex-specific transmission for efficient base editing in flies while providing new insights for future genetic biocontrol designs in insects.

Combining Molecular, Macroevolutionary, and Macroecological Perspectives on the Generation of Diversity

Charles Darwin presented a unified process of diversification driven by the gradual accumulation of heritable variation. The growth in DNA databases and the increase in genomic sequencing, combined with advances in molecular phylogenetic analyses, gives us an opportunity to realize Darwin's vision, connecting the generation of variation to the diversification of lineages. The rate of molecular evolution is correlated with the rate of diversification across animals and plants, but the relationship between genome change and speciation is complex: Mutation rates evolve in response to life history and niche; substitution rates are influenced by mutation, selection, and population size; rates of acquisition of reproductive isolation vary between populations; and traits, niches, and distribution can influence diversification rates. The connection between mutation rate and diversification rate is one part of the complex and varied story of speciation, which has theoretical importance for understanding the generation of biodiversity and also practical impacts on the use of DNA to understand the dynamics of speciation over macroevolutionary timescales.

A novel de novo synonymous variant in GREB1L impacts the mRNA splicing associated with aplasia of the urogenital system

GREB1-like retinoic acid receptor coactivator (GREB1L) gene is associated with autosomal dominant renal hypodysplasia/aplasia 3 (RHDA3) and deafness, autosomal dominant 80 (DFNA80). Among the GREB1L variants reported, most of them are missense or frameshift, while no pathogenic synonymous variants have been recorded. Classical theory paid little attention to synonymous variants and classified it as nonpathogenic; however, recent studies suggest that the variants might be equally important. Here, we report a 7-year-old girl with new symptoms of clitoromegaly, uterovaginal, and ovarian agenesis as well as right kidney missing. A novel de novo GREB1L synonymous variant (NM_001142966: c.4731C>T, p.G1577=) was identified via whole exome sequencing. The variant was predicted to be disease-causing through in silico analysis and was classified as likely pathogenic. Minigene splicing assays confirmed a 6 bp deletion in mutant cDNA comparing with the wild type, leading to two amino acids lost in GREB1L protein. Secondary and tertiary structure modeling showed alterations in protein structure. Our finding reveals a novel GREB1L variant with a new phenotype of urogenital system and is the first to report a pathogenic synonymous variant in GREB1L which affects mRNA splicing, suggesting synonymous variants cannot be ignored in prenatal diagnosis and genetic counseling.

Translation Rates and Protein Folding

The mRNA coding sequence defines not only the amino acid sequence of the protein, but also the speed at which the ribosomes move along the mRNA while making the protein. The non-uniform local kinetics - denoted as translational rhythm - is similar among mRNAs coding for related protein folds. Deviations from this conserved rhythm can result in protein misfolding. In this review we summarize the experimental evidence demonstrating how local translation rates affect cotranslational protein folding, with the focus on the synonymous codons and patches of charged residues in the nascent peptide as best-studied examples. Alterations in nascent protein conformations due to disturbed translational rhythm can persist off the ribosome, as demonstrated by the effects of synonymous codon variants of several disease-related proteins. Charged amino acid patches in nascent chains also modulate translation and cotranslational protein folding, and can abrogate translation when placed at the N-terminus of the nascent peptide. During cotranslational folding, incomplete nascent chains navigate through a unique conformational landscape in which earlier intermediate states become inaccessible as the nascent peptide grows. Precisely tuned local translation rates, as well as interactions with the ribosome, guide the folding pathway towards the native structure, whereas deviations from the natural translation rhythm may favor pathways leading to trapped misfolded states. Deciphering the 'folding code' of the mRNA will contribute to understanding the diseases caused by protein misfolding and to rational protein design.

Prediction of the effects of the top 10 synonymous mutations from 26645 SARS-CoV-2 genomes

Background

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) had led to a global pandemic since December 2019. SARS-CoV-2 is a single-stranded RNA virus, which mutates at a higher rate. Multiple works had been done to study nonsynonymous mutations, which change protein sequences. However, there is little study on the effects of SARS-CoV-2 synonymous mutations, which may affect viral fitness. This study aims to predict the effect of synonymous mutations on the SARS-CoV-2 genome.

Methods

A total of 26645 SARS-CoV-2 genomic sequences retrieved from Global Initiative on Sharing all Influenza Data (GISAID) database were aligned using MAFFT. Then, the mutations and their respective frequency were identified. Multiple RNA secondary structures prediction tools, namely RNAfold, IPknot++ and MXfold2 were applied to predict the effect of the mutations on RNA secondary structure and their base pair probabilities was estimated using MutaRNA. Relative synonymous codon usage (RSCU) analysis was also performed to measure the codon usage bias (CUB) of SARS-CoV-2.

Results

A total of 150 synonymous mutations were identified. The synonymous mutation identified with the highest frequency is C3037U mutation in the nsp3 of ORF1a. Of these top 10 highest frequency synonymous mutations, C913U, C3037U, U16176C and C18877U mutants show pronounced changes between wild type and mutant in all 3 RNA secondary structure prediction tools, suggesting these mutations may have some biological impact on viral fitness. These four mutations show changes in base pair probabilities. All mutations except U16176C change the codon to a more preferred codon, which may result in higher translation efficiency.

Conclusion

Synonymous mutations in SARS-CoV-2 genome may affect RNA secondary structure, changing base pair probabilities and possibly resulting in a higher translation rate. However, lab experiments are required to validate the results obtained from prediction analysis.

Exploring the contribution of ARMS2 and HTRA1 genetic risk factors in age-related macular degeneration

Age-related macular degeneration (AMD) is the leading cause of severe irreversible central vision loss in individuals over 65 years old. Genome-wide association studies (GWASs) have shown that the region at chromosome 10q26, where the age-related maculopathy susceptibility (ARMS2/LOC387715) and HtrA serine peptidase 1 (HTRA1) genes are located, represents one of the strongest associated loci for AMD. However, the underlying biological mechanism of this genetic association has remained elusive. In this article, we extensively review the literature by us and others regarding the ARMS2/HTRA1 risk alleles and their functional significance. We also review the literature regarding the presumed function of the ARMS2 protein and the molecular processes of the HTRA1 protein in AMD pathogenesis in vitro and in vivo, including those of transgenic mice overexpressing HtrA1/HTRA1 which developed Bruch's membrane (BM) damage, choroidal neovascularization (CNV), and polypoidal choroidal vasculopathy (PCV), similar to human AMD patients. The elucidation of the molecular mechanisms of the ARMS2 and HTRA1 susceptibility loci has begun to untangle the complex biological pathways underlying AMD pathophysiology, pointing to new testable paradigms for treatment.

Genes for highly abundant proteins in Escherichia coli avoid 5' codons that promote ribosomal initiation

In many species highly expressed genes (HEGs) over-employ the synonymous codons that match the more abundant iso-acceptor tRNAs. Bacterial transgene codon randomization experiments report, however, that enrichment with such "translationally optimal" codons has little to no effect on the resultant protein level. By contrast, consistent with the view that ribosomal initiation is rate limiting, synonymous codon usage following the 5' ATG greatly influences protein levels, at least in part by modifying RNA stability. For the design of bacterial transgenes, for simple codon based in silico inference of protein levels and for understanding selection on synonymous mutations, it would be valuable to computationally determine initiation optimality (IO) scores for codons for any given species. One attractive approach is to characterize the 5' codon enrichment of HEGs compared with the most lowly expressed genes, just as translational optimality scores of codons have been similarly defined employing the full gene body. Here we determine the viability of this approach employing a unique opportunity: for Escherichia coli there is both the most extensive protein abundance data for native genes and a unique large-scale transgene codon randomization experiment enabling objective definition of the 5' codons that cause, rather than just correlate with, high protein abundance (that we equate with initiation optimality, broadly defined). Surprisingly, the 5' ends of native genes that specify highly abundant proteins avoid such initiation optimal codons. We find that this is probably owing to conflicting selection pressures particular to native HEGs, including selection favouring low initiation rates, this potentially enabling high efficiency of ribosomal usage and low noise. While the classical HEG enrichment approach does not work, rendering simple prediction of native protein abundance from 5' codon content futile, we report evidence that initiation optimality scores derived from the transgene experiment may hold relevance for in silico transgene design for a broad spectrum of bacteria.

Subfunctionalisation of paralogous genes and evolution of differential codon usage preferences: The showcase of polypyrimidine tract binding proteins

Gene paralogs are copies of an ancestral gene that appear after gene or full genome duplication. When two sister gene copies are maintained in the genome, redundancy may release certain evolutionary pressures, allowing one of them to access novel functions. Here, we focused our study on gene paralogs on the evolutionary history of the three polypyrimidine tract binding protein genes (PTBP) and their concurrent evolution of differential codon usage preferences (CUPrefs) in vertebrate species. PTBP1-3 show high identity at the amino acid level (up to 80%) but display strongly different nucleotide composition, divergent CUPrefs and, in humans and in many other vertebrates, distinct tissue-specific expression levels. Our phylogenetic inference results show that the duplication events leading to the three extant PTBP1-3 lineages predate the basal diversification within vertebrates, and genomic context analysis illustrates that local synteny has been well preserved over time for the three paralogs. We identify a distinct evolutionary pattern towards GC3-enriching substitutions in PTBP1, concurrent with enrichment in frequently used codons and with a tissue-wide expression. In contrast, PTBP2s are enriched in AT-ending, rare codons, and display tissue-restricted expression. As a result of this substitution trend, CUPrefs sharply differ between mammalian PTBP1s and the rest of PTBPs. Genomic context analysis suggests that GC3-rich nucleotide composition in PTBP1s is driven by local substitution processes, while the evidence in this direction is thinner for PTBP2-3. An actual lack of co-variation between the observed GC composition of PTBP2-3 and that of the surrounding non-coding genomic environment would raise an interrogation on the origin of CUPrefs, warranting further research on a putative tissue-specific translational selection. Finally, we communicate an intriguing trend for the use of the UUG-Leu codon, which matches the trends of AT-ending codons. Our results are compatible with a scenario in which a combination of directional mutation-selection processes would have differentially shaped CUPrefs of PTBPs in vertebrates: the observed GC-enrichment of PTBP1 in placental mammals may be linked to genomic location and to the strong and broad tissue-expression, while AT-enrichment of PTBP2 and PTBP3 would be associated with rare CUPrefs and thus, possibly to specialized spatio-temporal expression. Our interpretation is coherent with a gene subfunctionalisation process by differential expression regulation associated with the evolution of specific CUPrefs.

Novel COL4A3 synonymous mutation causes Alport syndrome coexistent with immunoglobulin A nephropathy in a woman: A case report

BACKGROUND

Alport syndrome (AS) is an inherited disease of the glomerular basement membrane caused by mutations in genes encoding α3, α4, or α5 chains of type IV collagen. It manifests with hematuria or proteinuria, which is often accompanied by hearing impairments and ocular abnormalities. Histopathologically, AS shows mesangial proliferation and sometimes incidental immunoglobulin A (IgA) deposition. Hematuria or proteinuria is also a common presentation in patients with IgA nephropathy that makes it difficult to differentially diagnose AS and IgA nephropathy solely based on these clinical and pathological features.

CASE SUMMARY

Herein, we present the case of a 59-year-old female patient who was admitted to our hospital with persistent microscopic hematuria and occasional proteinuria that had lasted for > 2 years. This patient had a familial history of renal disease and was diagnosed with autosomal dominant AS (ADAS) and IgA nephropathy based on the findings of renal biopsy as well as genetic testing performed using whole-exome sequencing, which suggested that the patient carried a novel heterozygous variation (c.888G>A:p.Gln296Gln) in the COL4A3 gene that enriches the mutation spectrum of ADAS. The proband received an angiotensin receptor blocker therapy after a definitive diagnosis was established. After one year of therapy, a significant reduction in proteinuria was observed. The number of microscopic red blood cells per high-power field decreased to one-quarter of the baseline levels. Renal function also maintained well during the follow-up.

CONCLUSION

Our case highlights the significance of performing kidney biopsy and genetic testing in the diagnosis of AS and familial IgA nephropathy.

Computational Quantification of Cancer Immunoediting

The remarkable success of cancer immunotherapy has revolutionized cancer treatment, emphasizing the importance of tumor-immune interactions in cancer evolution and treatment. Cancer immunoediting describes the dual effect of tumor-immune interactions: inhibiting tumor growth by destroying tumor cells and facilitating tumor escape by shaping tumor immunogenicity. To better understand tumor-immune interactions, it is critical to develop computational methods to measure the extent of cancer immunoediting. In this review, we provide a comprehensive overview of the computational methods for quantifying cancer immunoediting. We focus on describing the basic ideas, computational processes, advantages, limitations, and influential factors. We also summarize recent advances in quantifying cancer immunoediting studies and highlight future research directions. As the methods for quantifying cancer immunoediting are continuously improved, future research will further help define the role of immunity in tumorigenesis and hopefully provide a basis for the design of new personalized cancer immunotherapy strategies.

Dynamic Landscapes of tRNA Transcriptomes and Translatomes in Diverse Mouse Tissues

Although the function of tRNAs in the translational process is well established, it remains controversial whether tRNA abundance is tightly associated with translational efficiency (TE) in mammals. Moreover, how critically the expression of tRNAs contributes to the establishment of tissue-specific proteomes in mammals has not been well addressed. Here, we measured both tRNA expression using demethylase-tRNA sequencing (DM-tRNA-seq) and TE of mRNAs using ribosome-tagging sequencing (RiboTag-seq) in the brain, heart, and testis of mice. Remarkable variation in the expression of tRNA isodecoders was observed among different tissues. When the statistical effect of isodecoder-grouping on reducing variations is considered through permutating the anticodons, we observed an expected reduction in the variation of anticodon expression across all samples, an unexpected smaller variation of anticodon usage bias, and an unexpected larger variation of tRNA isotype expression at amino acid level. Regardless of whether or not they share the same anticodons, the isodecoders encoding the same amino acids are co-expressed across different tissues. Based on the expression of tRNAs and the TE of mRNAs, we find that the tRNA adaptation index (tAI) and TE are significantly correlated in the same tissues but not between tissues; and tRNA expression and the amino acid composition of translating peptides are positively correlated in the same tissues but not between tissues. We therefore hypothesize that the tissue-specific expression of tRNAs might be due to post-transcriptional mechanisms. This study provides a resource for tRNA and translation studies, as well as novel insights into the dynamics of tRNAs and their roles in translational regulation.

Transcriptome analysis identification of A-to-I RNA editing in granulosa cells associated with PCOS

Background

Polycystic ovary syndrome (PCOS) is a complex, multifactor disorder in women of reproductive age worldwide. Although RNA editing may contribute to a variety of diseases, its role in PCOS remains unclear.

Methods

A discovery RNA-Seq dataset was obtained from the NCBI Gene Expression Omnibus database of granulosa cells from women with PCOS and women without PCOS (controls). A validation RNA-Seq dataset downloaded from the European Nucleotide Archive Databank was used to validate differential editing. Transcriptome-wide investigation was conducted to analyze adenosine-to-inosine (A-to-I) RNA editing in PCOS and control samples.

Results

A total of 17,395 high-confidence A-to-I RNA editing sites were identified in 3,644 genes in all GC samples. As for differential RNA editing, there were 545 differential RNA editing (DRE) sites in 259 genes with Nucleoporin 43 (NUP43), Retinoblastoma Binding Protein 4 (RBBP4), and leckstrin homology-like domain family A member 1 (PHLDA) showing the most significant three 3'-untranslated region (3'UTR) editing. Furthermore, we identified 20 DRE sites that demonstrated a significant correlation between editing levels and gene expression levels. Notably, MIR193b-365a Host Gene (MIR193BHG) and Hook Microtubule Tethering Protein 3 (HOOK3) exhibited significant differential expression between PCOS and controls. Functional enrichment analysis showed that these 259 differentially edited genes were mainly related to apoptosis and necroptosis pathways. RNA binding protein (RBP) analysis revealed that RNA Binding Motif Protein 45 (RBM45) was predicted as the most frequent RBP binding with RNA editing sites. Additionally, we observed a correlation between editing levels of differential editing sites and the expression level of the RNA editing enzyme Adenosine Deaminase RNA Specific B1 (ADARB1). Moreover, the existence of 55 common differentially edited genes and nine differential editing sites were confirmed in the validation dataset.

Conclusion

Our current study highlighted the potential role of RNA editing in the pathophysiology of PCOS as an epigenetic process. These findings could provide valuable insights into the development of more targeted and effective treatment options for PCOS.

gtAI: an improved species-specific tRNA adaptation index using the genetic algorithm

The tRNA adaptation index (tAI) is a translation efficiency metric that considers weighted values (S _ij values) for codon-tRNA wobble interaction efficiencies. The initial implementation of the tAI had significant flaws. For instance, generated S _ij weights were optimized based on gene expression in Saccharomyces cerevisiae, which is expected to vary among different species. Consequently, a species-specific approach (stAI) was developed to overcome those limitations. However, the stAI method employed a hill climbing algorithm to optimize the S _ij weights, which is not ideal for obtaining the best set of S _ij weights because it could struggle to find the global maximum given a complex search space, even after using different starting positions. In addition, it did not perform well in computing the tAI of fungal genomes in comparison with the original implementation. We developed a novel approach named genetic tAI (gtAI) implemented as a Python package (https://github.com/AliYoussef96/gtAI), which employs a genetic algorithm to obtain the best set of S _ij weights and follows a new codon usage-based workflow that better computes the tAI of genomes from the three domains of life. The gtAI has significantly improved the correlation with the codon adaptation index (CAI) and the prediction of protein abundance (empirical data) compared to the stAI.

PredDSMC: A predictor for driver synonymous mutations in human cancers

Introduction: Driver mutations play a critical role in the occurrence and development of human cancers. Most studies have focused on missense mutations that function as drivers in cancer. However, accumulating experimental evidence indicates that synonymous mutations can also act as driver mutations.

Methods: Here, we proposed a computational method called PredDSMC to accurately predict driver synonymous mutations in human cancers. We first systematically explored four categories of multimodal features, including sequence features, splicing features, conservation scores, and functional scores. Further feature selection was carried out to remove redundant features and improve the model performance. Finally, we utilized the random forest classifier to build PredDSMC.

Results: The results of two independent test sets indicated that PredDSMC outperformed the state-of-the-art methods in differentiating driver synonymous mutations from passenger mutations.

Discussion: In conclusion, we expect that PredDSMC, as a driver synonymous mutation prediction method, will be a valuable method for gaining a deeper understanding of synonymous mutations in human cancers.

Gene-by-Sex Interactions: Genome-Wide Association Study Reveals Five SNPs Associated with Obesity and Overweight in a Male Population

Obesity is a chronic health problem associated with severe complications and with an increasing prevalence in the Western world. Body-fat composition and distribution are closely associated with obesity, but the human body’s composition is a sexually dimorphic trait, as differences between the two sexes are evident even from fetal life. The effect of sex hormones contributes to this phenomenon. However, studies investigating gene-by-sex interactions for obesity are limited. Therefore, the aim of the present study was to identify single-nucleotide polymorphisms (SNPs) associated with obesity and overweight in a male population. A genome-wide association study (GWAS) that included 104 control, 125 overweight, and 61 obese subjects revealed four SNPs associated with overweight (rs7818910, rs7863750, rs1554116, and rs7500401) and one SNP (rs114252547) associated with obesity in males. An in silico functional annotation was subsequently used to further investigate their role. Most of the SNPs were found in genes regulating energy metabolism and homeostasis, and some of them were expression quantitative trait loci (eQTL). These findings contribute to the understanding of the molecular mechanisms underlying obesity-related traits, especially in males, and pave the road for future research toward the improvement of the diagnosis and therapy of obese individuals.

Association of SNPs in the PAI1 Gene with Disease Recurrence and Clinical Outcome in Bladder Cancer

PURPOSE

Bladder cancer (BCa) is one of the most common cancer types worldwide and is characterized by a high rate of recurrence. In previous studies, we and others have described the functional influence of plasminogen activator inhibitor-1 (PAI1) in bladder cancer development. While polymorphisms in PAI1 have been associated with increased risk and worsened prognosis in some cancers, the mutational status of PAI1 in human bladder tumors has not been well defined.

METHODS

In this study, we evaluated the mutational status of PAI1 in a series of independent cohorts, comprised of a total of 660 subjects.

RESULTS

Sequencing analyses identified two clinically relevant 3' untranslated region (UTR) single nucleotide polymorphisms (SNPs) in PAI1 (rs7242; rs1050813). Somatic SNP rs7242 was present in human BCa cohorts (overall incidence of 72%; 62% in Caucasians and 72% in Asians). In contrast, the overall incidence of germline SNP rs1050813 was 18% (39% in Caucasians and 6% in Asians). Furthermore, Caucasian patients with at least one of the described SNPs had worse recurrence-free survival and overall survival (p = 0.03 and p = 0.03, respectively). In vitro functional studies demonstrated that SNP rs7242 increased the anti-apoptotic effect of PAI1, and SNP rs1050813 was related to a loss of contact inhibition associated with cellular proliferation when compared to wild type.

CONCLUSION

Further investigation of the prevalence and potential downstream influence of these SNPs in bladder cancer is warranted.

Effects of synonymous mutations on kinetic properties and structure of firefly luciferase: Molecular dynamics simulation, molecular docking, RNA folding, and experimental study

Although synonymous mutations have long been thought to lack striking results, a growing body of research shows these mutations have highly variable effects. In this study, the impact of synonymous mutations in the development of thermostable luciferase was investigated using a combination of experimental and theoretical approaches. Using bioinformatics analysis, the codon usage features in the Lampyridae family's luciferases were studied and four synonymous mutations of Arg in luciferase were created. An exciting result was that the analysis of kinetic parameters showed a slight increase in the thermal stability of the mutant luciferase. AutoDock Vina, %MinMax algorithm, and UNAFold Server were used to perform molecular docking, folding rate, and RNA folding, respectively. Here, it was assumed that in the region (Arg337) with a moderate propensity for coil, synonymous mutation altered the rate of translation, which in turn may lead to a slight change in the structure of the enzyme. According to the molecular dynamics simulation data, local minor global flexibility is observed in the context of the protein conformation. A plausible explanation is that this flexibility may strengthen hydrophobic interactions due to its sensitivity to a molecular collision. Accordingly, thermostability originated mainly from hydrophobic interaction.

Transcriptomic, proteomic, and functional consequences of codon usage bias in human cells during heterologous gene expression

Differences in codon frequency between genomes, genes, or positions along a gene, modulate transcription and translation efficiency, leading to phenotypic and functional differences. Here, we present a multiscale analysis of the effects of synonymous codon recoding during heterologous gene expression in human cells, quantifying the phenotypic consequences of codon usage bias at different molecular and cellular levels, with an emphasis on translation elongation. Six synonymous versions of an antibiotic resistance gene were generated, fused to a fluorescent reporter, and independently expressed in HEK293 cells. Multiscale phenotype was analyzed by means of quantitative transcriptome and proteome assessment, as proxies for gene expression; cellular fluorescence, as a proxy for single-cell level expression; and real-time cell proliferation in absence or presence of antibiotic, as a proxy for the cell fitness. We show that differences in codon usage bias strongly impact the molecular and cellular phenotype: (i) they result in large differences in mRNA levels and protein levels, leading to differences of over 15 times in translation efficiency; (ii) they introduce unpredicted splicing events; (iii) they lead to reproducible phenotypic heterogeneity; and (iv) they lead to a trade-off between the benefit of antibiotic resistance and the burden of heterologous expression. In human cells in culture, codon usage bias modulates gene expression by modifying mRNA availability and suitability for translation, leading to differences in protein levels and eventually eliciting functional phenotypic changes.

Seasonal Adaptation: Geographic Photoperiod-Temperature Patterns Explain Genetic Variation in the Common Vole Tsh Receptor

The vertebrate photoperiodic neuroendocrine system uses the photoperiod as a proxy to time the annual rhythms in reproduction. The thyrotropin receptor (TSHR) is a key protein in the mammalian seasonal reproduction pathway. Its abundance and function can tune sensitivity to the photoperiod. To investigate seasonal adaptation in mammals, the hinge region and the first part of the transmembrane domain of the Tshr gene were sequenced for 278 common vole (Microtus arvalis) specimens from 15 localities in Western Europe and 28 localities in Eastern Europe. Forty-nine single nucleotide polymorphisms (SNPs; twenty-two intronic and twenty-seven exonic) were found, with a weak or lack of correlation with pairwise geographical distance, latitude, longitude, and altitude. By applying a temperature threshold to the local photoperiod-temperature ellipsoid, we obtained a predicted critical photoperiod (pCPP) as a proxy for the spring onset of local primary food production (grass). The obtained pCPP explains the distribution of the genetic variation in Tshr in Western Europe through highly significant correlations with five intronic and seven exonic SNPs. The relationship between pCPP and SNPs was lacking in Eastern Europe. Thus, Tshr, which plays a pivotal role in the sensitivity of the mammalian photoperiodic neuroendocrine system, was targeted by natural selection in Western European vole populations, resulting in the optimized timing of seasonal reproduction.

Do Noncoding and Coding Sites in Angiosperm Chloroplast DNA Have Different Mutation Processes?

Fourfold degenerate sites within coding regions and intergenic sites have both been used as estimates of neutral evolution. In chloroplast DNA, the pattern of substitution at intergenic sites is strongly dependent on the composition of the surrounding hexanucleotide composed of the three base pairs on each side, which suggests that the mutation process is highly context-dependent in this genome. This study examines the context-dependency of substitutions at fourfold degenerate sites in protein-coding regions and compares the pattern to what has been observed at intergenic sites. Overall, there is strong similarity between the two types of sites, but there are some intriguing differences. One of these is that substitutions of G and C are significantly higher at fourfold degenerate sites across a range of contexts. In fact, A → T and T → A substitutions are the only substitution types that occur at a lower rate at fourfold degenerate sites. The data are not consistent with selective constraints being responsible for the difference in substitution patterns between intergenic and fourfold degenerate sites. Rather, it is suggested that the difference may be a result of different epigenetic modifications that result in slightly different mutation patterns in coding and intergenic DNA.

Hemizygosity can reveal variant pathogenicity on the X-chromosome

Pathogenic variants on the X-chromosome can have more severe consequences for hemizygous males, while heterozygote females can avoid severe consequences due to diploidy and the capacity for nonrandom expression. Thus, when an allele is more common in females this could indicate that it increases the probability of early death in the male hemizygous state, which can be considered a measure of pathogenicity. Importantly, large-scale genomic data now makes it possible to compare allele proportions between the sexes. To discover pathogenic variants on the X-chromosome, we analyzed exome data from 125,748 ancestrally diverse participants in the Genome Aggregation Database (gnomAD). After filtering out duplicates and extremely rare variants, 44,606 of the original 348,221 remained for analysis. We divided the proportion of variant alleles in females by the proportion in males for all variant sites, and then placed each variant into one of three a priori categories: (1) Reference (Primarily synonymous and intronic), (2) Unlikely-to-be-tolerated (Primarily missense), and (3) Least-likely-to-be-tolerated (Primarily frameshift). To assess the impact of ploidy, we compared the distribution of these ratios between pseudoautosomal and non-pseudoautosomal regions. In the non-pseudoautosomal regions, mean female-to-male ratios were lowest among Reference (2.40), greater for Unlikely-to-be-tolerated (2.77) and highest for Least-likely-to-be-tolerated (3.28) variants. Corresponding ratios were lower in the pseudoautosomal regions (1.52, 1.57, and 1.68, respectively), with the most extreme ratio being just below 11. Because pathogenic effects in the pseudoautosomal regions should not drive ratio increases, this maximum ratio provides an upper bound for baseline noise. In the non-pseudoautosomal regions, 319 variants had a ratio over 11. In sum, we identified a measure with a dataset specific threshold for identifying pathogenicity in non-pseudoautosomal X-chromosome variants: the female-to-male allele proportion ratio.

Identification and analysis of putative tRNA genes in baculovirus genomes

Transfer RNAs (tRNAs) genes are both coded for and arranged along some viral genomes representing the entire virosphere and seem to play different biological functions during infection, other than transferring the correct amino acid to a growing peptide chain. Baculovirus genome description and annotation has focused mostly on protein-coding genes, microRNA, and homologous regions. Here we carried out a large-scale in silico search for putative tRNA genes in baculovirus genomes. Ninety-six of 257 baculovirus genomes analyzed was found to contain at least one putative tRNA gene. We found great diversity in primary and secondary structure, in location within the genome, in intron presence and size, and in anti-codon identity. In some cases, genes of tRNA-containing genomes were found to have a bias for the codons specified by the tRNAs present in such genomes. Moreover, analysis revealed that most of the putative tRNA genes possessed conserved motifs for tRNA type 2 promoters, including the A-box and B-box motifs with few mismatches from the eukaryotic canonical motifs. From publicly available small RNA deep sequencing datasets of baculovirus-infected insect cells, we found evidence that a putative Autographa californica multiple nucleopolyhedrovirus Gln-tRNA gene was transcribed and modified with the addition of the non-templated 3'-CCA tail found at the end of all tRNAs. Further research is needed to determine the expression and functionality of these viral tRNAs.

Useful genes for predicting the efficacy of transarterial chemoembolization in hepatocellular carcinoma

Transarterial chemoembolization (TACE) is generally used to treat patients with hepatocellular carcinoma (HCC), a common and deadly cancer; however, its efficacy varies according to factors such as tumor volume, stage, serum alpha-fetoprotein level, and chosen feeding artery. In addition, gene-related factors have been recently suggested to be involved in the regulation and prediction of TACE outcomes. Accordingly, genes could serve as effective biomarkers to select patients who can benefit from TACE. These gene-related factors can activate signaling pathways affecting cancer cell survival while regulating the epithelial-mesenchymal transition, angiogenesis, and the tumor microenvironment, all directly associated with tumor progression, thereby affecting TACE efficacy. Moreover, this disordered gene expression is associated with poor prognosis in patients with HCC, including TACE resistance, postoperative recurrence, and metastasis. To identify the exact relationship between various genes and TACE efficacy, this review summarizes the involvement of protein-coding and non-coding genes and single nucleotide polymorphisms in TACE efficacy for predicting the efficacy of TACE; the present findings may help improve the efficacy of TACE in clinical settings.

dN/dS-H, a New Test to Distinguish Different Selection Modes in Protein Evolution and Cancer Evolution

One of the most popular measures in the analysis of protein sequence evolution is the ratio of nonsynonymous distance (d_N) to synonymous distance (d_S). Under the assumption that synonymous substitutions in the coding region are selectively neutral, the d_N/d_S ratio can be used to statistically detect the adaptive evolution (or purifying selection) if d_N/d_S > 1 (or d_N/d_S < 1) significantly. However, due to strong structural constraints and/or variable functional constraints imposed on amino acid sites, most encoding genes in most species have demonstrated d_N/d_S < 1. Consequently, the statistical power for testing d_N/d_S = 1 may be insufficient to distinguish between different selection modes. In this paper, we propose a more powerful test, called d_N/d_S-H, in which a new parameter H, a relative measure of rate variation among sites, was introduced. Given the condition of strong purifying selections at some sites, the d_N/d_S-H model predicts d_N/d_S = 1-H for neutral evolution, d_N/d_S < 1-H for nearly neutral selection, and d_N/d_S > 1-H for adaptive evolution. The potential of this new method for resolving the neutral-adaptive debates is illustrated by the protein sequence evolution in vertebrates, Drosophila and yeasts, as well as somatic cancer evolution (specialized as the C_N/C_S-H test).

Sustainable solutions for indoor pollution abatement during COVID phase: A critical study on current technologies & challenges

The appearance of the contagious virus COVID-19, several revelations and environmental health experts punctually predicted the possibly disastrous public health complications of coexisting catching and airborne contamination-arbitrated disease. But much attention has been given on the outdoor-mediated interactions. Almost 3.8 million premature deaths occur every year globally due to the illness from indoor air pollution. Considering the human staying longer span indoors due to restricted human activities or work from home, the indoor air quality (IAQ) might show prominent role for individual health life. Currently, the Environmental Protection Agency (EPA) ensures no regulation of indoor airborne pollution. Herein, the paper underlines the common bases of indoor air pollution, poor IAQ, and impacts of the aerosolized airborne particles on the human health. In order to address these challenges and collective contagion events in indoor environment, several emerging control techniques and preventive sustainable solutions are suggested. By this, more innovations need to be investigated in future to measure the impact of indoor air pollution on individual health.

Assessment of phylogenetic approaches to study the timing of recombination cessation on sex chromosomes

The evolution of sex chromosomes is hypothesized to be punctuated by consecutive recombination cessation events, forming "evolutionary strata" that ceased to recombine at different time points. The demarcation of evolutionary strata is often assessed by estimates of the timing of recombination cessation (t_RC ) along the sex chromosomes, commonly inferred from the level of synonymous divergence or with species phylogenies at gametologous (X-Y or Z-W) sequence data. However, drift and selection affect sequences unpredictably and introduce uncertainty when inferring t_RC . Here, we assess two alternative phylogenetic approaches to estimate t_RC ; (i) the expected likelihood weight (ELW) approach that finds the most likely topology among a set of hypothetical topologies and (ii) the BEAST approach that estimates t_RC with specified calibration priors on a reference species topology. By using Z and W gametologs of an old and a young evolutionary stratum on the neo-sex chromosome of Sylvioidea songbirds, we show that the ELW and BEAST approaches yield similar t_RC estimates, and that both outperform two frequently applied approaches utilizing synonymous substitution rates (dS) and maximum likelihood (ML) trees, respectively. Moreover, we demonstrate that both ELW and BEAST provide more precise t_RC estimates when sequences of multiple species are included in the analyses.

Stop codon usage as a window into genome evolution: mutation, selection, biased gene conversion and the TAG paradox

Protein coding genes terminate with one of three stop codons (TAA, TGA, or TAG) that, like synonymous codons, are not employed equally. With TGA and TAG having identical nucleotide content, analysis of their differential usage provides an unusual window into the forces operating on what are ostensibly functionally identical residues. Across genomes and between isochores within the human genome, TGA usage increases with G + C content but, with a common G + C → A + T mutation bias, this cannot be explained by mutation bias-drift equilibrium. Increased usage of TGA in G + C-rich genomes or genomic regions is also unlikely to reflect selection for the optimal stop codon, as TAA appears to be universally optimal, probably because it has the lowest read-through rate. Despite TAA being favored by selection and mutation bias, as with codon usage bias G + C pressure is the prime determinant of between-species TGA usage trends. In species with strong G + C-biased gene conversion (gBGC), such as mammals and birds, the high usage and conservation of TGA is best explained by an A + T → G + C repair bias. How to explain TGA enrichment in other G + C-rich genomes is less clear. Enigmatically, across bacterial and archaeal species and between human isochores TAG usage is mostly unresponsive to G + C pressure. This unresponsiveness we dub the TAG paradox as currently no mutational, selective, or gBGC model provides a well-supported explanation. That TAG does increase with G + C usage across eukaryotes makes the usage elsewhere yet more enigmatic. We suggest resolution of the TAG paradox may provide insights into either an unknown but common selective preference (probably at the DNA/RNA level) or an unrecognized complexity to the action of gBGC.

Preferred synonymous codons are translated more accurately: Proteomic evidence, among-species variation, and mechanistic basis

A commonly stated cause of unequal uses of synonymous codons is their differential translational accuracies. A key prediction of this long-standing translational accuracy hypothesis (TAH) of codon usage bias is higher translational accuracies of more frequently used synonymous codons, which, however, has had no direct evidence beyond case studies. Analyzing proteomic data from Escherichia coli, we present direct, global evidence for this prediction. The experimentally measured codon-specific translational accuracies validate a sequence-based proxy; this proxy provides support for the TAH from the vast majority of over 1000 taxa surveyed in all domains of life. We find that the relative translational accuracies of synonymous codons vary substantially among taxa and are strongly correlated with the amounts of cognate transfer RNAs (tRNAs) relative to those of near-cognate tRNAs. These and other observations suggest a model in which selections for translational efficiency and accuracy drive codon usage bias and its coevolution with the tRNA pool.

Molecular Diversity and Phylogeny Reconstruction of Genus Colobanthus (Caryophyllaceae) Based on Mitochondrial Gene Sequences

Mitochondrial genomes have become an interesting object of evolutionary and systematic study both for animals and plants, including angiosperms. Although the framework of the angiosperm phylogeny was built on the information derived from chloroplast and nuclear genes, mitochondrial sequences also revealed their usefulness in solving the phylogenetic issues at different levels of plant systematics. Here, we report for the first time the complete sequences of 26 protein-coding genes of eight Colobanthus species (Caryophyllaceae). Of these, 23 of them represented core mitochondrial genes, which are directly associated with the primary function of that organelle, and the remaining three genes represented a facultative set of mitochondrial genes. Comparative analysis of the identified genes revealed a generally high degree of sequence conservation. The Ka/Ks ratio was <1 for most of the genes, which indicated purifying selection. Only for rps12 was Ka/Ks > 1 in all studied species, suggesting positive selection. We identified 146−165 potential RNA editing sites in genes of the studied species, which is lower than in most angiosperms. The reconstructed phylogeny based on mitochondrial genes was consistent with the taxonomic position of the studied species, showing the separate character of the family Caryophyllaceae and close relationships between all studied Colobanthus species, with C. lycopodioides sharing less similarity.

The Ramp Atlas: facilitating tissue and cell-specific ramp sequence analyses through an intuitive web interface

Ramp sequences occur when the average translational efficiency of codons near the 5′ end of highly expressed genes is significantly lower than the rest of the gene sequence, which counterintuitively increases translational efficiency by decreasing downstream ribosomal collisions. Here, we show that the relative codon adaptiveness within different tissues changes the existence of a ramp sequence without altering the underlying genetic code. We present the first comprehensive analysis of tissue and cell type-specific ramp sequences and report 3108 genes with ramp sequences that change between tissues and cell types, which corresponds with increased gene expression within those tissues and cells. The Ramp Atlas (https://ramps.byu.edu/) allows researchers to query precomputed ramp sequences in 18 388 genes across 62 tissues and 66 cell types and calculate tissue-specific ramp sequences from user-uploaded FASTA files through an intuitive web interface. We used The Ramp Atlas to identify seven SARS-CoV-2 genes and seven human SARS-CoV-2 entry factor genes with tissue-specific ramp sequences that may help explain viral proliferation within those tissues. We anticipate that The Ramp Atlas will facilitate personalized and creative tissue-specific ramp sequence analyses for both human and viral genes that will increase our ability to utilize this often-overlooked regulatory region.

Genomic insights into positive selection during barley domestication

BACKGROUND

Cultivated barley (Hordeum vulgare) is widely used in animal feed, beverages, and foods and has become a model crop for molecular evolutionary studies. Few studies have examined the evolutionary fates of different types of genes in barley during the domestication process.

RESULTS

The rates of nonsynonymous substitution (Ka) to synonymous substitution (Ks) were calculated by comparing orthologous genes in different barley groups (wild vs. landrace and landrace vs. improved cultivar). The rates of evolution, properties, expression patterns, and diversity of positively selected genes (PSGs) and negatively selected genes (NSGs) were compared. PSGs evolved more rapidly, possessed fewer exons, and had lower GC content than NSGs; they were also shorter and had shorter intron, exon, and first exon lengths. Expression levels were lower, the tissue specificity of expression was higher, and codon usage bias was weaker for PSGs than for NSGs. Nucleotide diversity analysis revealed that PSGs have undergone a more severe genetic bottleneck than NSGs. Several candidate PSGs were involved in plant growth and development, which might make them as excellent targets for the molecular breeding of barley.

CONCLUSIONS

Our comprehensive analysis of the evolutionary, structural, and functional divergence between PSGs and NSGs in barley provides new insight into the evolutionary trajectory of barley during domestication. Our findings also aid future functional studies of PSGs in barley.

Synonymous mutations in representative yeast genes are mostly strongly non-neutral

Synonymous mutations in protein-coding genes do not alter protein sequences and are thus generally presumed to be neutral or nearly neutral^1-5. Here, to experimentally verify this presumption, we constructed 8,341 yeast mutants each carrying a synonymous, nonsynonymous or nonsense mutation in one of 21 endogenous genes with diverse functions and expression levels and measured their fitness relative to the wild type in a rich medium. Three-quarters of synonymous mutations resulted in a significant reduction in fitness, and the distribution of fitness effects was overall similar-albeit nonidentical-between synonymous and nonsynonymous mutations. Both synonymous and nonsynonymous mutations frequently disturbed the level of mRNA expression of the mutated gene, and the extent of the disturbance partially predicted the fitness effect. Investigations in additional environments revealed greater across-environment fitness variations for nonsynonymous mutants than for synonymous mutants despite their similar fitness distributions in each environment, suggesting that a smaller proportion of nonsynonymous mutants than synonymous mutants are always non-deleterious in a changing environment to permit fixation, potentially explaining the common observation of substantially lower nonsynonymous than synonymous substitution rates. The strong non-neutrality of most synonymous mutations, if it holds true for other genes and in other organisms, would require re-examination of numerous biological conclusions about mutation, selection, effective population size, divergence time and disease mechanisms that rely on the assumption that synoymous mutations are neutral.

Rates of Mutations and Transcript Errors in the Foodborne Pathogen Salmonella enterica subsp. enterica

Because errors at the DNA level power pathogen evolution, a systematic understanding of the rate and molecular spectra of mutations could guide the avoidance and treatment of infectious diseases. We thus accumulated tens of thousands of spontaneous mutations in 768 repeatedly bottlenecked lineages of 18 strains from various geographical sites, temporal spread, and genetic backgrounds. Entailing over ∼1.36 million generations, the resultant data yield an average mutation rate of ∼0.0005 per genome per generation, with a significant within-species variation. This is one of the lowest bacterial mutation rates reported, giving direct support for a high genome stability in this pathogen resulting from high DNA-mismatch-repair efficiency and replication-machinery fidelity. Pathogenicity genes do not exhibit an accelerated mutation rate, and thus, elevated mutation rates may not be the major determinant for the diversification of toxin and secretion systems. Intriguingly, a low error rate at the transcript level is not observed, suggesting distinct fidelity of the replication and transcription machinery. This study urges more attention on the most basic evolutionary processes of even the best-known human pathogens and deepens the understanding of their genome evolution.

WT1, NR0B1, NR5A1, LHX9, ZFP92, ZNF275, INSL3, and NRIP1 Genetic Variants in Patients with Premature Ovarian Insufficiency in a Mexican Cohort

Premature ovarian insufficiency (POI) is one of the main causes of female premature infertility. POI is a genetically heterogeneous disorder with a complex molecular etiology; as such, the genetic causes remain unknown in the majority of patients. Therefore, this study aimed to identify mutations and characterize the associated molecular contribution of gonadogenesis-determinant genes to POI. Genomic assays, including PCR-SSCP and Sanger sequencing, followed by in silico analyses were used to investigate the underpinnings of ovarian deficiency in 11 women affected by POI. Large deletions and nucleotide insertions and duplications were excluded by PCR. Thirteen genetic variants were identified in the WT1 (c.213G>T, c.609T>C, c.873A>G, c.1122G>A), NR0B1 (c.353C>T, c.425G>A), NR5A1 (c.437G>C, IVS4-20C>T), LHX9 (IVS2-12G>C, IVS3+13C>T, c.741T>C), ZNF275 (c.969C>T), and NRIP1 (c.3403C>T) genes. Seven novel genetic variants and five unpublished substitutions were identified. No genetic aberrations were detected in the ZFP92 and INSL3 genes. Each variant was genotyped using PCR-SSCP in 100 POI-free subjects, and their allelic frequencies were similar to the patients. These analyses indicated that allelic variation in the WT1, NR0B1, NR5A1, LHX9, ZFP92, ZNF275, INSL3, and NRIP1 genes may be a non-disease-causing change or may not contribute significantly to the genetics underlying POI disorders. Findings support the polygenic nature of this clinical disorder, with the SNVs identified representing only a probable contribution to the variability of the human genome.

Next-generation sequencing of the whole mitochondrial genome identifies functionally deleterious mutations in patients with multiple sclerosis

Multiple sclerosis (MS) is an immune-mediated disease of the central nervous system with genetics and environmental determinants. Studies focused on the neurogenetics of MS showed that mitochondrial DNA (mtDNA) mutations that can ultimately lead to mitochondrial dysfunction, alter brain energy metabolism and cause neurodegeneration. We analyzed the whole mitochondrial genome using next-generation sequencing (NGS) from 47 Saudi individuals, 23 patients with relapsing-remitting MS and 24 healthy controls to identify mtDNA disease-related mutations/variants. A large number of variants were detected in the D-loop and coding genes of mtDNA. While distinct unique variants were only present in patients or only occur in controls, a number of common variants were shared among the two groups. The prevalence of some common variants differed significantly between patients and controls, thus could be implicated in susceptibility to MS. Of the unique variants only present in the patients, 34 were missense mutations, located in different mtDNA-encoded genes. Seven of these mutations were not previously reported in MS, and predicted to be deleterious with considerable impacts on the functions and structures of encoded-proteins and may play a role in the pathogenesis of MS. These include two heteroplasmic mutations namely 10237T>C in MT-ND3 gene and 15884G>C in MT-CYB gene; and three homoplasmic mutations namely 9288A>G in MT-CO3 gene, 14484T>C in MT-ND6 gene, 15431G>A in MT-CYB gene, 8490T>C in MT-ATP8 gene and 5437C>T in MT-ND2 gene. Notably some patients harboured multiple mutations while other patients carried the same mutations. This study is the first to sequence the entire mitochondrial genome in MS patients in an Arab population. Our results expanded the mutational spectrum of mtDNA variants in MS and highlighted the efficiency of NGS in population-specific mtDNA variant discovery. Further investigations in a larger cohort are warranted to confirm the role of mtDNA MS.

Elevated mutation rates underlie the evolution of the aquatic plant family Podostemaceae

Molecular evolutionary rates vary among lineages and influence the evolutionary process. Here, we report elevated genome-wide mutation rates in Podostemaceae, a family of aquatic plants with a unique body plan that allows members to live on submerged rocks in fast-flowing rivers. Molecular evolutionary analyses using 1640 orthologous gene groups revealed two historical increases in evolutionary rates: the first at the emergence of the family and the second at the emergence of Podostemoideae, which is the most diversified subfamily. In both branches, synonymous substitution rates were elevated, indicating higher mutation rates. On early branches, mutations were biased in favour of AT content, which is consistent with a role for ultraviolet light-induced mutation and habitat shift. In ancestors of Podostemoideae, DNA-repair genes were enriched in genes under positive selection, which may have responded to the meristem architectural changes.

Case Report: Identification and Functional Analysis of a Homozygous Synonymous Variant in the PLOD1 Gene in a Chinese Neonatal With the Ehlers-Danlos Syndrome

BACKGROUND

Kyphoscoliotic Ehlers-Danlos syndrome (kEDS; OMIM225400) is a rare autosomal recessive genetic disease caused by variants in the PLOD1 gene. This research was conducted to verify the disease-causing gene in a Chinese neonatal family with the EDS.

METHODS

We recruited a Han Chinese neonate with PLOD1-related kEDS without kyphoscoliosis. Detailed clinical examination and laboratory tests were performed and whole exome sequencing (WES) was used to detect the pathogenic genes of the proband. In vivo experiments (reverse-transcription PCR, quantitative real-time PCR) and in vitro experiments (minigene analysis) were used to verify the function of variants suspected of affecting the splicing process. The effect of the splice site variant on the PLOD1 transcript was analyzed using splice prediction programs NetGene2 and Alternative Splice Site Predictor (ASSP).

RESULTS

A homozygous synonymous variant c.1095C>T (p.Gly365, rs1032781250) in the PLOD1 gene was found and verified in the family with kEDS. This splicing variant resulted in a premature termination codon of exon 10 and affected the expression of the four bases GCGC.

CONCLUSION

Our research showed that the homozygous synonymous variant in PLOD1 was the pathogenic cause in the proband. The combined application of WES and functional studies verified the effect of uncertain gene variants on splicing, upgrading pathogenicity evidence, and determining the cause of disease. This is helpful for the early diagnosis and treatment of kEDS.

OUP accepted manuscript

Recognition of the important role of transposable elements (TEs) in eukaryotic genomes quickly led to a burgeoning literature modeling and estimating the effects of selection on TEs. Much of the empirical work on selection has focused on analyzing the site frequency spectrum (SFS) of TEs. But TE evolution differs from standard models in a number of ways that can impact the power and interpretation of the SFS. For example, rather than mutating under a clock-like model, transposition often occurs in bursts which can inflate particular frequency categories compared with expectations under a standard neutral model. If a TE burst has been recent, the excess of low-frequency polymorphisms can mimic the effect of purifying selection. Here, we investigate how transposition bursts affect the frequency distribution of TEs and the correlation between age and allele frequency. Using information on the TE age distribution, we propose an age-adjusted SFS to compare TEs and neutral polymorphisms to more effectively evaluate whether TEs are under selective constraints. We show that our approach can minimize instances of false inference of selective constraint, remains robust to simple demographic changes, and allows for a correct identification of even weak selection affecting TEs which experienced a transposition burst. The results presented here will help researchers working on TEs to more reliably identify the effects of selection on TEs without having to rely on the assumption of a constant transposition rate.

Vulnerable exons, like ACADM exon 5, are highly dependent on maintaining a correct balance between splicing enhancers and silencers

It is now widely accepted that aberrant splicing of constitutive exons is often caused by mutations affecting cis-acting splicing regulatory elements, but there is a misconception that all exons have an equal dependency on splicing regulatory elements and thus a similar susceptibility to aberrant splicing. We investigated exonic mutations in ACADM exon 5 to experimentally examine their effect on splicing and found that 7 out of 11 tested mutations affected exon inclusion, demonstrating that this constitutive exon is particularly vulnerable to exonic splicing mutations. Employing ACADM exon 5 and 6 as models, we demonstrate that the balance between splicing enhancers and silencers, flanking intron length, and flanking splice site strength are important factors that determine exon definition and splicing efficiency of the exon in question. Our study shows that two constitutive exons in ACADM have different inherent vulnerabilities to exonic splicing mutations. This suggests that in silico prediction of potential pathogenic effects on splicing from exonic mutations may be improved by also considering the inherent vulnerability of the exon. Moreover, we show that single nucleotide polymorphism that affect either of two different exonic splicing silencers, located far apart in exon 5, all protect against both immediately flanking and more distant exonic splicing mutations.

Decoding the effects of synonymous variants

Synonymous single nucleotide variants (sSNVs) are common in the human genome but are often overlooked. However, sSNVs can have significant biological impact and may lead to disease. Existing computational methods for evaluating the effect of sSNVs suffer from the lack of gold-standard training/evaluation data and exhibit over-reliance on sequence conservation signals. We developed synVep (synonymous Variant effect predictor), a machine learning-based method that overcomes both of these limitations. Our training data was a combination of variants reported by gnomAD (observed) and those unreported, but possible in the human genome (generated). We used positive-unlabeled learning to purify the generated variant set of any likely unobservable variants. We then trained two sequential extreme gradient boosting models to identify subsets of the remaining variants putatively enriched and depleted in effect. Our method attained 90% precision/recall on a previously unseen set of variants. Furthermore, although synVep does not explicitly use conservation, its scores correlated with evolutionary distances between orthologs in cross-species variation analysis. synVep was also able to differentiate pathogenic vs. benign variants, as well as splice-site disrupting variants (SDV) vs. non-SDVs. Thus, synVep provides an important improvement in annotation of sSNVs, allowing users to focus on variants that most likely harbor effects.

Selection on Accessible Chromatin Regions in Capsella grandiflora

Accurate estimates of genome-wide rates and fitness effects of new mutations are essential for an improved understanding of molecular evolutionary processes. Although eukaryotic genomes generally contain a large noncoding fraction, functional noncoding regions and fitness effects of mutations in such regions are still incompletely characterized. A promising approach to characterize functional noncoding regions relies on identifying accessible chromatin regions (ACRs) tightly associated with regulatory DNA. Here, we applied this approach to identify and estimate selection on ACRs in Capsella grandiflora, a crucifer species ideal for population genomic quantification of selection due to its favorable population demography. We describe a population-wide ACR distribution based on ATAC-seq data for leaf samples of 16 individuals from a natural population. We use population genomic methods to estimate fitness effects and proportions of positively selected fixations (α) in ACRs and find that intergenic ACRs harbor a considerable fraction of weakly deleterious new mutations, as well as a significantly higher proportion of strongly deleterious mutations than comparable inaccessible intergenic regions. ACRs are enriched for expression quantitative trait loci (eQTL) and depleted of transposable element insertions, as expected if intergenic ACRs are under selection because they harbor regulatory regions. By integrating empirical identification of intergenic ACRs with analyses of eQTL and population genomic analyses of selection, we demonstrate that intergenic regulatory regions are an important source of nearly neutral mutations. These results improve our understanding of selection on noncoding regions and the role of nearly neutral mutations for evolutionary processes in outcrossing Brassicaceae species.

A synonymous mutation in IGF-1 impacts the transcription and translation process of gene expression

Insulin-like growth factor 1 (IGF-1) is considered to be a crucial gene in the animal development of bone and body size. In this study, a unique synonymous mutation (c.258 A > G) of the IGF-1 gene was modified with an adenine base editor to observe the growth and developmental situation of mutant mice. Significant expression differences and molecular mechanisms among vectors with different alanine synonymous codons were explored. Although modification of a single synonymous codon rarely interferes with animal phenotypes, we observed that the expression and secretion of IGF-1 were different between 8-week-old homozygous (Ho) and wild-type (WT) mice. In addition, the IGF-1 with optimal codon combinations showed a higher expression content than other codon combination modes at both transcription and translation levels and performed proliferation promotion. The gene stability and translation initiation efficiency also changed significantly. Our findings illustrated that the synonymous mutation altered the IGF-1 gene expression in individual mice and suggested that the synonymous mutation affected the IGF-1 expression and biological function through the transcription and translation processes.

Prediction of the Effects of Synonymous Variants on SARS-CoV-2 Genome

Background: The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) had led to a global pandemic since December 2019. SARS-CoV-2 is a single-stranded RNA virus, which mutates at a higher rate. Multiple studies had been done to identify and study nonsynonymous mutations, which change amino acid residues of SARS-CoV-2 proteins. On the other hand, there is little study on the effects of SARS-CoV-2 synonymous mutations. Although these mutations do not alter amino acids, some studies suggest that they may affect viral fitness. This study aims to predict the effect of synonymous mutations on the SARS-CoV-2 genome.   Methods: A total of 30,229 SARS-CoV-2 genomic sequences were retrieved from Global Initiative on Sharing all Influenza Data (GISAID) database and aligned using MAFFT. Then, the mutations and their respective frequency were identified. A prediction of RNA secondary structures and their base pair probabilities was performed to study the effect of synonymous mutations on RNA structure and stability. Relative synonymous codon usage (RSCU) analysis was also performed to measure the codon usage bias (CUB) of SARS-CoV-2.  Results: A total of 150 synonymous mutations were identified. The synonymous mutation identified with the highest frequency is C3037U mutation in the nsp3 of ORF1a, followed by C313U and C9286U mutation in nsp1 and nsp4 of ORF1a, respectively.   Conclusion: Among the synonymous mutations identified, C913U mutation in ORF1a and C26735U in membrane (M) protein may affect RNA secondary structure, reducing the stability of RNA folding and possibly resulting in a higher translation rate. However, lab experiments are required to validate the results obtained from prediction analysis.

Frequent lineage-specific substitution rate changes support an episodic model for protein evolution

Since the inception of the molecular clock model for sequence evolution, the investigation of protein divergence has revolved around the question of a more or less constant change of amino acid sequences, with specific overall rates for each family. Although anomalies in clock-like divergence are well known, the assumption of a constant decay rate for a given protein family is usually taken as the null model for protein evolution. However, systematic tests of this null model at a genome-wide scale have lagged behind, despite the databases’ enormous growth. We focus here on divergence rate comparisons between very closely related lineages since this allows clear orthology assignments by synteny and reliable alignments, which are crucial for determining substitution rate changes. We generated a high-confidence dataset of syntenic orthologs from four ape species, including humans. We find that despite the appearance of an overall clock-like substitution pattern, several hundred protein families show lineage-specific acceleration and deceleration in divergence rates, or combinations of both in different lineages. Hence, our analysis uncovers a rather dynamic history of substitution rate changes, even between these closely related lineages, implying that one should expect that a large fraction of proteins will have had a history of episodic rate changes in deeper phylogenies. Furthermore, each of the lineages has a separate set of particularly fast diverging proteins. The genes with the highest percentage of branch-specific substitutions are ADCYAP1 in the human lineage (9.7%), CALU in chimpanzees (7.1%), SLC39A14 in the internal branch leading to humans and chimpanzees (4.1%), RNF128 in gorillas (9%), and S100Z in gibbons (15.2%). The mutational pattern in ADCYAP1 suggests a biased mutation process, possibly through asymmetric gene conversion effects. We conclude that a null model of constant change can be problematic for predicting the evolutionary trajectories of individual proteins.

Variability in Codon Usage in Coronaviruses Is Mainly Driven by Mutational Bias and Selective Constraints on CpG Dinucleotide

The Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the third human-emerged virus of the 21st century from the Coronaviridae family, causing the ongoing coronavirus disease 2019 (COVID-19) pandemic. Due to the high zoonotic potential of coronaviruses, it is critical to unravel their evolutionary history of host species breadth, host-switch potential, adaptation and emergence, to identify viruses posing a pandemic risk in humans. We present here a comprehensive analysis of the composition and codon usage bias of the 82 Orthocoronavirinae members, infecting 47 different avian and mammalian hosts. Our results clearly establish that synonymous codon usage varies widely among viruses, is only weakly dependent on their primary host, and is dominated by mutational bias towards AU-enrichment and by CpG avoidance. Indeed, variation in GC3 explains around 34%, while variation in CpG frequency explains around 14% of total variation in codon usage bias. Further insight on the mutational equilibrium within Orthocoronavirinae revealed that most coronavirus genomes are close to their neutral equilibrium, the exception being the three recently infecting human coronaviruses, which lie further away from the mutational equilibrium than their endemic human coronavirus counterparts. Finally, our results suggest that, while replicating in humans, SARS-CoV-2 is slowly becoming AU-richer, likely until attaining a new mutational equilibrium.

The rate and molecular spectrum of mutation are selectively maintained in yeast

What determines the rate (μ) and molecular spectrum of mutation is a fundamental question. The prevailing hypothesis asserts that natural selection against deleterious mutations has pushed μ to the minimum achievable in the presence of genetic drift, or the drift barrier. Here we show that, contrasting this hypothesis, μ substantially exceeds the drift barrier in diverse organisms. Random mutation accumulation (MA) in yeast frequently reduces μ, and deleting the newly discovered mutator gene PSP2 nearly halves μ. These results, along with a comparison between the MA and natural yeast strains, demonstrate that μ is maintained above the drift barrier by stabilizing selection. Similar comparisons show that the mutation spectrum such as the universal AT mutational bias is not intrinsic but has been selectively preserved. These findings blur the separation of mutation from selection as distinct evolutionary forces but open the door to alleviating mutagenesis in various organisms by genome editing.

How natural selection shapes the rate and molecular spectrum of mutations is debated. Yeast mutation accumulation experiments identify a gene promoting mutagenesis and show stabilizing selection maintaining the mutation rate above the drift barrier. Selection also preserves the mutation spectrum.

Evolutionary and phylogenetic insights from a nuclear genome sequence of the extinct, giant, "subfossil" koala lemur Megaladapis edwardsi

No endemic Madagascar animal with body mass >10 kg survived a relatively recent wave of extinction on the island. From morphological and isotopic analyses of skeletal "subfossil" remains we can reconstruct some of the biology and behavioral ecology of giant lemurs (primates; up to ∼160 kg) and other extraordinary Malagasy megafauna that survived into the past millennium. Yet, much about the evolutionary biology of these now-extinct species remains unknown, along with persistent phylogenetic uncertainty in some cases. Thankfully, despite the challenges of DNA preservation in tropical and subtropical environments, technical advances have enabled the recovery of ancient DNA from some Malagasy subfossil specimens. Here, we present a nuclear genome sequence (∼2× coverage) for one of the largest extinct lemurs, the koala lemur Megaladapis edwardsi (∼85 kg). To support the testing of key phylogenetic and evolutionary hypotheses, we also generated high-coverage nuclear genomes for two extant lemurs, Eulemur rufifrons and Lepilemur mustelinus, and we aligned these sequences with previously published genomes for three other extant lemurs and 47 nonlemur vertebrates. Our phylogenetic results confirm that Megaladapis is most closely related to the extant Lemuridae (typified in our analysis by E. rufifrons) to the exclusion of L. mustelinus, which contradicts morphology-based phylogenies. Our evolutionary analyses identified significant convergent evolution between M. edwardsi and an extant folivore (a colobine monkey) and an herbivore (horse) in genes encoding proteins that function in plant toxin biodegradation and nutrient absorption. These results suggest that koala lemurs were highly adapted to a leaf-based diet, which may also explain their convergent craniodental morphology with the small-bodied folivore Lepilemur.