Trends of mutation accumulation across global SARS-CoV-2 genomes: Implications for the evolution of the novel coronavirus

To understand SARS-CoV-2 microevolution, this study explored the genome-wide frequency, gene-wise distribution, and molecular nature of all point-mutations detected across its 71,703 RNA-genomes deposited in GISAID till 21 August 2020. Globally, nsp1/nsp2 and orf7a/orf3a were the most mutation-ridden non-structural and structural genes respectively. Phylogeny of 4618 spatiotemporally-representative genomes revealed that entities belonging to the early lineages are mostly spread over Asian countries, including India, whereas the recently-derived lineages are more globally distributed. Of the total 20,163 instances of polymorphism detected across global genomes, 12,594 and 7569 involved transitions and transversions, predominated by cytidine-to-uridine and guanosine-to-uridine conversions, respectively. Positive selection of nonsynonymous mutations (dN/dS >1) in most of the structural, but not the non-structural, genes indicated that SARS-CoV-2 has already harmonized its replication/transcription machineries with the host metabolism, while it is still redefining virulence/transmissibility strategies at the molecular level. Mechanistic bases and evolutionary/pathogenicity-related implications are discussed for the predominant mutation-types.

Clinical utility of KRAS and BRAF mutations in a cohort of patients with colorectal neoplasms submitted for microsatellite instability testing

BACKGROUND

Molecular analysis has become important in colorectal carcinoma (CRC) evaluation. Alterations in KRAS, BRAF, or mismatch repair (MMR) genes may determine therapeutic response or define a hereditary cancer syndrome. Correlation of DNA studies with clinical findings will further clarify the clinical utility of these markers.

PATIENTS AND METHODS

A retrospective study was performed on 111 paraffin-embedded tumor specimens submitted for microsatellite instability (MSI) testing based on clinical history or histologic examination, or both. DNA samples were screened for 7 KRAS mutations and the BRAF p.V600E mutation using fluorescent allele-specific polymerase-chain reaction (PCR) and capillary electrophoresis. Clinical data were collected through chart review.

RESULTS

Fifty-eight male and 53 female patients were studied. The incidence of KRAS and BRAF mutations was 49.5% and 7.2%, respectively. Dideoxy sequencing verified KRAS mutation status in 46 of 49 specimens tested. There was a trend toward significance of individual KRAS mutations on survival (P = .003). Dually positive KRAS and MSI tumors exclusively demonstrated p.G12D and p.G13D mutations (G>A transitions). BRAF-mutated tumors were predominantly right-sided and associated with a borderline worse prognosis. Forty-eight percent of tumors with MSI were present in the left colon or rectum.

CONCLUSION

Allele-specific PCR is an accurate and convenient method to assess KRAS and BRAF mutations and may detect mutations not identified by dideoxy sequencing. KRAS mutation status, in conjunction with morphologic or clinical parameters, may be useful in determining whether a tumor should be tested for MSI. MSI testing should not be considered exclusively in right-sided lesions. BRAF analysis may not be useful in rectal adenocarcinomas and should be evaluated in larger studies.

Genomic and evolutionary aspects of chloroplast tRNA in monocot plants

BACKGROUND

Chloroplasts are one of the most indispensable organelles that make life forms on the earth possible by their capacity to photosynthesize. These organelles possess a circular genome with a number of coding genes responsible for self-regulation. tRNAs are an important evolutionary-conserved gene family that are responsible for protein translation. However, within the chloroplast genome, tRNA machinery are poorly understood.

RESULTS

In the present study, the chloroplast genome of six monocot plants, Oryza nivara (NC_005973), Oryza sativa (NC_001320), Sachharum officinarum (NC_006084), Sorghum bicolor (NC_008602), Triticum aestivum (NC_002762), and Zea mays (NC_001666) were downloaded and analyzed to identify tRNA sequences. Further analysis of the tRNA sequences in the chloroplast genomes of the monocot plants resulted in the identification of several novel features. The length of tRNAs in the chloroplast genome of the monocot plants ranged from 59 to 155 nucleotides. Pair-wise sequence alignment revealed the presence of a conserved A-C-x-U-A-x-U-A-x-U-x5-U-A-A nucleotide consensus sequence. In addition, the tRNAs in chloroplast genomes of the monocot plants also contain 21-28 anti-codons against 61 sense codons in the genome. They also contain a group I intron and a C-A-U anti-codon for tRNAIle, which is a common anti-codon of tRNAMet. Evolutionary analysis indicates that tRNAs in the chloroplast genome have evolved from multiple common ancestors, and tRNAMet appears to be the ancestral tRNA that underwent duplication and diversification to give rise to other tRNAs.

CONCLUSION

The results obtained from the study of chloroplast tRNA will greatly help to increase our understanding of tRNA biology at a new level. Functional studies of the reported novel aspects of the chloroplast tRNA of the monocot plants will greatly help to decipher their roles in diverse cellular processes.

Genetic mutation screen in early non--small-cell lung cancer (NSCLC) specimens

BACKGROUND

Testing for genetic abnormalities in epithelial growth factor receptor (EGFR), anaplastic lymphoma receptor tyrosine kinase (ALK), and potentially additional genes is a critical tool in the care of advanced NSCLC. There is conflicting evidence for the role of such tests in early NSCLC. We report a single-institute Sequenom testing for a wide range of mutations and their clinical correlations in early-resected NSCLC specimens.

MATERIALS AND METHODS

Early NSCLC paraffin-embedded, formalin-fixed (FFPE) specimens were collected, DNA extracted, and using Sequenom-based matrix-assisted laser desorption/ionization-time of flight analysis, mutations in 22 oncogenes and tumor suppressor genes were evaluated. Clinical data was collected retrospectively.

RESULTS

The technique was found to be feasible. Thirty-six of 96 patients (37.5%) had any genetic abnormality identified, and 8 (8.3%) had 2 or more mutations. Kirsten rat sarcoma viral oncogene homolog (KRAS) and EGFR were the most common genes to appear mutated (15.6%); phosphatidylinositol-4,5-bisphosphate 3-kinase, catalytic subunit alpha (PIK3CA) was the gene to be found most commonly in tumors with co-mutations. Transversions were found mostly in KRAS gene mutations and to be nonprognostic. No difference in the spectrum of mutations was found between squamous-cell and non-squamous-cell lung cancers. Ever-smokers showed a trend for worse prognosis, with a similar spectrum of mutations.

CONCLUSION

Sequenom-based mutation screen is feasible using FFPE samples. More than a third of the patients were found to harbor some genetic abnormality, and 8% were found to have more than a single mutated gene. Wide-range gene screens using large sample depositories are required for further insight into the important genes at play in early NSCLC.

Quantitative characterization of hepatitis delta virus genome edition by next-generation sequencing

AIM

To determine the capacity of next-generation sequencing (NGS) for quantifying edited and unedited HDV populations, and to confirm if edition is a general phenomenon taking place along the entire HDV region analyzed, as we previously reported (Homs M et al. PLoS One 2016, 11, e0158557).

METHODS

Four serum samples from 4 patients with chronic HDV/HBV infection were included in the study. The region selected for analysis covered 360 nucleotides (nt), positions 910-1270 of the HDV genome, which included the HDAg ORF editing site (nt 1014 within codon 196). Quantification of edited and unedited genomes was performed by molecular cloning and Sanger sequencing and by NGS. To evaluate the reliability of the NGS values obtained, we combined a clone with an edited codon and one with an unedited codon in known percentages in a series of artificial mixtures, which were then analyzed by NGS. In addition, we determined the nt changes occurring over the complete amplified region after excluding the editing codon (196) to evaluate edition along it.

RESULTS

In total, 11,208 quality-filtered sequences were obtained in the 4 samples. The 95% confidence intervals for the proportions of unedited populations by molecular cloning and NGS were overlapping, and those of cloning were wider, indicating that they are comparable and that NGS is more precise than cloning. Unedited genomes predominated over edited ones in all 4 samples analyzed by NGS and in 3 of the 4 samples analyzed by molecular cloning. In total, 83,276 quality-filtered sequences were obtained from the artificial mixtures. Percentages of the two viral populations detected by NGS in these mixtures were comparable to the expected percentages. Evaluation of edition along the HDV coding region showed that transitions were more frequent than transversions, accounting for 63.09% and 36.91%, respectively. Interestingly, among the 4 possible transition-type changes, G:A and A:G accounted for 73.86% of the total.

CONCLUSION

Next-generation sequencing proved useful to quantify edited and unedited HDV genomes, and provided relevant information on the HDV quasispecies.

Single nucleotide polymorphisms at heat shock protein 90 gene and their association with thermo-tolerance potential in selected indigenous Nigerian cattle

Heat shock protein (HSP) 90 gene provides protection and adaptation to thermal assault and certain polymorphisms have been associated to heat tolerance in humans and animals. Single nucleotide polymorphisms (SNPs) of HSP 90 gene were used to evaluate the scientific basis of heat tolerance in four zebu breeds of Nigeria. The DNA was extracted from skin tissue of 90 adult bulls representing White Fulani (WF), Sokoto Gudali (SG), Red Bororo (RB), and Ambala (AM). The SNPs were determined in DNAs using PCR, sequencing, and visualization and bio-editing by chromatogram in SeqMan Ngen tool. Subsequently, respective genotypes were constructed and genotypic and allelic frequencies were computed. Also, body parameters related to heat stress (HS) including body temperature (BT), rectal temperature (RT), and respiratory rates (RR) were taken for each animal before biological sampling and heat tolerance coefficient (HTC) was calculated. We detected four SNPs distinct/specific for each breed as follows: change from thymine (T) to guanine (G) at position 116 (T116G) in RB, G to cytosine (C) at 220 (G220C) in SG, G to adenine (A) at two positions, 346 (G346A) and 390 (G390A) in AM and WF, respectively. Heterozygous SNPs showed significantly lower values (P < 0.0001) for BT, RT, RR, and HTC than homozygous genotypes at all positions. We hypothesize that animals with heterozygous SNPs in exon 3 of HSP 90 may be tolerant to HS. These SNPs can be used as bio-markers for screening large populations of cattle for tolerance to hot tropical conditions in Nigeria and other sub-humid places.

A compound heterozygote case of isolated sulfite oxidase deficiency

We report an isolated sulfite oxidase deficiency in the first child boy of a non-consanguineous Caucasian family. He's a compound heterozygote for the sulfite oxidase gene, presenting low cystine, undetectable homocysteine and normal uric acid blood concentrations and undetectable sulfite oxidase activity in his cultured fibroblasts. Both mutations are not reported yet. The clinical presentation was typical and severe, with generalized status epilepticus and premature death.

Stem Region of tRNA Genes Favors Transition Substitution Towards Keto Bases in Bacteria

Transversion and transition mutations have variable effects on the stability of RNA secondary structure considering that the former destabilizes the double helix geometry to a greater extent by introducing purine:purine (R:R) or pyrimidine:pyrimidine (Y:Y) base pairs. Therefore, transversion frequency is likely to be lower than that of transition in the secondary structure regions of RNA genes. Here, we performed an analysis of transition and transversion frequencies in tRNA genes defined well with secondary structure and compared with the intergenic regions in five bacterial species namely Escherichia coli, Klebsiella pneumoniae, Salmonella enterica, Staphylococcus aureus and Streptococcus pneumoniae using a large genome sequence data set. In general, the transversion frequency was observed to be lower than that of transition in both tRNA genes and intergenic regions. The transition to transversion ratio was observed to be greater in tRNA genes than that in the intergenic regions in all the five bacteria that we studied. Interestingly, the intraspecies base substitution analysis in tRNA genes revealed that non-compensatory substitutions were more frequent than compensatory substitutions in the stem region. Further, transition to transversion ratio in the loop region was observed to be significantly lesser than that among the non-compensatory substitutions in the stem region. This indicated that the transversion is more deleterious than transition in the stem regions. In addition, substitutions from amino bases (A/C) to keto bases (G/T) were also observed to be more than the reverse substitutions in the stem region. Substitution from amino bases to keto bases are likely to facilitate the stable G:U pairing unlike the reverse substitution that facilitates the unstable A:C pairing in the stem region of tRNA. This work provides additional support that the secondary structure of tRNA molecule is what drives the different substitutions in its gene sequence.

Selection signatures in melanocortin-1 receptor gene of turkeys (Meleagris gallopavo) raised in hot humid tropics

Feather colours are used by avian species for defense, adaptation and signaling. Melanocortin-1 receptor (MC1R) gene is one of the genes responsible for feather colour. This study identified selection signatures in MC1R gene of Nigerian indigenous turkeys (NIT) using British United turkeys (BUT) as control breed to investigate the evolutionary processes that have shaped NIT with various feather colours. Complete MC1R gene of 146 NIT (76 males and 70 females) and 32 BUT (18 males and 14 females) were sequenced. Transition/transversion and codon usage biases were predicted using MEGA v6 software. The selective force acting on the gene was predicted using HyPhy software. The FST values were estimated using Arlequin v3.5. The highest transition/transversion bias was predicted for white BUT (1.00) while the lowest was predicted for black NIT (0.50). Negative dN-dS values, indicative of purifying selection, were observed in MC1R gene of all the turkeys. The highest pairwise FST was observed between the MC1R gene of white BUT and black NIT while the least was observed between lavender NIT and white NIT. No recombination event was observed in black NIT and white BUT. The relative synonymous codon usage was the same among different colours for some codons. Presence of purifying selection in MC1R gene of all the turkeys with different feather colours confirms that the gene plays role in many biological processes such as feather colouration, behaviour, pain perception, immunity, growth and adaptation. The results also suggested that the genetic mechanisms generating different feather colours in turkeys are conserved.

High frequency of transition to transversion ratio in the stem region of RNA secondary structure of untranslated region of SARS-CoV-2

Introduction

The propensity of nucleotide bases to form pairs, causes folding and the formation of secondary structure in the RNA. Therefore, purine (R): pyrimidine (Y) base-pairing is vital to maintain uniform lateral dimension in RNA secondary structure. Transversions or base substitutions between R and Y bases, are more detrimental to the stability of RNA secondary structure, than transitions derived from substitutions between A and G or C and T. The study of transversion and transition base substitutions is important to understand evolutionary mechanisms of RNA secondary structure in the 5'  and 3'  untranslated (UTR) regions of SARS-CoV-2. In this work, we carried out comparative analysis of transition and transversion base substitutions in the stem and loop regions of RNA secondary structure of SARS-CoV-2.

Methods

We have considered the experimentally determined and well documented stem and loop regions of 5' and 3' UTR regions of SARS-CoV-2 for base substitution analysis. The secondary structure comprising of stem and loop regions were visualized using the RNAfold web server. The GISAID repository was used to extract base sequence alignment of the UTR regions. Python scripts were developed for comparative analysis of transversion and transition frequencies in the stem and the loop regions.

Results

The results of base substitution analysis revealed a higher transition (ti) to transversion (tv) ratio (ti/tv) in the stem region of UTR of RNA secondary structure of SARS-CoV-2 reported during the early stage of the pandemic. The higher ti/tv ratio in the stem region suggested the influence of secondary structure in selecting the pattern of base substitutions. This differential pattern of ti/tv values between stem and loop regions was not observed among the Delta and Omicron variants that dominated the later stage of the pandemic. It is noteworthy that the ti/tv values in the stem and loop regions were similar among the later dominant Delta and Omicron variant strains which is to be investigated to understand the rapid evolution and global adaptation of SARS-CoV-2.

Conclusion

Our findings implicate the lower frequency of transversions than the transitions in the stem regions of UTRs of SARS-CoV-2. The RNA secondary structures are associated with replication, translation, and packaging, further investigations are needed to understand these base substitutions across different variants of SARS-CoV-2.

Detection of Transversions and Transitions in HBG2 Cis-Elements Associated with Sickle Cell Allele in Ghanaians

Short tandem repeats located 5' prime to the β-globin gene, have been observed to be in linkage disequilibrium with the HbS allele, and thought to affect the severity of sickle cell disease. Here, we report on new mutants within the HBG2 region that may impact sickle cell disease. To determine the cis-acting elements microsatellites, indels and single nucleotide polymorphisms (SNPs), within the HBG2 region by sequencing, in subjects with sickle cell disease. The case-control study was located at the Center for Clinical Genetics, Sickle cell unit, Korle-Bu Teaching Hospital. A questionnaire was used for demographic data and clinical information. Hematological profile (red blood cell, white blood cell, platelet, hemoglobin and mean corpuscular volume) were assessed in 83 subjects. A set of 45 samples comprising amplified DNA on the HBG2 gene from HbSS (22), HbSC (17) and 6 controls (HbAA) were sequenced. Differences in the microsatellite region between sickle cell disease (SCD) (HbSS and HbSC) genotypes and control subjects were identified by counting and assessed by Chi-square analysis. Red blood cells, hematocrit, platelets, white blood cells and hemoglobin indices differed in genotypic groups. HbSS subjects were affirmed to have severer hemolytic anemia than HbSC subjects. Two indels (T1824 and C905) were seen in both SS and SC genotypes. Two peculiar SNPs: G:T1860 (transition) and A:G1872 transversions were found within the HBG2 gene that were significantly associated with the HbSS genotype (Fisher's exact test, p = 0.006) and HbS allele respectively (Fisher's exact test, p = 0.006). Cis-acting elements in HbSS and HbSC were different and may contribute to the phenotype seen in the disease state.