Single Non-synonymous Mutation Research Articles

The Role of BrKS in Leafy Head Formation Was Confirmed by Two Allelic Mutants of Chinese Cabbage (Brassica rapa L. ssp. pekinensis)

Considering that the presence of a leafy head is a key morphological characteristic that determines the yield and quality of Chinese cabbage, identifying the major genes regulating the formation of a leafy head is crucial for variety improvement. A heading-related gene, BrKS, was previously predicted from a non-heading mutant, nhm1, derived from a heading variety, which encodes a key enzyme involved in gibberellin biosynthesis. Here, another mutant, nhm2, was identified from the same EMS-mutagenized population, and the phenotype of nhm2 was consistent with that of nhm1. We crossed mutants nhm1 and nhm2, and their F1 plants exhibited the mutant phenotype, which indicated that their mutant genes were allelic. A single non-synonymous mutation in the fourth exon of BrKS in mutant nhm1 and another single non-synonymous mutation in the tenth exon of BrKS in mutant nhm2. The same gene mutation in mutants nhm1 and nhm2 produced a similar non-heading phenotype, which confirmed the role of BrKS in the leafy head formation of Chinese cabbage. RNA-Seq analysis indicated that a transcription factor gene, BrERF1A, which is associated with leaf development, significantly down-regulated expression in mutant nhm1, and after the mutant was treated with GA3, the expression level of BrERF1A was recovered, which indicated that BrKS might be involved in leafy head formation through regulating the expression level of BrERF1A. Our findings provide important clues for revealing the molecular mechanism of leafy head formation in Chinese cabbage.

Deciphering the Genomic Characterization of the GGP Gene Family and Expression Verification of CmGGP1 Modulating Ascorbic Acid Biosynthesis in Melon Plants.

Ascorbic acid (AsA), also known as vitamin C, is a well-known antioxidant found in living entities that plays an essential role in growth and development, as well as in defensive mechanisms. GDP-L-galactose phosphorylase (GGP) is a candidate gene regulating AsA biosynthesis at the translational and transcriptional levels in plants. In the current study, we conducted genome-wide bioinformatic analysis and pinpointed a single AsA synthesis rate-limiting enzyme gene in melon (CmGGP1). The protein prediction analysis depicted that the CmGGP1 protein does not have a signaling peptide or transmembrane structure and mainly functions in the chloroplast or nucleus. The constructed phylogenetic tree analysis in multispecies showed that the CmGGP1 protein has a highly conserved motif in cucurbit crops. The structural variation analysis of the CmGGP1 gene in different domesticated melon germplasms showed a single non-synonymous type-base mutation and indicated that this gene was selected by domestication during evolution. Wild-type (WT) and landrace (LDR) germplasms of melon depicted close relationships to each other, and improved-type (IMP) varieties showed modern domestication selection. The endogenous quantification of AsA content in both the young and old leaves of nine melon varieties exhibited the major differentiations for AsA synthesis and metabolism. The real-time quantitative polymerase chain reaction (qRT-PCR) analysis of gene co-expression showed that AsA biosynthesis in leaves was greater than AsA metabolic consumption, and four putative interactive genes (MELO3C025552.2, MELO3C007440.2, MELO3C023324.2, and MELO3C018576.2) associated with the CmGGP1 gene were revealed. Meanwhile, the CmGGP1 gene expression pattern was noticed to be up-regulated to varying degrees in different acclimated melons. We believe that the obtained results would provide useful insights for an in-depth genetic understanding of the AsA biosynthesis mechanism, aimed at the development of improving crop plants for melon.

In silico structural and mechanical insights into bedaquiline resistance associated with high-grade non-synonymous mutations in atpE, mmpR5, and pepQ

Clinical resistance against bedaquiline (BDQ) remains intractable to anti-tuberculosis therapies since its introduction to the market over a decade ago. Herein, we investigated the structural and mechanical aspects of BDQ resistance in AtpE, MmpR5, and PepQ. The known target-specific resistant single non-synonymous mutations were refined to high-grade candidates. Thus, 7 (AtpE), 5 (MmpR5), and 1 (PepQ) single nucleotide polymorphisms (SNPs) and one insertion frameshift mutation in MmpR5 were recreated at the molecular level, and these phenotypic models were then directed to stringent dynamics to define time-scaled changes. The AtpE variants destabilized the structure; mainly, L59V, E61D, and I66M were detrimental to the complex fitness, while L74V and L114P boosted the BDQ binding to MmpR5. The first three and last two alterations gave rise to loss- and gain-of-function to AtpE and MmpR5, respectively. Hence, these five mutants are functionally relevant and therapeutically targetable hotspots of BDQ resistance. There were no noticeable changes in PepQ data analysis. The present study revealed that MmpR5 mutations confer BDQ resistance, whereas AtpE and PepQ SNPs display low susceptibility. These results were tallied with the published findings, which testified to the pursued method’s reliability and accuracy. We hope these data and inferences could be helpful for the futuristic design of novel TB drugs. Communicated by Ramaswamy H. Sarma

The novel HLA-B*51:328 allele identified by sequencing-based typing in a Chinese individual.

HLA-B*51:328 differs from HLA-B*51:01:01:01 in exon 2 at position 241 by a single non-synonymous mutation.

A single nonsynonymous mutation on ZIKV E protein-coding sequences leads to markedly increased neurovirulence in vivo.

Zika virus (ZIKV) can infect a wide range of tissues including the developmental brain of human fetus. Whether specific viral genetic variants are linked to neuropathology is incompletely understood. To address this, we have intracranially serially passaged a clinical ZIKV isolate (SW01) in neonatal mice and discovered variants that exhibit markedly increased virulence and neurotropism. Deep sequencing analysis combining with molecular virology studies revealed that a single 67D (Aspartic acid) to N (Asparagine) substitution on E protein is sufficient to confer the increased virulence and neurotropism in vivo. Notably, virus clones with D67N mutation had higher viral production and caused more severe cytopathic effect (CPE) in human neural astrocytes U251 ​cells in vitro, indicating its potential neurological toxicity to human brain. These findings revealed that a single mutation D67N on ZIKV envelope may lead to severe neuro lesion that may help to explain the neurovirulence of ZIKV and suggest monitoring the occurrence of this mutation during nature infection may be important.

Identification and Fine Mapping of the Recessive Gene BK-5, Which Affects Cell Wall Biosynthesis and Plant Brittleness in Maize.

The cellulose of the plant cell wall indirectly affects the cell shape and straw stiffness of the plant. Here, the novel brittleness mutant brittle stalk-5 (bk-5) of the maize inbred line RP125 was characterized. We found that the mutant displayed brittleness of the stalk and even the whole plant, and that the brittleness phenotype existed during the whole growth period from germination to senescence. The compressive strength was reduced, the cell wall was thinner, and the cellulose content was decreased compared to that of the wild type. Genetic analysis and map-based cloning indicated that bk-5 was controlled by a single recessive nuclear gene and that it was located in a 90.2-Kb region on chromosome 3 that covers three open reading frames (ORFs). Sequence analysis revealed a single non-synonymous missense mutation, T-to-A, in the last exon of Zm00001d043477 (B73: version 4, named BK-5) that caused the 951th amino acid to go from leucine to histidine. BK-5 encodes a cellulose synthase catalytic subunit (CesA), which is involved with cellulose synthesis. We found that BK-5 was constitutively expressed in all tissues of the germinating stage and silking stage, and highly expressed in the leaf, auricula, and root of the silking stage and the 2-cm root and bud of the germinating stage. We found that BK-5 mainly localized to the Golgi apparatus, suggesting that the protein might move to the plasma membrane with the aid of Golgi in maize. According to RNA-seq data, bk-5 had more downregulated genes than upregulated genes, and many of the downregulated genes were enzymes and transcription factors related to cellulose, hemicellulose, and lignin biosynthesis of the secondary cell wall. The other differentially expressed genes were related to metabolic and cellular processes, and were significantly enriched in hormone signal transduction, starch and sucrose metabolism, and the plant–pathogen interaction pathway. Taken together, we propose that the mutation of gene BK-5 causes the brittle stalk phenotype and provides important insights into the regulatory mechanism of cellulose biosynthesis and cell wall development in maize.

Optimal flexibility of the linker region of Zika virus NS5 methyltransferase-polymerase is critical for virus replication

The flavivirus NS5 protein contains an N-terminal methyl-transferase (MTase) connected through a flexible linker with a C-terminal RNA-dependent RNA-polymerase (RdRp) domain, that work cooperatively to replicate and methylate the viral genome. In this study we probed the importance of an evolutionary-conserved hydrophobic residue (Val266) located at the start of the ten-residue interdomain linker of Zika virus (ZIKV) NS5. In flavivirus NS5 crystal structures, the start of the linker forms a 310 helix when NS5 adopts a compact conformation, but becomes disordered or extended in open conformations. Using reverse genetics system, we either introduced rigidity in the linker through mutation to a proline or flexibility through a glycine mutation at position 266. ZIKV NS5 Val 266 to Pro mutation was lethal for viral RNA replication while the Gly mutation was severely attenuated. Serial passaging of cell culture supernatant derived from C6/36 mosquito cells transfected with mutant ZIKV RNA showed that the attenuation can be rescued. Next generation deep sequencing revealed four single nucleotide polymorphisms that occur with an allele frequency >98%. The single non-synonymous NS5 mutation Glu419 to Lys is adjacent to RdRp motif G at the tip of the fingers subdomain, while the remaining three are synonymous variants at nucleotide positions 1403, 4403 and 6653 in the genome. Reverse engineering the changes into the ZIKV NS5/Val266Gly background followed by serial passaging revealed that residue 266 is under strong positive selection to revert back to Val. The interaction of the specific conformation of the NS5 linker with Val at position 266 and the RNA binding motif G region may present a potential strategy for allosteric antiviral drug development.

Genetic and Transcriptomic Analysis Reveal the Molecular Basis of Photoperiod-Regulated Flowering in Xishuangbanna Cucumber (Cucumis sativus L. var. xishuangbannesis Qi et Yuan).

Xishuangbanna (XIS) cucumber (Cucumis sativus L. var. xishuangbannesis Qi et Yuan), is a botanical variety of cucumber cultivars native to southwest China that possesses excellent agronomic traits for cucumber improvement. However, breeding utilization of XIS cucumber is limited due to the current poor understanding of its photoperiod-sensitive flowering characteristics. In this study, genetic and transcriptomic analysis were conducted to reveal the molecular basis of photoperiod-regulated flowering in XIS cucumber. A major-effect QTL locus DFF1.1 was identified that controls the days to first flowering (DFF) of XIS cucumbers with a span of 1.38 Mb. Whole-genome re-sequencing data of 9 cucumber varieties with different flowering characteristics in response to photoperiod suggested that CsaNFYA1 was the candidate gene of DFF1.1, which harbored a single non-synonymous mutation in its fifth exon. Transcriptomic analysis revealed the positive roles of auxin and ethylene in accelerating flowering under short-day (SD) light-dark cycles when compared with equal-day/night treatment. Carbohydrate storage and high expression levels of related genes were important reasons explaining early flowering of XIS cucumber under SD conditions. By combining with the RNA-Seq data, the co-expression network suggested that CsaNFYA1 integrated multiple types of genes to regulate the flowering of XIS cucumber. Our findings explain the internal regulatory mechanisms of a photoperiodic flowering pathway. These findings may guide the use of photoperiod shifts to promote flowering of photoperiod-sensitive crops.

A novel HLA-C allele, HLA-C*14:125.

HLA-C*14:125 allele is identical to HLA-C*14:02:01:01 except for a single nonsynonymous mutation C199T.

A novel HLA‐A*02 allele, A*02:828, identified by next‐generation sequencing in a Chinese individual

HLA-A*02:828 differs from A*02:01:01:01 by a single non-synonymous mutation in exon 2.

Altered expression of K13 disrupts DNA replication and repair in Plasmodium falciparum

BackgroundPlasmodium falciparum exhibits resistance to the artemisinin component of the frontline antimalarial treatment Artemisinin-based Combination Therapy in South East Asia. Millions of lives will be at risk if artemisinin resistance (ART-R) spreads to Africa. Single non-synonymous mutations in the propeller region of PF3D7_1343700,“K13” are implicated in resistance. In this work, we use transcriptional profiling to characterize a laboratory-generated k13 insertional mutant previously demonstrated to have increased sensitivity to artemisinins to explore the functional role of k13.ResultsA set of RNA-seq and microarray experiments confirmed that the expression profile of k13 is specifically altered during the early ring and early trophozoite stages of the mutant intraerythrocytic development cycle. The down-regulation of k13 transcripts in this mutant during the early ring stage is associated with a transcriptome advance towards a more trophozoite-like state. To discover the specific downstream effect of k13 dysregulation, we developed a new computational method to search for differential gene expression while accounting for the temporal sequence of transcription. We found that the strongest biological signature of the transcriptome shift is an up-regulation of DNA replication and repair genes during the early ring developmental stage and a down-regulation of DNA replication and repair genes during the early trophozoite stage; by contrast, the expressions of housekeeping genes are unchanged. This effect, due to k13 dysregulation, is antagonistic, such that k13 levels are negatively correlated with DNA replication and repair gene expression.ConclusionOur results support a role for k13 as a stress response regulator consistent with the hypothesis that artemisinins mode of action is oxidative stress and k13 as a functional homolog of Keap1 which in humans regulates DNA replication and repair genes in response to oxidative stress.

A mutation in class III homeodomain-leucine zipper (HD-ZIP III) transcription factor results in curly leaf (cul) in cucumber (Cucumis sativus L.).

We identified two curly-leaf (cul) mutants in cucumber. Map-based cloning revealed that both mutants are due to allelic mutations in the CsPHB gene, a homolog of the Arabidopsis PHABULOSA which encodes a class III homeodomain-leucine zipper (HD-ZIP III) transcription factor. Leaf rolling is an important agronomic trait in crop breeding. Moderate leaf rolling minimizes shadowing between leaves, leading to improved photosynthetic efficiency. Although a number of genes controlling rolled leaf have been identified from rice and other plant species, none have been mapped or cloned in cucurbit crops. In this study, we identified and characterized two curly leaf (cul) mutants, cul-1 and cul-2 in cucumber. With map-based cloning, we show that cul-1 and cul-2 are allelic mutations and CsPHB (Csa6G525430) was the candidate gene for both mutants. The CsPHB gene encoded a class III homeodomain-leucine zipper (HD-ZIP III) transcription factor. A single non-synonymous mutation in the fourth and fifth exons of the CsPHB was responsible for the cul-1 and cul-2 mutant phenotypes, respectively. The single-nucleotide substitutions in cul-1 and cul-2 were both located in cs-miRNA165/166 complementary sites of CsPHB. The expression level of CsPHB gene in multiple organs of cul-1 and cul-2 mutants was higher than that in the wild type, while the expression of cs-miRNA165/166 in the two genotypes showed the opposite trend. We speculate that disruption of the binding between the mutant allele of CsPHB and cs-miRNA165/166 leads to the curly-leaf phenotype. This is the first report to clone and characterize the CsPHB gene in the family Cucurbitaceae. Taken together, these results support CsPHB as an important player in the modulation of leaf shape development in cucumber.

Description of the novel KIR2DL4*039 allele identified in a southern Chinese Han individual.

The novel KIR2DL4*039 allele differs from the closest allele 2DL4*00602 by a single non-synonymous mutation.

Single-mutation fitness landscapes for an enzyme on multiple substrates reveal specificity is globally encoded

Our lack of total understanding of the intricacies of how enzymes behave has constrained our ability to robustly engineer substrate specificity. Furthermore, the mechanisms of natural evolution leading to improved or novel substrate specificities are not wholly defined. Here we generate near-comprehensive single-mutation fitness landscapes comprising >96.3% of all possible single nonsynonymous mutations for hydrolysis activity of an amidase expressed in E. coli with three different substrates. For all three selections, we find that the distribution of beneficial mutations can be described as exponential, supporting a current hypothesis for adaptive molecular evolution. Beneficial mutations in one selection have essentially no correlation with fitness for other selections and are dispersed throughout the protein sequence and structure. Our results further demonstrate the dependence of local fitness landscapes on substrate identity and provide an example of globally distributed sequence-specificity determinants for an enzyme.

Identification and characterization of the novel KIR2DL1*030 allele by sequence-based typing in a southern Chinese Han individual.

KIR2DL1*030 differs from KIR2DL1*00302 by a single non-synonymous mutation in exon 4.

Identification of the novel KIR2DL1*033 allele from a southern Chinese Han individual by sequencing-based typing.

The novel KIR2DL1*033 allele differs from the closest allele 2DL1*00302 by a single nonsynonymous mutation.

Inhibition of nervous necrosis virus by ribavirin in a zebrafish larvae model

The guanosine analog ribavirin is a broad-spectrum antiviral drug, mostly used in human clinical practice. It has in vitro and in vivo activity against a broad range of RNA and DNA viruses. Here, we report that treatment of zebrafish larvae with ribavirin prior to infection with nervous necrosis virus (NNV) significantly reduces the mortality caused by the virus during the first 10 days post-infection. The RNA genome of NNV harvested from ribavirin-treated infected larvae contains three synonymous and one single non-synonymous mutation, resulting in the replacement of a serine codon with a glycine codon in the RNA-dependent RNA polymerase gene. Adding increasing amounts of guanosine to ribavirin prior to larvae infection did not impede the antiviral activity. Ribavirin treatment of uninfected larvae reduces the basal level of IFNγ, but increases the level of IL-1β mRNA expression. Furthermore, infecting larvae with NNV following ribavirin treatment reduces the expression levels of IFNγ, IFN-I, Mx, and TNF-α genes, while the expression of IL-1β is increased. These results suggest that cytokine modulation plays an important role in the activity of ribavirin against NNV. Mortality of more than 40 species of teleost fish, mostly larvae and juveniles, from NNV is a major obstacle in hatcheries, and impedes the supply of young fish to farms. Hence, cost-effective ribavirin treatment should be considered as an efficient means to reduce the peril of NNV.

A CRISPR/Cas9 mediated point mutation in the alpha 6 subunit of the nicotinic acetylcholine receptor confers resistance to spinosad in Drosophila melanogaster

Spinosad, a widely used and economically important insecticide, targets the nicotinic acetylcholine receptor (nAChRs) of the insect nervous system. Several studies have associated loss of function mutations in the insect nAChR α6 subunit with resistance to spinosad, and in the process identified this particular subunit as the specific target site. More recently a single non-synonymous point mutation, that does not result in loss of function, was identified in spinosad resistant strains of three insect species that results in an amino acid substitution (G275E) of the nAChR α6 subunit. The causal role of this mutation has been called into question as, to date, functional evidence proving its involvement in resistance has been limited to the study of vertebrate receptors. Here we use the CRISPR/Cas9 gene editing platform to introduce the G275E mutation into the nAChR α6 subunit of Drosophila melanogaster. Reverse transcriptase-PCR and sequencing confirmed the presence of the mutation in Dα6 transcripts of mutant flies and verified that it does not disrupt the normal splicing of the two exons in close vicinity to the mutation site. A marked decrease in sensitivity to spinosad (66-fold) was observed in flies with the mutation compared to flies of the same genetic background minus the mutation, clearly demonstrating the functional role of this amino acid substitution in resistance to spinosad. Although the resistance levels observed are 4.7-fold lower than exhibited by a fly strain with a null mutation of Dα6, they are nevertheless predicated to be sufficient to result in resistance to spinosad at recommended field rates. Reciprocal crossings with susceptible fly strains followed by spinosad bioassays revealed G275E is inherited as an incompletely recessive trait, thus resembling the mode of inheritance described for this mutation in the western flower thrips, Frankliniella occidentalis. This study both resolves a debate on the functional significance of a target-site mutation and provides an example of how recent advances in genome editing can be harnessed to study insecticide resistance.

Molecular Basis for Resistance to Fluazifop-P-Butyl in Itchgrass (Rottboellia cochinchinensis) from Costa Rica

Rottboellia cochinchinensis is an annual grass weed species known as itchgrass, or "caminadora" in America´s Spanish speaking countries, and has become a major and troublesome weed in several crops. The application of fluazifop-P-butyl at recommended rates (125 g a.i. ha-1) was observed to be failing to control itchgrass in a field in San José, Upala county, Alajuela province, Costa Rica. Plants from the putative resistant R. cochinchinensis population survived fluazifop-P-butyl when treated with 250 g a.i. ha-1 (2X label rate) at the three- to four-leaf stage under greenhouse conditions. PCR amplification and sequencing of partial carboxyl transferase domain (CT) of the acetyl-CoA carboxylase (ACCase) gene were used to determine the molecular mechanism of resistance. A single non-synonymous point mutation from TGG (susceptible plants) to TGC (putative resistant plants) that leads to a Trp-2027-Cys substitution was found. This Trp-2027-Cys mutation is known to confer resistance to all aryloxyphenoxyproprionate (APP) herbicides to which fluazifop-P-butyl belongs. To the best of our knowledge, this is the first report of fluazifop-P-butyl resistance and a mutation at position 2027 for a Costa Rican R. cochinchinensis population.

Variability of non-structural proteins of equine arteritis virus during persistent infection of the stallion.

The genetic stability of ORF1a encoding non-structural proteins nsp1, nsp2, nsp3 and nsp4 of equine arteritis virus (EAV) has been analysed for nearly seven years in a persistently infected stallion of the Malopolska breed. Between November 2004 and June 2011, 11 semen samples were collected. Viral RNA extracted from semen of this carrier stallion was amplified, sequenced and compared with the sequences of the other known strains of EAV. Sequence analysis of ORF1a showed 84 synonymous and 16 non-synonymous mutations. The most variable part of ORF1a was the region encoding nsp2 protein with 13 non-synonymous substitutions. The degree of amino acid identity between isolates ranged from 98.91 to 100%. Only single non-synonymous mutations were detected in nsp1 (one substitution) and nsp4 (two substitutions). The most stable was nsp3 in which no amino acid substitutions were observed during the whole period of observation.