Exposing synonymous mutations

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.

Daily Practice Assessment of KRAS Status in NSCLC Patients: A New Challenge for the Thoracic Pathologist Is Right around the Corner

Simple Summary

RAS mutation is the most frequent oncogenic alteration in human cancers and KRAS is the most frequently mutated, notably in non-small cell lung carcinomas (NSCLC). Various attempts to inhibit KRAS in the past were unsuccessful in these latter tumors. However, recently, several small molecules (AMG510, MRTX849, JNJ-74699157, and LY3499446) have been developed to specifically target KRAS G12C-mutated tumors, which seems promising for patient treatment and should soon be administered in daily practice for non-squamous (NS)-NSCLC. In this context, it will be mandatory to systematically assess the KRAS status in routine clinical practice, at least in advanced NS-NSCLC, leading to new challenges for thoracic oncologists.

Abstract

KRAS mutations are among the most frequent genomic alterations identified in non-squamous non-small cell lung carcinomas (NS-NSCLC), notably in lung adenocarcinomas. In most cases, these mutations are mutually exclusive, with different genomic alterations currently known to be sensitive to therapies targeting EGFR, ALK, BRAF, ROS1, and NTRK. Recently, several promising clinical trials targeting KRAS mutations, particularly for KRAS G12C-mutated NSCLC, have established new hope for better treatment of patients. In parallel, other studies have shown that NSCLC harboring co-mutations in KRAS and STK11 or KEAP1 have demonstrated primary resistance to immune checkpoint inhibitors. Thus, the assessment of the KRAS status in advanced-stage NS-NSCLC has become essential to setting up an optimal therapeutic strategy in these patients. This stimulated the development of new algorithms for the management of NSCLC samples in pathology laboratories and conditioned reorganization of optimal health care of lung cancer patients by the thoracic pathologists. This review addresses the recent data concerning the detection of KRAS mutations in NSCLC and focuses on the new challenges facing pathologists in daily practice for KRAS status assessment.

Purging of deleterious burden in the endangered Iberian lynx

SignificanceThe dynamics of deleterious variation under contrasting demographic scenarios remain poorly understood in spite of their relevance in evolutionary and conservation terms. Here we apply a genomic approach to study differences in the burden of deleterious alleles between the endangered Iberian lynx (Lynx pardinus) and the widespread Eurasian lynx (Lynx lynx). Our analysis unveils a significantly lower deleterious burden in the former species that should be ascribed to genetic purging, that is, to the increased opportunities of selection against recessive homozygotes due to the inbreeding caused by its smaller population size, as illustrated by our analytical predictions. This research provides theoretical and empirical evidence on the evolutionary relevance of genetic purging under certain demographic conditions.

Ovine ELOVL5 and FASN genes polymorphisms and their correlations with sheep tail fat deposition

Reducing tail fat deposition can increase the economic value of a carcass and improve feed efficiency. This study aimed to explore ELOVL5 and FASN polymorphisms associated with tail fat deposition and their expression levels of sheep. Association analysis showed that ELOVL5 synonymous mutation g.62534 C > T was associated with tail width, tail fat weight, and relative tail fat weight (P < 0.05). FASN synonymous mutation g.12694 A > G was associated with tail length and width (P < 0.05). Combined effect analyses indicated significant differences between the combined genotypes and tail fat deposition. Quantitative real-time reverse transcription PCR indicated that the ELOVL5 and FASN expression levels were significantly higher in tail fat than in other tissues (P < 0.05). ELOVL5 expression levels in tail-fat tissue of big-tail sheep was significantly higher than that in small-tail sheep (P < 0.01). FASN expression levels were significantly higher in tail-fat tissue of small-tail sheep than in that of big-tail sheep (P < 0.05). During development, ELOVL5 tail fat expression increased significantly from 0 to 6 months old (P < 0.05), and FASN expression at 3 months old was significantly higher than that at 0 (minimum) and 6 months old (P < 0.05). Therefore, ELOVL5 and FASN polymorphisms could represent new candidate molecular markers and targets to reduce tail fat deposition in sheep.

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.

QTL mapping of a Brazilian bioethanol strain links the cell wall protein-encoding gene GAS1 to low pH tolerance in S. cerevisiae

BACKGROUND

Saccharomyces cerevisiae is largely applied in many biotechnological processes, from traditional food and beverage industries to modern biofuel and biochemicals factories. During the fermentation process, yeast cells are usually challenged in different harsh conditions, which often impact productivity. Regarding bioethanol production, cell exposure to acidic environments is related to productivity loss on both first- and second-generation ethanol. In this scenario, indigenous strains traditionally used in fermentation stand out as a source of complex genetic architecture, mainly due to their highly robust background-including low pH tolerance.

RESULTS

In this work, we pioneer the use of QTL mapping to uncover the genetic basis that confers to the industrial strain Pedra-2 (PE-2) acidic tolerance during growth at low pH. First, we developed a fluorescence-based high-throughput approach to collect a large number of haploid cells using flow cytometry. Then, we were able to apply a bulk segregant analysis to solve the genetic basis of low pH resistance in PE-2, which uncovered a region in chromosome X as the major QTL associated with the evaluated phenotype. A reciprocal hemizygosity analysis revealed the allele GAS1, encoding a β-1,3-glucanosyltransferase, as the casual variant in this region. The GAS1 sequence alignment of distinct S. cerevisiae strains pointed out a non-synonymous mutation (A631G) prevalence in wild-type isolates, which is absent in laboratory strains. We further showcase that GAS1 allele swap between PE-2 and a low pH-susceptible strain can improve cell viability on the latter of up to 12% after a sulfuric acid wash process.

CONCLUSION

This work revealed GAS1 as one of the main causative genes associated with tolerance to growth at low pH in PE-2. We also showcase how GAS1PE-2 can improve acid resistance of a susceptible strain, suggesting that these findings can be a powerful foundation for the development of more robust and acid-tolerant strains. Our results collectively show the importance of tailored industrial isolated strains in discovering the genetic architecture of relevant traits and its implications over productivity.

Genomic Variation Prediction: A Summary From Different Views

Structural variations in the genome are closely related to human health and the occurrence and development of various diseases. To understand the mechanisms of diseases, find pathogenic targets, and carry out personalized precision medicine, it is critical to detect such variations. The rapid development of high-throughput sequencing technologies has accelerated the accumulation of large amounts of genomic mutation data, including synonymous mutations. Identifying pathogenic synonymous mutations that play important roles in the occurrence and development of diseases from all the available mutation data is of great importance. In this paper, machine learning theories and methods are reviewed, efficient and accurate pathogenic synonymous mutation prediction methods are developed, and a standardized three-level variant analysis framework is constructed. In addition, multiple variation tolerance prediction models are studied and integrated, and new ideas for structural variation detection based on deep information mining are explored.

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.

A novel synonymous homozygous variant [c.2538G>A (p.Thr846Thr)] in TRPM6 in a patient with hypomagnesemia with secondary hypocalcemia

OBJECTIVES

Hypomagnesemia 1, intestinal (HOMG1) is characterized by neurological symptoms that occur due to hypocalcemia and hypomagnesemia and caused by mutations in the TRPM6. Most of the identified variants in TRPM6 lead to premature termination: nonsense, frameshift, deletion, and splice site mutations.

CASE PRESENTATION

Herein, we report a 1.5 month-old case who presented with convulsion due to hypocalcemia and hypomagnesemia in the early infancy. Sequencing of TRPM6 revealed a novel homozygous synonymous variant [c.2538G > A (p.Thr846Thr)] in the last codon of exon 19, which is most likely to affect the splicing. We report a novel homozygous synonymous variant in the TRPM6 leading to HOMG1, expanding the mutational spectrum.

CONCLUSIONS

Synonymous mutations that were previously considered as harmless should be evaluated at the nucleotide level, keeping in mind that they may affect splicing and cause to the disease.

Genetic spectrum and clinical features in a cohort of Chinese patients with autosomal recessive cerebellar ataxias

Background

Although many causative genes have been uncovered in recent years, genetic diagnosis is still missing for approximately 50% of autosomal recessive cerebellar ataxia (ARCA) patients. Few studies have been performed to determine the genetic spectrum and clinical profile of ARCA patients in the Chinese population.

Methods

Fifty-four Chinese index patients with unexplained autosomal recessive or sporadic ataxia were investigated by whole-exome sequencing (WES) and copy number variation (CNV) calling with ExomeDepth. Likely causal CNV predictions were validated by CNVseq.

Results

Thirty-eight mutations including 29 novel ones were identified in 25 out of the 54 patients, providing a 46.3% positive molecular diagnostic rate. Ten different genes were involved, of which four most common genes were SACS, SYNE1, ADCK3 and SETX, which accounted for 76.0% (19/25) of the positive cases. The de novo microdeletion in SACS was reported for the first time in China and the uniparental disomy of ADCK3 was reported for the first time worldwide. Clinical features of the patients carrying SACS, SYNE1 and ADCK3 mutations were summarized.

Conclusions

Our results expand the genetic spectrum and clinical profiles of ARCA patients, demonstrate the high efficiency and reliability of WES combined with CNV analysis in the diagnosis of suspected ARCA, and emphasize the importance of complete bioinformatics analysis of WES data for accurate diagnosis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40035-021-00264-z.

A Code Within a Code: How Codons Fine-Tune Protein Folding in the Cell

The genetic code sets the correspondence between the sequence of a given nucleotide triplet in an mRNA molecule, called a codon, and the amino acid that is added to the growing polypeptide chain during protein synthesis. With four bases (A, G, U, and C), there are 64 possible triplet codons: 61 sense codons (encoding amino acids) and 3 nonsense codons (so-called, stop codons that define termination of translation). In most organisms, there are 20 common/standard amino acids used in protein synthesis; thus, the genetic code is redundant with most amino acids (with the exception of Met and Trp) are being encoded by more than one (synonymous) codon. Synonymous codons were initially presumed to have entirely equivalent functions, however, the finding that synonymous codons are not present at equal frequencies in mRNA suggested that the specific codon choice might have functional implications beyond coding for amino acid. Observation of nonequivalent use of codons in mRNAs implied a possibility of the existence of auxiliary information in the genetic code. Indeed, it has been found that genetic code contains several layers of such additional information and that synonymous codons are strategically placed within mRNAs to ensure a particular translation kinetics facilitating and fine-tuning co-translational protein folding in the cell via step-wise/sequential structuring of distinct regions of the polypeptide chain emerging from the ribosome at different points in time. This review summarizes key findings in the field that have identified the role of synonymous codons and their usage in protein folding in the cell.

A Deep Learning-Based Method for Identification of Bacteriophage-Host Interaction

Multi-drug resistance (MDR) has become one of the greatest threats to human health worldwide, and novel treatment methods of infections caused by MDR bacteria are urgently needed. Phage therapy is a promising alternative to solve this problem, to which the key is correctly matching target pathogenic bacteria with the corresponding therapeutic phage. Deep learning is powerful for mining complex patterns to generate accurate predictions. In this study, we develop PredPHI (Predicting Phage-Host Interactions), a deep learning-based tool capable of predicting the host of phages from sequence data. We collect >3000 phage-host pairs along with their protein sequences from PhagesDB and GenBank databases and extract a set of features. Then we select high-quality negative samples based on the K-Means clustering method and construct a balanced training set. Finally, we employ a deep convolutional neural network to build the predictive model. The results indicate that PredPHI can achieve a predictive performance of 81 percent in terms of the area under the receiver operating characteristic curve on the test set, and the clustering-based method is significantly more robust than that based on randomly selecting negative samples. These results highlight that PredPHI is a useful and accurate tool for identifying phage-host interactions from sequence data.

Rare genetic coding variants associated with human longevity and protection against age-related diseases

Extreme longevity in humans has a strong genetic component, but whether this involves genetic variation in the same longevity pathways as found in model organisms is unclear. Using whole-exome sequences of a large cohort of Ashkenazi Jewish centenarians to examine enrichment for rare coding variants, we found most longevity-associated rare coding variants converge upon conserved insulin/insulin-like growth factor 1 signaling and AMP-activating protein kinase signaling pathways. Centenarians have a number of pathogenic rare coding variants similar to control individuals, suggesting that rare variants detected in the conserved longevity pathways are protective against age-related pathology. Indeed, we detected a pro-longevity effect of rare coding variants in the Wnt signaling pathway on individuals harboring the known common risk allele APOE4. The genetic component of extreme human longevity constitutes, at least in part, rare coding variants in pathways that protect against aging, including those that control longevity in model organisms.

Genes underlying the evolution of tetrapod testes size

BACKGROUND

Testes vary widely in mass relative to body mass across species, but we know very little about which genes underlie and contribute to such variation. This is partly because evidence for which genes are implicated in testis size variation tends to come from investigations involving just one or a few species. Contemporary comparative phylogenetic methods provide an opportunity to test candidate genes for their role in phenotypic change at a macro-evolutionary scale-across species and over millions of years. Previous attempts to detect genotype-phenotype associations across species have been limited in that they can only detect where genes have driven directional selection (e.g. brain size increase).

RESULTS

Here, we introduce an approach that uses rates of evolutionary change to overcome this limitation to test whether any of twelve candidate genes have driven testis size evolution across tetrapod vertebrates-regardless of directionality. We do this by seeking a relationship between the rates of genetic and phenotypic evolution. Our results reveal five genes (Alkbh5, Dmrtb1, Pld6, Nlrp3, Sp4) that each have played unique and complex roles in tetrapod testis size diversity. In all five genes, we find strong significant associations between the rate of protein-coding substitutions and the rate of testis size evolution. Such an association has never, to our knowledge, been tested before for any gene or phenotype.

CONCLUSIONS

We describe a new approach to tackle one of the most fundamental questions in biology: how do individual genes give rise to biological diversity? The ability to detect genotype-phenotype associations that have acted across species has the potential to build a picture of how natural selection has sculpted phenotypic change over millions of years.

Assessment of escape from X chromosome inactivation and gene expression in single human immune cells

X-chromosomal genes escaping from X chromosome inactivation (XCI) in immune cells can contribute to sex-specific differences in immune responses. This protocol describes the specific steps to determine escape from XCI and to simultaneously quantify mRNA expression of multiple genes at the single immune cell level using a single-nucleotide polymorphism approach. The protocol furthermore allows the analysis of allele-specific expression of X-chromosomal genes.

For complete details on the use and execution of this protocol, please refer to Hagen et al. (2020).

•

Approach to investigate escape from XCI and gene expression in single cells

•

Simultaneous gene expression measurement of over 100 genes in one cell

•

Assessment of allele-specific expression of genes with monoallelic expression pattern

Single Nucleotide Polymorphisms in Genes Encoding Toll-Like Receptors 7 and 8 and Their Association with Proviral Load of SRLVs in Goats of Polish Carpathian Breed

Toll-like receptors (TLRs) 7 and 8 are important in single-stranded viral RNA recognition, so genetic variation of these genes may play a role in SRLVs infection and disease progression. Present study aimed to identify SNPs in genes encoding TLR7 and TLR8 in goats of Carpathian breed and analyze their association with the SRLVs provirus concentration as index of disease progression. A total of 14 SNPs were detected, 6 SNPs in the TLR7 gene locus and 8 SNPs in the TLR8 gene. Nine of the 14 identified polymorphisms, 4 in the TLR7 gene and 5 in TLR8 gene, were significantly associated with the SRLVs proviral concentration. These SNPs were located in 3'UTR, 5'UTR and intron sequences as well as in the coding sequences, but they led to silent changes. Homozygous genotypes of three TLR7 SNPs (synonymous variant 1:50703293, 3'UTR variant 1:50701297 and 5'UTR variant 1:50718645) were observed in goats with lower provirus copy number as well as in seronegative animals. The results obtained in this study suggest that SNPs of TLR7/TLR8 genes may induce differential innate immune response towards SRLVs affecting proviral concentration and thereby disease pathogenesis and progression. These findings support a role for genetic variations of TLR7 and TLR8 in SRLVs infection and warrants further studies on the effect of TLR7/TLR8 polymorphisms on SRLVs infection in different populations.

Challenges in the diagnosis and discovery of rare genetic disorders using contemporary sequencing technologies

Next generation sequencing (NGS) has revolutionised rare disease diagnostics. Concomitant with advancing technologies has been a rise in the number of new gene disorders discovered and diagnoses made for patients and their families. However, despite the trend towards whole exome and whole genome sequencing, diagnostic rates remain suboptimal. On average, only ~30% of patients receive a molecular diagnosis. National sequencing projects launched in the last 5 years are integrating clinical diagnostic testing with research avenues to widen the spectrum of known genetic disorders. Consequently, efforts to diagnose genetic disorders in a clinical setting are now often shared with efforts to prioritise candidate variants for the detection of new disease genes. Herein we discuss some of the biggest obstacles precluding molecular diagnosis and discovery of new gene disorders. We consider bioinformatic and analytical challenges faced when interpreting next generation sequencing data and showcase some of the newest tools available to mitigate these issues. We consider how incomplete penetrance, non-coding variation and structural variants are likely to impact diagnostic rates, and we further discuss methods for uplifting novel gene discovery by adopting a gene-to-patient-based approach.

Effects of synonymous mutations beyond codon bias: The evidence for adaptive synonymous substitutions from microbial evolution experiments

Synonymous mutations are often assumed to be neutral with respect to fitness because they do not alter the encoded amino acid and so cannot be 'seen' by natural selection. Yet a growing body of evidence suggests that synonymous mutations can have fitness effects that drive adaptive evolution through their impacts on gene expression and protein folding. Here, we review what microbial experiments have taught us about the contribution of synonymous mutations to adaptation. A survey of site-directed mutagenesis experiments reveals the distributions of fitness effects for nonsynonymous and synonymous mutations are more similar, especially for beneficial mutations, than expected if all synonymous mutations were neutral, suggesting they should drive adaptive evolution more often than is typically observed. A review of experimental evolution studies where synonymous mutations have contributed to adaptation shows they can impact fitness through a range of mechanisms including the creation of illicit RNA polymerase binding sites impacting transcription and changes to mRNA folding stability that modulate translation. We suggest that clonal interference in evolving microbial populations may be the reason synonymous mutations play a smaller role in adaptive evolution than expected based on their observed fitness effects. We finish by discussing the impacts of falsely assuming synonymous mutations are neutral and discuss directions for future work exploring the role of synonymous mutations in adaptive evolution.

Computational Tools to Assess the Functional Consequences of Rare and Noncoding Pharmacogenetic Variability

Interindividual differences in drug response are a common concern in both drug development and across layers of care. While genetics clearly influences drug response and toxicity of many drugs, a substantial fraction of the heritable pharmacological and toxicological variability remains unexplained by known genetic polymorphisms. In recent years, population-scale sequencing projects have unveiled tens of thousands of coding and noncoding pharmacogenetic variants with unclear functional effects that might explain at least part of this missing heritability. However, translating these personalized variant signatures into drug response predictions and actionable advice remains challenging and constitutes one of the most important frontiers of contemporary pharmacogenomics. Conventional prediction methods are primarily based on evolutionary conservation, which drastically reduces their predictive accuracy when applied to poorly conserved pharmacogenes. Here, we review the current state-of-the-art of computational variant effect predictors across variant classes and critically discuss their utility for pharmacogenomics. Besides missense variants, we discuss recent progress in the evaluation of synonymous, splice, and noncoding variations. Furthermore, we discuss emerging possibilities to assess haplotypes and structural variations. We advocate for the development of algorithms trained on pharmacogenomic instead of pathogenic data sets to improve the predictive accuracy in order to facilitate the utilization of next-generation sequencing data for personalized clinical decision support and precision pharmacogenomics.

Genome-Wide Association Analysis for Triazole Resistance in Aspergillus fumigatus

Aspergillus fumigatus is a ubiquitous fungus and the main agent of aspergillosis, a common fungal infection in the immunocompromised population. Triazoles such as itraconazole and voriconazole are the common first-line drugs for treating aspergillosis. However, triazole resistance in A. fumigatus has been reported in an increasing number of countries. While most studies of triazole resistance have focused on mutations in the triazole target gene cyp51A, >70% of triazole-resistant strains in certain populations showed no mutations in cyp51A. To identify potential non-cyp51A mutations associated with triazole resistance in A. fumigatus, we analyzed the whole genome sequences and triazole susceptibilities of 195 strains from 12 countries. These strains belonged to three distinct clades. Our genome-wide association study (GWAS) identified a total of six missense mutations significantly associated with itraconazole resistance and 18 missense mutations with voriconazole resistance. In addition, to investigate itraconazole and pan-azole resistance, Fisher's exact tests revealed 26 additional missense variants tightly linked to the top 20 SNPs obtained by GWAS, of which two were consistently associated with triazole resistance. The large number of novel mutations related to triazole resistance should help further investigations into their molecular mechanisms, their clinical importance, and the development of a comprehensive molecular diagnosis toolbox for triazole resistance in A. fumigatus.

A novel synonymous variant in exon 1 of GNAS gene results in a cryptic splice site and causes pseudohypoparathyroidism type 1A and pseudo-pseudohypoparathyroidism in a French family

INTRODUCTION

Pseudohypoparathyroidism type 1A (PHP1A) and pseudopseudohypoparathyroidism (PPHP) (Inactivating PTH/PTHrP Signaling Disorders type 2, IPPSD2) are two rare autosomal disorders caused by loss-of-function mutations on either maternal or paternal allele, respectively, in the imprinted GNAS gene, which encodes the α subunit of the ubiquitously-expressed stimulatory G protein (Gαs).

CASE PRESENTATION

We investigated a synonymous GNAS variant NM_001077488.2: c.108C>A / p.(Val36=) identified in a family presenting with IPPSD2 phenotype. In silico splicing prediction algorithms were in favor of a deleterious effect of this variant, by creating a new donor splicing site. The GNAS expression studies in blood suggested haploinsufficiency and showed an alternate splice product demonstrating the unmasking of a cryptic site, leading to a 34 base pairs deletion and the creation of a probable unstable RNA.We present the first familial case of IPPSD2 caused by a pathogenic synonymous variant in GNAS gene.

A Genotype Signature for Predicting Pathologic Complete Response in Locally Advanced Rectal Cancer

PURPOSE

To construct and validate a predicting genotype signature for pathologic complete response (pCR) in locally advanced rectal cancer (PGS-LARC) after neoadjuvant chemoradiation.

METHODS AND MATERIALS

Whole exome sequencing was performed in 15 LARC tissues. Mutation sites were selected according to the whole exome sequencing data and literature. Target sequencing was performed in a training cohort (n = 202) to build the PGS-LARC model using regression analysis, and internal (n = 76) and external validation cohorts (n = 69) were used for validating the results. Predictive performance of the PGS-LARC model was compared with clinical factors and between subgroups. The PGS-LARC model comprised 15 genes.

RESULTS

The area under the curve (AUC) of the PGS model in the training, internal, and external validation cohorts was 0.776 (0.697-0.849), 0.760 (0.644-0.867), and 0.812 (0.690-0.915), respectively, and demonstrated higher AUC, accuracy, sensitivity, and specificity than cT stage, cN stage, carcinoembryonic antigen level, and CA19-9 level for pCR prediction. The predictive performance of the model was superior to clinical factors in all subgroups. For patients with clinical complete response (cCR), the positive prediction value was 94.7%.

CONCLUSIONS

The PGS-LARC is a reliable predictive tool for pCR in patients with LARC and might be helpful to enable nonoperative management strategy in those patients who refuse surgery. It has the potential to guide treatment decisions for patients with different probability of tumor regression after neoadjuvant therapy, especially when combining cCR criteria and PGS-LARC.

A synonymous variant in MYO15A enriched in the Ashkenazi Jewish population causes autosomal recessive hearing loss due to abnormal splicing

Nonsyndromic hearing loss is genetically heterogeneous. Despite comprehensive genetic testing, many cases remain unsolved because the clinical significance of identified variants is uncertain or because biallelic pathogenic variants are not identified for presumed autosomal recessive cases. Common synonymous variants are often disregarded. Determining the pathogenicity of synonymous variants may improve genetic diagnosis. We report a synonymous variant c.9861 C > T/p.(Gly3287=) in MYO15A in homozygosity or compound heterozygosity with another pathogenic or likely pathogenic MYO15A variant in 10 unrelated families with nonsyndromic sensorineural hearing loss. Biallelic variants in MYO15A were identified in 21 affected and were absent in 22 unaffected siblings. A mini-gene assay confirms that the synonymous variant leads to abnormal splicing. The variant is enriched in the Ashkenazi Jewish population. Individuals carrying biallelic variants involving c.9861 C > T often exhibit progressive post-lingual hearing loss distinct from the congenital profound deafness typically associated with biallelic loss-of-function MYO15A variants. This study establishes the pathogenicity of the c.9861 C > T variant in MYO15A and expands the phenotypic spectrum of MYO15A-related hearing loss. Our work also highlights the importance of multicenter collaboration and data sharing to establish the pathogenicity of a relatively common synonymous variant for improved diagnosis and management of hearing loss.

The Identification of the SARS-CoV-2 Whole Genome: Nine Cases Among Patients in Banten Province, Indonesia

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the strain of virus that causes coronavirus disease 2019 (COVID-19), the respiratory illness responsible for the current pandemic. Viral genome sequencing has been widely applied during outbreaks to study the relatedness of this virus to other viruses, its transmission mode, pace, evolution and geographical spread, and also its adaptation to human hosts. To date, more than 90,000 SARS-CoV-2 genome sequences have been uploaded to the GISAID database. The availability of sequencing data along with clinical and geographical data may be useful for epidemiological investigations. In this study, we aimed to analyse the genetic background of SARS-CoV-2 from patients in Indonesia by whole genome sequencing. We examined nine samples from COVID-19 patients with RT-PCR cycle threshold (Ct) of less than 25 using ARTIC Network protocols for Oxford Nanopore’s Gridi On sequencer. The analytical methods were based on the ARTIC multiplex PCR sequencing protocol for COVID-19. In this study, we found that several genetic variants within the nine COVID-19 patient samples. We identified a mutation at position 614 P323L mutation in the ORF1ab gene often found in our severe patient samples. The number of SNPs and their location within the SARS-CoV-2 genome seems to vary. This diversity might be responsible for the virulence of the virus and its clinical manifestation.

Role of Synonymous Mutations in the Evolution of TEM β-Lactamase Genes

Nonsynonymous mutations are well documented in TEM β-lactamases. The resulting amino acid changes often alter the conferred phenotype from broad spectrum (2b) conferred by TEM-1 to extended spectrum (2be), inhibitor resistant (2br), or both extended spectrum and inhibitor resistant (2ber). The encoding bla _TEM genes also deviate in numerous synonymous mutations, which are not well understood. bla _TEM-3 (2be), bla _TEM-33 (2br), and bla _TEM-109 (2ber) were studied in comparison to bla _TEM-1 bla _TEM-33 was chosen for more detailed studies because it deviates from bla _TEM-1 by a single nonsynonymous mutation and three additional synonymous mutations. Genes encoding the enzymes with only nonsynonymous or all (including synonymous) mutations plus all permutations between bla _TEM-1 and bla _TEM-33 were expressed in Escherichia coli cells. In disc diffusion assays, genes encoding TEM-3, TEM-33, and TEM-109 with all synonymous mutations resulted in higher resistance levels than genes without synonymous mutations. Disc diffusion assays with the 16 genes carrying all possible nucleotide change combinations between bla _TEM-1 and bla _TEM-33 indicated different susceptibilities for different variants. Nucleotide BLAST searches did not identify genes without synonymous mutations but did identify some without nonsynonymous mutations. Energies of possible secondary mRNA structures calculated with mfold are generally higher with synonymous mutations, suggesting that their role could be to destabilize the mRNA and facilitate its unfolding for efficient translation. In summary, our data indicate that transition from bla _TEM-1 to other variant genes by simply acquiring the nonsynonymous mutations is not favored. Instead, synonymous mutations seem to support the transition to other variant genes with nonsynonymous mutations leading to different phenotypes.

Genetic characterizations of Toll-like receptors in the brown rat and their associations with pathogen infections

Toll-like receptors (TLRs) are important initiators of innate immune responses that target host-pathogen interactions. However, further research into the molecular characteristics of TLRs in wild populations is required, as well as how TLRs genetically influenced pathogen infections in the brown rat (Rattus norvegicus). Here, we explored the genetic characterization and evolution of 2 sensing nucleic acid TLRs (TLR7 and TLR8) and 2 sensing non-nucleic acid TLRs (TLR2 and TLR4) in the wild brown rat, and assessed their associations with 2 RNA viruses (Seoul hantavirus and rat hepatitis E virus (HEV)) and 2 bacteria (Leptospira and Bartonella). In these 4 TLRs, we discovered a total of 16 variants. Furthermore, TLR8 had high genetic diversity among 7 variants, while TLR2 had low genetic diversity with only 1 variant. According to selective pressure analyses, TLR4, TLR7, and TLR8 genes evolved under purifying selection. Interestingly, significant associations were found between 3 TLR8 variants and HEV infection, as well as 1 TLR2 variant and Bartonella infection. Overall, our findings provided a glimpse into the genetic characterization of TLRs in the brown rat, and further demonstrated that TLR2 and TLR8 genetic variations were related to Bartonella and HEV infection, respectively. Especially, TLR8 may be a good candidate immune gene for future research on molecular ecology and functional adaptation in wild populations.

Synonymous point mutation of gtfB gene caused by therapeutic X-rays exposure reduced the biofilm formation and cariogenic abilities of Streptococcus mutans

BACKGROUND

The shift of oral microbiota is a critical factor of radiation caries in head and neck cancer patients after the radiotherapy. However, the direct effects of irradiation on the genome and virulence of cariogenic bacteria are poorly described. Here we investigated the genomic mutations and virulence change of Streptococcus mutans (S. mutans), the major cariogenic bacteria, exposed to the therapeutic doses of X-rays.

RESULTS

X-ray reduced the survival fraction of S. mutans and impacted its biofilm formation. We isolated a biofilm formation-deficient mutant #858 whose genome only possessed three synonymous mutations (c.2043 T > C, c.2100C > T, c.2109A > G) in gtfB gene. The "silent mutation" of c.2043 T > C in gtfB gene can cause the down-regulation of all of the gtfs genes' expression and decrease the GtfB enzyme secretion without the effect on the growth due to the codon bias. #858 and synonymous point mutation strain gtfB ^{2043 T>C}, similar to the gtfB gene null mutant Δ gtfB, can significantly decrease the extracellular polysaccharide production, biofilm formation and cariogenic capabilities both in vitro and in vivo compared with wild type.

CONCLUSION

The direct exposure of X-ray radiation can affect the genome and virulence of oral bacteria even at therapeutic doses. The synonymous mutations of genome are negligent factors for gene expression and related protein translation due to the codon usage frequency.

Association of whole mtDNA, an NADPH G11914A variant, and haplogroups with high physical performance in an elite military troop

Physical performance is a multifactorial and complex trait influenced by environmental and hereditary factors. Environmental factors alone have been insufficient to characterize all outstanding phenotypes. Recent advances in genomic technologies have enabled the investigation of whole nuclear and mitochondrial genome sequences, increasing our ability to understand interindividual variability in physical performance. Our objective was to evaluate the association of mitochondrial polymorphic loci with physical performance in Brazilian elite military personnel. Eighty-eight male military personnel who participated in the Command Actions Course of the Army were selected. Total DNA was obtained from blood samples and a complete mitochondrial genome (mtDNA) was sequenced using Illumina MiSeq platform. Twenty-nine subjects completed the training program (FINISHED, 'F'), and fifty-nine failed to complete (NOT_FINISHED, 'NF'). The mtDNA from NF was slightly more similar to genomes from African countries frequently related to endurance level. Twenty-two distinct mtDNA haplogroups were identified corroborating the intense genetic admixture of the Brazilian population, but their distribution was similar between the two groups (FST=0.0009). Of 745 polymorphisms detected in the mtDNA, the position G11914A within the NADPH gene component of the electron transport chain, was statistically different between F and NF groups (P=0.011; OR: 4.286; 95%CI: 1.198-16.719), with a higher frequency of the G allele in group F individuals). The high performance of military personnel may be mediated by performance-related genomic traits. Thus, mitochondrial genetic markers such as the ND4 gene may play an important role on physical performance variability.

CDH1 gene rs1801552 C/T polymorphism increases susceptibility to esophageal squamous cell carcinoma but not to gastric cardiac adenocarcinoma

PURPOSE

The present study aimed to investigate whether the single nucleotide polymorphism (SNP) rs1801552 C/T in CDH1 gene is correlated with the risk of esophageal squamous cell carcinoma (ESCC) and gastric cardiac adenocarcinoma (GCA), as a preliminary study.

METHODS

The rs1801552 C/T polymorphism was genotyped by the method of polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) in 1316 cancer patients (810 ESCC and 506 GCA) and 1966 controls in north China. We performed two case-control studies, each of which included a population-based set and a hospital-based set.

RESULTS

The data showed that the rs1801552 C/T polymorphism was associated with the risk of ESCC. Allelotype and genotype distributions of the rs1801552 C/T polymorphism in ESCC patients of high-incidence region and hospital were significantly different from that in their respective controls (P < 0.05). Compared with C/C genotype, T/T genotype increased the risk of ESCC in high-incidence region and hospital (age, sex, smoking status and family history of UGIC adjusted odds ratio (OR) = 1.79 and 2.10, 95% confidence interval (CI) = 1.23 - 2.60 and 1.10 - 4.04, respectively). Allelotype and genotype distributions of the rs1801552 C/T polymorphism in GCA patients were not significantly different from that in their controls, respectively (P > 0.05).

CONCLUSIONS

The findings in the present pilot study suggest that the rs1801552 C/T polymorphism was associated with the risk of ESCC, but was not associated with the risk of GCA in high-incidence region and hospital.

Genetic Determinants of Neurobehavioral Responses to Caffeine Administration during Sleep Deprivation: A Randomized, Cross Over Study (NCT03859882)

This study investigated whether four single nucleotide polymorphisms (SNPs) moderated caffeine effects on vigilance and performance in a double-blind and crossover total sleep deprivation (TSD) protocol in 37 subjects. In caffeine (2 × 2.5 mg/kg/24 h) or placebo-controlled condition, subjects performed a psychomotor vigilance test (PVT) and reported sleepiness every six hours (Karolinska sleepiness scale (KSS)) during TSD. EEG was also analyzed during the 09:15 PVT. Carriers of the TNF-α SNP A allele appear to be more sensitive than homozygote G/G genotype to an attenuating effect of caffeine on PVT lapses during sleep deprivation only because they seem more degraded, but they do not perform better as a result. The A allele carriers of COMT were also more degraded and sensitive to caffeine than G/G genotype after 20 h of sleep deprivation, but not after 26 and 32 h. Regarding PVT reaction time, ADORA2A influences the TSD effect but not caffeine, and PER3 modulates only the caffeine effect. Higher EEG theta activity related to sleep deprivation was observed in mutated TNF-α, PER3, and COMT carriers, in the placebo condition particularly. In conclusion, there are genetic influences on neurobehavioral impairments related to TSD that appear to be attenuated by caffeine administration. (NCT03859882).

Synonymous variants that disrupt messenger RNA structure are significantly constrained in the human population

Background

The role of synonymous single-nucleotide variants in human health and disease is poorly understood, yet evidence suggests that this class of “silent” genetic variation plays multiple regulatory roles in both transcription and translation. One mechanism by which synonymous codons direct and modulate the translational process is through alteration of the elaborate structure formed by single-stranded mRNA molecules. While tools to computationally predict the effect of non-synonymous variants on protein structure are plentiful, analogous tools to systematically assess how synonymous variants might disrupt mRNA structure are lacking.

Results

We developed novel software using a parallel processing framework for large-scale generation of secondary RNA structures and folding statistics for the transcriptome of any species. Focusing our analysis on the human transcriptome, we calculated 5 billion RNA-folding statistics for 469 million single-nucleotide variants in 45,800 transcripts. By considering the impact of all possible synonymous variants globally, we discover that synonymous variants predicted to disrupt mRNA structure have significantly lower rates of incidence in the human population.

Conclusions

These findings support the hypothesis that synonymous variants may play a role in genetic disorders due to their effects on mRNA structure. To evaluate the potential pathogenic impact of synonymous variants, we provide RNA stability, edge distance, and diversity metrics for every nucleotide in the human transcriptome and introduce a “Structural Predictivity Index” (SPI) to quantify structural constraint operating on any synonymous variant. Because no single RNA-folding metric can capture the diversity of mechanisms by which a variant could alter secondary mRNA structure, we generated a SUmmarized RNA Folding (SURF) metric to provide a single measurement to predict the impact of secondary structure altering variants in human genetic studies.

β-Arrestin 2 (ARRB2) Polymorphism is Associated With Adverse Consequences of Chronic Heroin Use

BACKGROUND AND OBJECTIVES

β-arrestin 2 is an intracellular protein recruited during the activation of G-protein-coupled receptors. In preclinical studies, β-arrestin 2 has been implicated in µ-opioid receptor desensitization and internalization and the development of opioid tolerance and dependence. The present study investigated relationships between variants in the gene encoding β-arrestin 2 (ARRB2) and clinically relevant phenotypes among individuals with opioid use disorder (OUD). We hypothesized that ARRB2 variants would be associated with the negative effects of long-term heroin use.

METHODS

Chronic heroin users (N = 201; n = 103 African American; n = 98 Caucasian) were genotyped for ARRB2 r1045280 (synonymous, also affecting binding motif of transcription factor GTF2IRD1), rs2036657 (3'UTR) and rs3786047 (intron) and assessed for the past-month frequency of use, injection use, and lifetime duration of heroin use, number of heroin quit-attempts, and heroin use-related consequences.

RESULTS

Lifetime heroin-use consequences (especially occupational and health-related) were significantly lower for African American ARRB2 r1045280 C-allele carriers compared with the TT genotype. There was no significant genotype difference in the Caucasian group. ARRB2 rs2036657 was in strong linkage disequilibrium with rs1045280.

DISCUSSION AND CONCLUSIONS

These results, consistent with extant data, illustrate a role for ancestry-dependent allelic variation in ARRB2 r1045280 on heroin-use consequences. The ARRB2 r1045280 C-allele played a protective role in African-descent participants.

SCIENTIFIC SIGNIFICANCE

These first-in-human findings, which should be replicated, provide support for mechanistic investigations of ARRB2 and related intracellular signaling molecules in OUD etiology, treatment, and relapse prevention. (Am J Addict 2021;00:00-00).

Clinical and immunological effects of mRNA vaccines in malignant diseases

In vitro-transcribed messenger RNA-based therapeutics represent a relatively novel and highly efficient class of drugs. Several recently published studies emphasize the potential efficacy of mRNA vaccines in treating different types of malignant and infectious diseases where conventional vaccine strategies and platforms fail to elicit protective immune responses. mRNA vaccines have lately raised high interest as potent vaccines against SARS-CoV2. Direct application of mRNA or its electroporation into dendritic cells was shown to induce polyclonal CD4+ and CD8+ mediated antigen-specific T cell responses as well as the production of protective antibodies with the ability to eliminate transformed or infected cells. More importantly, the vaccine composition may include two or more mRNAs coding for different proteins or long peptides. This enables the induction of polyclonal immune responses against a broad variety of epitopes within the encoded antigens that are presented on various MHC complexes, thus avoiding the restriction to a certain HLA molecule or possible immune escape due to antigen-loss. The development and design of mRNA therapies was recently boosted by several critical innovations including the development of technologies for the production and delivery of high quality and stable mRNA. Several technical obstacles such as stability, delivery and immunogenicity were addressed in the past and gradually solved in the recent years.This review will summarize the most recent technological developments and application of mRNA vaccines in clinical trials and discusses the results, challenges and future directions with a special focus on the induced innate and adaptive immune responses.

Is Gene-Size an Issue for the Diagnosis of Skeletal Muscle Disorders?

Human genes have a variable length. Those having a coding sequence of extraordinary length and a high number of exons were almost impossible to sequence using the traditional Sanger-based gene-by-gene approach. High-throughput sequencing has partly overcome the size-related technical issues, enabling a straightforward, rapid and relatively inexpensive analysis of large genes.

Several large genes (e.g. TTN, NEB, RYR1, DMD) are recognized as disease-causing in patients with skeletal muscle diseases. However, because of their sheer size, the clinical interpretation of variants in these genes is probably the most challenging aspect of the high-throughput genetic investigation in the field of skeletal muscle diseases.

The main aim of this review is to discuss the technical and interpretative issues related to the diagnostic investigation of large genes and to reflect upon the current state of the art and the future advancements in the field.

Association between genetic polymorphisms of telomere pathway genes and hydrogen peroxide level in omethoate exposure workers

OBJECTIVE

The aim of this study was to explore the association between genetic variations in telomere pathway genes and the level of hydrogen peroxide (H₂O₂) in omethoate exposure workers.

METHODS

A total of 180 omethoate exposure workers and 115 healthy controls were recruited. The level of H₂O₂ in plasma was determined with molybdenic acid colorimetry. Polymerase chain reaction and restriction fragment length was used to detect polymorphisms in POT1 rs1034794, POT1 rs10250202, TERF1 rs3863242, and TERT rs2736098.

RESULTS

The level of H₂O₂ in exposure group (4.26 ± 0.71) was significantly higher than that in control group (3.29 ± 0.46). Generalized linear models indicated that risk factors for the increase H₂O₂ level were exposure [β(95 % CI) = 0.951 (0.806, 1.096), P < 0.001] and AA + AT genotype in POT1 rs034794 [β(95 % CI) = 0.397 (0.049, 0.745), P = 0.025].

CONCLUSION

The increase H₂O₂ level was associated with omethoate exposure and AA + AT genotypes in POT1 gene rs1034794. It provided a new idea that polymorphisms in telomere pathway genes may indirectly regulate telomere length by influencing oxidative stress.

Understanding the codon usage patterns of mitochondrial CO genes among Amphibians

A universal phenomenon of using synonymous codons unequally in coding sequences known as codon usage bias (CUB) is observed in all forms of life. Mutation and natural selection drive CUB in many species but the relative role of evolutionary forces varies across species, genes and genomes. We studied the CUB in mitochondrial (mt) CO genes from three orders of Amphibia using bioinformatics approach as no work was reported yet. We observed that CUB of mt CO genes of Amphibians was weak across different orders. Order Caudata had higher CUB followed by Gymnophiona and Anura for all genes and CUB also varied across genes. Nucleotide composition analysis showed that CO genes were AT-rich. The AT content in Caudata was higher than that in Gymnophiona while Anura showed the least content. Multiple investigations namely nucleotide composition, correspondence analysis, parity plot analysis showed that the interplay of mutation pressure and natural selection caused CUB in these genes. Neutrality plot suggested the involvement of natural selection was more than the mutation pressure. The contribution of natural selection was higher in Anura than Gymnophiona and the lowest in Caudata. The codons CGA, TGA, AAA were found to be highly favoured by nature across all genes and orders.

Genome-Wide Association for Itraconazole Sensitivity in Non-resistant Clinical Isolates of Aspergillus fumigatus

Aspergillus fumigatus is a potentially lethal opportunistic pathogen that infects over ~200,000 people and causes ~100,000 deaths per year globally. Treating A. fumigatus infections is particularly challenging because of the recent emergence of azole-resistance. The majority of studies focusing on the molecular mechanisms underlying azole resistance have examined azole-resistant isolates. However, isolates that are susceptible to azoles also display variation in their sensitivity, presenting a unique opportunity to identify genes contributing to azole sensitivity. Here, we used genome-wide association (GWA) analysis to identify loci involved in azole sensitivity by analyzing the association between 68,853 SNPs and itraconazole (ITCZ) minimum inhibitory concentration (MIC) in 76 clinical isolates of A. fumigatus from Japan. Population structure analysis suggests the presence of four distinct populations, with ITCZ MICs distributed relatively evenly across populations. We independently conducted GWA when treating ITCZ MIC as a quantitative trait and a binary trait, and identified two SNPs with strong associations in both analyses. These SNPs fell within the coding regions of Afu2g02220 and Afu2g02140. We functionally validated Afu2g02220 by knocking it out using a CRISPR/Cas9 approach, because orthologs of this gene are involved in sterol modification and ITCZ targets the ergosterol biosynthesis pathway. Knockout strains displayed no difference in growth compared to the parent strain in minimal media, yet a minor but consistent inhibition of growth in the presence of 0.15 μg/ml ITCZ. Our results suggest that GWA paired with efficient gene deletion is a powerful and unbiased strategy for identifying the genetic basis of complex traits in A. fumigatus.

A Closer Look at Estrogen Receptor Mutations in Breast Cancer and Their Implications for Estrogen and Antiestrogen Responses

Breast cancer (BC) is the most common cancer among women worldwide. More than 70% of BC cases express estrogen receptor alpha (ERα), a central transcription factor that stimulates the proliferation of breast cancer cells, usually in the presence of estrogen. While most cases of ER-positive BC initially respond to antiestrogen therapies, a high percentage of cases develop resistance to treatment over time. The recent discovery of mutated forms of ERα that result in constitutively active forms of the receptor in the metastatic-resistance stage of BC has provided a strong rationale for the development of new antiestrogens. These molecules targeting clinically relevant ERα mutants and a combination with other pharmacological inhibitors of specific pathways may constitute alternative treatments to improve clinical practice in the fight against metastatic-resistant ER-positive BC. In this review, we summarize the latest advances regarding the particular involvement of point mutations of ERα in endocrine resistance. We also discuss the involvement of synonymous ERα mutations with respect to co-translational folding of the receptor and ribosome biogenesis in breast carcinogenesis.

Novel mutations in Uridyl-diphosphate-glucuronosyl-transferase 1A1 (UGT1A1) gene in Tunisian patients with unconjugated hyperbilirubinemia

INTRODUCTION

Unconjugated hyperbilirubinemia (UCB) is a feature of Gilbert's syndrome (GS) and Crigler-Najjar's syndrome (CNS), which are two hereditary defects in bilirubin metabolism. Both syndromes are linked to mutations in the UGT1A1 gene, which cause either the decrease or the absence of the UGT1A1 enzymatic activity. Here, we investigated the molecular basis of the UGT1A1 gene in Tunisian patients presenting with unconjugated hyperbilirubinemia.

METHODS

Twenty-four patients with UCB were investigated. The screening protocol for hemoglobinopathies, enzymopathies, and membrane defects was executed in all patients. Afterward, the molecular analysis of the entire UGT1A1 gene was performed by DNA Sanger sequencing. Several bioinformatic tools were used to explore the effects of novel mutations.

RESULTS

Fifteen different UGT1A1 variations were identified, among which four are described here for the first time. In exon 5, the c.1412C > G; p.(Ala471Gly) and c.1589C > T; p.(Ser530Phe) mutations were detected in patients presenting with CNS type I and GS, respectively. In the 3'UTR region of UGT1A1, the c.*90C > T mutation was detected in 3 patients with CNS type I. In the same region, the c.*388C > T defect was found in a GS patient. A deleterious and damaging effect on the UGT1A1 protein were predicted for both exonic mutations. Furthermore, novel microRNAs were identified as targetting the mutated sequences for the 3'UTR mutations.

CONCLUSION

Our study provides novel data on UCB among Tunisians. Furthermore, we report four novel mutations associated with both GS and CNS. The identification of these mutations increases the spectrum of the UGT1A1 mutations and contributes to an understanding of the molecular abnormalities associated with unconjugated hyperbilirubinemia.

Identification of New, Functionally Relevant Mutations in the Coding Regions of the Human Fos and Jun Proto-Oncogenes in Rheumatoid Arthritis Synovial Tissue

In rheumatoid arthritis (RA), the expression of many pro-destructive/pro-inflammatory proteins depends on the transcription factor AP-1. Therefore, our aim was to analyze the presence and functional relevance of mutations in the coding regions of the AP-1 subunits of the fos and jun family in peripheral blood (PB) and synovial membranes (SM) of RA and osteoarthritis patients (OA, disease control), as well as normal controls (NC). Using the non-isotopic RNAse cleavage assay, one known polymorphism (T252C: silent; rs1046117; present in RA, OA, and NC) and three novel germline mutations of the cfos gene were detected: (i) C361G/A367G: Gln121Glu/Ile123Val, denoted as “fos121/123”; present only in one OA sample; (ii) G374A: Arg125Lys, “fos125”; and (iii) C217A/G374A: Leu73Met/Arg125Lys, “fos73/125”, the latter two exclusively present in RA. In addition, three novel somatic cjun mutations (604–606ΔCAG: ΔGln202, “jun202”; C706T: Pro236Ser, “jun236”; G750A: silent) were found exclusively in the RA SM. Tansgenic expression of fos125 and fos73/125 mutants in NIH-3T3 cells induced an activation of reporter constructs containing either the MMP-1 (matrix metalloproteinase) promoter (3- and 4-fold, respectively) or a pentameric AP-1 site (approximately 5-fold). Combined expression of these two cfos mutants with cjun wildtype or mutants (jun202, jun236) further enhanced reporter expression of the pentameric AP-1 construct. Finally, genotyping for the novel functionally relevant germline mutations in 298 RA, 288 OA, and 484 NC samples revealed no association with RA. Thus, functional cfos/cjun mutants may contribute to local joint inflammation/destruction in selected patients with RA by altering the transactivation capacity of AP-1 complexes.

Genome-wide analysis of codon usage pattern in herpesviruses and its relation to evolution

The preferential use of a specific codon, out of a group of synonymous codons encoding the same amino acid, in a gene transcript results from the bias in codon choice. Various evolutionary forces namely mutation pressure and natural selection influence the pattern of codon usage i.e. distinct for each gene/genome. We investigated the pattern of codon usage of eight human herpesvirus genomes and compared them with two other herpesvirus genomes namely murine herpesvirus 68 and bovine herpesvirus type 1.1 to elucidate its compositional features, pattern of codon usage across the genomes and report the differences of codon usage pattern of human herpesviruses from that of other two other viruses. We also identified the similarity of the codon usage of human herpesviruses with its host (human). The genes were found to be CG rich in HHV2, HHV3, HHV4, HHV6, HHV7 and BH genomes while TA rich in HHV1, HHV5, HHV8 and MH genomes. The codon usage bias (CUB) of genes was low. A highly significant correlation was found among compositional contents depicting the role of mutational pressure along with natural selection in framing CUB. Several more frequently used codons as well as less frequently used codons were identified to be similar between each human virus and its host (human), while murine herpesvirus 68 and bovine herpesvirus type 1.1 genomes did not possess similar adaptation strategy as human herpesviruses to human (host), thus we could conclude that viral CUB might have been shaped as per their host's nature for better surveillance. Neutrality plot revealed mutational pressure mostly influenced the CUB of HHV1, HHV8 and MH viruses, while natural selection had a major impact in the CUB of HHV2, HHV3, HHV4, HHV5, HHV6, HHV7 and BH genomes.

Variants of CAPN3 cause limb-girdle muscular dystrophy type 2A in two Chinese families

Limb-girdle muscular dystrophies (LGMDs) are a group of neuromuscular diseases that are characterized by progressive muscle weakness. LGMD type 2A (LGMD2A), caused by variants in the calpain-3 (CAPN3) gene, is the most prevalent type. The present study aimed to analyze pathogenic CAPN3 gene variants in two pedigrees affected by LGMD2A. Each family contains three patients who are siblings and sought genetic counseling. Genomic DNA was extracted from the peripheral blood samples collected from the probands and family members and whole-exome sequencing (WES) was used to detect the pathogenic genes in the probands. Suspected variants were subsequently validated by Sanger sequencing. In family 1, WES revealed that the proband carried the compound heterogeneous variants c.1194-9A>G and c.1437C>T (p.Ser479=) in CAPN3 (NM_000070.2). In family 2, WES identified that the proband carried the compound heterogeneous variants c.632+4A>G and c.1468C>T (p.Arg490Trp) in CAPN3 (NM_000070.2). In conclusion, the present study indicated that the compound heterogeneous variants of the CAPN3 gene were most likely responsible for LGMD2A in the two Chinese families.

The Effects of Single Nucleotide Polymorphisms in Cancer RNAi Therapies

Tremendous progress in RNAi delivery methods and design has allowed for the effective development of siRNA-based therapeutics that are currently under clinical investigation for various cancer treatments. This approach has the potential to revolutionize cancer therapy by providing the ability to specifically downregulate or upregulate the mRNA of any protein of interest. This exquisite specificity, unfortunately, also has a downside. Genetic variations in the human population are common because of the presence of single nucleotide polymorphisms (SNPs). SNPs lead to synonymous and non-synonymous changes and they occur once in every 300 base pairs in both coding and non-coding regions in the human genome. Much less common are the somatic mosaicism variations associated with genetically distinct populations of cells within an individual that is derived from postzygotic mutations. These heterogeneities in the population can affect the RNAi's efficacy or more problematically, which can lead to unpredictable and sometimes adverse side effects. From a more positive viewpoint, both SNPs and somatic mosaicisms have also been implicated in human diseases, including cancer, and these specific changes could offer the ability to effectively and, more importantly, selectively target the cancer cells. In this review, we discuss how SNPs in the human population can influence the development and success of novel anticancer RNAi therapies and the importance of why SNPs should be carefully considered.

A de novo synonymous variant in EFTUD2 disrupts normal splicing and causes mandibulofacial dysostosis with microcephaly: case report

Background

Mandibulofacial dysostosis with microcephaly (MFDM) is a rare autosomal dominant genetic disease characterized by intellectual and growth retardations, as well as major microcephaly, induced by missense and splice site variants or microdeletions in the EFTUD2 gene.

Case presentation

Here, we investigate the case of a young girl with symptoms of MFDM and a normal karyotype. Whole-exome sequencing of the family was performed to identify genetic alterations responsible for this phenotype. We identified a de novo synonymous variant in the EFTUD2 gene. We demonstrated that this synonymous variant disrupts the donor splice-site in intron 9 resulting in the skipping of exon 9 and a frameshift that leads to a premature stop codon.

Conclusions

We present the first case of MFDM caused by a synonymous variant disrupting the donor splice site, leading to exon skipping.

SARS-CoV-2 Codon Usage Bias Downregulates Host Expressed Genes With Similar Codon Usage

Severe acute respiratory syndrome has spread quickly throughout the world and was declared a pandemic by the World Health Organization (WHO). The pathogenic agent is a new coronavirus (SARS-CoV-2) that infects pulmonary cells with great effectiveness. In this study we focus on the codon composition for the viral protein synthesis and its relationship with the protein synthesis of the host. Our analysis reveals that SARS-CoV-2 preferred codons have poor representation of G or C nucleotides in the third position, a characteristic which could result in an unbalance in the tRNAs pools of the infected cells with serious implications in host protein synthesis. By integrating this observation with proteomic data from infected cells, we observe a reduced translation rate of host proteins associated with highly expressed genes and that they share the codon usage bias of the virus. The functional analysis of these genes suggests that this mechanism of epistasis can contribute to understanding how this virus evades the immune response and the etiology of some deleterious collateral effect as a result of the viral replication. In this manner, our finding contributes to the understanding of the SARS-CoV-2 pathogeny and could be useful for the design of a vaccine based on the live attenuated strategy.

Associations between the Bovine Myostatin Gene and Milk Fatty Acid Composition in New Zealand Holstein-Friesian × Jersey-Cross Cows

Simple Summary

The gene that encodes myostatin influences more than one trait, and its expression has been observed in skeletal muscle, as well as the mammary gland. In this study, association analysis revealed that variation in the bovine myostatin gene affects milk fatty acid composition, raising the possibility that this genetic variation may be utilized to increase the amount of unsaturated fatty acid and decrease the amount of saturated fatty acid in milk.

Abstract

The myostatin gene (MSTN), which encodes the protein myostatin, is pleiotropic, and its expression has been associated with both increased and decreased adipogenesis and increased skeletal muscle mass in animals. In this study, the polymerase chain reaction, coupled with single strand conformation polymorphism analysis, was utilized to reveal nucleotide sequence variation in bovine MSTN in 410 New Zealand (NZ) Holstein-Friesian × Jersey (HF × J)-cross cows. These cows ranged from 3 to 9 years of age and over the time studied, produced an average 22.53 ± 2.18 L of milk per day, with an average milk fat content of 4.94 ± 0.17% and average milk protein content of 4.03 ± 0.10%. Analysis of a 406-bp amplicon from the intron 1 region, revealed five nucleotide sequence variants (A–E) that contained seven nucleotide substitutions. Using general linear mixed-effect model analyses the AD genotype was associated with reduced C10:0, C12:0, and C12:1 levels when compared to levels in cows with the AA genotype. These associations in NZ HF × J cross cows are novel, and they suggest that this variation in bovine MSTN could be explored for increasing the amount of milk unsaturated fatty acid and decreasing the amount of saturated fatty acid.

Critical Analysis of Genome-Wide Association Studies: Triple Negative Breast Cancer Quae Exempli Causa

Genome-wide association studies (GWAS) are useful in assessing and analyzing either differences or variations in DNA sequences across the human genome to detect genetic risk factors of diseases prevalent within a target population under study. The ultimate goal of GWAS is to predict either disease risk or disease progression by identifying genetic risk factors. These risk factors will define the biological basis of disease susceptibility for the purposes of developing innovative, preventative, and therapeutic strategies. As single nucleotide polymorphisms (SNPs) are often used in GWAS, their relevance for triple negative breast cancer (TNBC) will be assessed in this review. Furthermore, as there are different levels and patterns of linkage disequilibrium (LD) present within different human subpopulations, a plausible strategy to evaluate known SNPs associated with incidence of breast cancer in ethnically different patient cohorts will be presented and discussed. Additionally, a description of GWAS for TNBC will be presented, involving various identified SNPs correlated with miRNA sites to determine their efficacies on either prognosis or progression of TNBC in patients. Although GWAS have identified multiple common breast cancer susceptibility variants that individually would result in minor risks, it is their combined effects that would likely result in major risks. Thus, one approach to quantify synergistic effects of such common variants is to utilize polygenic risk scores. Therefore, studies utilizing predictive risk scores (PRSs) based on known breast cancer susceptibility SNPs will be evaluated. Such PRSs are potentially useful in improving stratification for screening, particularly when combining family history, other risk factors, and risk prediction models. In conclusion, although interpretation of the results from GWAS remains a challenge, the use of SNPs associated with TNBC may elucidate and better contextualize these studies.

Synonymous Mutations of Porcine Igf1r Extracellular Domain Affect Differentiation and Mineralization in MC3T3-E1 Cells

Owing to the wide application of miniature pigs in biomedicine, the formation mechanism of its short stature must be elucidated. The insulin-like growth factor 1 receptor (IGF-1R), which receives signals through the extracellular domain (ECD) binding with ligands, is crucial in regulating cell growth and bone matrix mineralization. In this study, two haplotypes of Igf1r with four synonymous mutations in the coding sequences of IGF-1R ECD between large pigs (LP) and Bama pigs (BM) were stably expressed in the Igf1r-knockout MC3T3-E1 cells and named as MC3T3-LP cells (LP group) and MC3T3-BM cells (BM group), respectively. IGF-1R expression was lower in the BM group than in the LP group both in terms of transcription and translation levels, and IGF-1R expression inhibited cell proliferation. In addition, IGF-1R expression in the BM group promoted early-stage differentiation but delayed late-stage differentiation, which not only suppressed the expression of bone-related factors but also reduced alkaline phosphatase activity and calcium deposition. Moreover, different haplotypes of Igf1r changed the stability and conformation of the protein, further affecting the binding with IGF-1. Our data indicated that the four synonymous mutations of IGF1R ECD encoded by affect gene transcription and translation, thereby further leading to differences in the downstream pathways and functional changes of osteoblasts.

The Silence Speaks, but We Do Not Listen: Synonymous c.1824C>T Gene Variant in the Last Exon of the Prothrombin Gene as a New Prothrombotic Risk Factor

BACKGROUND

Thrombosis is a major global disease burden with almost 60% of cases related to underlying heredity and most cases still idiopathic. Synonymous single nucleotide polymorphisms (sSNPs) are considered silent and phenotypically neutral. Our previous study revealed a novel synonymous FII c.1824C>T variant as a potential risk factor for pregnancy loss, but it has not yet been associated with thrombotic diseases.

METHODS

To determine the frequency of the FII c.1824C>T variant we have sequenced patients' DNA. Prothrombin RNA expression was measured by quantitative PCR. Functional analyses included routine hemostasis tests, western blotting and ELISA to determine prothrombin levels in plasma, and global hemostasis assays for thrombin and fibrin generation in carriers of the FII c.1824C>T variant. Scanning electron microscopy was used to examine the structure of fibrin clots.

RESULTS

Frequency of the FII c.1824C>T variant was significantly increased in patients with venous thromboembolism and cerebrovascular insult. Examination in vitro demonstrated increased expression of prothrombin mRNA in FII c.1824T transfected cells. Our ex vivo study of FII c.1824C>T carriers showed that the presence of this variant was associated with hyperprothrombinemia, hypofibrinolysis, and formation of densely packed fibrin clots resistant to fibrinolysis.

CONCLUSION

Our data indicate that FII c.1824C>T, although a synonymous variant, leads to the development of a prothrombotic phenotype and could represent a new prothrombotic risk factor. As a silent variant, FII c.1824C>T would probably be overlooked during genetic screening, and our results show that it could not be detected in routine laboratory tests.

A homozygote variant in the tRNA splicing endonuclease subunit 54 causes pontocerebellar hypoplasia in a consanguineous Iranian family

Background

Homozygous loss‐of‐function mutations in TSEN54 (tRNA splicing endonuclease subunit 54; OMIM: 608755) cause different types of pontocerebellar hypoplasias (PCH) including PCH2, PCH4, and PCH5. The study aimed to determine the possible genetic factors contributing to PCH phenotypes in two affected male infants in an Iranian family.

Methods

We subjected two affected individuals in a consanguineous Iranian family. To systematically investigate the susceptible gene(s), whole‐exome sequencing was performed on the proband and a novel identified variant was confirmed by Sanger sequencing. We also analyzed 26 relatives in three generations using PCR‐restriction fragment length polymorphism (PCR‐RFLP) followed and confirmed by Sanger sequencing.

Results

Physical and medical examinations confirmed PCH in the patients. Besides, the proband showed bilateral moderate sensorineural hearing loss and structural heart defects as the novel phenotypes. The molecular findings also verified that two affected individuals were homozygote for the novel synonymous variant, NM_207346.2: c.1170G>A; p.(Val390Val), in TSEN54. PCR‐RFLP and Sanger sequencing elucidated that the parents and 16 relatives were heterozygote for the novel variant.

Conclusion

We identified a novel synonymous variant, c.1170G>A, in TSEN54 associated with PCH in an Iranian family. Based on this study, we strongly suggest using “TSENopathies” to show the overlapped phenotypes among different types of PCH resulted from TSEN causative mutations.