Mitochondrial Molecular Adaptations and Life History Strategies Coevolve in Plants

Messenger RNA secondary structure prevents mutations at functionally important sites. Mutations at exposed sites would cause micro-adaptations, niche-specialization, and therefore, can be thought to promote K-strategists. Exposing, rather than protecting, conserved sites, is also potentially adaptive because they probably promote macro-adaptive changes. This presumably fits r-strategists: their population dynamics tolerate decreased survival. We found that helix-forming tendencies are greater at evolutionary conserved sites of plant mitochondrial mRNAs than at evolutionary variable sites in a majority (73%) of species–gene combinations. K-strategists preferentially protect conserved sites in short genes, r-strategists protect them most in larger genes. This adaptive scenario resembles our earlier findings in chloroplast genes. Protection levels at various codon positions also display disparity with respect to life history strategies of the plants. Conserved site protection increases overall mRNA folding stabilities for some genes, while decreases it for some others. This contrast exists between homologous genes of r- and K- strategists. Such compensating interactions between variability, mRNA size, codon position, and secondary structure factors within r- and K-strategists are most likely, molecular adaptations of plants belonging to the two extreme life history strategies. Our results suggest coevolution between molecular and ecological adaptive strategies

MutS recognition: multiple mismatches and sequence context effects

Escherichia coli MutS is a versatile repair protein that specifically recognizes not only various types of mismatches but also single stranded loops of up to 4 nucleotides in length. Specific binding, followed by the next step of tracking the DNA helix that locates hemi-methylated sites, is regulated by the conformational state of the protein as a function of ATP binding/hydrolysis. Here, we study how various molecular determinants of a heteroduplex regulate mismatch recognition by MutS, the critical first step of mismatch repair. Using classical DNase I footprinting assays, we demonstrate that the hierarchy of MutS binding to various types of mismatches is identical whether the mismatches are present singly or in multiples. Moreover, this unique hierarchy is indifferent both to the differential level of DNA helical flexibility and to the unpaired status of the mismatched bases in a heteroduplex. Surprisingly, multiple mismatches exhibit reduced affinity of binding to MutS, compared to that of a similar single mismatch. Such a reduction in the affinity might be due to sequence context effects, which we established more directly by studying two identical single mismatches in an altered sequence background. A mismatch, upon simply being flipped at the same location, elicits changes in MutS specific contacts, thereby underscoring the importance of sequence context in modulating MutS binding to mismatches

IMACULAT - an open access package for the quantitative analysis of chromosome localization in the nucleus

The alteration in the location of the chromosomes within the nucleus upon action of internal or external stimuli has been implicated in altering genome function. The effect of stimuli at a whole genome level is studied by using two-dimensional fluorescence in situ hybridization (FISH) to delineate whole chromosome territories within a cell nucleus, followed by a quantitative analysis of the spatial distribution of the chromosome. However, to the best of our knowledge, open access software capable of quantifying spatial distribution of whole chromosomes within cell nucleus is not available. In the current work, we present a software package that computes localization of whole chromosomes - Image Analysis of Chromosomes for computing localization (IMACULAT). We partition the nucleus into concentric elliptical compartments of equal area and the variance in the quantity of any chromosome in these shells is used to determine its localization in the nucleus. The images are pre-processed to remove the smudges outside the cell boundary. Automation allows high throughput analysis for deriving statistics. Proliferating normal human dermal fibroblasts were subjected to standard a two-dimensional FISH to delineate territories for all human chromosomes. Approximately 100 images from each chromosome were analyzed using IMACULAT. The analysis corroborated that these chromosome territories have non-random gene density based organization within the interphase nuclei of human fibroblasts. The ImageMagick Perl API has been used for pre-processing the images

The electrostatic profile of consecutive Cβ atoms applied to protein structure quality assessment.

The structure of a protein provides insight into its physiological interactions with other components of the cellular soup. Methods that predict putative structures from sequences typically yield multiple, closely-ranked possibilities. A critical component in the process is the model quality assessing program (MQAP), which selects the best candidate from this pool of structures. Here, we present a novel MQAP based on the physical properties of sidechain atoms. We propose a method for assessing the quality of protein structures based on the electrostatic potential difference (EPD) of Cβ atoms in consecutive residues. We demonstrate that the EPDs of Cβ atoms on consecutive residues provide unique signatures of the amino acid types. The EPD of Cβ atoms are learnt from a set of 1000 non-homologous protein structures with a resolution cuto of 1.6 Å obtained from the PISCES database. Based on the Boltzmann hypothesis that lower energy conformations are proportionately sampled more, and on Annsen's thermodynamic hypothesis that the native structure of a protein is the minimum free energy state, we hypothesize that the deviation of observed EPD values from the mean values obtained in the learning phase is minimized in the native structure. We achieved an average specificity of 0.91, 0.94 and 0.93 on hg_structal, 4state_reduced and ig_structal decoy sets, respectively, taken from the Decoys `R' Us database. The source code and manual is made available at https://github.com/sanchak/mqap and permanently available on 10.5281/zenodo.7134

Local repeat sequence organization of an intergenic spacer in the chloroplast genome of Chlamydomonas reinhardtii leads to DNA expansion and sequence scrambling: a complex mode of "copychoice replication"?

Parent-specific, randomly amplified polymorphic DNA (RAPD) markers were obtained from total genomic DNA ofChlamydomonas reinhardtii. Such parent-specific RAPD bands (genomic fingerprints) segregated uniparentally (through mt+) in a cross between a pair of polymorphic interfertile strains ofChlamydomonas (C. reinhardtii andC. minnesotti), suggesting that they originated from the chloroplast genome. Southern analysis mapped the RAPD-markers to the chloroplast genome. One of the RAPD-markers, "P2" (1.6 kb) was cloned, sequenced and was fine mapped to the 3 kb region encompassing 3′ end of 23S, full 5S and intergenic region between 5S and psbA. This region seems divergent enough between the two parents, such that a specific PCR designed for a parental specific chloroplast sequence within this region, amplified a marker in that parent only and not in the other, indicating the utility of RAPD-scan for locating the genomic regions of sequence divergence. Remarkably, the RAPD-product, "P2" seems to have originated from a PCR-amplification of a much smaller (about 600 bp), but highly repeat-rich (direct and inverted) domain of the 3 kb region in a manner that yielded no linear sequence alignment with its own template sequence. The amplification yielded the same uniquely "sequence-scrambled" product, whether the template used for PCR was total cellular DNA, chloroplast DNA or a plasmid clone DNA corresponding to that region. The PCR product, a "unique" new sequence, had lost the repetitive organization of the template genome where it had originated from and perhaps represented a "complex path" of copy-choice replication

Relationship between mRNA secondary structure and sequence variability in Chloroplast genes: possible life history implications

<p>Abstract</p> <p>Background</p> <p>Synonymous sites are freer to vary because of redundancy in genetic code. Messenger RNA secondary structure restricts this freedom, as revealed by previous findings in mitochondrial genes that mutations at third codon position nucleotides in helices are more selected against than those in loops. This motivated us to explore the constraints imposed by mRNA secondary structure on evolutionary variability at all codon positions in general, in chloroplast systems.</p> <p>Results</p> <p>We found that the evolutionary variability and intrinsic secondary structure stability of these sequences share an inverse relationship. Simulations of most likely single nucleotide evolution in <it>Psilotum nudum </it>and <it>Nephroselmis olivacea </it>mRNAs, indicate that helix-forming propensities of mutated mRNAs are greater than those of the natural mRNAs for short sequences and vice-versa for long sequences. Moreover, helix-forming propensity estimated by the percentage of total mRNA in helices increases gradually with mRNA length, saturating beyond 1000 nucleotides. Protection levels of functionally important sites vary across plants and proteins: <it>r</it>-strategists minimize mutation costs in large genes; <it>K</it>-strategists do the opposite.</p> <p>Conclusion</p> <p>Mrna length presumably predisposes shorter mRNAs to evolve under different constraints than longer mRNAs. The positive correlation between secondary structure protection and functional importance of sites suggests that some sites might be conserved due to packing-protection constraints at the nucleic acid level in addition to protein level constraints. Consequently, nucleic acid secondary structure <it>a priori </it>biases mutations. The converse (exposure of conserved sites) apparently occurs in a smaller number of cases, indicating a different evolutionary adaptive strategy in these plants. The differences between the protection levels of functionally important sites for <it>r</it>- and <it>K-</it>strategists reflect their respective molecular adaptive strategies. These converge with increasing domestication levels of <it>K</it>-strategists, perhaps because domestication increases reproductive output.</p

Structural phylogeny by profile extraction and multiple superimposition using electrostatic congruence as a discriminator

Phylogenetic analysis of proteins using multiple sequence alignment (MSA) assumes an underlying evolutionary relationship in these proteins which occasionally remains undetected due to considerable sequence divergence. Structural alignment programs have been developed to unravel such fuzzy relationships. However, none of these structure based methods have used electrostatic properties to discriminate between spatially equivalent residues. We present a methodology for MSA of a set of related proteins with known structures using electrostatic properties as an additional discriminator (STEEP). STEEP first extracts a profile, then generates a multiple structural superimposition providing a consolidated spatial framework for comparing residues and finally emits the MSA. Residues that are aligned differently by including or excluding electrostatic properties can be targeted by directed evolution experiments to transform the enzymatic properties of one protein into another. We have compared STEEP results to those obtained from a MSA program (ClustalW) and a structural alignment method (MUSTANG) for chymotrypsin serine proteases. Subsequently, we used PhyML to generate phylogenetic trees for the serine and metallo-&#946;-lactamase superfamilies from the STEEP generated MSA, and corroborated the accepted relationships in these superfamilies. We have observed that STEEP acts as a functional classifier when electrostatic congruence is used as a discriminator, and thus identifies potential targets for directed evolution experiments. In summary, STEEP is unique among phylogenetic methods for its ability to use electrostatic congruence to specify mutations that might be the source of the functional divergence in a protein family. Based on our results, we also hypothesize that the active site and its close vicinity contains enough information to infer the correct phylogeny for related proteins