Generalization of Entropy Based Divergence Measures for Symbolic Sequence Analysis

Entropy based measures have been frequently used in symbolic sequence analysis. A symmetrized and smoothed form of Kullback-Leibler divergence or relative entropy, the Jensen-Shannon divergence (JSD), is of particular interest because of its sharing properties with families of other divergence measures and its interpretability in different domains including statistical physics, information theory and mathematical statistics. The uniqueness and versatility of this measure arise because of a number of attributes including generalization to any number of probability distributions and association of weights to the distributions. Furthermore, its entropic formulation allows its generalization in different statistical frameworks, such as, non-extensive Tsallis statistics and higher order Markovian statistics. We revisit these generalizations and propose a new generalization of JSD in the integrated Tsallis and Markovian statistical framework. We show that this generalization can be interpreted in terms of mutual information. We also investigate the performance of different JSD generalizations in deconstructing chimeric DNA sequences assembled from bacterial genomes including that of E. coli, S. enterica typhi, Y. pestis and H. influenzae. Our results show that the JSD generalizations bring in more pronounced improvements when the sequences being compared are from phylogenetically proximal organisms, which are often difficult to distinguish because of their compositional similarity. While small but noticeable improvements were observed with the Tsallis statistical JSD generalization, relatively large improvements were observed with the Markovian generalization. In contrast, the proposed Tsallis-Markovian generalization yielded more pronounced improvements relative to the Tsallis and Markovian generalizations, specifically when the sequences being compared arose from phylogenetically proximal organisms.     

Accelerating public sector rice breeding with high-density KASP markers derived from whole genome sequencing of <Emphasis Type="Italic">indica</Emphasis> rice

Few public sector rice breeders have the capacity to use NGS-derived markers in their breeding programmes despite rapidly expanding repositories of rice genome sequence data. They rely on >?18,000 mapped microsatellites (SSRs) for marker-assisted selection (MAS) using gel analysis. Lack of knowledge about target SNP and InDel variant loci has hampered the uptake by many breeders of Kompetitive allele-specific PCR (KASP), a proprietary technology of LGC genomics that can distinguish alleles at variant loci. KASP is a cost-effective single-step genotyping technology, cheaper than SSRs and more flexible than genotyping by sequencing (GBS) or array-based genotyping when used in selection programmes. Before this study, there were 2015 rice KASP marker loci in the public domain, mainly identified by array-based screening, leaving large proportions of the rice genome with no KASP coverage. Here we have addressed the urgent need for a wide choice of appropriate rice KASP assays and demonstrated that NGS can detect many more KASP to give full genome coverage. Through re-sequencing of nine indica rice breeding lines or released varieties, this study has identified 2.5 million variant sites. Stringent filtering of variants generated 1.3 million potential KASP assay designs, including 92,500 potential functional markers. This strategy delivers a 650-fold increase in potential selectable KASP markers at a density of 3.1 per 1 kb in the indica crosses analysed and 377,178 polymorphic KASP design sites on average per cross. This knowledge is available to breeders and has been utilised to improve the efficiency of public sector breeding in Nepal, enabling identification of polymorphic KASP at any region or quantitative trait loci in relevant crosses. Validation of 39 new KASP was carried out by genotyping progeny from a range of crosses to show that they detected segregating alleles. The new KASP have replaced SSRs to aid trait selection during marker-assisted backcrossing in these crosses, where target traits include rice blast and BLB resistance loci. Furthermore, we provide the software for plant breeders to generate KASP designs from their own datasets.     

Multiple signaling systems target a core set of transition metal homeostasis genes using similar binding motifs

Bacterial response to metals can require complex regulation. We report an overlapping regulation for copper and zinc resistance genes in the denitrifying bacterium, Pseudomonas stutzeri RCH2, by three two‐component regulatory proteins CopR1, CopR2 and CzcR. We conducted genome‐wide evaluations to identify gene targets of two paralogous regulators, CopR1 and CopR2, annotated for copper signaling, and compared the results with the gene targets for CzcR, implicated in zinc signaling. We discovered that the CopRs and CzcR have largely common targets, and crossregulate a core set of P. stutzeri copper and zinc responsive genes. We established that this crossregulation is enabled by a conserved binding motif in the upstream regulatory regions of the target genes. The crossregulation is physiologically relevant as these regulators synergistically and antagonistically target multicopper oxidases, metal efflux and sequestration systems. CopR1 and CopR2 upregulate two cop operons encoding copper tolerance genes, while all three regulators downregulate a putative copper chaperone, Psest_1595. CzcR also upregulated the oprD gene and the CzcIABC Zn²⁺ efflux system, while CopR1 and CopR2 downregulated these genes. Our study suggests that crossregulation of copper and zinc homeostasis can be advantageous, and in P. stutzeri this is enabled by shared binding motifs for multiple response regulators.     

Detection of pea wilt pathogen <Emphasis Type="Italic">Fusarium oxysporum</Emphasis> f. sp. <Emphasis Type="Italic">pisi</Emphasis> using DNA-based markers

Identification of the fungus Fusarium oxysporum f. sp. pisi (Fop), the causal organism of wilt disease of pea, is a time consuming and arduous task. Diagnosis of Fop by traditional means requires more than 2 months and involves two steps, identification of species using morphological characters and formae specialis ‘pisi’ using pathogenicity assays. The ambiguous morphological differences between F. solani and F. oxysporum further complicate the diagnosis of F. oxysporum. A polymerase chain reaction–restriction fragment length polymorphism (PCR–RFLP) based method was developed to detect Fop from India. A PCR–RFLP marker, HPACAPS1₃₈₀, generated after restriction of 28S rDNA region with enzyme MvaI, detected accurately the Fop among several other fungi with detection sensitivity of 5 fg of Fop genomic DNA. In a mixture of Fop and pea DNA, the sensitivity was 500 pg of Fop DNA in 50 ng of pea DNA. The assay was further refined to detect the Fop from infected tissues and infested soil. The current assay can detect Fop from culture, plant tissues and soil in a considerably shorter period of time compared to traditional methods.     

Development and evaluation of high‐density Axiom®CicerSNP Array for high‐resolution genetic mapping and breeding applications in chickpea

To accelerate genomics research and molecular breeding applications in chickpea, a high‐throughput SNP genotyping platform ‘Axiom^®CicerSNP Array’ has been designed, developed and validated. Screening of whole‐genome resequencing data from 429 chickpea lines identified 4.9 million SNPs, from which a subset of 70 463 high‐quality nonredundant SNPs was selected using different stringent filter criteria. This was further narrowed down to 61 174 SNPs based on p‐convert score ≥0.3, of which 50 590 SNPs could be tiled on array. Among these tiled SNPs, a total of 11 245 SNPs (22.23%) were from the coding regions of 3673 different genes. The developed Axiom^®CicerSNP Array was used for genotyping two recombinant inbred line populations, namely ICCRIL03 (ICC 4958 × ICC 1882) and ICCRIL04 (ICC 283 × ICC 8261). Genotyping data reflected high success and polymorphic rate, with 15 140 (29.93%; ICCRIL03) and 20 018 (39.57%; ICCRIL04) polymorphic SNPs. High‐density genetic maps comprising 13 679 SNPs spanning 1033.67 cM and 7769 SNPs spanning 1076.35 cM were developed for ICCRIL03 and ICCRIL04 populations, respectively. QTL analysis using multilocation, multiseason phenotyping data on these RILs identified 70 (ICCRIL03) and 120 (ICCRIL04) main‐effect QTLs on genetic map. Higher precision and potential of this array is expected to advance chickpea genetics and breeding applications.     

Chill-responsive dehydrins in blueberry: Are they associated with cold hardiness or dormancy transitions?

To survive winters, woody perennials of temperate zones must enter into endodormancy. Resumption of spring growth requires sufficient exposure to low temperature (chill units, CUs) in winter (chilling requirement), which also plays a role in the development of cold hardiness (cold acclimation). Physiological studies on dormancy breaking have focused on identifying markers, such as appearance or disappearance of proteins in response to varying degrees of chill unit accumulation. However, whether these changes are associated with dormancy transitions or cold acclimation is not clear. In the present study, greenhouse-grown blueberry (Vaccinium section Cyanococcus) plants were used to address this question. Three blueberry cultivars, Bluecrop, Tifblue, and Gulfcoast having chilling requirement of approximately 1 200, 900 and 600 CUs, respectively, were first exposed to 4°C for long enough to provide chill units equivalent to one-half of their respective chilling requirement. This treatment was expected to result in cold acclimation. A fraction of plants was then subjected to a 15/12°C (light/dark) regime for 2 weeks, a treatment expected to be “dormancy-neutral” but cause deacclimation. Before and after each treatment, cold hardiness and dormancy status of floral buds were determined; proteins were extracted from the buds collected on the same sampling date, and separated by one-dimensional SDS-PAGE. Dehydrin-like proteins were identified by immunoblotting, using anti-dehydrin antiserum. Results indicate that the chilling treatment resulted in cold acclimation as indicated by increased bud hardiness in all three cultivars. Data also indicate a distinct accumulation of three dehydrin-like proteins of 65, 60, and 14 kDa during cold acclimation. The cold hardiness and levels of dehydrin proteins decreased during the exposure to 15/12°C for 2 weeks. Results also confirmed that this treatment had no negative effect on chill unit accumulation. Densitometric scans of protein gels indicated a close association between the abundance of dehydrins and degree of cold hardiness in these cultivars. In addition, levels of the dehydrin proteins and cold hardiness remained about the same between 100% and >100% satisfaction of chilling requirement. These results suggest that changes in dehydrin expression are more closely associated with cold hardiness than with dormancy transitions.     

The single-copy gene encoding high-mobility-group protein HMG-I/Y from pea contains a single intron and is expressed in all organs

The coding and 3-downstream regions of the gene encoding the high mobility group protein HMG-I/Y from pea have been isolated, sequenced and characterised. A 795 bp pea genomic fragment containing the coding region of the pea HMG-I/Y gene with a single intron of 201 bp was isolated by PCR. The gene encodes a protein of 197 amino acid residues with four copies of the AT-hook DNA-binding motif encoded by exon 2. Southern blot analysis on genomic DNA revealed the presence of a single copy of the HMG-I/Y gene in the haploid genome. The pea HMG-I/Y gene is expressed in all organs of pea including roots, stems, leaves, flowers, tendrils and developing seeds, as determined by northern blot analysis.     

Triplex formation at physiological pH by 5-Me-dC-N4-(spermine) [X] oligodeoxynucleotides: non protonation of N3 in X of X*G:C triad and effect of base mismatch/ionic strength on triplex stabilities.

Oligodeoxynucleotide (ODN) directed triplex formation has therapeutic importance and depends on Hoogsteen hydrogen bonds between a duplex DNA and a third DNA strand. T*A:T triplets are formed at neutral pH and C+*G:C are favoured at acidic pH. It is demonstrated that spermine conjugation at N4 of 5-Me-dC in ODNs 1-5 (sp-ODNs) imparts zwitterionic character, thus reducing the net negative charge of ODNs 1-5. sp-ODNs form triplexes with complementary 24mer duplex 8:9 show foremost stability at neutral pH 7.3 and decrease in stability towards lower pH, unlike the normal ODNs where optimal stability is found at an acidic pH 5.5. At pH 7.3, control ODNs 6 and 7 carrying dC or 5-Me-dC, respectively, do not show any triple helix formation. The stability order of triplex containing 5-Me-dC-N4-(spermine) with normal and mismatched duplex was found to be X*G:C approximately X*A:T > X*C:G > X*T:A. The hysteresis curve of sp-ODN triplex 3*8:9 indicated a better association with complementary duplex 8:9 as compared to unmodified ODN 6 in triplex 6*8:9. pH-dependent UV difference spectra suggest that N3 protonation is not a requirement for triplex formation by sp-ODN and interstrand interaction of conjugated spermine more than compensates for loss in stability due to absence of a single Hoogsteen hydrogen bond. These results may have importance in designing oligonucleotides for antigene applications.     

Triplex formation at physiological pH: comparative studies on DNA triplexes containing 5-Me-dC tethered at N4 with spermine and tetraethyleneoxyamine.

Oligodeoxynucleotides with spermine conjugation at C4 of 5-Me-dC ( sp -ODN) exhibit triple helix formation with complementary Watson-Crick duplexes, and were optimally stable at physiological pH 7.3 and low salt concentration. This was attributed to a favored reassociation of the polycationic third strand with the anionic DNA duplex. To gain further insights into the factors that contribute to the enhancement of triplex stability and for engineering improved triplex systems, the spermine appendage at C4 of 5-Me-dC was replaced with 1,11-diamino-3,6,9-trioxaundecane to create teg -ODNs. From the triple helix forming abilities of these modified ODNs studied by hysteresis behaviour and the effect of salts on triplex stability, it is demonstrated here that teg- ODNs stabilise triplexes through hydrophobic desolvation while sp -ODNs stabilise triplexes by charge effects. The results imply that factors in addition to base stacking effects and interstrand hydrogen bonds are significantly involved in modulation of triplex stability by base modified oligonucleotides.