Multivariate variance-components analysis in DTI

Twin studies are a major research direction in imaging genetics, a new field, which combines algorithms from quantitative genetics and neuroimaging to assess genetic effects on the brain. In twin imaging studies, it is common to estimate the intraclass correlation (ICC), which measures the resemblance between twin pairs for a given phenotype. In this paper, we extend the commonly used Pearson correlation to a more appropriate definition, which uses restricted maximum likelihood methods (REML). We computed proportion of phenotypic variance due to additive (A) genetic factors, common (C) and unique (E) environmental factors using a new definition of the variance components in the diffusion tensor-valued signals. We applied our analysis to a dataset of Diffusion Tensor Images (DTI) from 25 identical and 25 fraternal twin pairs. Differences between the REML and Pearson estimators were plotted for different sample sizes, showing that the REML approach avoids severe biases when samples are smaller. Measures of genetic effects were computed for scalar and multivariate diffusion tensor derived measures including the geodesic anisotropy (tGA) and the full diffusion tensors (DT), revealing voxel-wise genetic contributions to brain fiber microstructure.

Analysis of differential gene expression based on Bayesian estimation of variance

Gene expression is arguably the most important indicator of biological function. Thus identifying differentially expressed genes is one of the main aims of high throughout studies that use microarray and RNAseq platforms to study deregulated cellular pathways. There are many tools for analysing differentia gene expression from transciptomic datasets. The major challenge of this topic is to estimate gene expression variance due to the high amount of ‘background noise’ that is generated from biological equipment and the lack of biological replicates. Bayesian inference has been widely used in the bioinformatics field. In this work, we reveal that the prior knowledge employed in the Bayesian framework also helps to improve the accuracy of differential gene expression analysis when using a small number of replicates. We have developed a differential analysis tool that uses Bayesian estimation of the variance of gene expression for use with small numbers of biological replicates. Our method is more consistent when compared to the widely used cyber-t tool that successfully introduced the Bayesian framework to differential analysis. We also provide a user-friendly web based Graphic User Interface for biologists to use with microarray and RNAseq data. Bayesian inference can compensate for the instability of variance caused when using a small number of biological replicates by using pseudo replicates as prior knowledge. We also show that our new strategy to select pseudo replicates will improve the performance of the analysis. - See more at: http://www.eurekaselect.com/node/138761/article#sthash.VeK9xl5k.dpuf

Enhanced sensitivity in detection of Kerr rotation by double modulation and time averaging based on Allan variance

Experiments in spintronics necessarily involve the detection of spin polarization. The sensitivity of this detection becomes an important factor to consider when extending the low temperature studies on semiconductor spintronic devices to room temperature, where the spin signal is weaker. In pump-probe experiments, which optically inject and detect spins, the sensitivity is often improved by using a photoelastic modulator (PEM) for lock-in detection. However, spurious signals can arise if diode lasers are used as optical sources in such experiments, along with a PEM. In this work, we eliminated the spurious electromagnetic coupling of the PEM onto the probe diode laser, by the double modulation technique. We also developed a test for spurious modulated interference in the pump-probe signal, due to the PEM. Besides, an order of magnitude enhancement in the sensitivity of detection of spin polarization by Kerr rotation, to 3x10(-8) rad was obtained by using the concept of Allan variance to optimally average the time series data over a period of 416 s. With these improvements, we are able to experimentally demonstrate at room temperature, photoinduced steady-state spin polarization in bulk GaAs. Thus, the advances reported here facilitate the use of diode lasers with a PEM for sensitive pump-probe experiments. They also constitute a step toward detection of spin-injection in Si at room temperature.

Modeling strategies in longitudinal data analysis: Covariate, variance function and correlation structure selection

A modeling paradigm is proposed for covariate, variance and working correlation structure selection for longitudinal data analysis. Appropriate selection of covariates is pertinent to correct variance modeling and selecting the appropriate covariates and variance function is vital to correlation structure selection. This leads to a stepwise model selection procedure that deploys a combination of different model selection criteria. Although these criteria find a common theoretical root based on approximating the Kullback-Leibler distance, they are designed to address different aspects of model selection and have different merits and limitations. For example, the extended quasi-likelihood information criterion (EQIC) with a covariance penalty performs well for covariate selection even when the working variance function is misspecified, but EQIC contains little information on correlation structures. The proposed model selection strategies are outlined and a Monte Carlo assessment of their finite sample properties is reported. Two longitudinal studies are used for illustration.

Effects of variance-function misspecification in analysis of longitudinal data

The approach of generalized estimating equations (GEE) is based on the framework of generalized linear models but allows for specification of a working matrix for modeling within-subject correlations. The variance is often assumed to be a known function of the mean. This article investigates the impacts of misspecifying the variance function on estimators of the mean parameters for quantitative responses. Our numerical studies indicate that (1) correct specification of the variance function can improve the estimation efficiency even if the correlation structure is misspecified; (2) misspecification of the variance function impacts much more on estimators for within-cluster covariates than for cluster-level covariates; and (3) if the variance function is misspecified, correct choice of the correlation structure may not necessarily improve estimation efficiency. We illustrate impacts of different variance functions using a real data set from cow growth.

Method for comparing the capture efficiency of benthic sampling devices

Traditional comparisons between the capture efficiency of sampling devices have generally looked at the absolute differences between devices. We recommend that the signal-to-noise ratio be used when comparing the capture efficiency of benthic sampling devices. Using the signal-to-noise ratio rather than the absolute difference has the advantages that the variance is taken into account when determining how important the difference is, the hypothesis and minimum detectable difference can be made identical for all taxa, it is independent of the units used for measurement, and the sample-size calculation is independent of the variance. This new technique is illustrated by comparing the capture efficiency of a 0.05 m(2) van Veen grab and an airlift suction device, using samples taken from Heron and One Tree lagoons, Australia.

Genome-wide association study of major depressive disorder: New results, meta-analysis, and lessons learned

Major depressive disorder (MDD) is a common complex disorder with a partly genetic etiology. We conducted a genome-wide association study of the MDD2000+ sample (2431 cases, 3673 screened controls and >1 M imputed single-nucleotide polymorphisms (SNPs)). No SNPs achieved genome-wide significance either in the MDD2000+ study, or in meta-analysis with two other studies totaling 5763 cases and 6901 controls. These results imply that common variants of intermediate or large effect do not have main effects in the genetic architecture of MDD. Suggestive but notable results were: (a) gene-based tests suggesting roles for adenylate cyclase 3 (ADCY3, 2p23.3) and galanin (GAL, 11q13.3); published functional evidence relates both of these to MDD and serotonergic signaling; (b) support for the bipolar disorder risk variant SNP rs1006737 in CACNA1C (P=0.020, odds ratio=1.10), and; (c) lack of support for rs2251219, a SNP identified in a meta-analysis of affective disorder studies (P=0.51). We estimate that sample sizes 1.8- to 2.4-fold greater are needed for association studies of MDD compared with those for schizophrenia to detect variants that explain the same proportion of total variance in liability. Larger study cohorts characterized for genetic and environmental risk factors accumulated prospectively are likely to be needed to dissect more fully the etiology of MDD.

Six novel susceptibility Loci for early-onset androgenetic alopecia and their unexpected association with common diseases

Androgenetic alopecia (AGA) is a highly heritable condition and the most common form of hair loss in humans. Susceptibility loci have been described on the X chromosome and chromosome 20, but these loci explain a minority of its heritable variance. We conducted a large-scale meta-analysis of seven genome-wide association studies for early-onset AGA in 12,806 individuals of European ancestry. While replicating the two AGA loci on the X chromosome and chromosome 20, six novel susceptibility loci reached genome-wide significance (p = 2.62x10(-)(9)-1.01x10(-)(1)(2)). Unexpectedly, we identified a risk allele at 17q21.31 that was recently associated with Parkinson's disease (PD) at a genome-wide significant level. We then tested the association between early-onset AGA and the risk of PD in a cross-sectional analysis of 568 PD cases and 7,664 controls. Early-onset AGA cases had significantly increased odds of subsequent PD (OR = 1.28, 95% confidence interval: 1.06-1.55, p = 8.9x10(-)(3)). Further, the AGA susceptibility alleles at the 17q21.31 locus are on the H1 haplotype, which is under negative selection in Europeans and has been linked to decreased fertility. Combining the risk alleles of six novel and two established susceptibility loci, we created a genotype risk score and tested its association with AGA in an additional sample. Individuals in the highest risk quartile of a genotype score had an approximately six-fold increased risk of early-onset AGA [odds ratio (OR) = 5.78, p = 1.4x10(-)(8)(8)]. Our results highlight unexpected associations between early-onset AGA, Parkinson's disease, and decreased fertility, providing important insights into the pathophysiology of these conditions.

A variant in LIN28B is associated with 2D:4D finger-length ratio, a putative retrospective biomarker of prenatal testosterone exposure

The ratio of the lengths of an individual's second to fourth digit (2D:4D) is commonly used as a noninvasive retrospective biomarker for prenatal androgen exposure. In order to identify the genetic determinants of 2D:4D, we applied a genome-wide association approach to 1507 11-year-old children from the Avon Longitudinal Study of Parents and Children (ALSPAC) in whom 2D:4D ratio had been measured, as well as a sample of 1382 12- to 16-year-olds from the Brisbane Adolescent Twin Study. A meta-analysis of the two scans identified a single variant in the LIN28B gene that was strongly associated with 2D:4D (rs314277: p = 4.1 x 10(-8)) and was subsequently independently replicated in an additional 3659 children from the ALSPAC cohort (p = 1.53 x 10(-6)). The minor allele of the rs314277 variant has previously been linked to increased height and delayed age at menarche, but in our study it was associated with increased 2D:4D in the direction opposite to that of previous reports on the correlation between 2D:4D and age at menarche. Our findings call into question the validity of 2D:4D as a simplistic retrospective biomarker for prenatal testosterone exposure.

Dosage Compensation And Sex Determination In Drosophila - Mechanism Of Measurement Of The X/A Ratio

We propose a molecular mechanism for the intra-cellular measurement of the ratio of the number of X chromosomes to the number of sets of autosomes, a process central to both sex determination and dosage compensation in Drosophila melanogaster. In addition to the two loci, da and Sxl, which have been shown by Cline (Genetics, 90, 683, 1978)and others to be involved in these processes, we postulate two other loci, one autosomal (ω) and the other, X-linked (π). The product of the autosomal locus da stimulates ω and initiates synthesis of a limited quantity of repressor. Sxl and π ,both of which are X-linked, compete for this repressor as well as for RNA polymerase. It is assumed that Sxl has lower affinity than π for repressor as well as polymerase and that the binding of polymerase to one of these sites modulates the binding affinity of the other site for the enzyme. It can be shown that as a result of these postulated interactions transcription from the Sxl site is proportional to the X/A ratio such that the levels of Sxl+ product are low in males, high in females and intermediate in the intersexes. If, as proposed by Cline, the Sxl- product is an inhibitor of X chromosome activity, this would result in dosage compensation. The model leads to the conclusion that high levels of Sxl+ product promote a female phenotype and low levels, a male phenotype. One interesting consequence of the assumptions on which the model is based is that the level of Sxl+ product in the cell, when examined as a function of increasing repressor concentration, first goes up and then decreases, yielding a bell-shaped curve. This feature of the model provides an explanation for some of the remarkable interactions among mutants at the Sxl, da and mle loci and leads to several predictions. The proposed mechanism may also have relevance to certain other problems, such as size regulation during development, which seem to involve measurement of ratios at the cellular level.

Brownian particles in stationary and moving traps: The mean and variance of the heat distribution function

A recent theoretical model developed by Imparato et al. Phys of the experimentally measured heat and work effects produced by the thermal fluctuations of single micron-sized polystyrene beads in stationary and moving optical traps has proved to be quite successful in rationalizing the observed experimental data. The model, based on the overdamped Brownian dynamics of a particle in a harmonic potential that moves at a constant speed under a time-dependent force, is used to obtain an approximate expression for the distribution of the heat dissipated by the particle at long times. In this paper, we generalize the above model to consider particle dynamics in the presence of colored noise, without passing to the overdamped limit, as a way of modeling experimental situations in which the fluctuations of the medium exhibit long-lived temporal correlations, of the kind characteristic of polymeric solutions, for instance, or of similar viscoelastic fluids. Although we have not been able to find an expression for the heat distribution itself, we do obtain exact expressions for its mean and variance, both for the static and for the moving trap cases. These moments are valid for arbitrary times and they also hold in the inertial regime, but they reduce exactly to the results of Imparato et al. in appropriate limits. DOI: 10.1103/PhysRevE.80.011118 PACS.

Meta-analysis identifies 13 new loci associated with waist-hip ratio and reveals sexual dimorphism in the genetic basis of fat distribution

Waist-hip ratio (WHR) is a measure of body fat distribution and a predictor of metabolic consequences independent of overall adiposity. WHR is heritable, but few genetic variants influencing this trait have been identified. We conducted a meta-analysis of 32 genome-wide association studies for WHR adjusted for body mass index (comprising up to 77,167 participants), following up 16 loci in an additional 29 studies (comprising up to 113,636 subjects). We identified 13 new loci in or near RSPO3, VEGFA, TBX15-WARS2, NFE2L3, GRB14, DNM3-PIGC, ITPR2-SSPN, LY86, HOXC13, ADAMTS9, ZNRF3-KREMEN1, NISCH-STAB1 and CPEB4 (P = 1.9 × 10−9 to P = 1.8 × 10−40) and the known signal at LYPLAL1. Seven of these loci exhibited marked sexual dimorphism, all with a stronger effect on WHR in women than men (P for sex difference = 1.9 × 10−3 to P = 1.2 × 10−13). These findings provide evidence for multiple loci that modulate body fat distribution independent of overall adiposity and reveal strong gene-by-sex interactions.

Genetic effects and genetic relationships among shrunken (sh2) sweet corn lines and F1 hybrids

The objectives of this study were to quantify the components of genetic variance and the genetic effects, and to examine the genetic relationship of inbred lines extracted from various shrunken2(sh2) breeding populations. Ten diverse inbred lines developed from genetic background, were crossed in half diallel. Parents and their F1 hybrids were evaluated at three environments. The parents were genotyped using 20 polymorphic simple sequence repeats (SSR). Agronomic and quality traits were analysed by a mixed linear model according to additive-dominance genetic model. Genetic effects were estimated using an adjusted unbiased prediction method. Additive variance was more important than dominance variance in the expression of traits related to ear aspects (husk ratio and percentage of ear filled) and eating quality (flavour and total soluble solids). For agronomic traits, however, dominance variance was more important than additive variance. The additive genetic correlation between flavour and tenderness was strong (r = 0.84, P <0.01). Flavour, tenderness and kernel colour additive genetic effects were not correlated with yield related traits. Genetic distance (GD), estimated from SSR profiles on the basis of Jaccard's similarity coefficient varied from 0.10 to 0.77 with an average of 0.56. Cluster analysis classified parents according to their pedigree relationships. In most studied traits, F1 performance was not associated with GD.

Resilience of a eucalypt forest woody understorey to long-term (34–55 years) repeated burning in subtropical Australia

We investigated the effects of annual burning since 1952, triennial burning since 1973, fire exclusion since 1946 and infrequent wildfire (one fire in 61 years) on woody understorey vegetation in a dry sclerophyll eucalypt forest, south-eastern Queensland, Australia. We determined the influence of these treatments, and other site variables (rainfall, understorey density, topsoil C : N ratio, tree basal area, distance to watercourse and burn coverage) on plant taxa density, richness and composition. The richness of woody understorey taxa 0–1 m in height was not affected by burning treatments, but richness of woody plants 1–7.5 m in height was lower in the annually burnt treatment than in the triennially burnt treatment from 1989 to 2007. Fire frequency and other site variables explained 34% of the variation in taxa composition (three taxon groups and 10 species), of which 33% of the explained variance was explained by fire treatment and 46% was explained by other site variables. Annual burning between 1974 and 1993 was associated with lower understorey densities mainly due to reduced densities of eucalypts 1–7.5 m in height. Triennial burning during the same period was associated with higher densities of eucalypts 0–7.5 m in height relative to the annually burnt and unburnt treatments. Most woody taxa persisted in the frequently burnt treatments through resprouting mechanisms (e.g. lignotuberous regeneration), and fire patchiness associated with low-intensity burning was also found to be important. Persistence of plants <1 m tall demonstrates the resilience of woody taxa to repeated burning in this ecosystem, although they mainly exist in a suppressed growth state under annual burning.

Effects of dietary phosphorus and calcium-to-phosphorus ratio on calcium and bone metabolism in healthy 20- to 43-year-old Finnish women

Dietary habits have changed during the past decades towards an increasing consumption of processed foods, which has notably increased not only total dietary phosphorus (P) intake, but also intake of P from phosphate additives. While the intake of calcium (Ca) in many Western countries remains below recommended levels (800 mg/d), the usual daily P intake in a typical Western diet exceeds by 2- to 3-fold the dietary guidelines (600 mg/d). The effects of high P intake in healthy humans have been investigated seldom. In this thesis healthy 20- to 43-year-old women were studied. In the first controlled study (n = 14), we examined the effects of P doses, and in a cross-sectional study (n = 147) the associations of habitual P intakes with Ca and bone metabolism. In this same cross-sectional study, we also investigated whether differences exist between dietary P originating from natural P sources and phosphate additives. The second controlled study (n = 12) investigated whether by increasing the Ca intake, the effects of a high P intake could be reduced. The associations of habitual dietary calcium-to-phosphorus ratios (Ca:P ratio) with Ca and bone metabolism were determined in a cross-sectional study design (n = 147). In the controlled study, the oral intake of P doses (495, 745, 1245 and 1995 mg/d) with a low Ca intake (250 mg/d) increased serum parathyroid hormone (S-PTH) concentration in a dose-dependent manner. In addition, the highest P dose decreased serum ionized calcium (S-iCa) concentration and bone formation and increased bone resorption. In the second controlled study with a dietary P intake of 1850 mg/d, by increasing the Ca intake from 480 mg/d to 1080 mg/d and then to 1680 mg/d, the S-PTH concentration decreased, the S-iCa concentration increased and bone resorption decreased dose-dependently. However, not even the highest Ca intake could counteract the effect of high dietary P on bone formation, as indicated by unchanged bone formation activity. In the cross-sectional studies, a higher habitual dietary P intake (>1650 mg/d) was associated with lower S-iCa and higher S-PTH concentrations. The consumption of phosphate additive-containing foods was associated with a higher S-PTH concentration. Moreover, habitual low dietary Ca:P ratios (≤0.50, molar ratio) were associated with higher S-PTH concentrations and 24-h urinary Ca excretions, suggesting that low dietary Ca:P ratios may interfere with homeostasis of Ca metabolism and increase bone resorption. In summary, excessive dietary P intake in healthy Finnish women seems to be detrimental to Ca and bone metabolism, especially when dietary Ca intake is low. The results indicate that by increasing dietary Ca intake to the recommended level, the negative effects of high P intake could be diminished, but not totally prevented. These findings imply that phosphate additives may be more harmful than natural P. Thus, reduction of an excessively high dietary P intake is also beneficial for healthy individuals.