Large-Scale Genome Comparison Based on Cumulative Fourier Power and Phase Spectra: Central Moment and Covariance Vector.
Ontology highlight
ABSTRACT:
SUBMITTER: Pei S
PROVIDER: S-EPMC6661692 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:BackgroundFourier transforms and their associated power spectra are used for detecting periodicities and protein-coding genes and is generally regarded as a well established technique. Many of the periodicities which have been found with this method are quite well understood such as the periodicity of 3 nt which is associated to codon usage. But what is the origin of the peculiar frequency multiples k/21 which were reported for a tiny section of chromosome 2 in P. falciparum? Are these present in other chromosomes and perhaps in related organisms? And how should we interpret fractional periodicities in genomes?ResultsWe applied the binary indicator power spectrum to all chromosomes of P. falciparum, and found that the frequency overtones k/21 are present only in non-coding sections. We did not find such frequency overtones in any other related genomes. Furthermore, the frequency overtones were identified as artifacts of the way the genome is encoded into a numerical sequence, that is, they are frequency aliases. By choosing a different way to encode the sequence the overtones do not appear. In view of these results, we revisited early applications of this technique to proteins where frequency overtones were reported.ConclusionsSome authors hinted recently at the possibility of mapping artifacts and frequency aliases in power spectra. However, in the case of P. falciparum the frequency aliases are particularly strong and can mask the 1/3 frequency which is used for gene detecting. This shows that albeit being a well known technique, with a long history of application in proteins, few researchers seem to be aware of the problems represented by frequency aliases.
| S-EPMC3287587 | biostudies-literature
Project description:Nuclear magnetic resonance (NMR) studies of larger proteins are hampered by difficulties in assigning NMR resonances. Human intervention is typically required to identify NMR signals in 3D spectra, and subsequent procedures depend on the accuracy of this so-called peak picking. We present a method that provides sequential connectivities through correlation maps constructed with covariance NMR, bypassing the need for preliminary peak picking. We introduce two novel techniques to minimize false correlations and merge the information from all original 3D spectra. First, we take spectral derivatives prior to performing covariance to emphasize coincident peak maxima. Second, we multiply covariance maps calculated with different 3D spectra to destroy erroneous sequential correlations. The maps are easy to use and can readily be generated from conventional triple-resonance experiments. Advantages of the method are demonstrated on a 37 kDa nonribosomal peptide synthetase domain subject to spectral overlap.
| S-EPMC4183633 | biostudies-literature
Project description:MotivationStochastic molecular processes are a leading cause of cell-to-cell variability. Their dynamics are often described by continuous-time discrete-state Markov chains and simulated using stochastic simulation algorithms. As these stochastic simulations are computationally demanding, ordinary differential equation models for the dynamics of the statistical moments have been developed. The number of state variables of these approximating models, however, grows at least quadratically with the number of biochemical species. This limits their application to small- and medium-sized processes.ResultsIn this article, we present a scalable moment-closure approximation (sMA) for the simulation of statistical moments of large-scale stochastic processes. The sMA exploits the structure of the biochemical reaction network to reduce the covariance matrix. We prove that sMA yields approximating models whose number of state variables depends predominantly on local properties, i.e. the average node degree of the reaction network, instead of the overall network size. The resulting complexity reduction is assessed by studying a range of medium- and large-scale biochemical reaction networks. To evaluate the approximation accuracy and the improvement in computational efficiency, we study models for JAK2/STAT5 signalling and NF κ B signalling. Our method is applicable to generic biochemical reaction networks and we provide an implementation, including an SBML interface, which renders the sMA easily accessible.Availability and implementationThe sMA is implemented in the open-source MATLAB toolbox CERENA and is available from https://github.com/CERENADevelopers/CERENA .Contactjan.hasenauer@helmholtz-muenchen.de or atefeh.kazeroonian@tum.de.Supplementary informationSupplementary data are available at Bioinformatics online.
| S-EPMC5870627 | biostudies-literature
Project description:Huntington's disease (HD) is a progressive neurodegenerative disorder that can be diagnosed with certainty decades before symptom onset. Studies using structural MRI have identified grey matter (GM) loss predominantly in the striatum, but also involving various cortical areas. So far, voxel-based morphometric studies have examined each brain region in isolation and are thus unable to assess the changes in the interrelation of brain regions. Here, we examined the structural covariance in GM volumes in pre-specified motor, working memory, cognitive flexibility, and social-affective networks in 99 patients with manifest HD (mHD), 106 presymptomatic gene mutation carriers (pre-HD), and 108 healthy controls (HC). After correction for global differences in brain volume, we found that increased GM volume in one region was associated with increased GM volume in another. When statistically comparing the groups, no differences between HC and pre-HD were observed, but increased positive correlations were evident for mHD, relative to pre-HD and HC. These findings could be explained by a HD-related neuronal loss heterogeneously affecting the examined network at the pre-HD stage, which starts to dominate structural covariance globally at the manifest stage. Follow-up analyses identified structural connections between frontoparietal motor regions to be linearly modified by disease burden score (DBS). Moderator effects of disease load burden became significant at a DBS level typically associated with the onset of unequivocal HD motor signs. Together with existing findings from functional connectivity analyses, our data indicates a critical role of these frontoparietal regions for the onset of HD motor signs.
| S-EPMC6867397 | biostudies-literature
Project description:Motivated by applications in phenome-wide association studies (PheWAS), we consider in this paper simultaneous testing of columns of high-dimensional cross-covariance matrices and develop a multiple testing procedure with theoretical guarantees. It is shown that the proposed testing procedure maintains a desired false discovery rate (FDR) and false discovery proportion (FDP) under mild regularity conditions. We also provide results on the magnitudes of the signals that can be detected with high power. Simulation studies demonstrate that the proposed procedure can be substantially more powerful than existing FDR controlling procedures in the presence of correlation of unknown structure. The proposed multiple testing procedure is applied to a PheWAS of two auto-immune genetic markers using a rheumatoid arthritis patient cohort constructed from the electronic medical records of Partners Healthcare System.
| S-EPMC6937133 | biostudies-literature
Project description:International efforts to restore degraded ecosystems will continue to expand over the coming decades, yet the factors contributing to the effectiveness of long-term restoration across large areas remain largely unexplored. At large scales, outcomes are more complex and synergistic than the additive impacts of individual restoration projects. Here, we propose a cumulative-effects conceptual framework to inform restoration design and implementation and to comprehensively measure ecological outcomes. To evaluate and illustrate this approach, we reviewed long-term restoration in several large coastal and riverine areas across the US: the greater Florida Everglades; Gulf of Mexico coast; lower Columbia River and estuary; Puget Sound; San Francisco Bay and Sacramento-San Joaquin Delta; Missouri River; and northeastern coastal states. Evidence supported eight modes of cumulative effects of interacting restoration projects, which improved outcomes for species and ecosystems at landscape and regional scales. We conclude that cumulative effects, usually measured for ecosystem degradation, are also measurable for ecosystem restoration. The consideration of evidence-based cumulative effects will help managers of large-scale restoration capitalize on positive feedback and reduce countervailing effects.
| S-EPMC8597595 | biostudies-literature
Project description:Higher-rank correlation spectroscopy is introduced as an alternative to 3D Fourier-transform (FT) NMR spectroscopy for resonance assignment and molecular structure determination. The method combines standard 2D FT spectra that share a common frequency dimension, such as a 2D (13)C-(1)H HSQC and a 2D (1)H-(1)H TOCSY spectrum, and constructs higher-rank correlation spectra with ultra-high spectral resolution. Spectral overlap along a common dimension, in particular the (1)H dimension, is addressed by a spectral filtering method, which identifies mismatches between the 1(st) and 2(nd) moments of cross-peak profiles. The method, which provides a substantial speed-up over traditional 3D FT spectroscopy while effectively suppressing false peaks, is demonstrated for the triple-rank (13)C-(1)H HSQC-TOCSY spectrum of a cyclic decapeptide with different mixing times. Higher-rank correlation spectroscopy is usefully applicable to the analysis of a wide range of NMR spectra of synthetic and natural products.
| S-EPMC2860149 | biostudies-literature
Project description:To investigate the correlation and stationarity of noise in volumetric computed tomography (CT) using the local discrete noise-power spectrum (NPS) and off-diagonal elements of the covariance matrix of the discrete Fourier transform of noise-only images (denoted Σ(DFT)). Experimental conditions were varied to affect noise correlation and stationarity, the effects were quantified in terms of the NPS and Σ(DFT), and practical considerations in CT performance characterization were identified.Cone-beam CT (CBCT) images were acquired using a benchtop system comprising an x-ray tube and flat-panel detector for a range of acquisition techniques (e.g., dose and x-ray scatter) and three phantom configurations hypothesized to impart distinct effects on the NPS and Σ(DFT): (A) air, (B) a 20-cm-diameter water cylinder with a bowtie filter, and (C) the cylinder without a bowtie filter. The NPS and off-diagonal elements of the Σ(DFT) were analyzed as a function of position within the reconstructions.The local NPS varied systematically throughout the axial plane in a manner consistent with changes in fluence transmitted to the detector and view sampling effects. Variability in fluence was manifest in the NPS magnitude-e.g., a factor of ~2 variation in NPS magnitude within the axial plane for case C (cylinder without bowtie), compared to nearly constant NPS magnitude for case B (bowtie filter matched to the cylinder). View sampling effects were most prominent in case A (air) where the variance increased at greater distance from the center of reconstruction and in case C (cylinder) where the NPS exhibited correlations in the radial direction. The effects of detector lag were observed as azimuthal correlation. The cylinder (without bowtie) had the strongest nonstationarity because of the larger variability in fluence transmitted to the detector. The diagonal elements of the Σ(DFT) were equivalent to the NPS estimated from the periodogram, and the average off-diagonal elements of the Σ(DFT) exhibited amplitude of ~1% of the NPS for the experimental conditions investigated. Furthermore, the off-diagonal elements demonstrated fairly long tails of nearly constant amplitude, with magnitude somewhat reduced for experimental conditions associated with greater stationarity (viz., lower Σ(DFT) tails for cases A and B in comparison to case C).Volumetric CT exhibits nonstationarity in the NPS as hypothesized in relation to fluence uniformity and view sampling. Measurement of the NPS should seek to minimize such changes in noise correlations and include careful reporting of experimental conditions (e.g., phantom design and use of a bowtie filter) and spatial dependence (e.g., analysis at fixed radius within a phantom). Off-diagonal elements of the Σ(DFT) similarly depend on experimental conditions and can be readily computed from the same data as the NPS. This work begins to check assumptions in NPS analysis examine the extent to which NPS is an appropriate descriptor of noise correlations, and investigate the magnitude of off-diagonal elements of the Σ(DFT). While the magnitude of such off-diagonal elements appears to be low, their cumulative effect on space-variant detectability remains to be investigated-e.g., using task-specific figures of merit.
| S-EPMC3365918 | biostudies-other
Project description:BackgroundThe statistical power or multiple Type II error rate in large-scale multiple testing problems as, for example, in gene expression microarray experiments, depends on typically unknown parameters and is therefore difficult to assess a priori. However, it has been suggested to estimate the multiple Type II error rate post hoc, based on the observed data.MethodsWe consider a class of post hoc estimators that are functions of the estimated proportion of true null hypotheses among all hypotheses. Numerous estimators for this proportion have been proposed and we investigate the statistical properties of the derived multiple Type II error rate estimators in an extensive simulation study.ResultsThe performance of the estimators in terms of the mean squared error depends sensitively on the distributional scenario. Estimators based on empirical distributions of the null hypotheses are superior in the presence of strongly correlated test statistics.AvailabilityR-code to compute all considered estimators based on P-values and supplementary material is available on the authors web page http://statistics.msi.meduniwien.ac.at/index.php?page=pageszfnr.
| S-EPMC3500624 | biostudies-literature
Project description:When wind turbines are arranged in clusters, their performance is mutually affected, and their energy generation is reduced relative to what it would be if they were widely separated. Land-area power densities of small wind farms can exceed 10 W/m2, and wakes are several rotor diameters in length. In contrast, large-scale wind farms have an upper-limit power density in the order of 1 W/m2 and wakes that can extend several tens of kilometers. Here, we address two important questions: 1) How large can a wind farm be before its generation reaches energy replenishment limits and 2) How far apart must large wind farms be spaced to avoid inter-wind-farm interference? We characterize controls on these spatial and temporal scales by running a set of idealized atmospheric simulations using the Weather and Research Forecasting model. Power generation and wind speed within and over the wind farm show that a timescale inversely proportional to the Coriolis parameter governs such transition, and the corresponding length scale is obtained by multiplying the timescale by the geostrophic wind speed. A geostrophic wind of 8 m/s and a Coriolis parameter of 1.05 × 10-4 rad/s (latitude of ∼46°) would give a transitional scale of about 30 km. Wind farms smaller than this result in greater power densities and shorter wakes. Larger wind farms result instead in power densities that asymptotically reach their minimum and wakes that reach their maximum extent.
| S-EPMC8271749 | biostudies-literature