Lateral packing of mineral crystals in bone collagen fibrils

Combined small-angle x-ray scattering and transmission electron microscopy studies of intramuscular fish bone (shad and herring) indicate that the lateral packing of nanoscale calcium-phosphate crystals in collagen fibrils can be represented by irregular stacks of platelet-shaped crystals, intercalated with organic layers of collagen molecules. The scattering intensity distribution in this system can be described by a modified Zernike-Prins model, taking preferred orientation effects into account. Using the model, the diffuse fan-shaped small-angle x-ray scattering intensity profile, dominating the equatorial region of the scattering pattern, could be quantitatively analyzed as a function of the degree of mineralization. The mineral platelets were found to be very thin (1.5 nm ∼ 2.0 nm), having a narrow thickness distribution. The thickness of the organic layers between adjacent mineral platelets within a stack is more broadly distributed with the average value varying from 6 nm to 10 nm, depending on the extent of mineralization. The two-dimensional analytical scheme also leads to quantitative information about the preferred orientation of mineral stacks and the average height of crystals along the crystallographic c axis.     

Chromosomal regions 22q13 and 3p25 may harbor quantitative trait loci influencing both age at menarche and bone mineral density

Late age at menarche (AAM), an important type of endocrinopathy in females, is associated with lower bone mineral density (BMD), a major risk factor for osteoporosis. The correlation is mainly mediated through common genetic factors, which are largely unknown. A bivariate genome-wide linkage scan was conducted on 2,522 females from 414 Caucasian pedigrees to identify quantitative trait loci influencing both AAM and BMD. The strongest linkage signal was detected on chromosome 22q13. Other regions such as the 3q13, 3p25, 7p15, and 15q13 were also suggested. The inferred promising candidate genes in the linkage regions may contribute to our understanding of pathogenesis of endocrinopathy and osteoporosis in females.     

中国板栗居群间等位酶基因频率的空间分布

中国板栗21个自然群间等位酶遗传变异的空间自相关分析及F－统计分析结果表明：其多数等位基因频率在居群间呈随机分布模式,缺乏一定的空间结构;而部分等位基因表现为渐变或双向渐变的非随机分布模式,又具有特定空间结构。中国板栗遗传变异空间结构模式的形成可能是长距离基因流、自然气候、人类活动、地理距离隔离等诸因素综合作用的结果。文中还就居群等位基因分布格局的成因进行了讨论：在第四纪冰川后,中国板栗以长江流域中下游的孑遗中心为起点,等位基因分别沿着向北和向南的不同方向迁移形成现在的居群结构;季风气候和人类活动干扰是削弱居群分化的主要因素,而基于环境梯度的选择,是形成由北向南渐变分布的原因。     

Locomotion Guidance by Extracellular Matrix Is Adaptive and Can be Restored by a Transient Change in Ca2+ Level

Navigation of cell locomotion by gradients of soluble factors can be desensitized if the concentration of the chemo-attractant stays unchanged. It remains obscure if the guidance by immobilized extracellular matrix (ECM) as the substrate is also adaptive and if so, how can the desensitized ECM guidance be resensitized. When first interacting with a substrate containing micron-scale fibronectin (FBN) trails, highly motile fish keratocytes selectively adhere and migrate along the FBN paths. However, such guided motion become adaptive after about 10 min and the cells start to migrate out of the ECM trails. We found that a burst increase of intracellular calcium created by an uncaging technique immediately halts the undirected migration by disrupting the ECM-cytoskeleton coupling, as evidenced by the appearance of retrograde F-actin flow. When the motility later resumes, the activated integrin receptors render the cell selectively binding to the FBN path and reinitiates signaling events, including tyrosine phosphorylation of paxillin, that couple retrograde F-actin flow to the substrate. Thus, the calcium-resensitized cell can undergo a period of ECM-navigated movement, which later becomes desensitized. Our results also suggest that endogenous calcium transients as occur during spontaneous calcium oscillations may exert a cycling resensitization-desensitization control over cell''s sensing of substrate guiding cues.     

A Bayesian Analysis for Identifying DNA Copy Number Variations Using a Compound Poisson Process

To study chromosomal aberrations that may lead to cancer formation or genetic diseases, the array-based Comparative Genomic Hybridization (aCGH) technique is often used for detecting DNA copy number variants (CNVs). Various methods have been developed for gaining CNVs information based on aCGH data. However, most of these methods make use of the log-intensity ratios in aCGH data without taking advantage of other information such as the DNA probe (e.g., biomarker) positions/distances contained in the data. Motivated by the specific features of aCGH data, we developed a novel method that takes into account the estimation of a change point or locus of the CNV in aCGH data with its associated biomarker position on the chromosome using a compound Poisson process. We used a Bayesian approach to derive the posterior probability for the estimation of the CNV locus. To detect loci of multiple CNVs in the data, a sliding window process combined with our derived Bayesian posterior probability was proposed. To evaluate the performance of the method in the estimation of the CNV locus, we first performed simulation studies. Finally, we applied our approach to real data from aCGH experiments, demonstrating its applicability.     

影响中国女性峰值骨密度的潜在基因间相互作用

峰值骨密度是由遗传和环境因素及其相互作用共同决定的复杂性状。维生素D受体基因、雌激素α受体基因、白介素6基因、副甲状腺素基因、Ⅰ型胶原α2基因、骨钙素基因、α2巯基糖蛋白基因是与骨代谢相关的7个重要的候选基因。本研究旨在检测这7个候选基因之间的相互作用对中国女性峰值骨密度的影响。样本为中国上海的361个无关、健康的绝经前女性,均为汉族人,年龄为20—44岁。采用Hologic QDR2000＋双能X射线扫描仪对腰椎与髋部的骨密度进行了检测。采用聚合酶链式反应-限制性片段长度多态分析方法对每个个体的以下8个多态性标记位点进行基因分型：维生素D受体基因的Apa Ⅰ位点,雌激素α受体基因的PvuⅡ和XbaⅠ位点,白介素6基因的BsrB Ⅰ位点,副甲状腺素基因的BstB Ⅰ位点,Ⅰ型胶原α2基因的Msp Ⅰ位点,骨钙素基因的Hind Ⅲ位点,娃巯基糖蛋白基因的SacⅠ位点。采用二元方差分析对基因相互作用与骨密度的关系进行研究。结果表明,白介素基因和雌激素α受体基因（PvuⅡ）的相互作用对髋部（P=0．019）、转子间区（P=0．016）和股骨颈（P=0．019）的骨密度有显著作用。在这3个部位,GGPp基因型携带者比GGpp基因型携带者的骨密度值分别高出18．0％、19．5％和14．8％。另外观察到醒巯基糖蛋白基因与自介素6基因的相互作用对股骨颈骨密度有显著影响（P=0．046）。GGSS基因型携带者的股骨颈骨密度值比GGSs基因型携带者高出18．8％。该项群体水平的统计分析表明：对于中国女性峰值骨密度的遗传决定,白介素基因和雌激素α受体基因、α2巯基糖蛋白基因的相互作用显著。     

中国栗疫病菌群体遗传结构的空间自相关性分析

应用空间自相关分析方法对中国栗疫病菌17个居群RAPD遗传变异的空间结构进行研究,以探讨栗疫病菌居群遗传变异的分布特征及其形成机制。结果表明：中国栗疫病菌居群缺乏空间结构,绝大多数RAPD位点变异为随机分布的空间模式,但部分位点表现出渐变、斑块和双向渐变的非随机分布模式,又显示了一定的空间结构。推测其形成原因可能是长距离的基因流、人类活动、地理隔离以及栗疫病菌本身的繁殖特性综合作用的结果,并依据部分位点呈单向渐变的模式推测西南地区为中国栗疫病菌的起源中心。     

Univariate/Multivariate Genome-Wide Association Scans Using Data from Families and Unrelated Samples

As genome-wide association studies (GWAS) are becoming more popular, two approaches, among others, could be considered in order to improve statistical power for identifying genes contributing subtle to moderate effects to human diseases. The first approach is to increase sample size, which could be achieved by combining both unrelated and familial subjects together. The second approach is to jointly analyze multiple correlated traits. In this study, by extending generalized estimating equations (GEEs), we propose a simple approach for performing univariate or multivariate association tests for the combined data of unrelated subjects and nuclear families. In particular, we correct for population stratification by integrating principal component analysis and transmission disequilibrium test strategies. The proposed method allows for multiple siblings as well as missing parental information. Simulation studies show that the proposed test has improved power compared to two popular methods, EIGENSTRAT and FBAT, by analyzing the combined data, while correcting for population stratification. In addition, joint analysis of bivariate traits has improved power over univariate analysis when pleiotropic effects are present. Application to the Genetic Analysis Workshop 16 (GAW16) data sets attests to the feasibility and applicability of the proposed method.     

Critical assessment of coalescent simulators in modeling recombination hotspots in genomic sequences

Background

Coalescent simulation is pivotal for understanding population evolutionary models and demographic histories, as well as for developing novel analytical methods for genetic association studies for DNA sequence data. A plethora of coalescent simulators are developed, but selecting the most appropriate program remains challenging.

Results

We extensively compared performances of five widely used coalescent simulators – Hudson’s ms, msHOT, MaCS, Simcoal2, and fastsimcoal, to provide a practical guide considering three crucial factors, 1) speed, 2) scalability and 3) recombination hotspot position and intensity accuracy. Although ms represents a popular standard coalescent simulator, it lacks the ability to simulate sequences with recombination hotspots. An extended program msHOT has compensated for the deficiency of ms by incorporating recombination hotspots and gene conversion events at arbitrarily chosen locations and intensities, but remains limited in simulating long stretches of DNA sequences. Simcoal2, based on a discrete generation-by-generation approach, could simulate more complex demographic scenarios, but runs comparatively slow. MaCS and fastsimcoal, both built on fast, modified sequential Markov coalescent algorithms to approximate standard coalescent, are much more efficient whilst keeping salient features of msHOT and Simcoal2, respectively. Our simulations demonstrate that they are more advantageous over other programs for a spectrum of evolutionary models. To validate recombination hotspots, LDhat 2.2 rhomap package, sequenceLDhot and Haploview were compared for hotspot detection, and sequenceLDhot exhibited the best performance based on both real and simulated data.

Conclusions

While ms remains an excellent choice for general coalescent simulations of DNA sequences, MaCS and fastsimcoal are much more scalable and flexible in simulating a variety of demographic events under different recombination hotspot models. Furthermore, sequenceLDhot appears to give the most optimal performance in detecting and validating cross-over hotspots.