ABSTRACT. In this paper, we suggest a biomathematical numerical method for analysing mRNA nucleotides sequences based on UA/CG Fibonacci numbers proportions. This method is used to evaluate then compare the spike genes related to the main SARS-CoV2 VARIANTS currently circulating within the world population. The 10 main results proposed to be reproduced by peers are: 1/ SARS-CoV2 genome and spike evolution in one year 2020-2021. 2/ SARS-CoV2 Origins. 3/ Comparing 11 reference variants spikes. 4/ analysing 32 CAL.20C California variant patients spikes. 5/ Toward a meta mRNA Fibonacci gene end message code. 6/ Analysing S501 UK, S484 South Africa and « 2 mutations » INDIA variants. 7/ Suggesting a possible variants spike mRNA palindrome symmetry metastructure improving mRNA stability then infectiousness. 8/ Analysing Fibonacci Metastructures in the mRNA coding for the vaccines PFIZER and MODERNA. 9/ Does the CG-rich modification of the synonymous codons of the spikes of the 2 mRNA vaccines affect the expression and quantity of SARS-CoV2 antibodies? 10/ The exceptional case of the Brazilian variant P.1. Particularly, we suggest the following conjecture at mRNA folding level : CONJECTURE of SARS-CoV2 VARIANTS: The growth of long Fibonacci structures in the shape of "podiums" for almost all of the variants studied (UK, California, South Africa, India, etc.) suggests the probable folding of the Spike mRNA in the form of a "hairpin", which can strengthen the cohesion and the lifespan of this mRNA. Finally, we show that these kinds of Fibonacci matastructures disapear TOTALLY by analysing the published mRNA sequences of PFIZER and MODERNA vaccines. One fact is certain, the two mRNAs of the Moderna and Pfizer vaccines will result in a low functionality of the spike vaccine. This is because their designers by seeking greater stability, have doped to build CG rich sequences which, as soon as they are inserted into the human host, will, paradoxically, seek to mutate, like SARS-CoV2 variants, towards CG ==> UA forms in order to improve their STABILITY and LIFETIME. We conclude using new biomathematics theoretical methods (Master code and numerical standing waves), and comparing the Spikes of the two vaccines Moderna and Pfizer, that there will be very probable differences in stability and shelf life of the two respective mRNAs vaccines. However, “State of the Art” analyzes will disclose that their two protein sequences are strictly identical. By modified their synonymous codons using different strategies, no one can guarantee that the quantity of antibodies generated will be identical in the two cases. We wish to draw attention to the great ADAPTATION power - at the global scale of their genomes - of the most infectious VARIANTS, such as the BRAZIL 20J / 501Y.V3 variant (P.1). This is very worrying for the VACCINES <==> VARIANTS run: We demonstrate how the Brazilian variant P.1 which becomes uncontrollable in Brazil in April 2021 has a level of organization of long metastructures of 17,711 bases covering the genome which is 3.6 more important than that of the 2 reference genomes SARS-CoV2 and worldwide D614G. We suggest that this high level of overall structure of this variant contributes to the stability of this genome and, might explain its greater contagiousness.

This contribution is an attempt to reconstruct the favorable atmosphere that prevailed in Lau-sanne in 1986-87 and provided the backdrop of our invention of two-dimensional ion cyclotron resonance mass spectroscopy (2D ICR-MS). To avoid a self-centered histoire d’ancien combattant, we shall try to emphasize the context, the contributions of key players within our nascent research group at UNIL and the established group of Tino Gäumann at EPFL, the role of external speakers, and the open atmosphere that was not yet polluted by bibliometrics, obsessive concern with im-pact factors, and top-down management of research. We shall also explain why the idea of 2D ICR-MS has been ignored for many years and still has a limited impact: different scientific cul-tures in the ICR and NMR communities, different concerns with fundamental vs. applied research, different status of theory and numerical simulations, different levels of commitment of instrument manufacturers, not to mention many theoretical problems that appear to be at least as challenging in ICR as in NMR.

Some emergent SARS-CoV-2 variants raise concerns due to their altered biological properties. For both B.1.1.7 and B.1351 variants, named as variants of concern (VOC), increased transmissibility was reported, whereas B.1.351 was more resistant to multiple monoclonal antibodies (mAbs), as well as convalescent and vaccination sera. To test this hypothesis, we examined the proportion of VOC over time across different geographic areas where the two VOC, B.1.1.7 and B.1.351, co-circulate. Our comparative analysis was based on the number of SARS-CoV-2 sequences on GISAID database. We report that B.1.1.7 dominates over B.1.351 in geographic areas where both variants co-circulate and the B.1.1.7 was the first variant introduced in the population. The only areas where B.1.351 was detected at higher proportion were South Africa and Mayotte in Africa, where this strain was associated with increased community transmission before the detection of B.1.1.7. The dominance of B.1.1.7 over B.1.351 could be important since B.1.351 was more resistant to certain mAbs, as well as heterologous convalescent and vaccination sera, thus suggesting that it may be transmitted more effectively in people with pre-existing immunity to other VOC. This scenario would lessen the effectiveness of vaccine and urge the need to update them with new strains.

A signal analysis of the genoma sequenced of coronavirus variants: B.1.1.7, B.1.135, B.1.429-B.1.427, B.1.525 and P1 is presented. We deal with a certain type of finite alternating sum series having independently distributed terms associated with binary (0,1) indicators for the nucleotide bases A,C,G,T. This method provides additional information to conventional Similarity comparisons and Power Spectrum approaches. It leads to uncover distinctive patterns regarding the intrinsic data organization of genomic sequences according to its progression along the nucleotide bases position. Hence, the method could be useful for survelliance of genoma variants.

Background and Aims: Dyslipidemia and cardiovascular diseases (CAD) are comorbidities of nonalcoholic fatty liver disease (NAFLD), which ranges from steatosis to hepatocellular carcinoma (HCC). The rs599839 A>G variant, in the CELSR2-PSRC1-SORT1 cluster, has been associated CAD, but its impact on metabolic traits and liver damage in NAFLD has not been investigated yet. Methods: We evaluated the effect of the rs599839 variant in 1426 NAFLD patients (Overall cohort) of whom 131 have HCC (NAFLD-HCC), in 500,000 individuals from the UK Biobank Cohort (UKBBC) and in 366 HCC samples from The Cancer Genome Atlas (TCGA). Hepatic PSRC1, SORT1 and CELSR2 expressions were evaluated by RNAseq (n=125). Results: The rs599839 variant was associated with reduced circulating LDL, carotid intima-media thickness, carotid plaques and hypertension (p<0.05) in NAFLD patients and with protection against dyslipidemia in UKBBC. The G allele was associated with higher risk of HCC and advanced tumor stage (p<0.05) in the Overall cohort. Hepatic PSRC1, SORT1 and CELSR2 expressions were increased in NAFLD patients carrying the rs599839 variant (p<0.0001). SORT1 mRNA levels negatively correlated with circulating lipids and with those of genes involved in lipoprotein turnover (p<0.0001). Conversely, PSRC1 expression was positively related to that of genes implicated in cell proliferation (p<0.0001). In TCGA, PSRC1 over-expression promoted more aggressive HCC development (p<0.05). Conclusions: In sum, the rs599839 A>G variant improves dyslipidemia thus protecting against CAD in NAFLD patients, but as one it might promote HCC development by modulating SORT1 and PSRC1 expressions which impact on lipid metabolism and cell proliferation, respectively

Recently, much attention has been devoted to finding highly efficient and powerful activation functions for CNN layers. Because activation functions inject different nonlinearities between layers that affect performance, varying them is one method for building robust ensembles of CNNs. The objective of this study is to examine the performance of CNN ensembles made with different activation functions, including six new ones presented here: 2D Mexican ReLU, TanELU, MeLU+GaLU, Symmetric MeLU, Symmetric GaLU, and Flexible MeLU. The highest performing ensemble was built with CNNs having different activation layers that randomly replaced the standard ReLU. A comprehensive evaluation of the proposed approach was conducted across fifteen biomedical data sets representing various classification tasks. The proposed method was tested on two basic CNN architectures: Vgg16 and ResNet50. Results demonstrate the superiority in performance of this approach. The MATLAB source code for this study will be available at https://github.com/LorisNanni.

Germline mutations in predisposition genes account for only 20% of all familial colorectal cancer (CRC) and the remaining genetic burden may be due to rare high-to-moderate-penetrance germline variants that are not explored. With the aim of identifying such potential cancer predisposing variants, we performed whole exome sequencing on three CRC cases and three unaffected members of a Polish family and identified two novel heterozygous variants; a coding variant in APC down-regulated 1 gene (APCDD1, p.R299H) and a non-coding variant in the 5’ untranslated region (UTR) of histone deacetylase 5 gene (HDAC5). Sanger sequencing confirmed the variants segregating with the disease and Taqman assays revealed 8 additional APCDD1 variants in a cohort of 1705 familial CRC patients and no further HDAC5 variants. Proliferation assays indicated an insignificant proliferative impact for the APCDD1 variant. Luciferase reporter assays using the HDAC5 variant resulted in an enhanced promoter activity. Targeting of transcription factor binding sites of SNAI-2 and TCF4 interrupted by HDAC5 variant showed a significant impact of TCF4 on promoter activity of mutated HDAC5. Our findings contribute not only to the identification of unrecognized genetic causes of familial CRC but also underline the importance of 5´UTR variants affecting transcriptional regulation and the pathogenesis of complex disorders.

The oncogenic potential of high-risk HPVs is focused on producing the E6 and E7 oncoproteins responsible for disrupting the control of the cell cycle. Epidemiological studies propose the presence of the N29S and H51N variants of the HPV16 E7 protein as a significant association with cervical cancer. It has been suggested that changes in the amino acid sequence of E7 variants may affect the oncoprotein 3D structure; however, this remains unknown. Analysis of the structural differences of the HPV16 E7 protein and its variants (N29S and H51N) was performed through homology modeling and structural refinement by molecular dynamics simulation. We propose for the first time a 3D structure of the E7 reference protein and two of its variants (N29S and H51N) and conclude that the mutations induced by the variants in N29S and H51N have a significant influence on the 3D structure of the E7 protein of HPV16, which could be related to the oncogenic capacity of this protein.

Apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like (APOBEC) proteins are a diverse and evolutionarily conserved family of cytidine deaminases that provide a variety of functions from tissue-specific gene expression and immunoglobulin diversity to control of viruses and retrotransposons. APOBEC family expansion has been documented among mammalian species, suggesting a powerful selection for their activity. Enzymes with a duplicated zinc-binding domain often have catalytically active and inactive domains, yet both have antiviral function. Although APOBEC antiviral function was discovered through hypermutation of HIV-1 genomes lacking an active Vif protein, much evidence indicates that APOBECs also inhibit virus replication through mechanisms other than mutagenesis. Multiple steps of the viral replication cycle may be affected, although nucleic acid replication is a primary target. Packaging of APOBECs into virions was first noted with HIV-1, yet is not a prerequisite for viral inhibition. APOBEC antagonism may occur in viral producer and recipient cells. Signatures of APOBEC activity include G-to-A and C-to-T mutations in a particular sequence context. The importance of APOBEC activity for viral inhibition is reflected in the identification of numerous viral factors, including Vif, which are dedicated to antagonism of these deaminases. Such viral antagonists often are only partially successful, leading to selection for viral variants.

Parkinson’s disease (PD) is characterized by motor deficits and a wide variety of non-motor symptoms. The age of onset, rate of disease progression and the precise profile of motor and non-motor symptoms display considerable individual variation. Neuropathologically, the loss of substantia nigra dopaminergic neurons is a key feature of PD. The vast majority of PD patients exhibit alpha-synuclein aggregates in several brain regions, but there is also great variability in the neuropathology between individuals. While the dopamine replacement therapies can reduce motor symptoms, current therapies do not modify the disease progression. Numerous clinical trials using a wide variety of approaches have failed to achieve disease modification. It has been suggested that the heterogeneity of PD is a major contributing factor to the failure of disease modification trials, and that it is unlikely that a single treatment will be effective in all patients. Precision medicine, using drugs designed to target the pathophysiology in a manner that is specific to each individual with PD, has been suggested as a way forward. PD patients can be stratified according to whether they carry one of the risk variants associated with elevated PD risk. In this review we assess current clinical trials targeting two enzymes, leucine-rich repeat kinase 2 (LRRK2) and glucocerebrosidase (GBA), which are encoded by two most common PD risk genes. Because the details of the pathogenic processes coupled to the different LRRK2and GBA risk variants are not fully understood, we ask if these precision medicine-based intervention strategies will prove “precise“ or “personalized“ enough to modify the disease process in PD patients. We also consider at what phases of the disease that such strategies might be effective, in light of the genes being primarily associated with the risk of developing disease in the first place, and less clearly linked to the rate of disease progression. Finally, we critically evaluate the notion that therapies targeting LRRK2 and GBA might be relevant to a wider segment of PD patients, beyond those that actually carry risk variants of these genes.

Spike protein is the surface glycoprotein of the severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) necessary for the entry of the virus via the transmembrane receptors of the human endothelial cells of the respiratoty system for the virus to be engulfed causing COVID-19 disease after priming by type II transmembrane protease TMPRSS2 and then binding with the angiotensin-converting enzyme 2 (ACE2). Therefore, mutations and amino acid variants analysis are essential in understanding the mechanism of binding of spike protein with its receptor to have an insights on possibilities to design a peptide or nucleotide-based vaccine for COVID-19. Here, we employed Iterative Threading Assembly Refinement (I-TASSER) and Multiple Alignment using Fast Fourier Transform (MAFFT) to predict the three-dimensional monomer structure of spike protein of SARS-CoV-2 and to analyze the amino acid variants for protein sequences from GISAID database for samples collected from Jordan in a try to find an explanation for the low confirmed number of COVID-19 in Jordan. Our Protein Homology/analogY Recognition Engine V 2.0 (Phyre2) findings showed four single amino acid variants (SAV) found in 20 samples of SARS-CoV-2. What is equal to 5% of samples showed tyrosine deletion at Y144 located in the SARS-CoV-like_Spike_S1_NTD (N terminal domain), 62% showed aspartate substitution to glycine at D614G located in the SARS-CoV-2_Spike_S1_RBD (spike recognition binding site), 5% showed aspartate substitution to tyrosine at D1139Y and 5% showed glycine substitution to serine at G1167S both located in the Corona_S2 domain. The findings have shown lower mutational sensitivity in all variants that might not affect the function of spike glycoprotein except for D614G, which has the highest mutational sensitivity score (5 out of 9) indicating a higher likelihood to affect the function of the spike protein. This might suggest, in general, a reduced transmitability of SARS-CoV-2 in Jordan.

Citizen Science has come up to perform analytics over the SARS-CoV-2 genome. Public GALAXY servers provide an automated platform for genomics analysis. Study includes design of GALAXY workflows for RNASEQ assembly and annotation as well as genomic variant discovery and perform analysis across four samples of SARS-CoV-2 infected humans obtained from the local population of Wuhan, China. It provides information about transcriptomics and genomic variants across the SARS-CoV-2 genome. Study can be extended to perform evolutionary and comparative study across each species of coronaviruses. Augmented and integrated study with cheminformatics and immunoinformatics will be a way forward for drug discovery and vaccine development.

A striking aspect of COVID-19 is the difference in the outcome of the infection between different countries, and different ethnic groups within one country. We surveyed the literature on SARS-CoV-2 complemented with comparative publications on SARS-CoV and other coronaviruses to capture the current understanding of virus – host interactions including virus subtypes, transmission, zoonotic aspects, and potential host determinants.

Chromothripsis is a mutational mechanism leading to complex and relatively clustered chromosomal rearrangements resulting in diverse phenotypic outcomes depending on the involved genomic landscapes. It may occur both in the germ and the somatic cells resulting in congenital and developmental disorders and cancer, respectively. Asymptomatic individuals may be carriers of chromotriptic rearrangements and experience recurrent reproductive failures when two or more chromosomes are involved. Several mechanisms are postulated to underly chromothripsis. The most attractive hypothesis involves chromosome pulverization in micronuclei followed by incorrect reassembly of fragments through DNA repair to explain the clustered nature of the observed complex rearrangements. Moreover, exogenous or endogenous DNA damage induction and dicentric bridge formation may be involved. Chromosome instability is commonly observed in the cells of patients with DNA-repair disorders, such as ataxia telangiectasia, Nijmegen breakage syndrome and Bloom syndrome. In addition, germline variations of TP53 have been associated with chromothripsis in Sonic-Hedgehog medulloblastoma and acute myeloid leukemia. In the present review, we focus on the underlying mechanisms of chromothripsis and the involvement of defective DNA-repair genes resulting in chromosome instability and chromothripsis-like rearrangements.

Understanding the role of risk regions identified by genome-wide association studies (GWAS) have made considerable progress lately referred to the post-GWAS era. Annotation of the genes to the GWAS and fine-mapped functional variants, and understanding their biological pathway/gene networks enrichments is expected to give rich dividend by elucidating the mechanisms underlying prostate cancer. To this aim, we compiled and analysed currently available post-GWAS data on genes identified through GWAS and validated through experimental studies in prostate cancer to investigate molecular biological pathways enriched for assigned functional genes. The results highlight some well-known cancer signalling pathways, antigen presentation process and enrichment in cell growth and development gene networks suggesting prostate cancer may result from the accumulation of the effects of functional variants through multiple gene sets and pathways. The upstream analysis identifies critical transcription factors, which supplements the results regarding the regulatory role of the post-GWAS genes. We also identified the common genes between post-GWAS and three well-annotated prostate cancer Oncomine data in patient samples in order to uncover possible main genes in prostate cancer development/progression. Post-GWAS generated knowledge of gene networks and pathways, if analysed further and targeted appropriately, will have an important impact on clinical management of the disease.

The first detection of gypsum (CaSO4.2H2O) by the Mars Science Laboratory (MSL) rover Curiosity in the Gale Crater, Mars created a profound impact on planetary science and exploration. The unique capability of plasma spectroscopy involving in situ elemental analysis in extraterrestrial environments, suggesting the presence of water in the red planet based on phase characterization and providing a clue to Martian paleoclimate. The key to gypsum as an ideal paleoclimate proxy lies in its textural variants, and in this study terrestrial gypsum samples from varied locations and textural types have been analyzed by Laser Induced Breakdown Spectroscopy (LIBS) technique. Petrographic, sub-microscopic and powder X-ray diffraction characterizations confirm the presence of gypsum (hydrated calcium sulphate; CaSO4.2H2O), bassanite (semi-hydrated calcium sulphate; CaSO4.1/2H2O) and anhydrite (anhydrous calcium sulphate; CaSO4) along with accessory phases (quartz and jarosite). The principal component analysis of LIBS spectra from texturally varied gypsums can be differentiated from one another because of the chemical variability in their elemental concentrations. The concentration of gypsum is determined from the partial least-square regressions model. Rapid characterization of gypsum samples with LIBS is expected to work well in extraterrestrial environments.

Genetic testing allows for identification of germline DNA variations which are associated with a significant increase in risk of developing breast and ovarian cancer. Detection of a BRCA1 or BRCA2 pathogenic variant triggers several clinical management actions, which may include increased surveillance and prophylactic surgery for healthy carriers or treatment with PARP inhibitor therapy for carriers diagnosed with cancer. Thus, standardized validated criteria for annotation of BRCA1 and BRCA2 variants according to their pathogenicity are necessary to support clinical decision making and ensure improved outcomes. Upon detection, variants whose pathogenicity can be inferred by the genetic code are typically classified as pathogenic, likely pathogenic, likely benign, or benign. Variants whose impact on function cannot be directly inferred by the genetic code are labeled as Variants of Uncertain Clinical Significance (VUS) and are evaluated by multifactorial likelihood models that use personal and family history of cancer, segregation data, prediction tools, and co-occurrence with a pathogenic BRCA variant. Missense variants, coding alterations that replace a single amino acid residue with another, are a class of variants for which determination of clinical relevance is particularly challenging. Here, we discuss current issues in variant classification by following a typical life cycle of a BRCA1 missense variant through detection, annotation and information dissemination. Advances in massively parallel sequencing have led to a substantial increase in VUS findings. Although comprehensive assessment and classification of missense variants according to their pathogenicity remains the bottle neck, new developments in functional analysis, high throughput assays, data sharing, and statistical models are rapidly changing this scenario.

Current prenatal genetic evaluation showed a significantly increase in non-invasive screening and the reduction of invasive diagnostic procedures. To evaluate the diagnostic efficacy on detecting common aneuploidies, structural chromosomal rearrangements and pathogenic copy number variants (pCNV), we performed a retrospective analysis on a case series initially analyzed by aneuvysion fluorescence in situ hybridization (FISH) and karyotyping then followed by array comparative genomic hybridization (aCGH). Of the 386 cases retrieved from the past decade, common aneuploidies were detected in 137 cases (35.5%), other chromosomal structural rearrangements were detected in four cases (1%), and pCNV were detected in five cases (1.3%). The relative frequencies for common aneuploidies suggested a under detection of sex chromosome aneuploidies. Approximately 9.5% of cases with common aneuploidies showed a mosaic pattern. Inconsistent results between FISH and karyotyping were noted in cases with pseudo-mosaicism introduced by culture artifact or variable cellular proliferation from cells with mosaic karyotypic complements under in vitro cell culture. Based on findings from this case series, cell-based FISH and karyotyping should be performed to detect common aneuploidies, structural chromosomal abnormalities, and mosaic pattern. DNA-based aCGH and reflex FISH should be performed to detect and confirm genomic imbalances and pCNV. Practice points to ensure the diagnostic accuracy and efficacy were summarized.

The meticulous study of finite automata has produced many important and useful results. Automata are simple yet efficient finite state machines that can be utilized in a plethora of situations. It comes, therefore, as no surprise that they have been used in classic game theory in order to model players and their actions. Game theory has recently been influenced by ideas from the field of quantum computation. As a result, quantum versions of classic games have already been introduced and studied. The PQ penny flip game is a famous quantum game introduced by Meyer in 1999. In this paper we investigate all possible finite games that can be played between the two players Q and Picard of the original PQ game. For this purpose we establish a rigorous connection between finite automata and the PQ game along with all its possible variations. Starting from the automaton that corresponds to the original game, we construct more elaborate automata for certain extensions of the game, before finally presenting a semiautomaton that captures the intrinsic behavior of all possible variants of the PQ game. What this means is that from the semiautomaton in question, by setting appropriate initial and accepting states, one can construct deterministic automata able to capture every possible finite game that can be played between the two players Q and Picard. Moreover, we introduce the new concepts of a winning automaton and complete automaton for either player.

There are currently around 78 Nuclear Power Plants (NPP) in the world based on Boiling Water Reactors (BWR). The current parameter to assess BWR instability issues is the linear Decay Ratio (DR). However, it is well known that BWRs are complex non-linear dynamical systems that may even exhibit chaotic dynamics that normally preclude the use of the DR when the BWR is working at a specific operating point during instability. In this work a novel methodology based on an adaptive Shannon Entropy estimator and on Noise Assisted Empirical Mode Decomposition variants is presented. This methodology was developed for real-time implementation of a stability monitor. This methodology was applied to a set of signals stemming from several NPPs reactors (Ringhals-Sweden, Forsmark-Sweden and Laguna Verde-Mexico) under commercial operating conditions, that experienced instabilities events, each one of a different nature

Studies employing Arellano-Bond and Blundell-Bond GMM estimation for single linear dynamic panel data models are growing exponentially in number. However, for researchers it is hard to make a reasoned choice between many different possible implementations of these estimators and associated tests. By simulation the effects are examined of many options regarding: (i) reducing, extending or modifying the set of instruments; (ii) specifying the weighting matrix in relation to the type of heteroskedasticity; (iii) using (robustified) 1-step or (corrected) 2-step variance estimators; (iv) employing 1-step or 2-step residuals in Sargan-Hansen overall or incremental overidentification restrictions tests. This is all done for models in which some regressors may be either strictly exogenous, predetermined or endogenous. Surprisingly, particular asymptotically optimal and relatively robust weighting matrices are found to be superior in finite samples to ostensibly more appropriate versions. Most of the variants of tests for overidentification restrictions show serious deficiencies. A recently developed modification of GMM is found to have great potential when the cross-sectional heteroskedasticity is pronounced and the time-series dimension of the sample not too small. Finally all techniques are employed to actual data and lead to some profound insights.