Identification of Three Novel Frameshift Mutations in the PKD1 Gene in Iranian Families with Autosomal Dominant Polycystic Kidney Disease Using Efficient Targeted Next-Generation Sequencing.

BACKGROUND/AIMS: Autosomal dominant polycystic kidney disease (ADPKD) is one of the most common inherited cystic kidney diseases caused by mutations in two large multi-exon genes, PKD1 and PKD2. High allelic heterogeneity and duplication of PKD1 exons 1-32 as six pseudo genes on chromosome 16 complicate molecular analysis of this disease. METHODS: We applied targeted next-generation sequencing (NGS) in 9 non-consanguineous unrelated Iranian families with ADPKD to identify the genes hosting disease-causing mutations. This approach was confirmed by Sanger sequencing. RESULTS: Here, we determined three different novel frameshift mutations and four previously reported nonsense mutations in the PKD1 gene encoding polycystin1 in heterozygotes. CONCLUSION: This study demonstrates the effectiveness of NGS in significantly reducing the cost and time for simultaneous sequence analysis of PKD1 and PKD2, simplifying the genetic diagnostics of ADPKD. Although a probable correlation between the mutation types and phenotypic outcome is possible, however for more extensive studies in future, the consideration of renal hypouricemia (RHUC) and PKD1 coexistence may be helpful. The novel frameshift mutations reported by this study are p. Q1997X, P. D73X and p. V336X.

Resistance risk assessment for fludioxonil in Bipolaris maydis.

Bipolaris maydis (anamorph: Cochliobolus heterostrophus) is the causal agent of Southern Corn Leaf Blight (SCLB), leading to huge annually losses worldwide. Although fludioxonil, a phenylpyrrole fungicide with a broad spectrum of activity, was introduced in the 1990s, no baseline sensitivity has been established for B. maydis. One hundred field isolates were used to establish a baseline sensitivity of B. maydis against fludioxonil during 2015-2016. The results showed that the baseline sensitivity was distributed as a unimodal curve with a mean EC50 value of 0.044±0.022µgmL-1. With repeated exposure to fludioxonil, a total of five fludioxonil-resistant mutants (RF>100, RF=Resistance factor) were obtained in the laboratory. Compared with the parental isolates, the five fludioxonil-resistant mutants showed decreased fitness in sporulation and virulence, and exhibited different features of sensitivity to various stresses (oxidation and osmotic pressure, cell membrane and cell wall inhibitors), but not in mycelial growth on PDA without stress amendation. The five fludioxonil-resistant mutants showed a positive cross-resistance between fludioxonil and the dicarboximide fungicide procymidone, but not between fludioxonil and boscalid or fluazinam. All mutants exhibited stable resistance to fludioxonil after 10 transfers, as indicated by resistance factor values that ranged from 116.82 to 445.59. When treated with 1.0 M NaCl, all the fludioxonil-resistant mutants showed greater mycelial glycerol content than corresponding parental isolates. Sequencing alignment results of Bmos1 indicated that mutant R27-5 had a single point mutation (Z1125K), while the mutant R104 had a 34-bp deletion fragment between the codons of amino acid residues 1125 to 1236 and encodes a putative attenuated 1133-AA protein. The 34-bp deletion fragment led to not only a 11-AA deletion(DNAVNQKLAVR), but also the resulting frameshift mutation and early stop. The mutations of R27-5 and R104 were located in the Rec domain of the Bmos1 gene. No mutations at the Bmos1 were detected in the other three resistant mutants R27-1, R27-2 and R32.

Natural selection retains overrepresented out-of-frame stop codons against frameshift peptides in prokaryotes.

BACKGROUND: Out-of-frame stop codons (OSCs) occur naturally in coding sequences of all organisms, providing a mechanism of early termination of translation in incorrect reading frame so that the metabolic cost associated with frameshift events can be reduced. Given such a functional significance, we expect statistically overrepresented OSCs in coding sequences as a result of a widespread selection. Accordingly, we examined available prokaryotic genomes to look for evidence of this selection. RESULTS: The complete genome sequences of 990 prokaryotes were obtained from NCBI GenBank. We found that low G+C content coding sequences contain significantly more OSCs and G+C content at specific codon positions were the principal determinants of OSC usage bias in the different reading frames. To investigate if there is overrepresentation of OSCs, we modeled the trinucleotide and hexanucleotide biases of the coding sequences using Markov models, and calculated the expected OSC frequencies for each organism using a Monte Carlo approach. More than 93% of 342 phylogenetically representative prokaryotic genomes contain excess OSCs. Interestingly the degree of OSC overrepresentation correlates positively with G+C content, which may represent a compensatory mechanism for the negative correlation of OSC frequency with G+C content. We extended the analysis using additional compositional bias models and showed that lower-order bias like codon usage and dipeptide bias could not explain the OSC overrepresentation. The degree of OSC overrepresentation was found to correlate negatively with the optimal growth temperature of the organism after correcting for the G+C% and AT skew of the coding sequence. CONCLUSIONS: The present study uses approaches with statistical rigor to show that OSC overrepresentation is a widespread phenomenon among prokaryotes. Our results support the hypothesis that OSCs carry functional significance and have been selected in the course of genome evolution to act against unintended frameshift occurrences. Some results also hint that OSC overrepresentation being a compensatory mechanism to make up for the decrease in OSCs in high G+C organisms, thus revealing the interplay between two different determinants of OSC frequency.

Unique cost dynamics elucidate the role of frameshifting errors in promoting translational robustness.

There is now considerable evidence supporting the view that codon usage is frequently under selection for translational accuracy. There are, however, multiple forms of inaccuracy (missense, premature termination, and frameshifting errors) and pinpointing a particular error process behind apparently adaptive mRNA anatomy is rarely straightforward. Understanding differences in the fitness costs associated with different types of translational error can help us devise critical tests that can implicate one error process to the exclusion of others. To this end, we present a model that captures distinct features of frameshifting cost and apply this to 641 prokaryotic genomes. We demonstrate that, although it is commonly assumed that the ribosome encounters an off-frame stop codon soon after the frameshift and costs of mis-elongation are therefore limited, genomes with high GC content typically incur much larger per-error costs. We go on to derive the prediction, unique to frameshifting errors, that differences in translational robustness between the 5' and 3' ends of genes should be less pronounced in genomes with higher GC content. This prediction we show to be correct. Surprisingly, this does not mean that GC-rich organisms necessarily carry a greater fitness burden as a consequence of accidental frameshifting. Indeed, increased per-error costs are often more than counterbalanced by lower predicted error rates owing to more diverse anticodon repertoires in GC-rich genomes. We therefore propose that selection on tRNA repertoires may operate to reduce frameshifting errors.

Amino acid translation program for full-length cDNA sequences with frameshift errors.

Here we present an amino acid translation program designed to suggest the position of experimental frameshift errors and predict amino acid sequences for full-length cDNA sequences having phred scores. Our program generates artificial insertions into artificial deletions from low-accuracy positions of the original sequence, thereby generating many candidate sequences. The validity of the most probable sequence (the likelihood that it represents the actual protein) is evaluated by using a score (V(a)) that is calculated in light of the Kozak consensus, preferred codon usage, and position of the initiation codon. To evaluate the software, we have used a database in which, out of 612 cDNA sequences, 524 (86%) carried 773 frameshift errors in the coding sequence. Our software detected and corrected 48% of the total frameshift errors in 62% of the total cDNA sequences with frameshift errors. The false positive rate of frameshift correction was 9%, and 91% of the suggested frameshifts were true.