Quantitative genetic modeling and inference in the presence of nonignorable missing data.

Natural selection is typically exerted at some specific life stages. If natural selection takes place before a trait can be measured, using conventional models can cause wrong inference about population parameters. When the missing data process relates to the trait of interest, a valid inference requires explicit modeling of the missing process. We propose a joint modeling approach, a shared parameter model, to account for nonrandom missing data. It consists of an animal model for the phenotypic data and a logistic model for the missing process, linked by the additive genetic effects. A Bayesian approach is taken and inference is made using integrated nested Laplace approximations. From a simulation study we find that wrongly assuming that missing data are missing at random can result in severely biased estimates of additive genetic variance. Using real data from a wild population of Swiss barn owls Tyto alba, our model indicates that the missing individuals would display large black spots; and we conclude that genes affecting this trait are already under selection before it is expressed. Our model is a tool to correctly estimate the magnitude of both natural selection and additive genetic variance.

Arterial properties in relation to genetic variations in the adducin subunits in a white population.

BACKGROUND: Adducin is a membrane skeleton protein, which consists of either alpha- and beta- or alpha- and gamma-subunits. We investigated whether arterial characteristics might be related to the genes encoding ADD1 (Gly460Trp-rs4961), ADD2 (C1797T-rs4984), and ADD3 (IVS11+386A>G-rs3731566). METHODS: We randomly recruited 1,126 Flemish subjects (mean age, 43.8 years; 50.3% women). Using a wall-tracking ultrasound system, we measured the properties of the carotid, femoral, and brachial arteries. We studied multivariate-adjusted phenotype-genotype associations, using a population- and family-based approach. RESULTS: In single-gene analyses, brachial diameter was 0.15 mm (P = 0.0022) larger, and brachial distensibility and cross-sectional compliance were 1.55 x 10(-3)/kPa (P = 0.013) and 0.017 mm(2)/kPa (P = 0.0029) lower in ADD3 AA than ADD3 GG homozygotes with an additive effect of the G allele. In multiple-gene analyses, the association of brachial diameter and distensibility with the ADD3 G allele occurred only in ADD1 GlyGly homozygotes. Otherwise, the associations between the arterial phenotypes in the three vascular beds and the ADD1 or ADD2 polymorphisms were not significant. In family-based analyses, the multivariate-adjusted heritability was 0.52, 0.38, and 0.30 for brachial diameter, distensibility, and cross-sectional compliance, respectively (P < 0.001). There was no evidence for population stratification (0.07 < or = P < or = 0.96). Transmission of the mutated ADD3 G allele was associated with smaller brachial diameter in 342 informative offspring (-0.12 +/- 0.04 mm; P = 0.0085) and in 209 offspring, who were ADD1 GlyGly homozygotes (-0.14 +/- 0.06 mm; P = 0.018). CONCLUSIONS: In ADD1 GlyGly homozygotes, the properties of the brachial artery are related to the ADD3 (A386G) polymorphism, but the underlying mechanism needs further clarification.

Fine-scale spatial genetic structure and gene dispersal in Silene latifolia.

Plants are sessile organisms, often characterized by limited dispersal. Seeds and pollen are the critical stages for gene flow. Here we investigate spatial genetic structure, gene dispersal and the relative contribution of pollen vs seed in the movement of genes in a stable metapopulation of the white campion Silene latifolia within its native range. This short-lived perennial plant is dioecious, has gravity-dispersed seeds and moth-mediated pollination. Direct measures of pollen dispersal suggested that large populations receive more pollen than small isolated populations and that most gene flow occurs within tens of meters. However, these studies were performed in the newly colonized range (North America) where the specialist pollinator is absent. In the native range (Europe), gene dispersal could fall on a different spatial scale. We genotyped 258 individuals from large and small (15) subpopulations along a 60 km, elongated metapopulation in Europe using six highly variable microsatellite markers, two X-linked and four autosomal. We found substantial genetic differentiation among subpopulations (global F(ST)=0.11) and a general pattern of isolation by distance over the whole sampled area. Spatial autocorrelation revealed high relatedness among neighboring individuals over hundreds of meters. Estimates of gene dispersal revealed gene flow at the scale of tens of meters (5-30 m), similar to the newly colonized range. Contrary to expectations, estimates of dispersal based on X and autosomal markers showed very similar ranges, suggesting similar levels of pollen and seed dispersal. This may be explained by stochastic events of extensive seed dispersal in this area and limited pollen dispersal.

Isolation of genetic mutation leading to abnormal phenotype in human through ultra-high-throughput sequencing (UHTS)

The introduction of Next Generation Sequencing (NGS) facilitated the task of localizing DNA variation and identifying the genetic cause of yet unsolved Mendelian disorders. Using Whole Exome Capture method and NGS, we identified the causative genetic aberration responsible for a number of monogenic disorders previously undetermined. Due to the novelty of the NGS method we benchmarked different algorithms to assess their merits and defects. This allowed us to establish a pipeline that we successfully used to pinpoint genes responsible for a form of West's syndrome, a Complex Intellectual Disability syndrome associated with patellar dislocation and celiac disease, and correcting some erroneous molecular diagnosis of Alport's syndrome in a Saudi Arabian family.

Trabecular Bone Score: A Noninvasive Analytical Method Based Upon the DXA Image.

The trabecular bone score (TBS) is a gray-level textural metric that can be extracted from the two-dimensional lumbar spine dual-energy X-ray absorptiometry (DXA) image. TBS is related to bone microarchitecture and provides skeletal information that is not captured from the standard bone mineral density (BMD) measurement. Based on experimental variograms of the projected DXA image, TBS has the potential to discern differences between DXA scans that show similar BMD measurements. An elevated TBS value correlates with better skeletal microstructure; a low TBS value correlates with weaker skeletal microstructure. Lumbar spine TBS has been evaluated in cross-sectional and longitudinal studies. The following conclusions are based upon publications reviewed in this article: 1) TBS gives lower values in postmenopausal women and in men with previous fragility fractures than their nonfractured counterparts; 2) TBS is complementary to data available by lumbar spine DXA measurements; 3) TBS results are lower in women who have sustained a fragility fracture but in whom DXA does not indicate osteoporosis or even osteopenia; 4) TBS predicts fracture risk as well as lumbar spine BMD measurements in postmenopausal women; 5) efficacious therapies for osteoporosis differ in the extent to which they influence the TBS; 6) TBS is associated with fracture risk in individuals with conditions related to reduced bone mass or bone quality. Based on these data, lumbar spine TBS holds promise as an emerging technology that could well become a valuable clinical tool in the diagnosis of osteoporosis and in fracture risk assessment. © 2014 American Society for Bone and Mineral Research.

Genetic structure of Gymnures (genus Hylomys; Erinaceidae) on continental islands of Southeast Asia: Historical effects of fragmentation

Eustatic sea level changes during Pleistocene climatic fluctuations produced several cycles of connection-isolation among continental islands of the Sunda shelf. To explore the potential effects of these fluctuations, we reconstructed a model of the vicariant events that separated these islands, based on bathymetric information. Among many possible scenarios, two opposite phylogenetic patterns of evolution were predicted for terrestrial organisms living in this region: one is based on the classical allopatric speciation mode of evolution, while the other is the outcome of a sequential dispersal colonization of the archipelago. We tested the applicability of these predictions with an analysis of sequence variation of the cytochrome b gene from several taxa of Hylomys. They were sampled throughout SE-Asia and the Sunda islands. High levels of haplotype differentiation characterize the different island taxa. Such levels of differentiation support the existence of several allopatric species, as was suggested by previous allozyme and morphological data. Also in accordance with previous results, the occurrence of two sympatric species from Sumatra is suggested by their strongly divergent haplotypes. One species, Hylomys suillus maxi, is found both on Sumatra and in Peninsular Malaysia, while the other, H. parvus, is endemic to Sumatra. Its closest relative is H. suillus dorsalis from Borneo. Phylogenetic reconstructions also demonstrate the existence of a Sundaic clade composed of all island taxa, as opposed to those from the continent. Although there is no statistical support for either proposed biogeographic model of evolution, we argue that the sequential dispersal scenario is more appropriate to describe the genetic variation found among the Hylomys taxa. However, despite strong differentiation among island haplotypes, the cladistic relationships between some island taxa could not be resolved. We argue that this is evidence of a rapid radiation, suggesting that the separation of the islands may have been perceived as a simultaneous event rather than as a succession of vicariant events. Furthermore, the estimates of divergence times between the haplotypes of these taxa suggest that this radiation may actually have predated the climatic fluctuations of the Pleistocene. Further refinement of the initial palaeogeographic models of evolution are therefore needed to account for these results.

Corneal Opacities in Trisomy 8 Mosaic Syndrome : Clinical, Histologic and Genetic Features

Purpose: To describe the clinical, histologic and genetic findings of corneal opacities in the trisomy 8 mosaic syndrome. Methods: 3 children aged 8 years (Patients A), 6 years (Patients B) and 1 month (Patients C) respectively, were referred with corneal opacities for ophthalmologic evaluation. The 2 older patients had been previously diagnosed with trisomy 8 mosaicism, while the third was diagnosed after the ocular examination. Automated lamellar keratoplasty (ALTK) was performed on the most amblyopic eye. Histopathologic analysis with immunohistochemical markers and cytogenetic studies by FISH using haploid probes for chromosome 8 and chromosome 16 (control) were performed on the excised corneal lesion. Results: All patients presented vascularized corneal opacities involving the superficial stroma, and amblyopia with a bilateral involvement in two of them (Patients A and B). Post-operative follow-up (range 6-20 months) was satisfactory, with the graft remaining clear and improved visual acuity, allowing iso-acuity and stereoscopy in the one month old child (Patients C). The clinically observed corneal opacities corresponded histopathologically to the replacement of the normal anterior corneal stroma by a choristomatous loose richly vascularized connective tissue containing mucopolysacharides. Bowman's membrane was absent. There were no adnexal structures. The overlaying epithelium expressed keratin 3 in all three cases. Keratin 19 was found in the suprabasal epithelial cells in one case but was absent in the other cases. There were no expression of keratin 7 and 1 as well as MUC5AC in the epithelial cells. FISH analysis from 100 interphase cells of the affected tissue and normal conjontival probe revealed normal diploid cells. Conclusions: In this series, the corneal opacities associated with trisomy 8 mosaic syndrome share a common clinical, histopathological and genetic features. ALTK should be considered at diagnosis to prevent amblyopia in these children.

The prediction of cardiovascular events is improved by combination of Framingham risk score and subclinical atherosclerosis imaging

Purpose: In primary prevention of cardiovascular disease (CVD), it is accepted that the intensity of risk factor treatment should be guided by the magnitude of absolute risk. Risk factors tools like Framingham risk score (FHS) or noninvasive atherosclerosis imaging tests are available to detect high risk subjects. However, these methods are imperfect and may misclassify a large number of individuals. The purpose of this prospective study was to evaluate whether the prediction of future cardiovascular events (CVE) can be improved when subclinical imaging atherosclerosis (SCATS) is combined with the FRS in asymptomatic subjects. Methods: Overall, 1038 asymptomatic subjects (413 women, 625 men, mean age 49.1±12.8 years) were assessed for their cardiovascular risk using the FRS. B-mode ultrasonography on carotid and femoral arteries was performed by two investigators to detect atherosclerotic plaques (focal thickening of intima-media > 1.2 mm) and to measure carotid intima-media thickness (C-IMT). The severity of SCATS was expressed by an ATS-burden Score (ABS) reflecting the number of the arterial sites with >1 plaques (range 0-4). CVE were defined as fatal or non fatal acute coronary syndrome, stroke, or angioplasty for peripheral artery disease. Results: during a mean follow-up of 4.9±3.1 years, 61 CVE were recorded. Event rates the rate of CVE increased significantly from 2.7% to 39.1% according to the ABS (p<0.001) and from 4% to 24.6% according to the quartiles of C-IMT. Similarly, FRS predicted CVE (p<0.001). When computing the angiographic markers of SCATS in addition of FRS, we observed an improvement of net reclassification rate of 16.6% (p< 0.04) for ABS as compared to 5.5% (p = 0.26) for C-IMT. Conclusion: these results indicate that the detection of subjects requiring more attention to prevent CVE can be significantly improved when using both FRS and SCATS imaging.

Nosology and classification of genetic skeletal disorders: 2010 revision.

Genetic disorders involving the skeletal system arise through disturbances in the complex processes of skeletal development, growth and homeostasis and remain a diagnostic challenge because of their variety. The Nosology and Classification of Genetic Skeletal Disorders provides an overview of recognized diagnostic entities and groups them by clinical and radiographic features and molecular pathogenesis. The aim is to provide the Genetics, Pediatrics and Radiology community with a list of recognized genetic skeletal disorders that can be of help in the diagnosis of individual cases, in the delineation of novel disorders, and in building bridges between clinicians and scientists interested in skeletal biology. In the 2010 revision, 456 conditions were included and placed in 40 groups defined by molecular, biochemical, and/or radiographic criteria. Of these conditions, 316 were associated with mutations in one or more of 226 different genes, ranging from common, recurrent mutations to "private" found in single families or individuals. Thus, the Nosology is a hybrid between a list of clinically defined disorders, waiting for molecular clarification, and an annotated database documenting the phenotypic spectrum produced by mutations in a given gene. The Nosology should be useful for the diagnosis of patients with genetic skeletal diseases, particularly in view of the information flood expected with the novel sequencing technologies; in the delineation of clinical entities and novel disorders, by providing an overview of established nosologic entities; and for scientists looking for the clinical correlates of genes, proteins and pathways involved in skeletal biology. © 2011 Wiley-Liss, Inc.

Analysis of angiotensin II- and ACTH-driven mineralocorticoid functions and omental adiposity in a non-genetic, hyperadipose female rat phenotype

The hypothalamic damage induced by neonatal treatment with monosodium l-glutamate (MSG) induces several metabolic abnormalities, resulting in a rat hyperleptinemic-hyperadipose phenotype. This study was conducted to explore the impact of the neonatal MSG treatment, in the adult (120 days old) female rat on: (a) the in vivo and in vitro mineralocorticoid responses to ACTH and angiotensin II (AII); (b) the effect of leptin on ACTH- and AII-stimulated mineralocorticoid secretions by isolated corticoadrenal cells; and (c) abdominal adiposity characteristics. Our data indicate that, compared with age-matched controls, MSG rats displayed: (1) enhanced and reduced mineralocorticoid responses to ACTH and AII treatments, respectively, effects observed in both in vivo and in vitro conditions; (2) adrenal refractoriness to the inhibitory effect of exogenous leptin on ACTH-stimulated aldosterone output by isolated adrenocortical cells; and (3) distorted omental adiposity morphology and function. This study supports that the adult hyperleptinemic MSG female rat is characterized by enhanced ACTH-driven mineralocorticoid function, impaired adrenal leptin sensitivity, and disrupted abdominal adiposity function. MSG rats could counteract undesirable effects of glucocorticoid excess, by developing a reduced AII-driven mineralocorticoid function. Thus, chronic hyperleptinemia could play a protective role against ACTH-mediated allostatic loads in the adrenal leptin resistant, MSG female rat phenotype.

The shared allelic architecture of adiponectin levels and coronary artery disease.

OBJECTIVE: A large body of epidemiologic data strongly suggests an association between excess adiposity and coronary artery disease (CAD). Low adiponectin levels, a hormone secreted only from adipocytes, have been associated with an increased risk of CAD in observational studies. However, these associations cannot clarify whether this relationship is causal or due to a shared set of causal factors or even confounding. Genome-wide association studies have identified common variants that influence adiponectin levels, providing valuable tools to examine the genetic relationship between adiponectin and CAD. METHODS: Using 145 genome wide significant SNPs for adiponectin from the ADIPOGen consortium (n = 49,891), we tested whether adiponectin-decreasing alleles influenced risk of CAD in the CARDIoGRAM consortium (n = 85,274). RESULTS: In single-SNP analysis, 5 variants among 145 SNPs were associated with increased risk of CAD after correcting for multiple testing (P < 4.4 × 10(-4)). Using a multi-SNP genotypic risk score to test whether adiponectin levels and CAD have a shared genetic etiology, we found that adiponectin-decreasing alleles increased risk of CAD (P = 5.4 × 10(-7)). CONCLUSION: These findings demonstrate that adiponectin levels and CAD have a shared allelic architecture and provide rationale to undertake a Mendelian randomization studies to understand if this relationship is causal.

The additive genetic variance after bottlenecks is affected by the number of loci involved in epistatic interactions.

We investigated the role of the number of loci coding for a neutral trait on the release of additive variance for this trait after population bottlenecks. Different bottleneck sizes and durations were tested for various matrices of genotypic values, with initial conditions covering the allele frequency space. We used three different types of matrices. First, we extended Cheverud and Routman's model by defining matrices of "pure" epistasis for three and four independent loci; second, we used genotypic values drawn randomly from uniform, normal, and exponential distributions; and third we used two models of simple metabolic pathways leading to physiological epistasis. For all these matrices of genotypic values except the dominant metabolic pathway, we find that, as the number of loci increases from two to three and four, an increase in the release of additive variance is occurring. The amount of additive variance released for a given set of genotypic values is a function of the inbreeding coefficient, independently of the size and duration of the bottleneck. The level of inbreeding necessary to achieve maximum release in additive variance increases with the number of loci. We find that additive-by-additive epistasis is the type of epistasis most easily converted into additive variance. For a wide range of models, our results show that epistasis, rather than dominance, plays a significant role in the increase of additive variance following bottlenecks.

Genetic variation near IRS1 associates with reduced adiposity and an impaired metabolic profile.

Genome-wide association studies have identified 32 loci influencing body mass index, but this measure does not distinguish lean from fat mass. To identify adiposity loci, we meta-analyzed associations between ∼2.5 million SNPs and body fat percentage from 36,626 individuals and followed up the 14 most significant (P < 10(-6)) independent loci in 39,576 individuals. We confirmed a previously established adiposity locus in FTO (P = 3 × 10(-26)) and identified two new loci associated with body fat percentage, one near IRS1 (P = 4 × 10(-11)) and one near SPRY2 (P = 3 × 10(-8)). Both loci contain genes with potential links to adipocyte physiology. Notably, the body-fat-decreasing allele near IRS1 is associated with decreased IRS1 expression and with an impaired metabolic profile, including an increased visceral to subcutaneous fat ratio, insulin resistance, dyslipidemia, risk of diabetes and coronary artery disease and decreased adiponectin levels. Our findings provide new insights into adiposity and insulin resistance.

Beclin 1 mitigates motor and neuropathological deficits in genetic mouse models of Machado-Joseph disease.

Machado-Joseph disease or spinocerebellar ataxia type 3, the most common dominantly-inherited spinocerebellar ataxia, results from translation of the polyglutamine-expanded and aggregation prone ataxin 3 protein. Clinical manifestations include cerebellar ataxia and pyramidal signs and there is no therapy to delay disease progression. Beclin 1, an autophagy-related protein and essential gene for cell survival, is decreased in several neurodegenerative disorders. This study aimed at evaluating if lentiviral-mediated beclin 1 overexpression would rescue motor and neuropathological impairments when administered to pre- and post-symptomatic lentiviral-based and transgenic mouse models of Machado-Joseph disease. Beclin 1-mediated significant improvements in motor coordination, balance and gait with beclin 1-treated mice equilibrating longer periods in the Rotarod and presenting longer and narrower footprints. Furthermore, in agreement with the improvements observed in motor function beclin 1 overexpression prevented neuronal dysfunction and neurodegeneration, decreasing formation of polyglutamine-expanded aggregates, preserving Purkinje cell arborization and immunoreactivity for neuronal markers. These data show that overexpression of beclin 1 in the mouse cerebellum is able to rescue and hinder the progression of motor deficits when administered to pre- and post-symptomatic stages of the disease.

Trabecular bone score improves fracture risk prediction in non-osteoporotic women: the OFELY study.

The use of areal bone mineral density (aBMD) for fracture prediction may be enhanced by considering bone microarchitectural deterioration. Trabecular bone score (TBS) helped in redefining a significant subset of non-osteoporotic women as a higher risk group. INTRODUCTION: TBS is an index of bone microarchitecture. Our goal was to assess the ability of TBS to predict incident fracture. METHODS: TBS was assessed in 560 postmenopausal women from the Os des Femmes de Lyon cohort, who had a lumbar spine (LS) DXA scan (QDR 4500A, Hologic) between years 2000 and 2001. During a mean follow-up of 7.8 ± 1.3 years, 94 women sustained 112 fragility fractures. RESULTS: At the time of baseline DXA scan, women with incident fracture were significantly older (70 ± 9 vs. 65 ± 8 years) and had a lower LS_aBMD and LS_TBS (both -0.4SD, p < 0.001) than women without fracture. The magnitude of fracture prediction was similar for LS_aBMD and LS_TBS (odds ratio [95 % confidence interval] = 1.4 [1.2;1.7] and 1.6 [1.2;2.0]). After adjustment for age and prevalent fracture, LS_TBS remained predictive of an increased risk of fracture. Yet, its addition to age, prevalent fracture, and LS_aBMD did not reach the level of significance to improve the fracture prediction. When using the WHO classification, 39 % of fractures occurred in osteoporotic women, 46 % in osteopenic women, and 15 % in women with T-score > -1. Thirty-seven percent of fractures occurred in the lowest quartile of LS_TBS, regardless of BMD. Moreover, 35 % of fractures that occurred in osteopenic women were classified below this LS_TBS threshold. CONCLUSION: In conclusion, LS_aBMD and LS_TBS predicted fractures equally well. In our cohort, the addition of LS_TBS to age and LS_aBMD added only limited information on fracture risk prediction. However, using the lowest quartile of LS_TBS helped in redefining a significant subset of non-osteoporotic women as a higher risk group which is important for patient management.