Analyzing Clustered Data: Why and How to Account for Multiple Observations Nested within a Study Participant?

A conventional study design among medical and biological experimentalists involves col- lecting multiple measurements from a study subject. For example, experiments utilizing mouse models in neuroscience often involve collecting multiple neuron measurements per mouse to increase the number of observations without requiring a large number of mice. This leads to a form of statistical dependence referred to as clustering. Inappropriate analy- ses of clustered data have resulted in several recent critiques of neuroscience research that suggest the bar for statistical analyses within the field is set too low. We compare naïve ana- lytical approaches to marginal, fixed-effect, and mixed-effect models and provide guidelines for when each of these models is most appropriate based on study design. We demonstrate the influence of clustering on a between-mouse treatment effect, a within-mouse treatment effect, and an interaction effect between the two. Our analyses demonstrate that these sta- tistical approaches can give substantially different results, primarily when the analyses include a between-mouse treatment effect. In a novel analysis from a neuroscience per- spective, we also refine the mixed-effect approach through the inclusion of an aggregate mouse-level counterpart to a within-mouse (neuron level) treatment as an additional predic- tor by adapting an advanced modeling technique that has been used in social science research and show that this yields more informative results. Based on these findings, we emphasize the importance of appropriate analyses of clustered data, and we aim for this work to serve as a resource for when one is deciding which approach will work best for a given study

Gender and the systemic hypertension-snoring association: a questionnaire-based case-control study.

Since the role of gender in the association between hypertension and snoring is unknown, we studied it while accounting for age and body mass index (BMI) as confounding variables. A questionnaire on snoring was administered to 90 hypertensive (HT) subjects (45 men and 45 women) and to 90 normotensive (NT) subjects matched for gender, age and BMI. As expected, snoring was more commonly reported by men than by women, but no significant difference was found between HT and NT men, irrespective of age. Conversely, heavy snoring was more frequently reported by HT than NT women; habitual snoring was more common among young (age50 years) HT than NT women; and heavy snoring was more common among older (age50 years) HT than NT women. These data suggest an effect of gender on the hypertension-snoring association: in men, snoring may be accounted for by age and BMI whether or not hypertension is present, whereas in women the natural history of snoring appears different and more severe in HT than in NT. Although the mechanism(s) responsible for the differences between men and women are obscure at present, gender may be an important variable in the systemic hypertension-snoring association

Validity of needle core biopsy in the histological characterisation of mammary lesions

Over the last few years, there has been an enormous increase in the use of needle core biopsy (CB) for the histopathological characterisation of suspect lesions of the breast. The aim of this study was to verify the diagnostic reliability of CB by comparing the histological results obtained with the use of this technique with those obtained from the whole of the surgically resected specimen. We studied 198 out of 426 patients with clinically and/or radiologically suspect breast lesions. We found correspondence between the histological examination of the whole of the excised specimen and that of the CB in 94.9% of the cases of infiltrating carcinoma and in 71.4% of those involving ductal carcinomas in situ. The predictive value of CB was 98.9%, sensitivity was 96.1% and specificity 93.3%. These results confirmed that CB is an extremely reliable diagnostic tool in the definition of breast lesions

Microwave Properties of Ba(0.6)K(0.4)BiO(3) Crystals

We report on field-induced variations of the microwave surface resistance at 9.6 GHz of Ba(0.6)K(0.4)BiO(3) crystals. Energy losses have been investigated as a function of the static magnetic field in the range of temperatures 4.2 K - Tc. By analyzing the experimental results in the framework of the Coffey and Clem model we determine the temperature dependence of the first-penetration field, upper critical field and depinning frequency. The results show that the pinning energy of this bismuthate superconductor is weaker than those of cuprates.Comment: 6 pages, 8 embedded figure

Genomic prediction of grain yield in a barley MAGIC population modelling genotype per environment interaction

Multi-parent Advanced Generation Inter-crosses (MAGIC) lines have mosaic genomes that are generated shuffling the genetic material of the founder parents following predefined crossing schemes. In cereal crops, these experimental populations have been extensively used to investigate the genetic bases of several traits and dissect the genetic bases of epistasis. In plants, genomic prediction models are usually fitted using either diverse panels of mostly unrelated accessions or individuals of biparental families and several empirical analyses have been conducted to evaluate the predictive ability of models fitted to these populations using different traits. In this paper, we constructed, genotyped and evaluated a barley MAGIC population of 352 individuals developed with a diverse set of eight founder parents showing contrasting phenotypes for grain yield. We combined phenotypic and genotypic information of this MAGIC population to fit several genomic prediction models which were cross-validated to conduct empirical analyses aimed at examining the predictive ability of these models varying the sizes of training populations. Moreover, several methods to optimize the composition of the training population were also applied to this MAGIC population and cross-validated to estimate the resulting predictive ability. Finally, extensive phenotypic data generated in field trials organized across an ample range of water regimes and climatic conditions in the Mediterranean were used to fit and cross-validate multi-environment genomic prediction models including GE interaction, using both genomic best linear unbiased prediction and reproducing kernel Hilbert space along with a non-linear Gaussian Kernel. Overall, our empirical analyses showed that genomic prediction models trained with a limited number of MAGIC lines can be used to predict grain yield with values of predictive ability that vary from 0.25 to 0.60 and that beyond QTL mapping and analysis of epistatic effects, MAGIC population might be used to successfully fit genomic prediction models. We concluded that for grain yield, the single-environment genomic prediction models examined in this study are equivalent in terms of predictive ability while, in general, multi-environment models that explicitly split marker effects in main and environmentalspecific effects outperform simpler multi-environment models

Physical Mapping of Bread Wheat Chromosome 5A: An Integrated Approach

The huge size, redundancy, and highly repetitive nature of the bread wheat [Triticum aestivum (L.)] genome, makes it among the most difficult species to be sequenced. To overcome these limitations, a strategy based on the separation of individual chromosomes or chromosome arms and the subsequent production of physical maps was established within the frame of the International Wheat Genome Sequence Consortium (IWGSC). A total of 95,812 bacterial artificial chromosome (BAC) clones of short-arm chromosome 5A (5AS) and long-arm chromosome 5A (5AL) arm-specific BAC libraries were fingerprinted and assembled into contigs by complementary analytical approaches based on the FingerPrinted Contig (FPC) and Linear Topological Contig (LTC) tools. Combined anchoring approaches based on polymerase chain reaction (PCR) marker screening, microarray, and sequence homology searches applied to several genomic tools (i. e., genetic maps, deletion bin map, neighbor maps, BAC end sequences (BESs), genome zipper, and chromosome survey sequences) allowed the development of a high-quality physical map with an anchored physical coverage of 75% for 5AS and 53% for 5AL with high portions (64 and 48%, respectively) of contigs ordered along the chromosome. In the genome of grasses, Brachypodium [Brachypodium distachyon (L.) Beauv.], rice (Oryza sativa L.), and sorghum [Sorghum bicolor (L.) Moench] homologs of genes on wheat chromosome 5A were separated into syntenic blocks on different chromosomes as a result of translocations and inversions during evolution. The physical map presented represents an essential resource for fine genetic mapping and map-based cloning of agronomically relevant traits and a reference for the 5A sequencing projects

Differential responses of zooplankton assemblages to environmental variation in temporary and permanent ponds

Permanent and temporary wetlands in Mediterranean shrublands represent unique repositories of biodiversity, which are increasingly threatened by human-induced habitat loss. The zooplankton of a permanent (P1) and a temporary pond (T35) in the Natural Reserve of Castelporziano, a rare residual stretch of such a shrubland in Central Italy (Latium), was investigated to: (1) expand and deepen knowledge of these endangered freshwater habitats, which represent a crucial component of Mediterranean biodiversity; (2) identify environmental controls regulating the development of zooplankton communities of each environment; and (3) highlight differences in the adaptive responses of the zooplankton community in relation to the different ecological conditions experienced by permanent and temporary habitats. Despite summer desiccation in T35, the two ponds exhibited a relative homogeneity in hydrological and physico-chemical dynamics. Zooplankton assemblages contained 41 total taxa, of which 32 were found in P1 and 28 in T35. Out of the 41 taxa identified, 22 (> 50%) were exclusively present in one of the two ponds. On a yearly basis, the community dynamics of P1 seemed to be conditioned by physical and chemical factors and by hydrological cycle characteristics, while the community of T35 responded to algal blooms, food competition and predator/prey equilibria rather than correlating to abiotic factors. The main differences amongst zooplankton assemblages were observed over short time scales and occurred both within and between seasons, highlighting the role of some structural taxa that dominated the average composition of the community throughout the year, and the importance of "quick-response" taxa in determining the short-term composition and structure variation of pond zooplankton. A year-round cyclic community succession peculiar to each pond is described

Lipoprotein Lipase Links Dietary Fat to Solid Tumor Cell Proliferation

Many types of cancer cells require a supply of fatty acids (FA) for growth and survival, and interrupting de novo FA synthesis in model systems causes potent anticancer effects. We hypothesized that, in addition to synthesis, cancer cells may obtain pre-formed, diet-derived fatty acids by uptake from the bloodstream. This would require hydrolytic release of FA from triglyceride in circulating lipoprotein particles by the secreted enzyme lipoprotein lipase (LPL), and the expression of CD36, the channel for cellular FA uptake. We find that selected breast cancer and sarcoma cells express and secrete active LPL, and all express CD36. We further demonstrate that LPL, in the presence of triglyceride-rich lipoproteins, accelerates the growth of these cells. Providing LPL to prostate cancer cells, which express low levels of the enzyme, did not augment growth, but did prevent the cytotoxic effect of FA synthesis inhibition. Moreover, LPL knockdown inhibited HeLa cell growth. In contrast to the cell lines, immunohistochemical analysis confirmed the presence of LPL and CD36 in the majority of breast, liposarcoma, and prostate tumor tissues examined (n = 181). These findings suggest that, in addition to de novo lipogenesis, cancer cells can use LPL and CD36 to acquire FA from the circulation by lipolysis, and this can fuel their growth. Interfering with dietary fat intake, lipolysis, and/or fatty acid uptake will be necessary to target the requirement of cancer cells for FA

Genome-wide analysis of japonica rice performance under limited water and permanent flooding conditions

A rice GWAS panel of 281 accessions of japonica rice was phenotypically characterized for 26 traits related to phenology, plant and seed morphology, physiology and yield for 2 years in field conditions under permanent flooding (PF) and limited water (LW). A genome-wide analysis uncovered a total of 160 significant marker-trait associations (MTAs), of which 32 were LW-specific, 59 were PF-specific, and 69 were in common between the two water management systems. LW-specific associations were identified for several agronomic traits including days to maturation, days from flowering to maturation, leaf traits, plant height, panicle and seed traits, hundred grain weight, yield and tillering. Significant MTAs were detected across all the 12 rice chromosomes, while clusters of effects influencing different traits under LW or in both watering conditions were, respectively, observed on chromosomes 4, 8, and 12 and on chromosomes 1, 3, 4, 5, and 8. The analysis of genes annotated in the Nipponbare reference sequence and included in the regions associated to traits related to plant morphology, grain yield, and physiological parameters allowed the identification of genes that were demonstrated to affect the respective traits. Among these, three (OsOFP2, Dlf1, OsMADS56) and seven (SUI1, Sd1, OsCOL4, Nal1, OsphyB, GW5, Ehd1) candidate genes were, respectively, identified to co-localize with LW-specific associations and associations in common between the two water treatments. For several LW-specific MTAs, or in common among the two treatments, positional co-localizations with previously identified QTLs for rice adaptation to water shortages were observed, a result that further supports the role of the loci identified in this work in conferring adaptation to LW. The most robust associations identified here could represent suitable targets for genomic selection approaches to improve yield-related traits under LW

Genetic variants of HvCbf14 are statistically associated with frost tolerance in a European germplasm collection of Hordeum vulgare

Two quantitative trait loci (Fr-H1 and Fr-H2) for frost tolerance (FT) have been discovered on the long arm of chromosome 5H in barley. Two tightly linked groups of CBF genes, known to play a key role in the FT regulatory network in A. thaliana, have been found to co-segregate with Fr-H2. Here, we investigate the allelic variations of four barley CBF genes (HvCbf3, HvCbf6, HvCbf9 and HvCbf14) in a panel of European cultivars, landraces and H. spontaneum accessions. In the cultivars a reduction of nucleotide and haplotype diversities in CBFs compared with the landraces and the wild ancestor H. spontaneum, was evident. In particular, in cultivars the loss of HvCbf9 genetic variants was higher compared to other sequences. In order to verify if the pattern of CBF genetic variants correlated with the level of FT, an association procedure was adopted. The pairwise analysis of linkage disequilibrium (LD) among the genetic variants in four CBF genes was computed to evaluate the resolution of the association procedure. The pairwise plotting revealed a low level of LD in cultivated varieties, despite the tight physical linkage of CBF genes analysed. A structured association procedure based on a general liner model was implemented, including the variants in CBFs, of Vrn-H1, and of two reference genes not involved in FT (α-Amy1 and Gapdh) and considering the phenotypic data for FT. Association analysis recovered two nucleotide variants of HvCbf14 and one nucleotide variant of Vrn-H1 as statistically associated to FT