A bootstrap-based regression method for comprehensive discovery of differential gene expressions: An application to the osteoporosis study

A common purpose of microarray experiments is to study the variation in gene expression across the categories of an experimental factor such as tissue types and drug treatments. However, it is not uncommon that the studied experimental factor is a quantitative variable rather than categorical variable. Loss of information would occur by comparing gene-expression levels between groups that are factitiously defined according to the quantitative threshold values of an experimental factor. Additionally, lack of control for some sensitive clinical factors may bring serious false positive or negative findings.In the present study, we described a bootstrap-based regression method for analyzing gene-expression data from the non-categorical microarray experiments. To illustrate the utility of this method, we applied it to our recent gene-expression study of circulating monocytes in subjects with a wide range of variations in bone mineral density (BMD). This method allows a comprehensive discovery of gene expressions associated with osteoporosis-related traits while controlling other common confounding factors such as height, weight and age. Several genes identified in our study are involved in osteoblast and osteoclast functions and bone remodeling and/or menopause-associated estrogen-dependent pathways, which provide important clues to understand the etiology of osteoporosis. Availability: SAS code is available from the authors upon request. © 2011 Elsevier Masson SAS.postprin

Copy Number Variation in CNP267 Region May Be Associated with Hip Bone Size

Osteoporotic hip fracture (HF) is a serious global public health problem associated with high morbidity and mortality. Hip bone size (BS) has been identified as one of key measurable risk factors for HF, independent of bone mineral density (BMD). Hip BS is highly genetically determined, but genetic factors underlying BS variation are still poorly defined. Here, we performed an initial genome-wide copy number variation (CNV) association analysis for hip BS in 1,627 Chinese Han subjects using Affymetrix GeneChip Human Mapping SNP 6.0 Array and a follow-up replicate study in 2,286 unrelated US Caucasians sample. We found that a copy number polymorphism (CNP267) located at chromosome 2q12.2 was significantly associated with hip BS in both initial Chinese and replicate Caucasian samples with p values of 4.73E-03 and 5.66E-03, respectively. An important candidate gene, four and a half LIM domains 2 (FHL2), was detected at the downstream of CNP267, which plays important roles in bone metabolism by binding to several bone formation regulator, such as insulin-like growth factor-binding protein 5 (IGFBP-5) and androgen receptor (AR). Our findings suggest that CNP267 region may be associated with hip BS which might influence the FHL2 gene downstream

Gene expression analysis in human osteoblasts exposed to dexamethasone identifies altered developmental pathways as putative drivers of osteoporosis

BACKGROUND: Osteoporosis, a disease of decreased bone mineral density represents a significant and growing burden in the western world. Aging population structure and therapeutic use of glucocorticoids have contributed in no small way to the increase in the incidence of this disease. Despite substantial investigative efforts over the last number of years the exact molecular mechanism underpinning the initiation and progression of osteoporosis remain to be elucidated. This has meant that no significant advances in therapeutic strategies have emerged, with joint replacement surgery being the mainstay of treatment. METHODS: In this study we have used an integrated genomics profiling and computational biology based strategy to identify the key osteoblast genes and gene clusters whose expression is altered in response to dexamethasone exposure. Primary human osteoblasts were exposed to dexamethasone in vitro and microarray based transcriptome profiling completed. RESULTS: These studies identified approximately 500 osteoblast genes whose expression was altered. Functional characterization of the transcriptome identified developmental networks as being reactivated with 106 development associated genes found to be differentially regulated. Pathway reconstruction revealed coordinate alteration of members of the WNT signaling pathway, including frizzled-2, frizzled-7, DKK1 and WNT5B, whose differential expression in this setting was confirmed by real time PCR. CONCLUSION: The WNT pathway is a key regulator of skeletogenesis as well as differentiation of bone cells. Reactivation of this pathway may lead to altered osteoblast activity resulting in decreased bone mineral density, the pathological hallmark of osteoporosis. The data herein lend weight to the hypothesis that alterations in developmental pathways drive the initiation and progression of osteoporosis

Diurnally Entrained Anticipatory Behavior in Archaea

By sensing changes in one or few environmental factors biological systems can anticipate future changes in multiple factors over a wide range of time scales (daily to seasonal). This anticipatory behavior is important to the fitness of diverse species, and in context of the diurnal cycle it is overall typical of eukaryotes and some photoautotrophic bacteria but is yet to be observed in archaea. Here, we report the first observation of light-dark (LD)-entrained diurnal oscillatory transcription in up to 12% of all genes of a halophilic archaeon Halobacterium salinarum NRC-1. Significantly, the diurnally entrained transcription was observed under constant darkness after removal of the LD stimulus (free-running rhythms). The memory of diurnal entrainment was also associated with the synchronization of oxic and anoxic physiologies to the LD cycle. Our results suggest that under nutrient limited conditions halophilic archaea take advantage of the causal influence of sunlight (via temperature) on O2 diffusivity in a closed hypersaline environment to streamline their physiology and operate oxically during nighttime and anoxically during daytime

Genome-wide scans identify known and novel regions associated with prolificacy and reproduction traits in a sub-Saharan African indigenous sheep (Ovis aries)

Maximizing the number of offspring born per female is a key functionality trait in commercial- and/or subsistence-oriented livestock enterprises. Although the number of offspring born is closely associated with female fertility and reproductive success, the genetic control of these traits remains poorly understood in sub-Saharan Africa livestock. Using selection signature analysis performed on Ovine HD BeadChip data from the prolific Bonga sheep in Ethiopia, 41 candidate regions under selection were identified. The analysis revealed one strong selection signature on a candidate region on chromosome X spanning BMP15, suggesting this to be the primary candidate prolificacy gene in the breed. The analysis also identified several candidate regions spanning genes not reported before in prolific sheep but underlying fertility and reproduction in other species. The genes associated with female reproduction traits included SPOCK1 (age at first oestrus), GPR173 (mediator of ovarian cyclicity), HB-EGF (signalling early pregnancy success) and SMARCAL1 and HMGN3a (regulate gene expression during embryogenesis). The genes involved in male reproduction were FOXJ1 (sperm function and successful fertilization) and NME5 (spermatogenesis). We also observed genes such as PKD2L2, MAGED1 and KDM3B, which have been associated with diverse fertility traits in both sexes of other species. The results confirm the complexity of the genetic mechanisms underlying reproduction while suggesting that prolificacy in the Bonga sheep, and possibly African indigenous sheep is partly under the control of BMP15 while other genes that enhance male and female fertility are essential for reproductive fitness

Single nucleotide polymorphisms in bone turnover-related genes in Koreans: ethnic differences in linkage disequilibrium and haplotype

<p>Abstract</p> <p>Background</p> <p>Osteoporosis is defined as the loss of bone mineral density that leads to bone fragility with aging. Population-based case-control studies have identified polymorphisms in many candidate genes that have been associated with bone mass maintenance or osteoporotic fracture. To investigate single nucleotide polymorphisms (SNPs) that are associated with osteoporosis, we examined the genetic variation among Koreans by analyzing 81 genes according to their function in bone formation and resorption during bone remodeling.</p> <p>Methods</p> <p>We resequenced all the exons, splice junctions and promoter regions of candidate osteoporosis genes using 24 unrelated Korean individuals. Using the common SNPs from our study and the HapMap database, a statistical analysis of deviation in heterozygosity depicted.</p> <p>Results</p> <p>We identified 942 variants, including 888 SNPs, 43 insertion/deletion polymorphisms, and 11 microsatellite markers. Of the SNPs, 557 (63%) had been previously identified and 331 (37%) were newly discovered in the Korean population. When compared SNPs in the Korean population with those in HapMap database, 1% (or less) of SNPs in the Japanese and Chinese subpopulations and 20% of those in Caucasian and African subpopulations were significantly differentiated from the Hardy-Weinberg expectations. In addition, an analysis of the genetic diversity showed that there were no significant differences among Korean, Han Chinese and Japanese populations, but African and Caucasian populations were significantly differentiated in selected genes. Nevertheless, in the detailed analysis of genetic properties, the LD and Haplotype block patterns among the five sub-populations were substantially different from one another.</p> <p>Conclusion</p> <p>Through the resequencing of 81 osteoporosis candidate genes, 118 unknown SNPs with a minor allele frequency (MAF) > 0.05 were discovered in the Korean population. In addition, using the common SNPs between our study and HapMap, an analysis of genetic diversity and deviation in heterozygosity was performed and the polymorphisms of the above genes among the five populations were substantially differentiated from one another. Further studies of osteoporosis could utilize the polymorphisms identified in our data since they may have important implications for the selection of highly informative SNPs for future association studies.</p

The Cyst-Dividing Bacterium Ramlibacter tataouinensis TTB310 Genome Reveals a Well-Stocked Toolbox for Adaptation to a Desert Environment

Ramlibacter tataouinensis TTB310T (strain TTB310), a betaproteobacterium isolated from a semi-arid region of South Tunisia (Tataouine), is characterized by the presence of both spherical and rod-shaped cells in pure culture. Cell division of strain TTB310 occurs by the binary fission of spherical “cyst-like” cells (“cyst-cyst” division). The rod-shaped cells formed at the periphery of a colony (consisting mainly of cysts) are highly motile and colonize a new environment, where they form a new colony by reversion to cyst-like cells. This unique cell cycle of strain TTB310, with desiccation tolerant cyst-like cells capable of division and desiccation sensitive motile rods capable of dissemination, appears to be a novel adaptation for life in a hot and dry desert environment. In order to gain insights into strain TTB310's underlying genetic repertoire and possible mechanisms responsible for its unusual lifestyle, the genome of strain TTB310 was completely sequenced and subsequently annotated. The complete genome consists of a single circular chromosome of 4,070,194 bp with an average G+C content of 70.0%, the highest among the Betaproteobacteria sequenced to date, with total of 3,899 predicted coding sequences covering 92% of the genome. We found that strain TTB310 has developed a highly complex network of two-component systems, which may utilize responses to light and perhaps a rudimentary circadian hourglass to anticipate water availability at the dew time in the middle/end of the desert winter nights and thus direct the growth window to cyclic water availability times. Other interesting features of the strain TTB310 genome that appear to be important for desiccation tolerance, including intermediary metabolism compounds such as trehalose or polyhydroxyalkanoate, and signal transduction pathways, are presented and discussed

TNFRSF11A and TNFSF11 are associated with age at menarche and natural menopause in white women

Objective: Menarche and menopause mark the lower and upper limits of the female reproductive period. The timing of these events influences women's health in later life. The onsets of menarche and menopause have a strong genetic basis. We tested two genes, TNFRSF11A (RANK) and TNFSF11 (RANKL), for their association with age at menarche (AM) and age at natural menopause (ANM). Methods: Nineteen single nucleotide polymorphisms (SNPs) of TNFRSF11A and 12 SNPs of TNFSF11 were genotyped in a random sample of 306 unrelated white women. This sample was analyzed for the association of the SNPs and common haplotypes with AM. Then, a subsample of 211 women with natural menopause was analyzed for the association of both genes with ANM. Smoking, alcohol intake, and duration of lactation were applied as covariates in the association analyses. Results: Three polymorphisms of TNFSF11 were associated with AM: rs2200287 (P = 0.005), rs9525641 (P = 0.039), and rs1054016 (P = 0.047). Two SNPs of this gene, rs346578 and rs9525641, showed an association with ANM (P = 0.007 and P = 0.011, respectively). Two SNPs of TNFRSF11A were associated with AM (rs3826620; P = 0.022) and ANM (rs8086340; P = 0.015). Multiple SNP-SNP and SNP-environment interaction effects on AM and ANM were detected for both genes. One polymorphism of TNFRSF11A, rs4436867, was not directly associated with either trait but indicated significant interactions with four TNFSF11 polymorphisms on ANM. Two other TNFRSF11A polymorphisms, rs4941125 and rs7235803, showed interaction effects with several TNFSF11 polymorphisms on AM. Both genes manifested significant interaction with the duration of breast-feeding in their effect on ANM. Conclusions: The TNFRSF11A and TNFSF11 genes are associated with the onset of AM and ANM in white women. © 2010 by The North American Menopause Society.link_to_subscribed_fulltex

Association analyses suggest multiple interaction effects of the methylenetetrahydrofolate reductase polymorphisms on timing of menarche and natural menopause in white women

Objective: This study aimed to investigate whether polymorphisms of the methylenetetrahydrofolate reductase (MTHFR) gene are associated with age at menarche and age at natural menopause in white women. Methods: In a cross-sectional study, a total of 305 randomly selected unrelated white women were genotyped for six single nucleotide polymorphisms (SNPs) of the MTHFR gene (including one common replacement, rs1801133). This sample was comprehensively analyzed for the association of the SNPs with age at menarche. Then a subsample of 210 women who experienced natural menopause was analyzed for the association of the MTHFR gene with age at natural menopause. Results: Duration of breast-feeding was a significant predictor of earlier natural menopause (P < 0.05). No individual SNPs were associated with either age at menarche or age at natural menopause. However, three significant (P < 0.05) SNP-SNP interaction effects (rs2066470/rs1476413, rs2066470/rs4846049, and rs17037390/rs4846049) on the onset of menarche were determined. Three haplotypes were significantly associated with age at menopause (P < 0.05). Four SNPs (rs2066470, rs17037390, rs1801133, and rs4846048) indicated significant interaction effects with various lifestyle factors on age at natural menopause. Conclusions: The results of our study suggest that the MTHFR gene may influence the onset of menarche and natural menopause. This effect is probably due to the multiple SNP-SNP and SNP-environment interactions. More independent studies are needed to further clarify the possible contribution of this gene to the timing of menarche and menopause. © 2010 by The North American Menopause Society.link_to_OA_fulltex