Normalization of miRNA qPCR high-throughput data: a comparison of methods

Low-density quantitative real-time PCR (qPCR) arrays are often used to profile expression patterns of microRNAs in various biological milieus. To achieve accurate analysis of expression of miRNAs, non-biological sources of variation in data should be removed through precise normalization of data. We have systematically compared the performance of 19 normalization methods on different subsets of a real miRNA qPCR dataset that covers 40 human tissues. After robustly modeling the mean squared error (MSE) in normalized data, we demonstrate lower variability between replicates is achieved using various methods not applied to high-throughput miRNA qPCR data yet. Normalization methods that use splines or wavelets smoothing to estimate and remove Cq dependent non-linearity between pairs of samples best reduced the MSE of differences in Cq values of replicate samples. These methods also retained between-group variability in different subsets of the dataset.     

Mixed Mode of Ligand-DNA Binding Results in S-Shaped Binding Curves

Abstract

S-shaped binding curves often characterize interactions of ligands with nucleic acid molecules as analyzed by different physicochemical and biophysical techniques. S-shaped experimental binding curves are usually interpreted as indicative of the positive cooperative interactions between the bound ligand molecules. This paper demonstrates that S-shaped binding curves may occur as a result of the “mixed mode” of DNA binding by the same ligand molecule. Mixed mode of the ligand-DNA binding can occur, for example, due to 1) isomerization or dimerization of the ligands in solution or on the DNA lattice, 2) their ability to intercalate the DNA and to bind it within the minor groove in different orientations. DNA- ligand complexes are characterized by the length of the ligand binding site on the DNA lattice (so-called “multiple-contact” model). We show here that if two or more complexes with different lengths of the ligand binding sites could be produced by the same ligand, the dependence of the concentration of the complex with the shorter length of binding site on the total concentration of ligand should be S-shaped. Our theoretical model is confirmed by comparison of the calculated and experimental CD binding curves for bis-netropsin binding to poly(dA-dT) poly(dA-dT). Bis-netropsin forms two types of DNA complexes due to its ability to interact with the DNA as monomers and trimers. Experimental S-shaped bis-netropsin-DNA binding curve is shown to be in good correlation with those calculated on the basis of our theoretical model. The present work provides new insight into the analysis of ligand-DNA binding curves.     

Comprehensive phylogenetic analyses of the Ruppia maritima complex focusing on taxa from the Mediterranean

Recent molecular phylogenetic studies reported high diversity of Ruppia species in the Mediterranean. Multiple taxa, including apparent endemics, are known from that region, however, they have thus far not been exposed to phylogenetic analyses aimed at studying their relationships to taxa from other parts of the world. Here we present a comprehensive phylogenetic analyses of the R. maritima complex using data sets composed of DNA sequences of the plastid genome, the multi-copy nuclear ITS region, and the low-copy nuclear phyB gene with a primary focus on the Mediterranean representatives of the complex. As a result, a new lineage, “Drepanensis”, was identified as the seventh entity of the complex. This lineage is endemic to the Mediterranean. The accessions included in the former “Tetraploid” entity were reclassified into two entities: an Asia–Australia–Europe disjunct “Tetraploid_α” with a paternal “Diploid” origin, and a European “Tetraploid_γ” originating from a maternal “Drepanensis” lineage. Another entity, “Tetraploid_β”, is likely to have been originated as a result of chloroplast capture through backcrossing hybridization between paternal “Tetraploid_α” and maternal “Tetraploid_γ”. Additional discovery of multiple tetraploidizations as well as hybridization and chloroplast capture at the tetraploid level indicated that hybridization has been a significant factor in the diversification of Ruppia.     

Needle δ13C and mobile carbohydrates in Pinus koraiensis in relation to decreased temperature and increased moisture along an elevational gradient in NE China

A tree’s crown interacts with atmospheric variables such as CO₂, temperature, and humidity. Physioecology of leaves/needles (e.g. δ¹³C, mobile carbohydrates, and nitrogen) is, therefore, strongly affected by microclimate in and surrounding a tree crown. To understand the physiological responses of leaves to changes in air temperature and moisture, we measured δ¹³C, soluble sugars, starch, and total nitrogen (N) concentrations in current year and 1-yr-old needles of Pinus koraiensis trees, and compared the growing season air temperature and relative humidity within and outside P. koraiensis crowns along an elevational gradient from 760 to 1,420 m a.s.l. on Changbai Mountain, NE China. Our results indicated that needle N and mobile carbohydrates concentrations, as well as needle δ¹³C values changed continuously with increasing elevation, corresponding to a continuous decrease in air temperature and an increase in relative humidity. Needle carbon and nitrogen status is highly significantly negatively correlated with temperature, but positively correlated with relative humidity. These results indicate that increases in air temperature in combination with decreases in relative humidity may result in lower levels of N and mobile carbohydrates in P. koraiensis trees, suggesting that future climate changes such as global warming and changes in precipitation patterns will directly influence the N and carbon physiology at P. koraiensis individual level, and indirectly affect the competitive ability, species composition, productivity and functioning at the stand and ecosystem level in NE China. Due to the relatively limited range of the transect (760–1,420 m) studied, further research is needed to explain whether the present results are applicable to scales across large elevational gradients.     

Epistasis in iron metabolism: complex interactions between Cp, Mon1a, and Slc40a1 loci and tissue iron in mice

Disorders of iron metabolism are among the most common acquired and constitutive diseases. Hemochromatosis has a solid genetic basis and in Northern European populations it is usually associated with homozygosity for the C282Y mutation in the HFE protein. However, the penetrance of this mutation is incomplete and the clinical presentation is highly variable. The rare and common variants identified so far as genetic modifiers of HFE-related hemochromatosis are unable to account for the phenotypic heterogeneity of this disorder. There are wide variations in the basal iron status of common inbred mouse strains, and this diversity may reflect the genetic background of the phenotypic diversity under pathological conditions. We therefore examined the genetic basis of iron homeostasis using quantitative trait loci mapping applied to the HcB-15 recombinant congenic strains for tissue and serum iron indices. Two highly significant QTL containing either the N374S Mon1a mutation or the Ferroportin locus were found to be major determinants in spleen and liver iron loading. Interestingly, when considering possible epistatic interactions, the effects of Mon1a on macrophage iron export are conditioned by the genotype at the Slc40a1 locus. Only mice that are C57BL/10ScSnA homozygous at both loci display a lower spleen iron burden. Furthermore, the liver-iron lowering effect of the N374S Mon1a mutation is observed only in mice that display a nonsense mutation in the Ceruloplasmin (Cp) gene. This study highlights the existence of genetic interactions between Cp, Mon1a, and the Slc40a1 locus in iron metabolism, suggesting that epistasis may be a crucial determinant of the variable biological and clinical presentations in iron disorders.     

Nanoscale Distribution of Ryanodine Receptors and Caveolin-3 in Mouse Ventricular Myocytes: Dilation of T-Tubules near Junctions

We conducted super-resolution light microscopy (LM) imaging of the distribution of ryanodine receptors (RyRs) and caveolin-3 (CAV3) in mouse ventricular myocytes. Quantitative analysis of data at the surface sarcolemma showed that 4.8% of RyR labeling colocalized with CAV3 whereas 3.5% of CAV3 was in areas with RyR labeling. These values increased to 9.2 and 9.0%, respectively, in the interior of myocytes where CAV3 was widely expressed in the t-system but reduced in regions associated with junctional couplings. Electron microscopic (EM) tomography independently showed only few couplings with caveolae and little evidence for caveolar shapes on the t-system. Unexpectedly, both super-resolution LM and three-dimensional EM data (including serial block-face scanning EM) revealed significant increases in local t-system diameters in many regions associated with junctions. We suggest that this regional specialization helps reduce ionic accumulation and depletion in t-system lumen during excitation-contraction coupling to ensure effective local Ca²⁺ release. Our data demonstrate that super-resolution LM and volume EM techniques complementarily enhance information on subcellular structure at the nanoscale.The contraction of cardiac ventricular myocytes depends on the rapid cell-wide transient increase in intracellular [Ca²⁺] upon depolarization of the cell-membrane potential. The cardiac ryanodine receptor (RyR) (1), which is the intracellular Ca²⁺ release channel in the sarcoplasmic reticulum (SR), plays a central role in shaping Ca²⁺ transients. RyRs form clusters of various sizes (2,3) with the majority located within junctions between the SR and the surface membrane and its cytoplasmic extension, the transverse tubular (t-) system. It has been suggested that some RyR clusters are associated with caveolae, a specialized signaling microdomain of the surface membrane. Previous studies were complicated by the limited resolution of optical imaging methods of ∼250 nm, much larger than the nanometer scale of RyRs and caveolae. Accordingly, these studies report varying colocalization between RyRs and caveolin-3 (CAV3), a caveolar marker also expressed in the t-system (4,5).In this work, we investigated the relative distribution of CAV3 and RyRs in mouse ventricular myocytes both in the cytosol and near the cell surface with super-resolution fluorescence microscopy that achieves a resolution approaching 30 nm. Our data revealed unexpected local t-system swellings near junctional couplings, which was supported by two different three-dimensional electron microscopy (EM) modalities with <10-nm resolution: EM tomography and serial block-face scanning EM (SBFSEM).Super-resolution images of CAV3 and RyR labeling at the surface sarcolemma of mouse myocytes showed little overlap, suggesting that few RyRs were in couplings with caveolae (Fig. 1 A, for detailed methods, see the Supporting Material). Only ∼4.8% of RyR labeling was associated with CAV3 positive areas and ∼3.5% of CAV3 associated with RyR positive areas (n = 6 cells from three animals, Fig. 1 B, see also Table S1 in the Supporting Material), broadly consistent with previous data in rats (6). To support this finding, EM tomography was applied to mouse ventricular tissue that included a part of the surface sarcolemma, to our knowledge for the first time. Segmentation of peripheral couplings (containing RyR foot structures) and surface caveolae (∼60 nm in diameter and often interconnected) confirmed that the great majority of peripheral couplings were in regions devoid of caveolae (Fig. 1 C). A few junctional couplings containing feet were between caveolae and subsarcolemmal SR (Fig. 1 D, see also Fig. S1 and Movie S1 in the Supporting Material). We conducted a similar analysis in the cytosol where CAV3 expression occurs in the t-system (5) and RyRs are abundant in dyadic junctions between the t-system and SR terminal cisterns.Open in a separate windowFigure 1Colocalization of CAV3 and RyRs at the surface sarcolemma. (A) Super-resolution micrograph of the distribution of CAV3 (green) and RyRs (red) at the surface of a mouse cardiac myocyte. (B) Analysis of the association of CAV3 with RyRs. The fraction of RyR labeling within CAV3 positive areas was ∼4.8% (front data) whereas ∼3.5% of CAV3 was found in RyR-positive membrane areas. (C) Segmented EM tomogram containing a patch of surface sarcolemma (light blue) and associated caveolae (green) as well as peripheral couplings (red). (D) Detailed view of a region with abundant caveolae. (Arrows) Couplings with caveolae.As shown in Fig. 2 A, the spatial distribution of CAV3 and RyR clusters in super-resolution micrographs taken several microns below the surface sarcolemma is consistent with this view. The association of the two labels is slightly increased (as compared to the surface), according to distance analysis with 9% of CAV3 and 9.2% of RyR labeling associating with each other (Fig. 2 B, n = 6 cells from three animals). The similarity of manually traced t-system in EM tomograms (Fig. 2 C) and super-resolved CAV3 labeling suggested that CAV3 is widely distributed in the t-system except for regions where dyadic membrane junctions occur as CAV3 labeling was much weaker in regions with strong RyR labeling. It was notable that the t-system diameter appeared to increase at regions of strong RyR labeling (Fig. 2 D), broadly consistent with the behavior seen in tomograms (Fig. 2 C). This was confirmed by a quantitative analysis of t-tubule diameters in dyadic versus extradyadic regions on the basis of CAV3 and RyR labeling, with full-width at quarter-maximum mean diameters increasing from ∼150 nm distal to dyads, to ∼190 nm (using CAV3 signal only) or ∼280 nm (using CAV3 and RyR signal) near dyads (Fig. 2, G and H, see also Methods in the Supporting Material). The combined RyR and CAV3 signals seemed to be a better representation of the entire t-system lumen near junctions (see Fig. S2).Open in a separate windowFigure 2Distribution of CAV3 and RyRs in the cell interior. (A) Super-resolution micrograph of CAV3 (green) and RyR (red) distribution at t-system. (Arrow) Direction of longitudinal cell axis. (B) Distance analysis of the CAV3 and RyR association (N = 6 cells per group). (C) Segmented EM tomogram of a similar region with three-dimensional mesh models of t-system membrane (green) and dyadic couplings (red). (D) This image illustrates the tracing (white path) of t-tubules. The label distribution was extracted and linearized along the path (E) to calculate a mask that shows the full width at quarter-maximum diameter along tubules, CAV3 (green) and RyR (red) (F). (G) Histograms of local diameters extracted from traced t-tubules. (H) Mean diameters in junctional (dyad) and nonjunctional (ex-dyad) regions. See main text and the Supporting Material for details. ^**p < 0.01.Taken together, super-resolution imaging and EM tomography strongly support the presence of local t-system dilations in regions where the t-system opposes SR at dyads and such t-system bulges are connected by narrower tubule segments. Further support was provided by SBFSEM, another volume EM technique to study larger cell volumes (albeit at the expense of a slightly lower resolution). SBFSEM clearly showed local t-system dilations were regularly involved in the architecture of most (but not all) dyads (Fig. 3, see also Fig. S3 and Movie S2), as also observed in full three-dimensional super-resolution images (see Fig. S3 C).Open in a separate windowFigure 3Segmented SBFSEM data showing t-system dilations near dyadic junctions. (A) The overview shows t-system membranes (green) and jSR (red) in a mouse myocyte. (B, enlarged inset from panel A) Thin connecting tubules (arrows) and regular swellings in junctional regions at z-lines.Our data identify local dilations of the t-system associated with dyads in mouse cardiac myocytes. Frequent tubule distensions had been observed especially at the intersections of transverse and axial tubules (7), and constrictions were seen in rabbit myocytes although their relationship to dyads was unknown (8). The increased local t-system lumen near junctions may help reduce the predicted ionic accumulation/depletion during excitation-contraction coupling (9). Alternatively, it might simply be secondary to increasing local membrane area and allow the formation of large area junctions that harbor many RyRs. In connection with this point, it would be interesting to investigate the t-system near junctions in species that have larger average tubule diameters (e.g., human and rabbit (10)), or if this architecture changes in mouse heart failure models where t-tubule diameters are often increased.Most peripheral couplings were in regions void of surface caveolae, although a small number of RyR clusters were in junctional couplings between subsarcolemmal SR and caveolae as shown both by the low colocalization between CAV3 and RyRs as well as direct evidence from EM tomography. Similarly, a relatively small fraction of CAV3 colocalized with RyR clusters in the t-system although CAV3 was expressed widely in the t-system. A structural role of CAV3 in the t-system is still unclear—t-tubules in tomogram data did not reveal distinct caveolae shapes on the t-system membrane (see Fig. S4), although this might change in pathology (11). In any case, the t-system exhibits high curvature orthogonal to the tubule axis, which may be supported by CAV3 oligomerization. In addition, the presence of CAV3 in the t-system may be important for regulating other signaling systems (e.g., adrenergic signaling).Finally, our data demonstrate that complementary data from optical super-resolution and three-dimensional EM images assists data interpretation and reliability. We suggest that truly correlative optical and EM imaging approaches should provide further information and improve our knowledge of the basis of cardiac excitation-contraction coupling.     

Inhibiting neuropeptide Y Y1 receptor modulates melanocortin receptor- and NF-κB-mediated feeding behavior in phenylpropanolamine-treated rats

Neuropeptide Y (NPY) and nuclear factor-kappa B (NF-κB) are involved in regulating anorexia elicited by phenylpropanolamine (PPA), a sympathomimetic drug. This study explored whether NPY Y1 receptor (Y1R) is involved in this process, and a potential role for the proopiomelanocortin system was identified. Rats were given PPA once a day for 4 days. Changes in the hypothalamic expression of the NPY, Y1R, NF-κB, and melanocortin receptor 4 (MC4R) levels were assessed and compared. The results indicated that food intake and NPY expression decreased, with the largest reductions observed on Day 2 (approximately 50% and 45%, respectively), whereas NF-κB, MC4R, and Y1R increased, achieving maximums on Day 2 (160%, 200%, and 280%, respectively). To determine the role of Y1R, rats were pretreated with Y1R antisense or a Y1R antagonist via intracerebroventricular injection 1 h before the daily PPA dose. Y1R knockdown and inhibition reduced PPA anorexia and partially restored the normal expression of NPY, MC4R, and NF-κB. The data suggest that hypothalamic Y1R participates in the appetite-suppression from PPA by regulating MC4R and NF-κB. The results of this study increase our understanding of the molecular mechanisms in PPA-induced anorexia.     

Effect of plant stilbene precursors on the biosynthesis of resveratrol in Vitis amurensis Rupr. Cell cultures

The biosynthesis of resveratrol after the application of a precursor for biosynthesis, i.e., phenylalanine (Phe), has been studied. The application of Phe has been shown to increase significantly the expression of the phenylalanine-ammonia-lyase (PAL) and stilbene synthase (STS) genes and enhance the production of resveratrol by 8.5 times. Data on resveratrol production after the addition of Phe and coumaric acid (CA) were compared with known analogs.