Estimation of metabolite networks with regard to a specific covariable: applications to plant and human data

Introduction

In systems biology, where a main goal is acquiring knowledge of biological systems, one of the challenges is inferring biochemical interactions from different molecular entities such as metabolites. In this area, the metabolome possesses a unique place for reflecting “true exposure” by being sensitive to variation coming from genetics, time, and environmental stimuli. While influenced by many different reactions, often the research interest needs to be focused on variation coming from a certain source, i.e. a certain covariable \(\mathbf {X}_m\).

Objective

Here, we use network analysis methods to recover a set of metabolite relationships, by finding metabolites sharing a similar relation to \(\mathbf {X}_m\). Metabolite values are based on information coming from individuals’ \(\mathbf {X}_m\) status which might interact with other covariables.

Methods

Alternative to using the original metabolite values, the total information is decomposed by utilizing a linear regression model and the part relevant to \(\mathbf {X}_m\) is further used. For two datasets, two different network estimation methods are considered. The first is weighted gene co-expression network analysis based on correlation coefficients. The second method is graphical LASSO based on partial correlations.

Results

We observed that when using the parts related to the specific covariable of interest, resulting estimated networks display higher interconnectedness. Additionally, several groups of biologically associated metabolites (very large density lipoproteins, lipoproteins, etc.) were identified in the human data example.

Conclusions

This work demonstrates how information on the study design can be incorporated to estimate metabolite networks. As a result, sets of interconnected metabolites can be clustered together with respect to their relation to a covariable of interest.

     

Generalized enhanced suffix array construction in external memory

Background

Suffix arrays, augmented by additional data structures, allow solving efficiently many string processing problems. The external memory construction of the generalized suffix array for a string collection is a fundamental task when the size of the input collection or the data structure exceeds the available internal memory.

Results

In this article we present and analyze \(\mathsf {eGSA}\) [introduced in CPM (External memory generalized suffix and \(\mathsf {LCP}\) arrays construction. In: Proceedings of CPM. pp 201–10, 2013)], the first external memory algorithm to construct generalized suffix arrays augmented with the longest common prefix array for a string collection. Our algorithm relies on a combination of buffers, induced sorting and a heap to avoid direct string comparisons. We performed experiments that covered different aspects of our algorithm, including running time, efficiency, external memory access, internal phases and the influence of different optimization strategies. On real datasets of size up to 24 GB and using 2 GB of internal memory, \(\mathsf {eGSA}\) showed a competitive performance when compared to \(\mathsf {eSAIS}\) and \(\mathsf {SAscan}\), which are efficient algorithms for a single string according to the related literature. We also show the effect of disk caching managed by the operating system on our algorithm.

Conclusions

The proposed algorithm was validated through performance tests using real datasets from different domains, in various combinations, and showed a competitive performance. Our algorithm can also construct the generalized Burrows-Wheeler transform of a string collection with no additional cost except by the output time.

     

Dynamic Optimization with Particle Swarms (DOPS): a meta-heuristic for parameter estimation in biochemical models

Background

Mathematical modeling is a powerful tool to analyze, and ultimately design biochemical networks. However, the estimation of the parameters that appear in biochemical models is a significant challenge. Parameter estimation typically involves expensive function evaluations and noisy data, making it difficult to quickly obtain optimal solutions. Further, biochemical models often have many local extrema which further complicates parameter estimation. Toward these challenges, we developed Dynamic Optimization with Particle Swarms (DOPS), a novel hybrid meta-heuristic that combined multi-swarm particle swarm optimization with dynamically dimensioned search (DDS). DOPS uses a multi-swarm particle swarm optimization technique to generate candidate solution vectors, the best of which is then greedily updated using dynamically dimensioned search.

Results

We tested DOPS using classic optimization test functions, biochemical benchmark problems and real-world biochemical models. We performed \(\mathcal {T}\) = 25 trials with \(\mathcal {N}\) = 4000 function evaluations per trial, and compared the performance of DOPS with other commonly used meta-heuristics such as differential evolution (DE), simulated annealing (SA) and dynamically dimensioned search (DDS). On average, DOPS outperformed other common meta-heuristics on the optimization test functions, benchmark problems and a real-world model of the human coagulation cascade.

Conclusions

DOPS is a promising meta-heuristic approach for the estimation of biochemical model parameters in relatively few function evaluations. DOPS source code is available for download under a MIT license at http://www.varnerlab.org.

     

An improved Four-Russians method and sparsified Four-Russians algorithm for RNA folding

Background

The basic RNA secondary structure prediction problem or single sequence folding problem (SSF) was solved 35 years ago by a now well-known \(O(n^3)\)-time dynamic programming method. Recently three methodologies—Valiant, Four-Russians, and Sparsification—have been applied to speedup RNA secondary structure prediction. The sparsification method exploits two properties of the input: the number of subsequence Z with the endpoints belonging to the optimal folding set and the maximum number base-pairs L. These sparsity properties satisfy \(0 \le L \le n / 2\) and \(n \le Z \le n^2 / 2\), and the method reduces the algorithmic running time to O(LZ). While the Four-Russians method utilizes tabling partial results.

Results

In this paper, we explore three different algorithmic speedups. We first expand the reformulate the single sequence folding Four-Russians \(\Theta \left(\frac{n^3}{\log ^2 n}\right)\)-time algorithm, to utilize an on-demand lookup table. Second, we create a framework that combines the fastest Sparsification and new fastest on-demand Four-Russians methods. This combined method has worst-case running time of \(O(\tilde{L}\tilde{Z})\), where \(\frac{{L}}{\log n} \le \tilde{L}\le min\left({L},\frac{n}{\log n}\right)\) and \(\frac{{Z}}{\log n}\le \tilde{Z} \le min\left({Z},\frac{n^2}{\log n}\right)\). Third we update the Four-Russians formulation to achieve an on-demand \(O( n^2/ \log ^2n )\)-time parallel algorithm. This then leads to an asymptotic speedup of \(O(\tilde{L}\tilde{Z_j})\) where \(\frac{{Z_j}}{\log n}\le \tilde{Z_j} \le min\left({Z_j},\frac{n}{\log n}\right)\) and \(Z_j\) the number of subsequence with the endpoint j belonging to the optimal folding set.

Conclusions

The on-demand formulation not only removes all extraneous computation and allows us to incorporate more realistic scoring schemes, but leads us to take advantage of the sparsity properties. Through asymptotic analysis and empirical testing on the base-pair maximization variant and a more biologically informative scoring scheme, we show that this Sparse Four-Russians framework is able to achieve a speedup on every problem instance, that is asymptotically never worse, and empirically better than achieved by the minimum of the two methods alone.

     

Self-preserving lognormal volume-size distributions of starch granules in developing sweetpotatoes and modulation of their scale parameters by a starch synthase II (SSII)

Main conclusion

Starch granule size distributions in plant tissues, when determined in high resolution and specifiedproperly as a frequency function, could provide useful information on the granule formation and growth.

Abstract

To better understand genetic control of physical properties of starch granules, we attempted a new approach to analyze developmental and genotypic effects on morphology and size distributions of starch granules in sweetpotato storage roots. Starch granules in sweetpotatoes exhibited low sphericity, many shapes that appeared to be independent of genotypes or developmental stages, and non-randomly distributed sizes. Granule size distributions of sweetpotato starches were determined in high resolution as differential volume-percentage distributions of volume-equivalent spherical diameters, rigorously curve-fitted to be lognormal, and specified using their geometric means \(\bar{x}^{*}\) and multiplicative standard deviations \(s^{*}\) in a \(\bar{x}^{*} \times /({\text{multiply/divide}})s^{*}\) form. The scale (\(\bar{x}^{*}\)) and shape (\(\bar{s}^{*}\)) of these distributions were independently variable, ranging from 14.02 to 19.36 μm and 1.403 to 1.567, respectively, among 22 cultivars/clones. The shape (\(s^{*}\)) of granule lognormal volume-size distributions of sweetpotato starch were found to be highly significantly and inversely correlated with their apparent amylose contents. More importantly, granule lognormal volume-size distributions of starches in developing sweetpotatoes displayed the same self-preserving kinetics, i.e., preserving the shape but shifting upward the scale, as those of particles undergoing agglomeration, which strongly indicated involvement of agglomeration in the formation and growth of starch granules. Furthermore, QTL analysis of a segregating null allele at one of three homoeologous starch synthase II loci in a reciprocal-cross population, which was identified through profiling starch granule-bound proteins in sweetpotatoes of diverse genotypes, showed that the locus is a QTL modulating the scale of granule volume-size distributions of starch in sweetpotatoes.

     

Tapping the genetic diversity of landraces in allogamous crops with doubled haploid lines: a case study from European flint maize

Key message

Using landraces for broadening the genetic base of elite maize germplasm is hampered by heterogeneity and high genetic load. Production of DH line libraries can help to overcome these problems.

Abstract

Landraces of maize (Zea mays L.) represent a huge reservoir of genetic diversity largely untapped by breeders. Genetic heterogeneity and a high genetic load hamper their use in hybrid breeding. Production of doubled haploid line libraries (DHL) by the in vivo haploid induction method promises to overcome these problems. To test this hypothesis, we compared the line per se performance of 389 doubled haploid (DH) lines across six DHL produced from European flint landraces with that of four flint founder lines (FFL) and 53 elite flint lines (EFL) for 16 agronomic traits evaluated in four locations. The genotypic variance (\(\sigma _{G}^{2}\)) within DHL was generally much larger than that among DHL and exceeded also \(\sigma _{G}^{2}\) of the EFL. For most traits, the means and \(\sigma _{G}^{2}\) differed considerably among the DHL, resulting in different expected selection gains. Mean grain yield of the EFL was 25 and 62% higher than for the FFL and DHL, respectively, indicating considerable breeding progress in the EFL and a remnant genetic load in the DHL. Usefulness of the best 20% lines was for individual DHL comparable to the EFL and grain yield (GY) in the top lines from both groups was similar. Our results corroborate the tremendous potential of landraces for broadening the narrow genetic base of elite germplasm. To make best use of these “gold reserves”, we propose a multi-stage selection approach with optimal allocation of resources to (1) choose the most promising landraces for DHL production and (2) identify the top DH lines for further breeding.

     

Cellular mechanosensitivity to substrate stiffness decreases with increasing dissimilarity to cell stiffness

Computational modelling has received increasing attention to investigate multi-scale coupled problems in micro-heterogeneous biological structures such as cells. In the current study, we investigated for a single cell the effects of (1) different cell-substrate attachment (2) and different substrate modulus \(\textit{E}_\mathrm{s}\) on intracellular deformations. A fibroblast was geometrically reconstructed from confocal micrographs. Finite element models of the cell on a planar substrate were developed. Intracellular deformations due to substrate stretch of \(\lambda =1.1\), were assessed for: (1) cell-substrate attachment implemented as full basal contact (FC) and 124 focal adhesions (FA), respectively, and \(\textit{E}_\mathrm{s}\,=\,\)140 KPa and (2) \(\textit{E}_\mathrm{s}\,=\,10\), 140, 1000, and 10,000 KPa, respectively, and FA attachment. The largest strains in cytosol, nucleus and cell membrane were higher for FC (1.35\(\text {e}^{-2}\), 0.235\(\text {e}^{-2}\) and 0.6\(\text {e}^{-2}\)) than for FA attachment (0.0952\(\text {e}^{-2}\), 0.0472\(\text {e}^{-2}\) and 0.05\(\text {e}^{-2}\)). For increasing \(\textit{E}_\mathrm{s}\), the largest maximum principal strain was 4.4\(\text {e}^{-4}\), 5\(\text {e}^{-4}\), 5.3\(\text {e}^{-4}\) and 5.3\(\text {e}^{-4}\) in the membrane, 9.5\(\text {e}^{-4}\), 1.1\(\text {e}^{-4}\), 1.2\(\text {e}^{-3}\) and 1.2\(\text {e}^{-3}\) in the cytosol, and 4.5\(\text {e}^{-4}\), 5.3\(\text {e}^{-4}\), 5.7\(\text {e}^{-4}\) and 5.7\(\text {e}^{-4}\) in the nucleus. The results show (1) the importance of representing FA in cell models and (2) higher cellular mechanical sensitivity for substrate stiffness changes in the range of cell stiffness. The latter indicates that matching substrate stiffness to cell stiffness, and moderate variation of the former is very effective for controlled variation of cell deformation. The developed methodology is useful for parametric studies on cellular mechanics to obtain quantitative data of subcellular strains and stresses that cannot easily be measured experimentally.     

Superbubbles revisited

Background

Superbubbles are distinctive subgraphs in direct graphs that play an important role in assembly algorithms for high-throughput sequencing (HTS) data. Their practical importance derives from the fact they are connected to their host graph by a single entrance and a single exit vertex, thus allowing them to be handled independently. Efficient algorithms for the enumeration of superbubbles are therefore of important for the processing of HTS data. Superbubbles can be identified within the strongly connected components of the input digraph after transforming them into directed acyclic graphs. The algorithm by Sung et al. (IEEE ACM Trans Comput Biol Bioinform 12:770–777, 2015) achieves this task in \(\mathcal {O}(m~log(m))\)-time. The extraction of superbubbles from the transformed components was later improved to by Brankovic et al. (Theor Comput Sci 609:374–383, 2016) resulting in an overall \(\mathcal {O}(m+n)\)-time algorithm.

Results

A re-analysis of the mathematical structure of superbubbles showed that the construction of auxiliary DAGs from the strongly connected components in the work of Sung et al. missed some details that can lead to the reporting of false positive superbubbles. We propose an alternative, even simpler auxiliary graph that solved the problem and retains the linear running time for general digraph. Furthermore, we describe a simpler, space-efficient \(\mathcal {O}(m+n)\)-time algorithm for detecting superbubbles in DAGs that uses only simple data structures.

Implementation

We present a reference implementation of the algorithm that accepts many commonly used formats for the input graph and provides convenient access to the improved algorithm. https://github.com/Fabianexe/Superbubble.

     

Sorting signed circular permutations by super short operations

Background

One way to estimate the evolutionary distance between two given genomes is to determine the minimum number of large-scale mutations, or genome rearrangements, that are necessary to transform one into the other. In this context, genomes can be represented as ordered sequences of genes, each gene being represented by a signed integer. If no gene is repeated, genomes are thus modeled as signed permutations of the form \(\pi =(\pi _1 \pi _2 \ldots \pi _n)\), and in that case we can consider without loss of generality that one of them is the identity permutation \(\iota _n =(1 2 \ldots n)\), and that we just need to sort the other (i.e., transform it into \(\iota _n\)). The most studied genome rearrangement events are reversals, where a segment of the genome is reversed and reincorporated at the same location; and transpositions, where two consecutive segments are exchanged. Many variants, e.g., combining different types of (possibly constrained) rearrangements, have been proposed in the literature. One of them considers that the number of genes involved, in a reversal or a transposition, is never greater than two, which is known as the problem of sorting by super short operations (or SSOs).

Results and conclusions

All problems considering SSOs in permutations have been shown to be in \(\mathsf {P}\), except for one, namely sorting signed circular permutations by super short reversals and super short transpositions. Here we fill this gap by introducing a new graph structure called cyclic permutation graph and providing a series of intermediate results, which allows us to design a polynomial algorithm for sorting signed circular permutations by super short reversals and super short transpositions.

     

Glutathione <Emphasis Type="BoldItalic">S</Emphasis>-transferase P1 gene polymorphisms and susceptibility to coronary artery disease in a subgroup of north Indian population

The present study aimed to investigate the association of \(\hbox {g}.313\hbox {A}{>}\hbox {G}\) and \(\hbox {g}.341\hbox {C}{>}\hbox {T}\) polymorphisms of GSTP1 with coronary artery disease (CAD) in a subgroup of north Indian population. In the present case–control study, CAD patients (\(n = 200\)) and age-matched, sex-matched and ethnicity-matched healthy controls (\(n = 200\)) were genotyped for polymorphisms in GSTP1 using polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method. Genotype distribution of \(\hbox {g}.313\hbox {A}{>}\hbox {G}\) and \(\hbox {g}.341\hbox {C}{>}\hbox {T}\) polymorphisms of GSTP1 gene was significantly different between cases and controls (\(P = 0.005\) and 0.024, respectively). Binary logistic regression analysis showed significant association of A/G (odds ratio (OR): 1.6, 95% CI: 1.08–2.49, \(P = 0.020\)) and G/G (OR: 3.1, 95% CI: 1.41–6.71, P \(=\) 0.005) genotypes of GSTP1 \(\hbox {g}.313\hbox {A}{\!>\!}\hbox {G}\), and C/T (OR: 5.8, 95% CI: 1.26–26.34, \(P = 0.024\)) genotype of GSTP1 \(\hbox {g}.341\hbox {C}{>}\hbox {T}\) with CAD. The A/G and G/G genotypes of \(\hbox {g}.313\hbox {A}{>}\hbox {G}\) and C/T genotype of \(\hbox {g}.341\hbox {C}{>}\hbox {T}\) conferred 6.5-fold increased risk for CAD (OR: 6.5, 95% CI: 1.37–31.27, \(P = 0.018\)). Moreover, the recessive model of GSTP1 \(\hbox {g}.313\hbox {A}{>}\hbox {G}\) is the best fit inheritance model to predict the susceptible gene effect (OR: 2.3, 95% CI: 1.11–4.92, \(P = 0.020\)). In conclusion, statistically significant associations of GSTP1 \(\hbox {g}.313\hbox {A}{>}\hbox {G}\) (A/G, G/G) and \(\hbox {g}.341\hbox {C}{>}\hbox {T}\) (C/T) genotypes with CAD were observed.     

A Mathematical Model for the Macrophage Response to Respiratory Viral Infection in Normal and Asthmatic Conditions

Respiratory viral infections are common in the general population and one of the most important causes of asthma aggravation and exacerbation. Despite many studies, it is not well understood how viral infections cause more severe symptoms and exacerbations in asthmatics. We develop a mathematical model of two types of macrophages that play complementary roles in fighting viral infection: classically \((\hbox {CA}\)-\(\hbox {M}\Phi )\) and alternatively activated macrophages \((\hbox {AA}\)-\(\hbox {M}\Phi )\). \(\hbox {CA}\)-\(\hbox {M}\Phi \) destroy infected cells and tissues to remove viruses, while \(\hbox {AA}\)-\(\hbox {M}\Phi \) repair damaged tissues. We show that a higher viral load or longer duration of infection provokes a stronger immune response from the macrophage system. By adjusting the parameters, we model the differences in response to respiratory viral infection in normal and asthmatic subjects and show how this skews the system toward a response that generates more severe symptoms in asthmatic patients.     

FSH: fast spaced seed hashing exploiting adjacent hashes

Background

Patterns with wildcards in specified positions, namely spaced seeds, are increasingly used instead of k-mers in many bioinformatics applications that require indexing, querying and rapid similarity search, as they can provide better sensitivity. Many of these applications require to compute the hashing of each position in the input sequences with respect to the given spaced seed, or to multiple spaced seeds. While the hashing of k-mers can be rapidly computed by exploiting the large overlap between consecutive k-mers, spaced seeds hashing is usually computed from scratch for each position in the input sequence, thus resulting in slower processing.

Results

The method proposed in this paper, fast spaced-seed hashing (FSH), exploits the similarity of the hash values of spaced seeds computed at adjacent positions in the input sequence. In our experiments we compute the hash for each positions of metagenomics reads from several datasets, with respect to different spaced seeds. We also propose a generalized version of the algorithm for the simultaneous computation of multiple spaced seeds hashing. In the experiments, our algorithm can compute the hashing values of spaced seeds with a speedup, with respect to the traditional approach, between 1.6\(\times\) to 5.3\(\times\), depending on the structure of the spaced seed.

Conclusions

Spaced seed hashing is a routine task for several bioinformatics application. FSH allows to perform this task efficiently and raise the question of whether other hashing can be exploited to further improve the speed up. This has the potential of major impact in the field, making spaced seed applications not only accurate, but also faster and more efficient.

Availability

The software FSH is freely available for academic use at: https://bitbucket.org/samu661/fsh/overview.

     

Bayesian estimation of the number of protonation sites for urinary metabolites from NMR spectroscopic data

Introduction

To aid the development of better algorithms for \(^1\)H NMR data analysis, such as alignment or peak-fitting, it is important to characterise and model chemical shift changes caused by variation in pH. The number of protonation sites, a key parameter in the theoretical relationship between pH and chemical shift, is traditionally estimated from the molecular structure, which is often unknown in untargeted metabolomics applications.

Objective

We aim to use observed NMR chemical shift titration data to estimate the number of protonation sites for a range of urinary metabolites.

Methods

A pool of urine from healthy subjects was titrated in the range pH 2–12, standard \(^1\)H NMR spectra were acquired and positions of 51 peaks (corresponding to 32 identified metabolites) were recorded. A theoretical model of chemical shift was fit to the data using a Bayesian statistical framework, using model selection procedures in a Markov Chain Monte Carlo algorithm to estimate the number of protonation sites for each molecule.

Results

The estimated number of protonation sites was found to be correct for 41 out of 51 peaks. In some cases, the number of sites was incorrectly estimated, due to very close pKa values or a limited amount of data in the required pH range.

Conclusions

Given appropriate data, it is possible to estimate the number of protonation sites for many metabolites typically observed in \(^1\)H NMR metabolomics without knowledge of the molecular structure. This approach may be a valuable resource for the development of future automated metabolite alignment, annotation and peak fitting algorithms.

     

Identifying sequence features that drive ribosomal association for lncRNA

Background

With the increasing number of annotated long noncoding RNAs (lncRNAs) from the genome, researchers are continually updating their understanding of lncRNAs. Recently, thousands of lncRNAs have been reported to be associated with ribosomes in mammals. However, their biological functions or mechanisms are still unclear.

Results

In this study, we tried to investigate the sequence features involved in the ribosomal association of lncRNA. We have extracted ninety-nine sequence features corresponding to different biological mechanisms (i.e., RNA splicing, putative ORF, k-mer frequency, RNA modification, RNA secondary structure, and repeat element). An \(\mathcal {L}1\)-regularized logistic regression model was applied to screen these features. Finally, we obtained fifteen and nine important features for the ribosomal association of human and mouse lncRNAs, respectively.

Conclusion

To our knowledge, this is the first study to characterize ribosome-associated lncRNAs and ribosome-free lncRNAs from the perspective of sequence features. These sequence features that were identified in this study may shed light on the biological mechanism of the ribosomal association and provide important clues for functional analysis of lncRNAs.

     

Trend analysis of variations in carbon stock using stock big data

Changes in land use affect the terrestrial carbon stock through changes in the land cover. Research on land use and analysis of variations in carbon stock have practical applications in the optimization of land use and the mitigation of climate change effects. This study was conducted in Baixiang and Julu counties in the Taihang Piedmont by employing the trend analysis method to characterize the variation in county land use and carbon stock. The findings show that in both counties, agricultural and unused land areas decreased while built-up land area increased, and the reduction in cropland was the main reason behind the agricultural land reduction. An inflection point appeared on the cropland curves of Julu, because the cropland area decreased by 1576.97 hm\(^{2}\) from 2004 to 2006. Cropland area in Baixiang decreased from 1996 to 1998 by a total of 129.89 hm\(^{2}\) and then remained relatively stable after 1998. The total carbon storage and variation in land use in the two counties displayed similar trends. Total carbon reserves in Julu increased by 2.76 \(\times \) 10\(^{4}\) tC (carbon equivalent), while those in Baixiang decreased by 0.63 \(\times \) 10\(^{4}\) tC. Carbon stock of built-up land in Julu and Baixiang increased by 2.44 \(\times \) 10\(^{4}\) and 1.22 \(\times \) 10\(^{4}\) tC, respectively.     

Isolation and characterization of microsatellite markers in the spiny lobster, <Emphasis Type="Italic">Panulirus echinatus</Emphasis> Smith, 1869 (Decapoda: Palinuridae) by Illumina MiSeq sequencing

Okra’s (Abelmoschus esculentus (L.) Moench) commercial cultivation is threatened in the tropics due to high incidence of yellow vein mosaic virus (YVMV) disease. Okra geneticists across the world tried to understand the inheritance pattern of YVMV disease tolerance without much success. Therefore, the inheritance pattern of YVMV disease in okra was revisited by employing six generations (\(\hbox {P}_{1}\), \(\hbox {P}_{2}\), \(\hbox {F}_{1}\), \(\hbox {F}_{2}\), \(\hbox {BC}_{1}\) and \(\hbox {BC}_{2}\)) of four selected crosses (one tolerant \(\times \) tolerant, two tolerant \(\times \) susceptible and one susceptible \(\times \) susceptible) using two tolerant (BCO-1 and Lal Bhendi) and two susceptible (Japanese Jhar Bhendi and PAN 2127) genotypes. Qualitative genetic analysis was done on the basis of segregation pattern of tolerant and susceptible plants in \(\hbox {F}_{2}\) and backcross generations of all the four crosses. It revealed that a single dominant gene along with some minor factors governed the disease tolerant trait in both the tolerant parents used. However, it was observed that genes governing disease tolerance identified in both the tolerant variety used was different. It could be concluded that the gene governing YVMV disease tolerance in okra was genotype specific. Further, duplicate gene action as evident from an approximate ratio of 15 : 1 (tolerant : susceptible) in the \(\hbox {F}_{2}\) population in the cross of two tolerant varieties gave a scope of increasing the tolerance level of the hybrid plants when both the tolerant genes are brought together. However, generation mean analysis revealed involvement of both additive and nonadditive effects in the inheritance of disease tolerance. Thus, the present study confirms that a complicated genetic inheritance pattern is involved in the disease tolerance against YVMV trait. The major tolerance genes could be transferred to other okra varieties, but the tolerance breaking virus strains might not allow them to achieve tolerance in stable condition. Therefore, accumulation of additional genes may be needed for a sustainable tolerance phenotype in okra.     

Backbone resonance assignments for the SET domain of the human methyltransferase NSD2

Aberrant NSD2 methyltransferase activity is implicated as the oncogenic driver in multiple myeloma, suggesting opportunities for novel therapeutic intervention. The methyltransferase activity of NSD2 resides in its catalytic SET domain, which is conserved among most lysine methyltransferases. Here we report the backbone \(\hbox {H}^{\mathrm{N}}\), N, C\(^{\prime }\), \(\hbox {C}^\alpha\) and side-chain \(\hbox {C}^\beta\) assignments of a 25 kDa NSD2 SET domain construct, spanning residues 991–1203. A chemical shift analysis of C\(^{\prime }\), \(\hbox {C}^\alpha\) and \(\hbox {C}^\beta\) resonances predicts a secondary structural pattern that is in agreement with homology models.     

Loss of ncm<Superscript>5</Superscript> and mcm<Superscript>5</Superscript> wobble uridine side chains results in an altered metabolic profile

Introduction

The Elongator complex, comprising six subunits (Elp1p-Elp6p), is required for formation of 5-carbamoylmethyl (ncm⁵) and 5-methoxycarbonylmethyl (mcm⁵) side chains on wobble uridines in 11 out of 42 tRNA species in Saccharomyces cerevisiae. Loss of these side chains reduces the efficiency of tRNA decoding during translation, resulting in pleiotropic phenotypes. Overexpression of hypomodified \( {\text {tRNA}_{{\rm s^{2} {\rm UUU}}}^{{\rm Lys}} , {\rm tRNA}_{{\rm s^{2} {\rm UUG}}}^{{\rm Gln }} \;{\rm and}\;{\rm tRNA}_{{\rm s^{2} {\rm UUC}}}^{{\rm Glu}}} \), which in wild-type strains are modified with mcm⁵s²U, partially suppress phenotypes of an elp3Δ strain.

Objectives

To identify metabolic alterations in an elp3Δ strain and elucidate whether these metabolic alterations are suppressed by overexpression of hypomodified \( {\text {tRNA}_{{\rm s^{2} {\rm UUU}}}^{{\rm Lys}} , {\rm tRNA}_{{\rm s^{2} {\rm UUG}}}^{{\rm Gln }} \;{\rm and}\;{\rm tRNA}_{{\rm s^{2} {\rm UUC}}}^{{\rm Glu}}} \).

Method

Metabolic profiles were obtained using untargeted GC-TOF-MS of a temperature-sensitive elp3Δ strain carrying either an empty low-copy vector, an empty high-copy vector, a low-copy vector harboring the wild-type ELP3 gene, or a high-copy vector overexpressing \( {\text {tRNA}_{{\rm s^{2} {\rm UUU}}}^{{\rm Lys}} , {\rm tRNA}_{{\rm s^{2} {\rm UUG}}}^{{\rm Gln }} \;{\rm and}\;{\rm tRNA}_{{\rm s^{2} {\rm UUC}}}^{{\rm Glu}}} \). The temperature sensitive elp3Δ strain derivatives were cultivated at permissive (30 °C) or semi-permissive (34 °C) growth conditions.

Results

Culturing an elp3Δ strain at 30 or 34 °C resulted in altered metabolism of 36 and 46 %, respectively, of all metabolites detected when compared to an elp3Δ strain carrying the wild-type ELP3 gene. Overexpression of hypomodified \( {\text {tRNA}_{{\rm s^{2} {\rm UUU}}}^{{\rm Lys}} , {\rm tRNA}_{{\rm s^{2} {\rm UUG}}}^{{\rm Gln }} \;{\rm and}\;{\rm tRNA}_{{\rm s^{2} {\rm UUC}}}^{{\rm Glu}}} \) suppressed a subset of the metabolic alterations observed in the elp3Δ strain.

Conclusion

Our results suggest that the presence of ncm⁵- and mcm⁵-side chains on wobble uridines in tRNA are important for metabolic homeostasis.

     

StreAM-$$T_g$$: algorithms for analyzing coarse grained RNA dynamics based on Markov models of connectivity-graphs

Background

In this work, we present a new coarse grained representation of RNA dynamics. It is based on adjacency matrices and their interactions patterns obtained from molecular dynamics simulations. RNA molecules are well-suited for this representation due to their composition which is mainly modular and assessable by the secondary structure alone. These interactions can be represented as adjacency matrices of k nucleotides. Based on those, we define transitions between states as changes in the adjacency matrices which form Markovian dynamics. The intense computational demand for deriving the transition probability matrices prompted us to develop StreAM-\(T_g\), a stream-based algorithm for generating such Markov models of k-vertex adjacency matrices representing the RNA.

Results

We benchmark StreAM-\(T_g\) (a) for random and RNA unit sphere dynamic graphs (b) for the robustness of our method against different parameters. Moreover, we address a riboswitch design problem by applying StreAM-\(T_g\) on six long term molecular dynamics simulation of a synthetic tetracycline dependent riboswitch (500 ns) in combination with five different antibiotics.

Conclusions

The proposed algorithm performs well on large simulated as well as real world dynamic graphs. Additionally, StreAM-\(T_g\) provides insights into nucleotide based RNA dynamics in comparison to conventional metrics like the root-mean square fluctuation. In the light of experimental data our results show important design opportunities for the riboswitch.