PubMed Health⌕ Search

Biomedical subjects

Gary Churchill

Publications and source records attributed to Gary Churchill.

10 recordsLinked to original sources

Estimating p-values in small microarray experiments.

MOTIVATION: Microarray data typically have small numbers of observations per gene, which can result in low power for statistical tests. Test statistics that borrow information from data across all of the genes can improve power, but these statistics have non-standard distributions, and their significance must be assessed using permutation analysis. When sample sizes are small, the number of distinct permutations can be severely limited, and pooling the permutation-derived test statistics across all genes has been proposed. However, the null distribution of the test statistics under permutation is not the same for equally and differentially expressed genes. This can have a negative impact on both p-value estimation and the power of information borrowing statistics. RESULTS: We investigate permutation based methods for estimating p-values. One of methods that uses pooling from a selected subset of the data are shown to have the correct type I error rate and to provide accurate estimates of the false discovery rate (FDR). We provide guidelines to select an appropriate subset. We also demonstrate that information borrowing statistics have substantially increased power compared to the t-test in small experiments.

Algorithms↗

Complex genetic architecture revealed by analysis of high-density lipoprotein cholesterol in chromosome substitution strains and F2 crosses.

Intercrosses between inbred lines provide a traditional approach to analysis of polygenic inheritance in model organisms. Chromosome substitution strains (CSSs) have been developed as an alternative to accelerate the pace of gene identification in quantitative trait mapping. We compared a classical intercross and three CSS intercrosses to examine the genetic architecture underlying plasma high-density lipoprotein cholesterol (HDL) levels in the C57BL/6J (B) and A/J (A) mouse strains. The B x A intercross revealed significant quantitative trait loci (QTL) for HDL on chromosomes 1, 4, 8, 15, 17, 18, and 19. A CSS survey revealed that many have significantly different HDL levels compared to the background strain B, including chromosomes with no significant QTL in the intercross and, in some cases (CSS-1, CSS-17), effects that are opposite to those observed in the B x A intercross population. Intercrosses between B and three CSSs (CSS-3, CSS-11, and CSS-8) revealed significant QTL but with some unexpected differences from the B x A intercross. Our inability to predict the results of CSS intercrosses suggests that additional complexity will be revealed by further crosses and that the CSS mapping strategy should be viewed as a complement to, rather than a replacement for, classical intercross mapping.

Animals↗

Experimental design for three-color and four-color gene expression microarrays.

MOTIVATION: Three-color microarrays, compared with two-color microarrays, can increase design efficiency and power to detect differential expression without additional samples and arrays. Furthermore, three-color microarray technology is currently available at a reasonable cost. Despite the potential advantages, clear guidelines for designing and analyzing three-color experiments do not exist. RESULTS: We propose a three- and a four-color cyclic design (loop) and a complementary graphical representation to help design experiments that are balanced, efficient and robust to hybridization failures. In theory, three-color loop designs are more efficient than two-color loop designs. Experiments using both two- and three-color platforms were performed in parallel and their outputs were analyzed using linear mixed model analysis in R/MAANOVA. These results demonstrate that three-color experiments using the same number of samples (and fewer arrays) will perform as efficiently as two-color experiments. The improved efficiency of the design is somewhat offset by a reduced dynamic range and increased variability in the three-color experimental system. This result suggests that, with minor technological improvements, three-color microarrays using loop designs could detect differential expression more efficiently than two-color loop designs. AVAILABILITY: http://www.jax.org/staff/churchill/labsite/software SUPPLEMENTARY INFORMATION: Multicolor cyclic design construction methods and examples along with additional results of the experiment are provided at http://www.jax.org/staff/churchill/labsite/pubs/yong.

Algorithms↗

Quantitative trait loci associated with elevated thyroid-stimulating hormone in the Wistar-Kyoto rat.

Thyroid hormones are essential for the regulation of developmental and physiological processes. The genetic factors underlying naturally occurring variability in mammalian thyroid function are, however, only partially understood. Genetic control of thyroid function can be studied with animal models such as the inbred Wistar-Kyoto (WKY) rat strain. Previous studies established that WKY rats have elevated TSH, slightly elevated total T3, and normal total T4 levels compared with Wistar controls. The present study confirmed a persistent 24-h elevation of TSH in WKY rats compared with the Fisher 344 (F344) rat, another inbred strain. Acute T3 challenge (25 microg/100 g body weight ip) suppressed serum TSH and T4 levels in both strains. Quantitative trait locus analysis of elevated TSH in a reciprocally bred WKY x F344 F2 population identified one highly significant locus on chromosome 6 (LOD=11.7, TSH-1) and one suggestive locus on chromosome 5 (LOD=2.3, TSH-2). The confidence interval of TSH-1 contains the TSH receptor and type 2 deiodinase genes, and TSH-2 contains the type 1 deiodinase gene. The WKY alleles of each gene contain sequence alterations, but additional studies are indicated to identify the specific gene or genes responsible for altered regulation of the thyroid axis. These findings suggest that one or more genetic alterations within the TSH-1 locus significantly contribute to the altered thyroid function tests of the WKY rat.

3' Untranslated Regions↗

Influence of sex and diet on quantitative trait loci for HDL cholesterol levels in an SM/J by NZB/BlNJ intercross population.

To investigate the dependence of HDL quantitative trait loci (QTL) on sex and diet, we generated a large intercross population of mice from parental strains SM/J and NZB/BlNJ. We measured HDL levels in progeny fed a chow diet and measured them again after 6, 12, and 16 weeks of feeding a high-fat, high-cholesterol diet. QTL analysis was performed on the 260 female and 253 male F(2) progeny. A total of 13 significant QTL were found. Four QTL were specific to female mice: Hdlq23 (Chr 6, 26 cM), Hdlq26 (Chr 10, 70 cM), Hdlq27 (Chr 15, 48 cM), and Hdlq32 (Chr 19, 40 cM). One significant QTL was specific to male mice: Hdlq29 (Chr 17, 36 cM). In addition, several QTL were found to have effects that were dependent on diet. Sex- and diet-dependent effects were characterized using a linear model-based genome scan method that avoids the potential pitfalls of subdivided data analysis. The dependence of QTL effects on sex suggests an important role for the sex hormones in HDL regulation. We recommend that sex should be explicitly accounted for in future studies in the genetics of HDL regulation in both mice and humans.

Animals↗

A comparison of cDNA, oligonucleotide, and Affymetrix GeneChip gene expression microarray platforms.

We have conducted a study to compare the variability in measured gene expression levels associated with three types of microarray platforms. Total RNA samples were obtained from liver tissue of four male mice, two each from inbred strains A/J and C57BL/6J. The same four samples were assayed on Affymetrix Mouse Genome Expression Set 430 GeneChips (MOE430A and MOE430B), spotted cDNA microarrays, and spotted oligonucleotide microarrays using eight arrays of each type. Variances associated with measurement error were observed to be comparable across all microarray platforms. The MOE430A GeneChips and cDNA arrays had higher precision across technical replicates than the MOE430B GeneChips and oligonucleotide arrays. The Affymetrix platform showed the greatest range in the magnitude of expression levels followed by the oligonucleotide arrays. We observed good concordance in both estimated expression level and statistical significance of common genes between the Affymetrix MOE430A GeneChip and the oligonucleotide arrays. Despite their apparently high precision, cDNA arrays showed poor concordance with other platforms.

Analysis of Variance↗

Genetic modulation of tau phosphorylation in the mouse.

The axonal microtubule stabilizing protein tau is hyperphosphorylated in several neurodegenerative conditions, including Alzheimer's disease, yet the genes that regulate tau phosphorylation are largely unknown. Disabled-1 (Dab1) is a cytoplasmic adapter protein that interacts with apolipoprotein E (ApoE) receptors and controls neuronal positioning during embryonic brain development. We have investigated the role of Dab1 in tau phosphorylation. We found that wild-type Dab1, but not a mutant lacking tyrosine phosphorylation sites, protects mice from the hyperphosphorylation of tau. However, the absence of Dab1 is not sufficient to cause tau hyperphosphorylation, because hyperphosphorylation is manifested only when Dab1 is mutated in specific mouse strain backgrounds. Tau hyperphosphorylation correlates with early death in susceptible mouse strains, and it occurs in the neurons of the hippocampus and dentate gyrus. By quantitative trait locus (QTL) analysis of Dab1-deficient mice on a hybrid strain background, we uncovered one significant and three suggestive chromosomal loci that modulate tau phosphorylation. Two of these QTL regions contain genes that are defective in early onset Alzheimer's disease. Our findings suggest that Dab1 gene disruption sensitizes mice to tau hyperphosphorylation contingent on specific haplotypes that are linked to Alzheimer's disease loci. Dab1 mutant mice provide an animal model for studying the relationships between ApoE receptors, tau hyperphosphorylation, and Alzheimer's disease.

Alzheimer Disease↗

Major locus on mouse chromosome 17 and minor locus on chromosome 9 are linked with alopecia areata in C3H/HeJ mice.

Alopecia areata is an autoimmune disease that targets actively growing (anagen) hair follicles in humans, mice, rats, dogs, horses, and cattle. C3H/HeJ mice spontaneously develop alopecia areata from 5 mo of age and older in females and later in males. Frequency of disease approached 20% in a colony by 18 mo of age. C57BL/6J mice do not develop alopecia areata. A segregating F2 population of female mice (n=1096) was generated from crossing these two strains. Alopecia areata (n=138) and clinically normal (n=214) mice were genotyped at 12 mo of age using 211 microsatellite probes. The peak logarithm of odds ratio score on mouse chromosome 17 (10.9) was around marker D17Mit134 at 16.9 cM from the centromere. The mouse histocompatibility locus, H2, the mouse equivalent of human leukocyte antigen in humans, was a likely candidate. Twelve-month-old C3H.SW-H2b/SnJ mice (C3H/HeJ congenic mice in which the H2k purported susceptibility locus was replaced with the H2b purported resistance locus) did not develop alopecia areata, supporting this locus as being important in alopecia areata. A suggestive linkage was also found on mouse Chromosome 9 (logarithm of odds ratio score 2.0) around D9Mit206, 20 cM from the centromere. The interval on mouse Chromosome 17 contains several orthologous genes potentially associated with human alopecia areata.

Alopecia Areata↗

Pleiotropy, homeostasis, and functional networks based on assays of cardiovascular traits in genetically randomized populations.

A major problem in studying biological traits is understanding how genes work together to provide organismal structures and functions. Conventional reductionist paradigms attribute functions to particular proteins, motifs, and amino acids. An equally important but harder problem involves the synthesis of data at fundamental levels of biological systems to understand functionality at higher levels. We used subtle, naturally occurring, multigenic variation of cardiovascular (CV) properties in a panel of genetically randomized strains that are derived from the A/J and C57BL/6J strains of mice to perturb CV functions in nonpathologic ways. In this proof-of-concept study, computational analysis correctly identified the known relations among CV properties and revealed functionality at higher levels of the CV system. The network was then used to account for pleiotropies and homeostatic responses in single gene mutant mice and in mice treated with a pharmacologic agent (anesthesia). The CV network accounted for functional dependencies in complementary ways to the insights obtained from genetic networks and biochemical pathways. These networks are therefore an important approach for defining and characterizing functional relations in complex biological systems in health and disease.

Animals↗

Generation of a new congenic mouse strain to test the relationships among serum insulin-like growth factor I, bone mineral density, and skeletal morphology in vivo.

Insulin-like growth factor (IGF) I is a critical peptide for skeletal growth and consolidation. However, its regulation is complex and, in part, heritable. We previously indicated that changes in both serum and skeletal IGF-I were related to strain-specific differences in total femoral bone mineral density (BMD) in mice. In addition, we defined four quantitative trait loci (QTLs) that contribute to the heritable determinants of the serum IGF-I phenotype in F2 mice derived from progenitor crosses between C3H/HeJ (C3H; high total femoral BMD and high IGF-I) and C57BL/6J (B6; low total femoral BMD and low IGF-I) strains. The strongest QTL, IGF-I serum level 1 (Igflsl-1; log10 of the odds ratio [LOD] score, approximately 9.0), is located on the middle portion of chromosome (Chr) 6. For this locus, C3H alleles are associated with a significant reduction in serum IGF-I. To test the effect of this QTL in vivo, we generated a new congenic strain (B6.C3H-6T [6T]) by placing the Chr 6 QTL region (D6Mit93 to D6Mit150) from C3H onto the B6 background. We then compared serum and skeletal IGF-I levels, body weight, and several skeletal phenotypes from the N9 generation of 6T congenic mice against B6 control mice. Female 6T congenic mice had 11-21% lower serum IGF-I levels at 6, 8, and 16 weeks of age compared with B6 (p < 0.05 for all). In males, serum IGF-I levels were similar in 6T congenics and B6 controls at 6 weeks and 8 weeks but were lower in 6T congenic mice at 16 weeks (p < 0.02). In vitro, there was a 40% reduction in secreted IGF-I in the conditioned media (CMs) from 6T calvaria osteoblasts compared with B6 cells (p < 0.01). Total femoral BMD as measured by peripheral quantitative computed tomography (pQCT) was lower in both 6T male (-4.8%, p < 0.01) and 6T female (-2.3%, p = 0.06) congenic mice. Geometric features of middiaphyseal cortical bone were reduced in 6T congenic mice compared with control mice. Femoral cancellous bone volume (BV) density and trabecular number (Tb.N) were 50% lower, whereas trabecular separation (Tb.Sp) was 90% higher in 8-week-old female 6T congenic mice compared with B6 control mice (p < 0.01 for all). Similarly, vertebral cancellous BV density and Tb.N were lower (-29% and -19%, respectively), whereas Tb.Sp was higher (+29%) in 16-week-old female 6T congenic mice compared with B6 control mice (p < 0.001 for all). Histomorphometric evaluation of the proximal tibia indicated that 6T congenics had reduced BV fraction, labeled surface, and bone formation rates compared with B6 congenic mice. In summary, we have developed a new congenic mouse strain that confirms the Chr 6 QTL as a major genetic regulatory determinant for serum IGF-I. This locus also influences bone density and morphology, with more dramatic effects in cancellous bone than in cortical bone.

Animals↗