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

KRT14 is a promising prognostic biomarker of breast cancer related to immune infiltration

BACKGROUND

Breast cancer (BC) is the leading cancer among women globally, which has the highest incidence and mortality rate in over a hundred countries. This study was intended to discover a new prognostic biomarker, facilitating personalized treatment approaches.

METHODS

RNA sequencing data from The Cancer Genome Atlas database and Gene Expression Omnibus database were utilized to download to evaluate expression levels and prognostic significance of Keratin 14 (KRT14). Methylation of KRT14 was also assessed. The CIBERSORT and single-sample gene set enrichment analysis algorithms were applied to explore the connection between KRT14 and the tumor microenvironment. Primary drugs' sensitivity and potential small molecule therapeutic compounds were analyzed through the "pRRophetic" R package and the Connectivity Map. The prognostic value of KRT14 was additionally corroborated through a comparison of protein levels in peritumoral and cancerous tissues via immunohistochemistry. Moreover, an immune-related prognostic model based on KRT14 was designed to enhance the prediction accuracy for the prognosis of BC patients.

RESULTS

The study found that KRT14 expression was generally downregulated in BC, correlating strongly with poor prognosis. Compared to normal tissues, the methylation level of KRT14 was higher in BC tissues. Lower expression of KRT14 was linked to decreased anti-tumoral immune cells infiltration and increased immunosuppressive cells infiltration. Sensitivity to various key therapeutic drugs was lower in groups with diminished KRT14 expression. In addition, several potential anti-BC small molecule compounds were identified. The model designed in this study significantly enhanced the predictive capability for BC patients compared to predictions based solely on KRT14 expression levels.

CONCLUSION

Overall, KRT14 was closely correlated with the prognosis in BC, making it a reliable biomarker.

A Novel Homozygous Synonymous Variant in CCDC134 as a Cause of Osteogenesis Imperfecta Type XXII

Osteogenesis imperfecta (OI) is a heterogeneous group of rare, inherited connective tissue disorders. It includes over 20 defined subtypes, each of which is associated with distinct causative genes that are listed in the Online Mendelian Inheritance in Man (OMIM) database. Type XXII OI (OI 22) is caused by a homozygous variant in the coiled-coil domain containing 134 (CCDC134) gene, which is located on chromosome 22q13. OI, which is associated with CCDC134, is extremely rare with only five cases reported worldwide. All known cases involve the c.2 T > C (p. Met1Thr) homozygous missense variant in the CCDC134 gene. We present the case of a 13-year-old Chinese girl with non-union fracture, short stature and specific radiographic findings, which include scoliosis, pelvic tilt, thin clavicles, ribs, and limbs. Whole exome sequencing revealed a novel, homozygous c.492G > C (p. Leu164=) variation in the CCDC134 gene. RNA sequencing (RNA-seq) analysis identified this variant as an abnormal splicing variant that causes the deletion of Exon 5, which result in the observed disease phenotype. This case demonstrates the clinical phenotype of OI 22 associated with the c.492G > C (p. Leu164=) novel synonymous variation in the coding region of the CCDC134 gene in a female patient. This is the first reported case of OI 22 in the Chinese population, the sixth reported worldwide and the fourth reported genotype for diseases associated with a CCDC134 variant. It also enriches the global clinical phenotype spectrum of OI 22 patients.

Engineering of Peptide-Inserted Base Editors with Enhanced Accuracy and Security

Base editors are effective tools for introducing base conversions without double-strand breaks, showing broad applications in biotechnological and clinical areas. However, their non-negligible bystander mutations and off-target effects have raised extensive safety concerns. To address these issues, a novel method is developed by inserting specific peptide fragments into the substrate binding pocket of deaminases in base editors to modify these outcomes. It is validated that the composition and position of the inserted peptide can significantly impact the performance of A3A-based cytosine base editor and TadA-8e-based adenine base editor, leading to improved editing activity and precision in human HEK293T cells. Importantly, the TadA-8e variant with DPLVLRRRQ peptide inserted behind S116 residue showed a strong motif preference of Y4A5N6, which can accurately edit the A5 base in targeted protospacer with minimized bystander and off-target effects in DNA and RNA-level. By summarizing the regularity during engineering, a set of systematic procedures is established, which can potentially be used to modify other types of base editors and make them more accurate and secure. In addition, the peptide insertion strategy is also proven to be compatible with traditional amino acid changes which have been reported, exhibiting excellent compatibility.

Genome-wide identification and comparative evolution of 14-3-3 gene family members in five Brassicaceae species

BACKGROUND

The 14-3-3 proteins are highly conserved regulatory eukaryotic proteins, which are crucial in growth, development, and stress responses. However, systematic characterization of the 14-3-3 gene family in Brassicaceae species and their evolutionary relationships have not been comprehensively reported.

RESULTS

This study conducted genome-wide identification, structural characteristics, and comparative evolutionary analysis of 14-3-3 gene family members in Arabidopsis thaliana, A. lyrata, A. pumila, Camelina sativa, and Brassica oleracea using comparative genomics. Overall, a total of 108 14-3-3 genes, which were phylogenetically classified into ε and non-ε groups were identified in the five species, with the non-ε members exhibiting more similar exon-intron structures and conserved motif patterns. Collinearity analysis revealed that the Brassicaceae 14-3-3 gene family members underwent varying degrees of expansion following whole-genome duplication (WGD) events. Notably, the number of 14-3-3 gene family members between A. lyrata and A. thaliana remained similar despite the former having approximately 1.66-fold larger genome size. In contrast, the number of 14-3-3 gene family members in A. pumila and C. sativa increased in proportionately to their genome size, while gene members in the more distantly related species to A. thaliana, B. oleracea, showed irregular expansion patterns. Selection pressure analysis revealed that 14-3-3 homologs in all the five species underwent purifying selection, with the group ε members experiencing relatively weaker purifying selection. Cloning of ApGRF6-2 gene from A. pumila indicated that the ApGRF6-2 protein was localized in the cell membrane and cytoplasm, while ectopic overexpression of ApGRF6-2 in A. thaliana could promote early flowering by upregulating the expression of floral meristem identity genes.

CONCLUSION

This study provides a comprehensive and systematic identification of the 14-3-3 gene family members in five Brassicaceae species using updated genome sequences, and the results could form a basis for further validation of functional and molecular mechanisms of 14-3-3 genes in plant growth, development, abiotic stress responses, as well as flowering regulation.

Role of ERG3 mutation and expression in azole resistant Candida albicans isolated from vulvovaginal candidiasis patients

Objective

The study aimed to investigate mutations and mRNA expression of the ERG3 gene in resistant Candida albicans isolates.

Methods

This cross-sectional study was conducted from October 2018 to June 2019. High vaginal swab samples were collected from Hayatabad Medical Complex and transported to Khyber Medical University. Samples were inoculated on different media and identified by 20C AUX strips. Antifungal susceptibility was determined using the disc diffusion and broth microdilution methods. The ERG3 gene was amplified and sequenced to find amino acid polymorphisms. Real-time PCR was performed to study level of ERG3 expression.

Results

A total of seventy-three (n=73) Candida albicans out of 369 samples were isolated. Among the isolates 49.3%, 54.8%, 53.4%, 47.9%, 30.1% were resistant to fluconazole, Clotrimazole, Miconazole, Voriconazole and Itraconazole, respectively. Sanger sequencing of ERG3 gene of isolates revealed six synonymous mutations. Expression level of mRNA of ERG3 gene in azoles sensitive stains (3.72±2.22) was higher than those in the resistant Candida albican strains (1.74±0.96).

Conclusion

This study revealed synonymous mutations and low expression of ERG3 gene in azole-resistant C. albican.

Codon Usage Bias: A Potential Factor Affecting VGLUT Developmental Expression and Protein Evolution

More and more attention has been paid to the role of synonymous substitution in evolution, in which codon usage preference can affect gene expression distribution and protein structure and function. Vesicular glutamate transporter (VGLUT) consists of three isoforms, among which VGLUT3 is significantly different from other VGLUTs in functional importance, expression level, and distribution range, whose reason is still unclear. This study sought to analyze the role of codon preference in VGLUT differentiation. To conduct an evolutionary analysis of the three VGLUTs, this paper uses bioinformatics research methods to analyze the coding sequences of the three VGLUTs in different species and compare the codon usage patterns. Furthermore, the differences among the three VGLUTs were analyzed by combining functional importance, expression level, distribution range, gene structure, protein relationship network, expression at specific developmental stages, and phylogenetic tree, and the influence of codon usage pattern was explored. The results showed that the VGLUT with greater codon preference had less functional importance, lower expression levels, more peripheral distribution away from the CNS, smaller exon density of gene, less conserved and farther away from the CDS region miRNA regulatory sites, simpler and less tight protein interaction networks, delayed developmental expression, and more distant evolutionary relationships. Codon usage preference is a potential factor affecting VGLUT developmental expression and protein evolution.

Capsule mutations serve as a key strategy of phage resistance evolution of K54 hypervirulent Klebsiella pneumoniae

Phage therapy is a promising antibacterial strategy against the antibiotic resistance crisis. The evolved phage resistance could pose a big challenge to clinical phage therapy. Therefore, it is necessary to conduct a comprehensive analysis of phage resistance mechanisms during treatment. Here, we characterize 37 phage-resistant mutants of hypervirulent K. pneumoniae strain SCNJ1 under phage-imposed selection in both in vitro and in vivo experiments. We show that 97.3% (36/37) of phage-resistant clones possessed at least one mutation in genes related to the CPS biosynthesis. Notably, the wcaJ gene emerges as a mutation hotspot, as mutations in this gene are detected at a high frequency under both conditions. In contrast, mutations in wzc exhibit more association with in vivo samples. These CPS-related mutants all exhibit compromised bacterial fitness and attenuated virulence in mice. Strain CM8 is the only non-CPS-related mutant, which has a bglA mutation that confers phage resistance and retains full fitness and virulence. This study highlights that laboratory characterization of phage resistance evolution can give useful insights for clinical phage therapy.

Role of silent mutations in KRAS -mutant tumors

ABSTRACT

Silent mutations within the RAS  gene have garnered increasing attention for their potential roles in tumorigenesis and therapeutic strategies. Kirsten-RAS ( KRAS ) mutations, predominantly oncogenic, are pivotal drivers in various cancers. While extensive research has elucidated the molecular mechanisms and biological consequences of active KRAS  mutations, the functional significance of silent mutations remains relatively understudied. This review synthesizes current knowledge on KRAS  silent mutations, highlighting their impact on cancer development. Silent mutations, which do not alter protein sequences but can affect RNA stability and translational efficiency, pose intriguing questions regarding their contribution to tumor biology. Understanding these mutations is crucial for comprehensively unraveling KRAS -driven oncogenesis and exploring novel therapeutic avenues. Moreover, investigations into the clinical implications of silent mutations in KRAS -mutant tumors suggest potential diagnostic and therapeutic strategies. Despite being in early stages, research on KRAS  silent mutations holds promise for uncovering novel insights that could inform personalized cancer treatments. In conclusion, this review underscores the evolving landscape of KRAS  silent mutations, advocating for further exploration to bridge fundamental biology with clinical applications in oncology.

Ovalbumin gene polymorphism: Implications for hatchability and egg quality changes during storage in Japanese quail

The aim of the study was to identify polymorphisms in the ovalbumin gene - SERPINB14 gene and evaluate their effect on hatchability traits and egg quality changes during storage in two strains of Japanese quails: meat-type (F33) and laying-type (S22). To individually determine hatchability traits for each female, eggs were collected and incubated. To determine egg quality traits, 10 eggs were collected from each female and stored for 14 weeks. Egg quality was analyzed 10 times during storage. All exons and the 3'UTR of the SERPINB14 gene were sequenced. A total of 17 SNPs were identified in both strains: 4 in exons, 5 in the 3'UTR, and 8 in intron regions. Association analysis showed significant effects of SNP14 and SNP16 on the percentage of late died embryos. Fresh egg weight in F33 females was influenced by eight SNPs: SNP6, SNP7, SNP9, SNP11, SNP14, SNP15, SNP16, and SNP17, with significant diplotype effects observed. Individuals with H3H3 and H7H7 diplotypes showed the highest egg weight. SNPs 6, 7, and 11 influenced eggshell thickness on the laying day and at 2, 4, and 14 weeks of storage. The effects of haplotypes on this trait were also observed. Significant SNP effects were also found on albumen weight and albumen percentage at different storage times. Moreover, diplotypes from block 1 influenced albumen traits during storage. These studies provide new information on the SERPINB14 gene polymorphism in Japanese quail, and some of the markers merit further validation as useful tools for selection to improve hatchability and egg quality in poultry breeding programs.

Control of Botrytis cinerea from Chilean Grapevines by Pydiflumetofen: Baseline and Carboxamide-Mutant Sensitivity

The gray mold (Botrytis cinerea; Botrytis) is the main disease affecting grapevine production in Chile. Succinate dehydrogenase inhibitors belonging to the carboxamide fungicide family are a key tool for the control of Botrytis in grapevines from the Chilean Central Valley. This study aimed to determine the sensitivity of the Chilean Botrytis population to the new generation carboxamide pydiflumetofen. Conidial germination (CG) and germ-tube elongation (GTE) sensitivity assays were conducted on 200 single-spore isolates collected during the 2016 to 2017 season. The mean effective concentration that inhibited 50% (EC50) of CG in the Botrytis population was 0.0545 μg/ml, with mean values of 0.066 and 0.042 μg/ml for table and wine grapes, respectively. The mean EC50 value of GTE was 0.000245, 0.0003, and 0.0019 μg/ml for the total, table grape, and wine grape populations, respectively. The comparison between pydiflumetofen and fludioxonil, a highly efficient fungicide carrying a different mode of action, showed 87.5 and 97.5% of Botrytis control with an EC50 threshold of 0.1 μg/ml in table grape and wine grape populations, respectively. No cross-resistance between pydiflumetofen and fludioxonil was detected. For nine isolates with reduced pydiflumetofen sensitivity, we evaluated SdhB mutations using a qPCR-HRM diagnostic system. Two isolates carried the sdhBP225/H272R genotype, and two isolates carried the sdhBP225/H272Y genotype. Additional analysis of SdhB mutant isolates determined that pydiflumetofen controls wild-type as well as sdhBP225/H272R and sdhBP225H/H272 mutants. Pydiflumetofen does not control CG in the sdhBP225/H272Y mutant but is effective in the GTE control. Pydiflumetofen significantly controls Botrytis independently of the SdhB genotype in wounded berry assays. This condition resembles the berry cracking due to heavy rainfall right before harvest, as seen in recent years in the Chilean Central Valley. The findings demonstrate that pydiflumetofen effectively controls the grapevine Botrytis population, suggest a moderate risk of pydiflumetofen resistance, and highlight the significance of incorporating genetic data into the design of control programs.

Dynamic Outlier Slicing Allows Broader Exploration of Adaptive Divergence: A Comparison of Individual Genome and Pool-Seq Data Linked to Humic Adaptation in Perch

How genetic variation contributes to adaptation at different environments is a central focus in evolutionary biology. However, most free-living species still lack a comprehensive understanding of the primary molecular mechanisms of adaptation. Here, we characterised the targets of selection associated with drastically different aquatic environments-humic and clear water-in the common freshwater fish, Eurasian perch (Perca fluviatilis). By using whole-genome sequencing (WGS) on a large population dataset (n = 42 populations) and analysing 873,788 SNPs, our primary aim was to uncover novel and confirm known footprints of selection. We compared individual and pooled WGS, and developed a novel approach, termed dynamic outlier slicing, to assess how the choice of outlier-calling stringency influences functional and Gene Ontology (GO) enrichment. By integrating genome-environment association (GEA) analysis with allele frequency-based approaches, we estimated composite selection signals (CSS) and identified 2679 outlier SNPs distributed across 324 genomic regions, involving 468 genes. Dynamic outlier slicing identified robust enrichment signals in five annotation categories (upstream, downstream, synonymous, 5'UTR and 3'UTR) highlighting the crucial role of regulatory elements in adaptive evolution. Furthermore, GO analyses revealed strong enrichment of molecular functions associated with gated channel activity, transmembrane transporter activity and ion channel activity, emphasising the importance of osmoregulation and ion balance maintenance. Our findings demonstrate that despite substantial random drift and divergence, WGS of high number of population pools enabled the identification of strong selection signals associated with adaptation to both humic and clear water environments, providing robust evidence of widespread adaptation. We anticipate that the dynamic outlier slicing method we developed will enable a more thorough exploration of adaptive divergence across a diverse range of species.

Successful Traceability of Wildlife Samples Contributes to Wildlife Conservation: A Case Study of Tracing the Snub-Nosed Monkey (Rhinopithecus spp.)

Rapid and effective methods for tracing the geographic origin of wildlife samples are essential for tackling the illegal wildlife trade. Traditional morphological categorization methods are often inadequate as relying on the mitochondrial COXI barcode is insufficient for determining geographic populations. To address these limitations, we developed a bioinformatics-based pipeline for the rapid identification of traceable nuclear genome loci. This pipeline has been applied to the whole-genome sequence (WGS) data of China's flagship species, the snub-nosed monkey (Rhinopithecus spp.). These species are known for sex-biased dispersal and hybrid speciation, which complicates genealogy tracing. Using phylogenetic principles, we employed the Robinson and Foulds (RF) distance and scanned over 1,850,726 population-specific loci, identifying five pairs that can trace genealogy origins rapidly and cost-effectively using PCR. Additionally, we found that relying only on mitochondrial genetic information is insufficient for rapid and accurate traceability to subspecies-level geographic populations. Our pipeline efficiently identifies loci and traces the geographic origin of snub-nosed monkey individuals, providing a valuable tool for species preservation and combating the wildlife trade. This approach can be extended to other species, aiding in the conservation of endangered wildlife and tracing criminal evidence.

Fixation of Expression Divergences by Natural Selection in Arabidopsis Coding Genes

Functional divergences of coding genes can be caused by divergences in their coding sequences and expression. However, whether and how expression divergences and coding sequence divergences coevolve is not clear. Gene expression divergences in differentiated cells and tissues recapitulate developmental models within a species, while gene expression divergences between analogous cells and tissues resemble traditional phylogenies in different species, suggesting that gene expression divergences are molecular traits that can be used for evolutionary studies. Using transcriptomes and evolutionary proxies to study gene expression divergences among differentiated cells and tissues in Arabidopsis, expression divergences of coding genes are shown to be strongly anti-correlated with phylostrata (gene ages), indicators of selective constraint Ka/Ks (nonsynonymous replacement rate/synonymous substitution rate) and indicators of positive selection (frequency of loci with Ka/Ks > 1), but only weakly or not correlated with indicators of neutral selection (Ks). Our results thus suggest that expression divergences largely coevolve with coding sequence divergences, suggesting that expression divergences of coding genes are selectively fixed by natural selection but not neutral selection, which provides a molecular framework for trait diversification, functional adaptation and speciation. Our findings therefore support that positive selection rather than negative or neutral selection is a major driver for the origin and evolution of Arabidopsis genes, supporting the Darwinian theory at molecular levels.

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.

Bona Fide Plant Steroid Receptors are Innovated in Seed Plants and Angiosperms through Successive Whole-Genome Duplication Events

Whole-genome duplication (WGD) events are widespread in plants and animals, thus their long-term evolutionary contribution has long been speculated, yet a specific contribution is difficult to verify. Here, we show that ɛ-WGD and ζ-WGD contribute to the origin and evolution of bona fide brassinosteroid (BR) signaling through the innovation of active BR biosynthetic enzymes and active BR receptors from their respective ancestors. We found that BR receptors BRI1 (BR INSENSITIVE 1) and BRL1/3 (BRI1-LIKES 1/3) derived by ɛ-WGD and ζ-WGD, which occurred in the common ancestor of angiosperms and seed plants, respectively, while orphan BR receptor BRL2 first appeared in stomatophytes. Additionally, CYP85A enzymes synthesizing the bioactive BRs derived from a common ancestor of seed plants, while its sister enzymes CYP90 synthesizing BR precursors presented in all land plants, implying possible ligand-receptor coevolution. Consistently, the island domains (IDs) responsible for BR perception in BR receptors were most divergent among different receptor branches, supporting ligand-driven evolution. As a result, BRI1 was the most diversified BR receptor in angiosperms. Importantly, relative to the BR biosynthetic DET2 gene presented in all land plants, BRL2, BRL1/3 and BRI1 had high expression in vascular plants ferns, gymnosperms and angiosperms, respectively. Notably, BRI1 is the most diversified BR receptor with the most abundant expression in angiosperms, suggesting potential positive selection. Therefore, WGDs initiate a neofunctionalization process diverged by ligand-perception and transcriptional expression, which might optimize both BR biosynthetic enzymes and BR receptors, likely contributing to the evolution of land plants, especially seed plants and angiosperms.

The complete mitochondrial genome assembly of Capsicum pubescens reveals key evolutionary characteristics of mitochondrial genes of two Capsicum subspecies

BACKGROUND

Pepper (Capsicum pubescens), one of five domesticated pepper species, has unique characteristics, such as numerous hairs on the epidermis of its leaves and stems, black seeds, and vibrant purple flowers. To date, no studies have reported on the complete assembly of the mitochondrial genome (mitogenome) of C. pubescens. Understanding the mitogenome is crucial for further research on C. pubescens.

RESULTS

In our study, we successfully assembled the first mitogenome of C. pubescens, which was assigned the GenBank accession number OP957066. This mitogenome has a length of 454,165 bp and exhibits the typical circular structure observed in most mitogenomes. We annotated a total of 70 genes, including 35 protein-coding genes (PCGs), 30 tRNA genes, 3 rRNA genes, and 2 pseudogenes. Compared to the other three pepper mitogenomes (KJ865409, KJ865410, and MN196478), C. pubescens OP957066 exhibited four unique PCGs (atp4, atp8, mttB, and rps1), while two PCGs (rpl10 and rps3) were absent. Notably, each of the three pepper mitogenomes from C. annuum (KJ865409, KJ865410, and MN196478) experienced the loss of four PCGs (atp4, atp8, mttB, and rps1). To further explore the evolutionary relationships, we reconstructed a phylogenetic tree using the mitogenomes of C. pubescens and fourteen other species. Structural comparison and synteny analysis of the above four pepper mitogenomes revealed that C. pubescens shares high sequence similarity with KJ865409 and that C. pubescens has rearranged with the other three pepper mitogenomes. Interestingly, we observed 72 similar sequences between the mitochondrial and chloroplast genomes, which accounted for 12.60% of the mitogenome, with a total length of 57,207 bp. These sequences encompassed 12 tRNA genes and the rRNA gene (rrn18). Remarkably, selective pressure analysis suggested that the nad5 gene underwent obvious positive selection. Furthermore, a single-base mutation in three genes (nad1, nad2, and nad4) resulted in an amino acid change.

CONCLUSION

This study provides a high-quality mitogenome of pepper, providing valuable molecular data for future investigations into the exchange of genetic information between pepper organelle genomes.

Novel compound heterozygous variants in LTBP2 associated with relative anterior microphthalmos

PURPOSE

Relative anterior microphthalmos (RAM) is a rare congenital defect associated with severe vision impairment that is primarily caused by genetic alterations. The purpose of this study was to identify the causative genetic variants in two Chinese families with RAM with an autosomal recessive inheritance pattern.

METHODS

DNA samples were obtained from two probands and their family members. Targeted next-generation sequencing (NGS) was used to screen 425 genes associated with inherited eye diseases to identify possible disease-causing variants in the two patients. Sanger sequencing was subsequently used to validate the results in both families.

RESULTS

The targeted NGS panel identified potentially causative novel variants of the latent transforming growth factor beta binding protein 2 (LTBP2) gene in the two RAM families: a missense variant (c.2771C > T; p.Ala924Val) and an intronic variant (c.4582 + 9A > G) in Family A and a different missense variant (c.5239C > A; p.Arg1747Ser) and a synonymous variant (c.951G > A; p.Pro317Pro) in Family B. These four novel variants all cosegregated with the disease phenotype.

CONCLUSION

To our knowledge, this is the first study to report novel LTBP2 gene variants related to RAM. Considering the importance of LTBP2 in ocular development, we provide initial insights into the potential pathogenic mechanisms of LTBP2 in RAM.

Further Evidence for Strong Nonneutrality of Yeast Synonymous Mutations

Although synonymous mutations are commonly assumed neutral or nearly so, recent years have seen reports of fitness effects of synonymous mutations detected under laboratory conditions. In a previous study, we used genome editing to construct thousands of yeast mutants each carrying a synonymous or nonsynonymous mutation in one of 21 genes, and discovered that most synonymous and most nonsynonymous mutations are deleterious. A concern was raised that this observation could be caused by the fitness effects of potential CRISPR/Cas9 off-target edits and/or secondary mutations, and an experiment that would be refractory to such effects was proposed. Using genome sequencing, we here show that no CRISPR/Cas9 off-target editing occurred, although some mutants did carry secondary mutations. Analysis of mutants with negligible effects from secondary mutations and new data collected from the proposed experiment confirms the original conclusion. These findings, along with other reports of fitness effects of synonymous mutations from both case and systematic studies, necessitate a paradigm shift from assuming (near) neutrality of synonymous mutations.

Direct inference of the distribution of fitness effects of spontaneous mutations from recombinant inbred Caenorhabditis elegans mutation accumulation lines

The distribution of fitness effects of new mutations plays a central role in evolutionary biology. Estimates of the distribution of fitness effect from experimental mutation accumulation lines are compromised by the complete linkage disequilibrium between mutations in different lines. To reduce the linkage disequilibrium, we constructed 2 sets of recombinant inbred lines from a cross of 2 Caenorhabditis elegans mutation accumulation lines. One set of lines ("RIAILs") was intercrossed for 10 generations prior to 10 generations of selfing; the second set of lines ("RILs") omitted the intercrossing. Residual linkage disequilibrium in the RIAILs is much less than in the RILs, which affects the inferred distribution of fitness effect when the sets of lines are analyzed separately. The best-fit model estimated from all lines (RIAILs + RILs) infers a large fraction of mutations with positive effects (∼40%); models that constrain mutations to have negative effects fit much worse. The conclusion is the same using only the RILs. For the RIAILs, however, models that constrain mutations to have negative effects fit nearly as well as models that allow positive effects. When mutations in high linkage disequilibrium are pooled into haplotypes, the inferred distribution of fitness effect becomes increasingly negative-skewed and leptokurtic. We conclude that the conventional wisdom-most mutations have effects near 0, a handful of mutations have effects that are substantially negative, and mutations with positive effects are very rare-is likely correct, and that unless it can be shown otherwise, estimates of the distribution of fitness effect that infer a substantial fraction of mutations with positive effects are likely confounded by linkage disequilibrium.

Genome-wide impact of codon usage bias on translation optimization in Drosophila melanogaster

Accuracy and efficiency are fundamental to mRNA translation. Codon usage bias is widespread across species. Despite the long-standing association between optimized codon usage and improved translation, our understanding of its evolutionary basis and functional effects remains limited. Drosophila is widely used to study codon usage bias, but genome-scale experimental data are scarce. Using high-resolution mass spectrometry data from Drosophila melanogaster, we show that optimal codons have lower translation errors than nonoptimal codons after accounting for these biases. Genomic-scale analysis of ribosome profiling data shows that optimal codons are translated more rapidly than nonoptimal codons. Although we find no long-term selection favoring synonymous mutations in D. melanogaster after diverging from D. simulans, we identify signatures of positive selection driving codon optimization in the D. melanogaster population. These findings expand our understanding of the functional consequences of codon optimization and serve as a foundation for future investigations.

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

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.

Synonymous rpsH variants: the common denominator in Escherichia coli adapting to ionizing radiation

Ionizing radiation (IR) in high doses is generally lethal to most organisms. Investigating mechanisms of radiation resistance is crucial for gaining insights into the underlying cellular responses and understanding the damaging effects of IR. In this study, we conducted a comprehensive analysis of sequencing data from an evolutionary experiment aimed at understanding the genetic adaptations to ionizing radiation in Escherichia coli. By including previously neglected synonymous mutations, we identified the rpsH c.294T > G variant, which emerged in all 17 examined isolates across four subpopulations. The identified variant is a synonymous mutation affecting the 30S ribosomal protein S8, and consistently exhibited high detection and low allele frequencies in all subpopulations. This variant, along with two additional rpsH variants, potentially influences translational control of the ribosomal spc operon. The early emergence and stability of these variants suggest their role in adapting to environmental stress, possibly contributing to radiation resistance. Our findings shed light on the dynamics of ribosomal variants during the evolutionary process and their potential role in stress adaptation, providing valuable implications for understanding clinical radiation sensitivity and improving radiotherapy.

Association of mutations in Mycobacterium tuberculosis complex (MTBC) respiration chain genes with hyper-transmission

BACKGROUND

The respiratory chain plays a key role in the growth of Mycobacterium tuberculosis complex (MTBC). However, the exact regulatory mechanisms of this system still need to be elucidated, and only a few studies have investigated the impact of genetic mutations within the respiratory chain on MTBC transmission. This study aims to explore the impact of respiratory chain gene mutations on the global spread of MTBC.

RESULTS

A total of 13,402 isolates of MTBC were included in this study. The majority of the isolates (n = 6,382, 47.62%) belonged to lineage 4, followed by lineage 2 (n = 5,123, 38.23%). Our findings revealed significant associations between Single Nucleotide Polymorphisms (SNPs) of specific genes and transmission clusters. These SNPs include Rv0087 (hycE, G178T), Rv1307 (atpH, C650T), Rv2195 (qcrA, G181C), Rv2196 (qcrB, G1250T), Rv3145 (nuoA, C35T), Rv3149 (nuoE, G121C), Rv3150 (nuoF, G700A), Rv3151 (nuoG, A1810G), Rv3152 (nuoH, G493A), and Rv3157 (nuoM, A1243G). Furthermore, our results showed that the SNPs of atpH C73G, atpA G271C, qcrA G181C, nuoJ G115A, nuoM G772A, and nuoN G1084T were positively correlated with cross-country transmission clades and cross-regional transmission clades.

CONCLUSIONS

Our study uncovered an association between mutations in respiratory chain genes and the transmission of MTBC. This important finding provides new insights for future research and will help to further explore new mechanisms of MTBC pathogenicity. By uncovering this association, we gain a more complete understanding of the processes by which MTBC increases virulence and spread, providing potential targets and strategies for preventing and treating tuberculosis.

Synonymous variant of TLR7 at restriction site rs864058 identified in Covid 19 Pakistani patients

BACKGROUND

TLR7, the receptor accountable for immune response to RNA viruses, has been studied extensively to identify its variants related to the severity of Covid-19 in different populations worldwide. However, the genotype of Pakistani population is still unknown. This study aimed to determine the TLR7 genotypes and their relation with severity in our population.

METHODS

This cross sectional study collected data on 151 Covid-19 positive patients (aged 18-80 years), from June 2022 to May 2023, after an informed consent, from Ziauddin University and Hospital. Prior to that approval from ethics review committee was taken. The demographic variables and comorbidities were recorded along with health status till LAMA (Leave Against Medical Advise), recovery or death. The DNA was extracted from collected blood samples, PCR and Sanger sequencing was done for identification of TLR7 variants. SPSS was used for data analyses and Chi-Square for categorical variables. P-values of <0.05 was considered significant.

RESULTS

Out of 151 patients' sequencing was done for 59 samples. The restriction site, rs864058 of TLR7 gene, identified G/A and G/G variants. This missense variant of TLR7 identified at rs864058 of TLR7 gene, has not been previously reported in population control databases. The genotype G/G was main variant of 49 (83%) patients, whereas, G/A was found in 10 (17%). Majority, 25 (51%) of patients with mild covid-19 had GG genotype but results were not significant (P=0.684). Among female patients the main genotype was GA 8 (80%) while male had G/G 29 (59.2%) with significant results (P=0.024). Since G/G genotype was the major genotype, high percentage was found in hypertensives [20 (40.8%)], Diabetics [13 (26.5%)], depression [24 (49%)] and pneumonia patients [20 (40.8%)]. However, significant association (P=0.023) was only found with pneumonia. Males, in majority had severe [17 (68%)] infection and death [40 (26.4%)], whereas, females had mild [14 (25%)] with [12 (7.9%)] deaths.

CONCLUSION

A variant rs864058 "G/A" of TLR7, in relation to covid-19 were found in our population. Males were found more at risk of morbidity and mortality due to covid-19. Larger studies are required to further confirm these results.

Effect of Plasma-Activated Water on the Cellulase-Producing Strain Aspergillus niger A32

To investigate the effect and mechanism of plasma-activated water (PAW) on Aspergillus niger, PAW was prepared using a needle array-plate dielectric barrier discharge plasma system. The concentrations of long-lived reactive oxygen and nitrogen species (RONS), namely, H2O2, NO2-, and NO3-, in the PAW were 48.76 mg/L, 0.046 mg/L, and 172.36 mg/L, respectively. Chemically activated water (CAW) with the same concentration of long-lived RONS was also prepared for comparison. A. niger A32 was treated with PAW and CAW. After treatment, the treated strains were observed and analyzed with scanning electron microscopy (SEM) and transmission electron microscopy (TEM) to screen probable mutants. The results indicated that the pH, conductivity, and ORP values of PAW were 2.42, 1935 μS/cm, and 517.07 mV, respectively. In contrast, the pH and ORP values of CAW were 6.15 and 301.73 mV, respectively, which differed significantly from those of PAW. In addition, the conductivity of CAW showed no change. SEM and TEM analyses revealed that A. niger A32 treated with CAW exhibited less damage compared with the control. In contrast, A. niger A32 treated with PAW showed significant shrinkage, deformation, and exudate attachment over time. Following PAW treatment, after four passages, a high cellulase-producing stable mutant strain A-WW5 was screened, exhibiting a filter paper enzyme activity of 29.66 U/mL, a cellulose endonuclease activity of 13.79 U/mL, and a β-glucosidase activity of 27.13 U/mL. These values were found to be 33%, 38%, and 2.1% higher than those of the original fungus sample, respectively. In total, 116 SNPs and 61 InDels were present in the genome of the mutant strain A-WW5. The above findings indicate that the impact of PAW on A. niger is not only attributed to long-lasting H2O2, NO2-, and NO3- particles but also to other short-lived active particles; PAW is expected to become a new microbial breeding mutagen.

FLT4 gene polymorphisms influence isolated ventricular septal defect predisposition in a Southwest China population

BACKGROUND

Ventricular septal defect (VSD) is the most common congenital heart disease. Although a small number of genes associated with VSD have been found, the genetic factors of VSD remain unclear. In this study, we evaluated the association of 10 candidate single nucleotide polymorphisms (SNPs) with isolated VSD in a population from Southwest China.

METHODS

Based on the results of 34 congenital heart disease whole-exome sequencing and 1000 Genomes databases, 10 candidate SNPs were selected. A total of 618 samples were collected from the population of Southwest China, including 285 VSD samples and 333 normal samples. Ten SNPs in the case group and the control group were identified by SNaPshot genotyping. The chi-square (χ2) test was used to evaluate the relationship between VSD and each candidate SNP. The SNPs that had significant P value in the initial stage were further analysed using linkage disequilibrium, and haplotypes were assessed in 34 congenital heart disease whole-exome sequencing samples using Haploview software. The bins of SNPs that were in very strong linkage disequilibrium were further used to predict haplotypes by Arlequin software. ViennaRNA v2.5.1 predicted the haplotype mRNA secondary structure. We evaluated the correlation between mRNA secondary structure changes and ventricular septal defects.

RESULTS

The χ2 results showed that the allele frequency of FLT4 rs383985 (P = 0.040) was different between the control group and the case group (P < 0.05). FLT4 rs3736061 (r2 = 1), rs3736062 (r2 = 0.84), rs3736063 (r2 = 0.84) and FLT4 rs383985 were in high linkage disequilibrium (r2 > 0.8). Among them, rs3736061 and rs3736062 SNPs in the FLT4 gene led to synonymous variations of amino acids, but predicting the secondary structure of mRNA might change the secondary structure of mRNA and reduce the free energy.

CONCLUSIONS

These findings suggest a possible molecular pathogenesis associated with isolated VSD, which warrants investigation in future studies.

Polymorphisms in the cysteine dioxygenase gene and their association with taurine content in the Pacific oyster Crassostrea gigas

The Pacific oyster Crassostrea gigas is rich in taurine, which is crucial for its adaptation to the fluctuating intertidal environment and presents significant potential in improving taurine nutrition and boosting immunity in humans. Cysteine dioxygenase (CDO) is a key enzyme involved in the initial step of taurine biosynthesis and plays a crucial role in regulating taurine content in the body. In the present study, polymorphisms of CDO gene in C. gigas (CgCDO) and their association with taurine content were evaluated in 198 individuals. A total of 24 single nucleotide polymorphism (SNP) loci were identified in the exonic region of CgCDO gene by direct sequencing. Among these SNPs, c.279G>A and c.287C>A were found to be significantly associated with taurine content, with the GG and AA genotype at the two loci exhibiting enhanced taurine accumulation (p < 0.05). Haplotype analysis revealed that the 279GG/287AA haplotype had the highest taurine content of 29.24 mg/g, while the 279AA/287CC haplotype showed the lowest taurine content of 21.19 mg/g. These results indicated that the SNPs of CgCDO gene could influence the taurine content in C. gigas and have potential applications in the selective breeding of high-taurine varieties.

Comparative Analysis of Codon Usage Bias in Six Eimeria Genomes

The codon usage bias (CUB) of genes encoded by different species' genomes varies greatly. The analysis of codon usage patterns enriches our comprehension of genetic and evolutionary characteristics across diverse species. In this study, we performed a genome-wide analysis of CUB and its influencing factors in six sequenced Eimeria species that cause coccidiosis in poultry: Eimeria acervulina, Eimeria necatrix, Eimeria brunetti, Eimeria tenella, Eimeria praecox, and Eimeria maxima. The GC content of protein-coding genes varies between 52.67% and 58.24% among the six Eimeria species. The distribution trend of GC content at different codon positions follows GC1 > GC3 > GC2. Most high-frequency codons tend to end with C/G, except in E. maxima. Additionally, there is a positive correlation between GC3 content and GC3s/C3s, but a significantly negative correlation with A3s. Analysis of the ENC-Plot, neutrality plot, and PR2-bias plot suggests that selection pressure has a stronger influence than mutational pressure on CUB in the six Eimeria genomes. Finally, we identified from 11 to 15 optimal codons, with GCA, CAG, and AGC being the most commonly used optimal codons across these species. This study offers a thorough exploration of the relationships between CUB and selection pressures within the protein-coding genes of Eimeria species. Genetic evolution in these species appears to be influenced by mutations and selection pressures. Additionally, the findings shed light on unique characteristics and evolutionary traits specific to the six Eimeria species.

DArTseq-Based, High-Throughput Identification of Novel Molecular Markers for the Detection of Blackleg (Leptosphaeria Spp.) Resistance in Rapeseed

Blackleg disease, caused by Leptosphaeria spp. fungi, is one of the most important diseases of Brassica napus, responsible for severe yield losses worldwide. Blackleg resistance is controlled by major R genes and minor quantitative trait loci (QTL). Due to the high adaptation ability of the pathogen, R-mediated resistance can be easily broken, while the resistance mediated via QTL is believed to be more durable. Thus, the identification of novel molecular markers linked to blackleg resistance for B. napus breeding programs is essential. In this study, 183 doubled haploid (DH) rapeseed lines were assessed in field conditions for resistance to Leptosphaeria spp. Subsequently, DArTseq-based Genome-Wide Association Study (GWAS) was performed to identify molecular markers linked to blackleg resistance. A total of 133,764 markers (96,121 SilicoDArT and 37,643 SNP) were obtained. Finally, nine SilicoDArT and six SNP molecular markers were associated with plant resistance to Leptosphaeria spp. at the highest significance level, p < 0.001. Importantly, eleven of these fifteen markers were found within ten genes located on chromosomes A06, A07, A08, C02, C03, C06 and C08. Given the immune-related functions of the orthologues of these genes in Arabidopsis thaliana, the identified markers hold great promise for application in rapeseed breeding programs.

Toxin-antitoxin system gene mutations driving Mycobacterium tuberculosis transmission revealed by whole genome sequencing

Background

The toxin-antitoxin (TA) system plays a vital role in the virulence and pathogenicity of Mycobacterium tuberculosis (M. tuberculosis). However, the regulatory mechanisms and the impact of gene mutations on M. tuberculosis transmission remain poorly understood.

Objective

To investigate the influence of gene mutations in the toxin-antitoxin system on M. tuberculosis transmission dynamics.

Method

We performed whole-genome sequencing on the analyzed strains of M. tuberculosis. The genes associated with the toxin-antitoxin system were obtained from the National Center for Biotechnology Information (NCBI) Gene database. Mutations correlating with enhanced transmission within the genes were identified by using random forest, gradient boosting decision tree, and generalized linear mixed models.

Results

A total of 13,518 M. tuberculosis isolates were analyzed, with 42.29% (n = 5,717) found to be part of genomic clusters. Lineage 4 accounted for the majority of isolates (n = 6488, 48%), followed by lineage 2 (n = 5133, 37.97%). 23 single nucleotide polymorphisms (SNPs) showed a positive correlation with clustering, including vapB1 G34A, vapB24 A76C, vapB2 T171C, mazF2 C85T, mazE2 G104A, vapB31 T112C, relB T226A, vapB11 C54T, mazE5 T344C, vapB14 A29G, parE1 (C103T, C88T), and parD1 C134T. Six SNPs, including vapB6 A29C, vapB31 T112C, parD1 C134T, vapB37 G205C, Rv2653c A80C, and vapB22 C167T, were associated with transmission clades across different countries. Notably, our findings highlighted the positive association of vapB6 A29C, vapB31 T112C, parD1 C134T, vapB37 G205C, vapB19 C188T, and Rv2653c A80C with transmission clades across diverse regions. Furthermore, our analysis identified 32 SNPs that exhibited significant associations with clade size.

Conclusion

Our study presents potential associations between mutations in genes related to the toxin-antitoxin system and the transmission dynamics of M. tuberculosis. However, it is important to acknowledge the presence of confounding factors and limitations in our study. Further research is required to establish causation and assess the functional significance of these mutations. These findings provide a foundation for future investigations and the formulation of strategies aimed at controlling TB transmission.

Annotation Vocabulary (Might Be) All You Need

Protein Language Models (pLMs) have revolutionized the computational modeling of protein systems, building numerical embeddings that are centered around structural features. To enhance the breadth of biochemically relevant properties available in protein embeddings, we engineered the Annotation Vocabulary, a transformer readable language of protein properties defined by structured ontologies. We trained Annotation Transformers (AT) from the ground up to recover masked protein property inputs without reference to amino acid sequences, building a new numerical feature space on protein descriptions alone. We leverage AT representations in various model architectures, for both protein representation and generation. To showcase the merit of Annotation Vocabulary integration, we performed 515 diverse downstream experiments. Using a novel loss function and only $3 in commercial compute, our premier representation model CAMP produces state-of-the-art embeddings for five out of 15 common datasets with competitive performance on the rest; highlighting the computational efficiency of latent space curation with Annotation Vocabulary. To standardize the comparison of de novo generated protein sequences, we suggest a new sequence alignment-based score that is more flexible and biologically relevant than traditional language modeling metrics. Our generative model, GSM, produces high alignment scores from annotation-only prompts with a BERT-like generation scheme. Of particular note, many GSM hallucinations return statistically significant BLAST hits, where enrichment analysis shows properties matching the annotation prompt - even when the ground truth has low sequence identity to the entire training set. Overall, the Annotation Vocabulary toolbox presents a promising pathway to replace traditional tokens with members of ontologies and knowledge graphs, enhancing transformer models in specific domains. The concise, accurate, and efficient descriptions of proteins by the Annotation Vocabulary offers a novel way to build numerical representations of proteins for protein annotation and design.

A new era of mutation rate analyses: Concepts and methods

The mutation rate is a pivotal biological characteristic, intricately governed by natural selection and historically garnering considerable attention. Recent advances in high-throughput sequencing and analytical methodologies have profoundly transformed our understanding in this domain, ushering in an unprecedented era of mutation rate research. This paper aims to provide a comprehensive overview of the key concepts and methodologies frequently employed in the study of mutation rates. It examines various types of mutations, explores the evolutionary dynamics and associated theories, and synthesizes both classical and contemporary hypotheses. Furthermore, this review comprehensively explores recent advances in understanding germline and somatic mutations in animals and offers an overview of experimental methodologies, mutational patterns, molecular mechanisms, and driving forces influencing variations in mutation rates across species and tissues. Finally, it proposes several potential research directions and pressing questions for future investigations.

The loci of environmental adaptation in a model eukaryote

While the underlying genetic changes have been uncovered in some cases of adaptive evolution, the lack of a systematic study prevents a general understanding of the genomic basis of adaptation. For example, it is unclear whether protein-coding or noncoding mutations are more important to adaptive evolution and whether adaptations to different environments are brought by genetic changes distributed in diverse genes and biological processes or concentrated in a core set. We here perform laboratory evolution of 3360 Saccharomyces cerevisiae populations in 252 environments of varying levels of stress. We find the yeast adaptations to be primarily fueled by large-effect coding mutations overrepresented in a relatively small gene set, despite prevalent antagonistic pleiotropy across environments. Populations generally adapt faster in more stressful environments, partly because of greater benefits of the same mutations in more stressful environments. These and other findings from this model eukaryote help unravel the genomic principles of environmental adaptation.

Downstream-of-gene (DoG) transcripts contribute to an imbalance in the cancer cell transcriptome

Downstream-of-gene (DoG) transcripts are an emerging class of noncoding RNAs. However, it remains largely unknown how DoG RNA production is regulated and whether alterations in DoG RNA signatures exist in major cancers. Here, through transcriptomic analyses of matched tumors and nonneoplastic tissues and cancer cell lines, we reveal a comprehensive catalog of DoG RNA signatures. Through separate lines of evidence, we support the biological importance of DoG RNAs in carcinogenesis. First, we show tissue-specific and stage-specific differential expression of DoG RNAs in tumors versus paired normal tissues with their respective host genes involved in tumor-promoting versus tumor-suppressor pathways. Second, we identify that differential DoG RNA expression is associated with poor patient survival. Third, we identify that DoG RNA induction is a consequence of treating colon cancer cells with the topoisomerase I (TOP1) poison camptothecin and following TOP1 depletion. Our results underlie the significance of DoG RNAs and TOP1-dependent regulation of DoG RNAs in diversifying and modulating the cancer transcriptome.

Mycobacterium tuberculosis is less likely to acquire pathogenic mutations during latent infection than during active disease

Most people infected with Mycobacterium tuberculosis (Mtb) are believed to be in a state of latent tuberculosis (TB) infection (LTBI). Although LTBI is asymptomatic and not infectious, there is a risk of developing active disease even decades after infection. Here, to characterize mutations acquired during LTBI, we collected and analyzed Mtb genomes from seven Japanese patient pairs, each pair consisting of two active TB patients whose starting dates of developing active disease were >3 years apart; one had a high suspicion of LTBI before developing active disease, whereas the other did not. Thereafter, we compared these genomes with those of longitudinal sample pairs within a host of chronic active TB infections combined with public data. The bacterial populations in patients with LTBI were genetically more homogeneous and accumulated single nucleotide polymorphisms (SNPs) slower than those from active disease. Moreover, the lower proportion of nonsynonymous SNPs indicated weaker selective pressures during LTBI than active disease. Finally, the different mutation spectrums indicated different mutators between LTBI and active disease. These results suggest that the likelihood of the acquisition of mutations responsible for antibiotic resistance and increased virulence was lower in the Mtb population from LTBI than active disease.IMPORTANCEControlling latent tuberculosis (TB) infection (LTBI) activation is an effective strategy for TB elimination, where understanding Mycobacterium tuberculosis (Mtb) dynamics within the host plays an important role. Previous studies on chronic active disease reported that Mtb accumulated genomic mutations within the host, possibly resulting in acquired drug resistance and increased virulence. However, several reports suggest that fewer mutations accumulate during LTBI than during the active disease, but the associated risk is largely unknown. Here, we analyzed the genomic dynamics of Mtb within the host during LTBI. Our results statistically suggest that Mtb accumulates mutations during LTBI, but most mutations are under low selective pressures, which induce mutations responsible for drug resistance and virulence. Thus, we propose that LTBI acts as a source for new TB disease rather than as a period for in-host genome evolution.

The relationship between lysozyme gene polymorphism and quality changes during the storage of eggs derived from 2 commercial strains of Japanese quail

During the storage irreversible changes occur in eggs that result in a deterioration of their quality. The most significant changes affect the albumen. One of the major proteins of albumen present in egg white is lysozyme, which protects the embryo from microorganisms. This enzyme also contributes to the qualitative characteristics of albumen. It is possible that its polymorphism also affects the quality and stability of the obtained raw material that is, table eggs. Therefore, the aim of this study was to assess the potential effect of polymorphism in the lysozyme gene and protein on the quality changes during the storage of eggs derived from 2 genetic strains of Japanese quail belonging to various utility types. Eggs from selected females of laying and meat-type breeds were stored for 14 wk. During this period the egg quality traits were evaluated 10 times. DNA was isolated from each female and all exons of the lysozyme gene had been sequenced. In total, fourteen SNPs' and one 4-bp indel mutation were identified in exons and adjacent intronic sequences, among which SNP1 (1:32140723) resulted in a substitution of lysine with glutamine (Q21K). The results showed that SNP1 (strain S22), as well as the SNP2, SNP5, SNP7, SNP8, SNP10, SNP11, SNP12 and SNP13 were significantly associated with breaking strength during egg storage in both investigated Japanese quail strains. Furthermore, a 3 haplotype blocks containing nine SNPs (2, 5, 6, 7, 8, 10, 11, 12 and 13) were identified. These blocks displayed 8 distinct haplotypes that had significant association with breaking strength at all storage time points where egg quality analyses were performed. The study also revealed significant effects of breed and storage time on the egg quality traits. These results provide new insights into the genetic basis of egg quality during storage and could be incorporated into the breeding programs involving these strains.

The Effects of Codon Usage on Protein Structure and Folding

The rate of protein synthesis is slower than many folding reactions and varies depending on the synonymous codons encoding the protein sequence. Synonymous codon substitutions thus have the potential to regulate cotranslational protein folding mechanisms, and a growing number of proteins have been identified with folding mechanisms sensitive to codon usage. Typically, these proteins have complex folding pathways and kinetically stable native structures. Kinetically stable proteins may fold only once over their lifetime, and thus, codon-mediated regulation of the pioneer round of protein folding can have a lasting impact. Supporting an important role for codon usage in folding, conserved patterns of codon usage appear in homologous gene families, hinting at selection. Despite these exciting developments, there remains few experimental methods capable of quantifying translation elongation rates and cotranslational folding mechanisms in the cell, which challenges the development of a predictive understanding of how biology uses codons to regulate protein folding.

Exploring MAP2K3 as a prognostic biomarker and potential immunotherapy target in glioma treatment

Glioma, the most prevalent primary brain tumor in adults, is characterized by significant invasiveness and resistance. Current glioma treatments include surgery, radiation, chemotherapy, and targeted therapy, but these methods often fail to eliminate the tumor completely, leading to recurrence and poor prognosis. Immune checkpoint inhibitors, a class of commonly used immunotherapeutic drugs, have demonstrated excellent efficacy in treating various solid malignancies. Recent research has indicated that unconventional levels of expression of the MAP2K3 gene closely correlates with glioma malignancy, hinting it could be a potential immunotherapy target. Our study unveiled substantial involvement of MAP2K3 in gliomas, indicating the potential of the enzyme to serve as a prognostic biomarker related to immunity. Through the regulation of the infiltration of immune cells, MAP2K3 can affect the prognosis of patients with glioma. These discoveries establish a theoretical foundation for exploring the biological mechanisms underlying MAP2K3 and its potential applications in glioma treatment.

PE/PPE mutations in the transmission of Mycobacterium tuberculosis in China revealed by whole genome sequencing

OBJECTIVE

This study aims to examine the impact of PE/PPE gene mutations on the transmission of Mycobacterium tuberculosis (M. tuberculosis) in China.

METHODS

We collected the whole genome sequencing (WGS) data of 3202 M. tuberculosis isolates in China from 2007 to 2018 and investigated the clustering of strains from different lineages. To evaluate the potential role of PE/PPE gene mutations in the dissemination of the pathogen, we employed homoplastic analysis to detect homoplastic single nucleotide polymorphisms (SNPs) within these gene regions. Subsequently, logistic regression analysis was conducted to analyze the statistical association.

RESULTS

Based on nationwide M. tuberculosis WGS data, it has been observed that the majority of the M. tuberculosis burden in China is caused by lineage 2 strains, followed by lineage 4. Lineage 2 exhibited a higher number of transmission clusters, totaling 446 clusters, of which 77 were cross-regional clusters. Conversely, there were only 52 transmission clusters in lineage 4, of which 9 were cross-regional clusters. In the analysis of lineage 2 isolates, regression results showed that 4 specific gene mutations, PE4 (position 190,394; c.46G > A), PE_PGRS10 (839,194; c.744 A > G), PE16 (1,607,005; c.620T > G) and PE_PGRS44 (2,921,883; c.333 C > A), were significantly associated with the transmission of M. tuberculosis. Mutations of PE_PGRS10 (839,334; c.884 A > G), PE_PGRS11 (847,613; c.1455G > C), PE_PGRS47 (3,054,724; c.811 A > G) and PPE66 (4,189,930; c.303G > C) exhibited significant associations with the cross-regional clusters. A total of 13 mutation positions showed a positive correlation with clustering size, indicating a positive association. For lineage 4 strains, no mutations were found to enhance transmission, but 2 mutation sites were identified as risk factors for cross-regional clusters. These included PE_PGRS4 (338,100; c.974 A > G) and PPE13 (976,897; c.1307 A > C).

CONCLUSION

Our results indicate that some PE/PPE gene mutations can increase the risk of M. tuberculosis transmission, which might provide a basis for controlling the spread of tuberculosis.

Cetuximab chemotherapy resistance: Insight into the homeostatic evolution of head and neck cancer (Review)

The complex evolution of genetic alterations in cancer that occurs in vivo is a selective process involving numerous factors and mechanisms. Chemotherapeutic agents that prevent the growth and spread of cancer cells induce selective pressure, leading to rapid artificial selection of resistant subclones. This rapid evolution is possible because antineoplastic drugs promote alterations in tumor‑cell metabolism, thus creating a bottleneck event. The few resistant cells that survive in this new environment obtain differential reproductive success that enables them to pass down the newly selected resistant gene pool. The present review aims to summarize key findings of tumor evolution, epithelial‑mesenchymal transition and resistance to cetuximab therapy in head and neck squamous cell carcinoma.

Mutational analysis of pulmonary large cell neuroendocrine carcinoma: APC gene mutations identify a good prognostic factor

Pulmonary large cell neuroendocrine carcinoma (LCNEC) is a highly aggressive neoplasm with biological heterogeneity. Mutations in multiple genes have been identified in LCNEC. However, associations between gene alterations, histopathological characteristics, and prognosis remain ambiguous. Here, we investigated the clinicopathologic, immunohistochemical, and genomic characteristics of 19 patients with LCNEC and 9 patients with atypical carcinoid (AC). We revealed high mutation frequencies of TP53 (89.5 %), RB1 (42.1 %), APC (31.6 %), and MCL1 (31.6 %) in LCNEC, while genetic alterations were rarely found in AC. APC alterations mainly occurred to the exon 16 and were only identified in LCNEC with wild-type RB1. The 19 LCNEC were further subgrouped into APC wild-type (LCNEC-APCMT, 6/19) and APC-mutated (LCNEC-APCWT, 13/19) subgroups. In comparison with LCNEC-APCWT, LCNEC-APCMT displayed lower TMB (median: 12.64 vs 4.20, P = 0.045), and relatively mild cytologic atypia. In addition, LCNEC-APCMT distinguished itself from AC and LCNEC-APCWT by obviously downregulated expression of neuroendocrine markers (CD56 and Syn, P < 0.01) and significantly altered expression of genes downstream of APC (β-catenin migrating into the cytoplasm and nucleus, P < 0.001; c-Myc upregulating, P = 0.005). The OS of LCNEC-APCMT was numerically intermediate between AC and LCNEC-APCWT. We first proposed that APC alterations were common in LCNEC with wild-type RB1 and that LCNEC-APCMT was associated with lower TMB and better OS in comparison with LCNEC-APCWT.

Clinical relevance of SCN and CyN induced by ELANE mutations: a systematic review

Introduction

According to the PRISMA criteria, a systematic review has been conducted to investigate the clinical relevance between patients with severe congenital neutropenia (SCN) and cyclic congenital neutropenia (CyN) induced by ELANE mutations.

Methods

We have searched PubMed, EMBASE, Web of Science, Scopus, Cochrane, CNKI, Wanfang Medicine, and VIP for ELANE mutation related literature published from 1997 to 2022. Using Microsoft Excel collect and organize data, SPSS 25, GraphPad Prism 8.0.1, and Omap analyze and plot statistical. Compare the gender, age, geography, mutation sites, infection characteristics, treatment, and other factors of SCN and CyN patients induced by ELANE mutations, with a focus on exploring the relationship between genotype and clinical characteristics, genotype and prognosis.

Results

This study has included a total of 467 patients with SCN and 90 patients with CyN. The onset age of SCN and CyN are both less than 1 year old, and the onset and diagnosis age of SCN are both younger than CyN. The mutation of ELANE gene is mainly missense mutation, and hot spot mutations include S126L, P139L, G214R, c.597+1G>A. The high-frequency mutations with severe outcomes are A57V, L121H, L121P, c.597+1G>A, c.597+1G>T, S126L, C151Y, C151S, G214R, C223X. Respiratory tract, skin and mucosa are the most common infection sites, Staphylococcus aureus, Pseudomonas aeruginosa and Escherichia coli are the most common.

Discussion

Patients with refractory G-CSF are more likely to develop severe outcomes. The commonly used pre-treatment schemes for transplantation are Bu-Cy-ATG and Flu-Bu-ATG. The prognosis of transplantation is mostly good, but the risk of GVHD is high.

Systematic Review Registration

https://www.crd.york.ac.uk/PROSPERO/. PROSPERO, identifier CRD42023434656.

Two synonymous single-nucleotide polymorphisms promoting fluoroquinolone resistance of Escherichia coli in the environment

Single-nucleotide polymorphism (SNP) is one of the core mechanisms that respond to antibiotic resistance of Escherichia coli (E. coli), which is a major issue in environmental pollution. A specific type of SNPs, synonymous SNPs, have been generally considered as the "silent" SNPs since they do not change the encoded amino acid. However, the impact of synonymous SNPs on mRNA splicing, nucleo-cytoplasmic export, stability, and translation was gradually discovered in the last decades. Figuring out the mechanism of synonymous SNPs in regulating antibiotic resistance is critical to improve antimicrobial therapy strategies in clinics and biological treatment strategies of antibiotic-resistant E. coli-polluted materials. With our newly designed antibiotic resistant SNPs prediction algorithm, Multilocus Sequence Type based Identification for Phenotype-single nucleotide polymorphism Analysis (MIPHA), and in vivo validation, we identified 2 important synonymous SNPs 522 G>A and 972 C>T, located at hisD gene, which was previously predicted as a fluoroquinolone resistance-related gene without a detailed mechanism in the E. coli samples with environmental backgrounds. We first discovered that hisD causes gyrA mutation via the upregulation of sbmC and its downstream gene umuD. Moreover, those 2 synonymous SNPs of hisD cause its own translational slowdown and further reduce the expression levels of sbmC and its downstream gene umuD, making the fluoroquinolone resistance determining region of gyrA remains unmutated, ultimately causing the bacteria to lose their ability to resist drugs. This study provided valuable insight into the role of synonymous SNPs in mediating antibiotic resistance of bacteria and a new perspective for the treatment of environmental pollution caused by drug-resistant bacteria.

From cyanobacteria and cyanophages to chloroplasts: the fate of the genomes of oxyphototrophs and the genes encoding photosystem II proteins

Chloroplasts are the result of endosymbiosis of cyanobacterial organisms with proto-eukaryotes. The psbA, psbD and psbO genes are present in all oxyphototrophs and encode the D1/D2 proteins of photosystem II (PSII) and PsbO, respectively. PsbO is a peripheral protein that stabilizes the O2-evolving complex in PSII. Of these genes, psbA and psbD remained in the chloroplastic genome, while psbO was transferred to the nucleus. The genomes of selected cyanobacteria, chloroplasts and cyanophages carrying psbA and psbD, respectively, were analysed. The highest density of genes and coding sequences (CDSs) was estimated for the genomes of cyanophages, cyanobacteria and chloroplasts. The synonymous mutation rate (rS) of psbA and psbD in chloroplasts remained almost unchanged and is lower than that of psbO. The results indicate that the decreasing genome size in chloroplasts is more similar to the genome reduction observed in contemporary endosymbiotic organisms than in streamlined genomes of free-living cyanobacteria. The rS of atpA, which encodes the α-subunit of ATP synthase in chloroplasts, suggests that psbA and psbD, and to a lesser extent psbO, are ancient and conservative and arose early in the evolution of oxygenic photosynthesis. The role of cyanophages in the evolution of oxyphototrophs and chloroplastic genomes is discussed.

Freshwater genome-reduced bacteria exhibit pervasive episodes of adaptive stasis

The emergence of bacterial species is rooted in their inherent potential for continuous evolution and adaptation to an ever-changing ecological landscape. The adaptive capacity of most species frequently resides within the repertoire of genes encoding the secreted proteome (SP), as it serves as a primary interface used to regulate survival/reproduction strategies. Here, by applying evolutionary genomics approaches to metagenomics data, we show that abundant freshwater bacteria exhibit biphasic adaptation states linked to the eco-evolutionary processes governing their genome sizes. While species with average to large genomes adhere to the dominant paradigm of evolution through niche adaptation by reducing the evolutionary pressure on their SPs (via the augmentation of functionally redundant genes that buffer mutational fitness loss) and increasing the phylogenetic distance of recombination events, most of the genome-reduced species exhibit a nonconforming state. In contrast, their SPs reflect a combination of low functional redundancy and high selection pressure, resulting in significantly higher levels of conservation and invariance. Our findings indicate that although niche adaptation is the principal mechanism driving speciation, freshwater genome-reduced bacteria often experience extended periods of adaptive stasis. Understanding the adaptive state of microbial species will lead to a better comprehension of their spatiotemporal dynamics, biogeography, and resilience to global change.

Genetic association and functional validation of ZFP36L2 in non-syndromic orofacial cleft subtypes

BACKGROUND

Non-syndromic orofacial cleft (NSOC) is one of the most common craniofacial malformations with complex etiology. This study aimed to explore the role of specific SNPs in ZFP36L2 and its functional relevance in zebrafish models.

METHODS

We analyzed genetic data of the Chinese Han population from two previous GWAS, comprising of 2512 cases and 2255 controls. Based on the Hardy-Weinberg Equilibrium (HWE) and minor allele frequency (MAF), SNPs in the ZFP36L2 were selected for association analysis. In addition, zebrafish models were used to clarify the in-situ expression pattern of zfp36l2 and the impact of its Morpholino-induced knockdown.

RESULTS

Via association analysis, rs7933 in ZFP36L2 was significantly associated with various non-syndromic cleft lip-only subtypes, potentially conferring a protective effect. Zebrafish embryos showed elevated expression of zfp36l2 in the craniofacial region during critical stages of oral cavity formation. Furthermore, Morpholino-induced knockdown of zfp36l2 led to craniofacial abnormalities, including cleft lip, which was partially rescued by the addition of zfp36l2 mRNA.

CONCLUSION

Our findings highlight the significance of ZFP36L2 in the etiology of NSOC, supported by both human genetic association data and functional studies in zebrafish. These results pave the way for further exploration of targeted interventions for craniofacial malformations.

Translational adaptation in breast cancer metastasis and emerging therapeutic opportunities

Breast cancer's tendency to metastasize poses a critical barrier to effective treatment, making it a leading cause of mortality among women worldwide. A growing body of evidence is showing that translational adaptation is emerging as a key mechanism enabling cancer cells to thrive in the dynamic tumor microenvironment (TME). Here, we systematically summarize how breast cancer cells utilize translational adaptation to drive metastasis, highlighting the intricate regulation by specific translation machinery and mRNA attributes such as sequences and structures, along with the involvement of tRNAs and other trans-acting RNAs. We provide an overview of the latest findings and emerging concepts in this area, discussing their potential implications for therapeutic strategies in breast cancer.

Contiguous identity between entire coding regions of transgenic and native genes rather than special regions is essential for a strong co-suppression

The discovery of co-suppression in plants has greatly boosted the study of gene silencing mechanisms, but its triggering mechanism has remained a mystery. In this study, we explored its possible trigger mechanism by using Fatty acid desaturase 2 (FAD2) and Fatty acid elongase 1 (FAE1) strong co-suppression systems. Analysis of small RNAs in FAD2 co-suppression lines showed that siRNAs distributed throughout the coding region of FAD2 with an accumulated peak. However, mutations of the peak siRNA-matched site and siRNA derived site had not alleviated the co-suppression of its transgenic lines. Synthetic FAD2 (AtFAD2sm), which has synonymous mutations in the entire coding region, failed to trigger any co-suppression. Furthermore, 5' and 3' portions of AtFAD2 and AtFAD2sm were swapped to form two hybrid genes, AtFAD2-3sm and AtFAD2-5sm. 80 % and 92 % of their transgenic lines exhibited co-suppression, respectively. Finally, FAE1s with different degrees of the continuous sequence identity compared with AtFAE1 were tested in their Arabidopsis transgenic lines, and the results showed the co-suppression frequency was reduced as their continuous sequence identity stepped down. This work suggests that contiguous identity between the entire coding regions of transgenic and native genes rather than a special region is essential for a strong co-suppression.

A cross-trait study of lung cancer and its related respiratory diseases based on large-scale exome sequencing population

Background

Genome-wide association studies (GWASs) explain the genetic susceptibility between diseases and common variants. Nevertheless, with the appearance of large-scale sequencing profiles, we could explore the rare coding variants in disease pathogenesis.

Methods

We estimated the genetic correlation of nine respiratory diseases and lung cancer in the UK Biobank (UKB) by linkage disequilibrium score regression (LDSC). Then, we performed exome-wide association studies at single-variant level and gene-level for lung cancer and lung cancer-related respiratory diseases using the whole-exome sequencing (WES) data of 427,934 European participants. Cross-trait meta-analysis was conducted by association analysis based on subsets (ASSET) to identify the pleiotropic variants, while in-silico functional analysis was performed to explore their function. Causal mediation analysis was used to explore whether these pleiotropic variants lead to lung cancer is mediated by affecting the chronic respiratory diseases.

Results

Five respiratory diseases [emphysema, pneumonia, asthma, chronic obstructive pulmonary disease (COPD), and fibrosis] were genetically correlated with lung cancer. We identified 102 significant independent variants at single-variant levels for lung cancer and five lung cancer-related diseases. 15:78590583:G>A (missense variant in CHRNA5) was shared in lung cancer, emphysema, and COPD. Meanwhile, 14 significant genes and 87 suggestive genes were identified in gene-based association tests, including HSD3B7 (lung cancer), SRSF2 (pneumonia), TNXB (asthma), TERT (fibrosis), MOSPD3 (emphysema). Based on the cross-trait meta-analysis, we detected 145 independent pleiotropic variants. We further identified abundant pathways with significant enrichment effects, demonstrating that these pleiotropic genes were functional. Meanwhile, the proportion of mediation effects of these variants ranged from 6 to 23 (emphysema: 23%; COPD: 20%; pneumonia: 20%; fibrosis: 7%; asthma: 6%) through these five respiratory diseases to the incidence of lung cancer.

Conclusions

The identified shared genetic variants, genes, biological pathways, and potential intermediate causal pathways provide a basis for further exploration of the relationship between lung cancer and respiratory diseases.