wnt16 regulates spine and muscle morphogenesis through parallel signals from notochord and dermomyotome.

Bone and muscle are coupled through developmental, mechanical, paracrine, and autocrine signals. Genetic variants at the CPED1-WNT16 locus are dually associated with bone- and muscle-related traits. While Wnt16 is necessary for bone mass and strength, this fails to explain pleiotropy at this locus. Here, we show wnt16 is required for spine and muscle morphogenesis in zebrafish. In embryos, wnt16 is expressed in dermomyotome and developing notochord, and contributes to larval myotome morphology and notochord elongation. Later, wnt16 is expressed at the ventral midline of the notochord sheath, and contributes to spine mineralization and osteoblast recruitment. Morphological changes in wnt16 mutant larvae are mirrored in adults, indicating that wnt16 impacts bone and muscle morphology throughout the lifespan. Finally, we show that wnt16 is a gene of major effect on lean mass at the CPED1-WNT16 locus. Our findings indicate that Wnt16 is secreted in structures adjacent to developing bone (notochord) and muscle (dermomyotome) where it affects the morphogenesis of each tissue, thereby rendering wnt16 expression into dual effects on bone and muscle morphology. This work expands our understanding of wnt16 in musculoskeletal development and supports the potential for variants to act through WNT16 to influence bone and muscle via parallel morphogenetic processes.

Proceedings of the Post-Genome Analysis for Musculoskeletal Biology Workshop.

PURPOSE OF THE REVIEW: Herein, we report on the proceedings of the workshop entitled "Post-Genome analysis for musculoskeletal biology" that was held in July of 2022 in Safed, Galilee, Israel. Supported by the Israel Science Foundation, the goal of this workshop was to bring together established investigators and their trainees who were interested in understanding the etiology of musculoskeletal disease, from Israel and from around the world. RECENT FINDINGS: Presentations at this workshop spanned the spectrum from basic science to clinical studies. A major emphasis of the discussion centered on genetic studies in humans, and the limitations and advantages of such studies. The power of coupling studies using human data with functional follow-up studies in pre-clinical models such as mice, rats, and zebrafish was discussed in depth. The advantages and limitations of mice and zebrafish for faithfully modelling aspects of human disease were debated, specifically in the context of age-related diseases such as osteoporosis, osteoarthritis, adult-onset auto-immune disease, and osteosarcopenia. There remain significant gaps in our understanding of the nature and etiology of human musculoskeletal disease. While therapies and medications exist, much work is still needed to find safe and effective interventions for all patients suffering from diseases associated with age-related deterioration of musculoskeletal tissues. The potential of forward and reverse genetic studies has not been exhausted for diseases of muscles, joints, and bones.

Zebrafish as a Model for Osteoporosis: Functional Validations of Genome-Wide Association Studies.

PURPOSE OF REVIEW: GWAS, as a largely correlational analysis, requires in vitro or in vivo validation. Zebrafish (Danio rerio) have many advantages for studying the genetics of human diseases. Since gene editing in zebrafish has been highly valuable for studying embryonic skeletal developmental processes that are prenatally or perinatally lethal in mammalian models, we are reviewing pros and cons of this model. RECENT FINDINGS: The true power for the use of zebrafish is the ease by which the genome can be edited, especially using the CRISPR/Cas9 system. Gene editing, followed by phenotyping, for complex traits such as BMD, is beneficial, but the major physiological differences between the fish and mammals must be considered. Like mammals, zebrafish do have main bone cells; thus, both in vivo stem cell analyses and in vivo imaging are doable. Yet, the "long" bones of fish are peculiar, and their bone cavities do not contain bone marrow. Partial duplication of the zebrafish genome should be taken into account. Overall, small fish toolkit can provide unmatched opportunities for genetic modifications and morphological investigation as a follow-up to human-first discovery.

Both indirect maternal and direct fetal genetic effects reflect the observational relationship between higher birth weight and lower adult bone mass.

BACKGROUND: Birth weight is considered not only to undermine future growth, but also to induce lifelong diseases; the aim of this study is to explore the relationship between birth weight and adult bone mass. METHODS: We performed multivariable regression analyses to assess the association of birth weight with bone parameters measured by dual-energy X-ray absorptiometry (DXA) and by quantitative ultrasound (QUS), independently. We also implemented a systemic Mendelian randomization (MR) analysis to explore the causal association between them with both fetal-specific and maternal-specific instrumental variables. RESULTS: In the observational analyses, we found that higher birth weight could increase the adult bone area (lumbar spine, ß-coefficient= 0.17, P < 2.00 × 10-16; lateral spine, ß-coefficient = 0.02, P = 0.04), decrease bone mineral content-adjusted bone area (BMCadjArea) (lumbar spine, ß-coefficient= - 0.01, P = 2.27 × 10-14; lateral spine, ß-coefficient = - 0.05, P = 0.001), and decrease adult bone mineral density (BMD) (lumbar spine, ß-coefficient = - 0.04, P = 0.007; lateral spine; ß-coefficient = - 0.03, P = 0.02; heel, ß-coefficient = - 0.06, P < 2.00 × 10-16), and we observed that the effect of birth weight on bone size was larger than that on BMC. In MR analyses, the higher fetal-specific genetically determined birth weight was identified to be associated with higher bone area (lumbar spine; ß-coefficient = 0.15, P = 1.26 × 10-6, total hip, ß-coefficient = 0.15, P = 0.005; intertrochanteric area, ß-coefficient = 0.13, P = 0.0009; trochanter area, ß-coefficient = 0.11, P = 0.03) but lower BMD (lumbar spine, ß-coefficient = - 0.10, P = 0.01; lateral spine, ß-coefficient = - 0.12, P = 0.0003, and heel ß-coefficient = - 0.11, P = 3.33 × 10-13). In addition, we found that the higher maternal-specific genetically determined offspring birth weight was associated with lower offspring adult heel BMD (ß-coefficient = - 0.001, P = 0.04). CONCLUSIONS: The observational analyses suggested that higher birth weight was associated with the increased adult bone area but decreased BMD. By leveraging the genetic instrumental variables with maternal- and fetal-specific effects on birth weight, the observed relationship could be reflected by both the direct fetal and indirect maternal genetic effects.

Acute hypoxia elevates arginase 2 and induces polyamine stress response in zebrafish via evolutionarily conserved mechanism.

Living organisms repeatedly encounter stressful events and apply various strategies to survive. Polyamines are omnipresent bioactive molecules with multiple functions. Their transient synthesis, inducible by numerous stressful stimuli, is termed the polyamine stress response. Animals developed evolutionarily conserved strategies to cope with stresses. The urea cycle is an ancient attribute that deals with ammonia excess in terrestrial species. Remarkably, most fish retain the urea cycle genes fully expressed during the early stages of development and silenced in adult animals. Environmental challenges instigate urea synthesis in fish despite substantial energetic costs, which poses the question of the urea cycle's evolutionary significance. Arginase plays a critical role in oxidative stress-dependent reactions being the final urea cycle enzyme. Its unique subcellular localization, high inducibility, and several regulation levels provide a supreme ability to control the polyamine synthesis rate. Notably, oxidative stress instigates the arginase-1 activity in mammals. Arginase is also dysregulated in aging organisms' brain and muscle tissues, indicating its role in the pathogenesis of age-associated diseases. We designed a study to investigate the levels of the urea cycle and polyamine synthesis-related enzymes in a fish model of acute hypoxia. We evidence synchronized elevation of arginase-2 and ornithine decarboxylase following oxidative stress in adult fish and aging animals signifying the specific function of arginase-2 in fish. Moreover, we demonstrate oxidative stress-associated polyamine synthesis' induction and urea cycle' arrest in adult fish. The subcellular arginase localization found in the fish seems to correspond to its possible evolutionary roles.

Genome-wide identification of novel long non-coding RNAs and their possible roles in hypoxic zebrafish brain.

Long non-coding RNAs (lncRNAs) are the master regulators of numerous biological processes. Hypoxia causes oxidative stress with severe and detrimental effects on brain function and acts as a critical initiating factor in the pathogenesis of Alzheimer's disease (AD). From the RNA-Seq in the forebrain (Fb), midbrain (Mb), and hindbrain (Hb) regions of hypoxic and normoxic zebrafish, we identified novel lncRNAs, whose potential cis targets showed involvement in neuronal development and differentiation pathways. Under hypoxia, several lncRNAs and mRNAs were differentially expressed. Co-expression studies indicated that the Fb and Hb regions' potential lncRNA target genes were involved in the AD pathogenesis. In contrast, those in Mb (cry1b, per1a, cipca) was responsible for regulating circadian rhythm. We identified specific lncRNAs present in the syntenic regions between zebrafish and humans, possibly functionally conserved. We thus identified several conserved lncRNAs as the probable regulators of AD genes (adrb3b, cav1, stat3, bace2, apoeb, psen1, s100b).

Life-Course Genome-wide Association Study Meta-analysis of Total Body BMD and Assessment of Age-Specific Effects.

Bone mineral density (BMD) assessed by DXA is used to evaluate bone health. In children, total body (TB) measurements are commonly used; in older individuals, BMD at the lumbar spine (LS) and femoral neck (FN) is used to diagnose osteoporosis. To date, genetic variants in more than 60 loci have been identified as associated with BMD. To investigate the genetic determinants of TB-BMD variation along the life course and test for age-specific effects, we performed a meta-analysis of 30 genome-wide association studies (GWASs) of TB-BMD including 66,628 individuals overall and divided across five age strata, each spanning 15 years. We identified variants associated with TB-BMD at 80 loci, of which 36 have not been previously identified; overall, they explain approximately 10% of the TB-BMD variance when combining all age groups and influence the risk of fracture. Pathway and enrichment analysis of the association signals showed clustering within gene sets implicated in the regulation of cell growth and SMAD proteins, overexpressed in the musculoskeletal system, and enriched in enhancer and promoter regions. These findings reveal TB-BMD as a relevant trait for genetic studies of osteoporosis, enabling the identification of variants and pathways influencing different bone compartments. Only variants in ESR1 and close proximity to RANKL showed a clear effect dependency on age. This most likely indicates that the majority of genetic variants identified influence BMD early in life and that their effect can be captured throughout the life course.

Novel model of restricted mobility induced osteopenia in zebrafish.

Immobilization, such as prolonged bed rest, is a risk factor for bone loss in humans. Motivated by the emerging utility of zebrafish (Danio rerio) as an animal of choice for the study of musculoskeletal disease, here we report a model of restricted mobility induced osteopenia in adult zebrafish. Aquatic tanks with small cubical compartments to restrict the movement and locomotion of single fish were designed and fabricated for this study. Adult zebrafish were divided into two groups: a normal control (CONT) and a restricted mobility group (RMG) (18 fish/group). Six fish from each group were euthanized on days 14, 21 and 35 of the movement restriction. By using microcomputed tomography (micro-CT), we assessed bone volume/tissue volume (BV/TV) and bone density in the whole skeleton of the fish. Furthermore, we assessed skeletal shape in the vertebrae (radius, length, volume, neural and haemal arch aperture areas, neural and haemal arch angle, and thickness of the intervertebral space), single vertebra bone volume and bone density. Movement restriction significantly decreased vertebral skeletal parameters such as radius, length, volume, arch aperture areas and angles as well as the thickness of the intervertebral space in RMG. Furthermore, restricted mobility significantly (P < 0.001) decreased BV/TV and bone density as compared to the CONT group, starting as early as 14 days. By analysing zebrafish from CONT and RMG, we show that micro-CT imaging is a sensitive method to quantify distinct skeletal properties in zebrafish. We further defined the micro-CT parameters which can be used to examine the effects of restricted mobility on the skeleton of the fish. Our findings propose a rapid and effective osteopenia "stabulation" model, which could be used widely for osteoporosis drug screening.

Low-Frequency Synonymous Coding Variation in CYP2R1 Has Large Effects on Vitamin D Levels and Risk of Multiple Sclerosis.

Vitamin D insufficiency is common, correctable, and influenced by genetic factors, and it has been associated with risk of several diseases. We sought to identify low-frequency genetic variants that strongly increase the risk of vitamin D insufficiency and tested their effect on risk of multiple sclerosis, a disease influenced by low vitamin D concentrations. We used whole-genome sequencing data from 2,619 individuals through the UK10K program and deep-imputation data from 39,655 individuals genotyped genome-wide. Meta-analysis of the summary statistics from 19 cohorts identified in CYP2R1 the low-frequency (minor allele frequency = 2.5%) synonymous coding variant g.14900931G>A (p.Asp120Asp) (rs117913124[A]), which conferred a large effect on 25-hydroxyvitamin D (25OHD) levels (-0.43 SD of standardized natural log-transformed 25OHD per A allele; p value = 1.5 × 10-88). The effect on 25OHD was four times larger and independent of the effect of a previously described common variant near CYP2R1. By analyzing 8,711 individuals, we showed that heterozygote carriers of this low-frequency variant have an increased risk of vitamin D insufficiency (odds ratio [OR] = 2.2, 95% confidence interval [CI] = 1.78-2.78, p = 1.26 × 10-12). Individuals carrying one copy of this variant also had increased odds of multiple sclerosis (OR = 1.4, 95% CI = 1.19-1.64, p = 2.63 × 10-5) in a sample of 5,927 case and 5,599 control subjects. In conclusion, we describe a low-frequency CYP2R1 coding variant that exerts the largest effect upon 25OHD levels identified to date in the general European population and implicates vitamin D in the etiology of multiple sclerosis.

Plasma exosomes stimulate breast cancer metastasis through surface interactions and activation of FAK signaling.

PURPOSE: The interaction between malignant cells and surrounding healthy tissues is a critical factor in the metastatic progression of breast cancer (BC). Extracellular vesicles, especially exosomes, are known to be involved in inter-cellular communication during cancer progression. In the study presented herein, we aimed to evaluate the role of circulating plasma exosomes in the metastatic dissemination of BC and to investigate the underlying molecular mechanisms of this phenomenon. METHODS: Exosomes isolated from plasma of healthy female donors were applied in various concentrations into the medium of MDA-MB-231 and MCF-7 cell lines. Motility and invasive properties of BC cells were examined by random migration and Transwell invasion assays, and the effect of plasma exosomes on the metastatic dissemination of BC cells was demonstrated in an in vivo zebrafish model. To reveal the molecular mechanism of interaction between plasma exosomes and BC cells, a comparison between un-treated and enzymatically modified exosomes was performed, followed by mass spectrometry, gene ontology, and pathway analysis. RESULTS: Plasma exosomes stimulated the adhesive properties, two-dimensional random migration, and transwell invasion of BC cells in vitro as well as their in vivo metastatic dissemination in a dose-dependent manner. This stimulatory effect was mediated by interactions of surface exosome proteins with BC cells and consequent activation of focal adhesion kinase (FAK) signaling in the tumor cells. CONCLUSIONS: Plasma exosomes have a potency to stimulate the metastasis-promoting properties of BC cells. This pro-metastatic property of normal plasma exosomes may have impact on the course of the disease and on its prognosis.

Genetics of Bone and Muscle Interactions in Humans.

PURPOSE OF REVIEW: To summarize the evidence from recent studies on the shared genetics between bone and muscle in humans. RECENT FINDINGS: Genome-wide association studies (GWAS) have successfully identified a multitude of loci influencing the variability of different bone or muscle parameters, with multiple loci overlapping between the traits. In addition, joint analyses of multiple correlated musculoskeletal traits (i.e., multivariate GWAS) have underscored several genes with possible pleiotropic effects on both bone and muscle including MEF2C and SREBF1. Notably, several of the proposed pleiotropic genes have been validated using human cells or animal models. It is clear that the study of pleiotropy may provide novel insights into disease pathophysiology potentially leading to the identification of new treatment strategies that simultaneously prevent or treat both osteoporosis and sarcopenia. However, the role of muscle factors (myokines) that stimulate bone metabolism, as well as osteokines that affect muscles, is in its earliest stage of understanding.

Targeted sequencing of genome wide significant loci associated with bone mineral density (BMD) reveals significant novel and rare variants: the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) targeted sequencing study.

Background: Bone mineral density (BMD) is a heritable phenotype that predicts fracture risk. We performed fine-mapping by targeted sequencing at WLS, MEF2C, ARHGAP1/F2 and JAG1 loci prioritized by eQTL and bioinformatic approaches among 56 BMD loci from our previous GWAS meta-analysis. Targeted sequencing was conducted in 1,291 Caucasians from the Framingham Heart Study ( n = 925) and Cardiovascular Health Study ( n = 366), including 206 women and men with extreme low femoral neck (FN) BMD. A total of 4,964 sequence variants (SNVs) were observed and 80% were rare with MAF <1%. The associations between previously identified SNPs in these loci and BMD, while nominally significant in sequenced participants, were no longer significant after multiple testing corrections. Conditional analyses did not find protein-coding variants that may be responsible for GWAS signals. On the other hand, in the sequenced subjects, we identified novel associations in WLS , ARHGAP1 , and 5' of MEF2C ( P- values < 8x10 - 5 ; false discovery rate (FDR) q-values < 0.01) that were much more strongly associated with BMD compared to the GWAS SNPs. These associated SNVs are less-common; independent from previous GWAS signals in the same loci; and located in gene regulatory elements. Our findings suggest that protein-coding variants in selected GWAS loci did not contribute to GWAS signals. By performing targeted sequencing in GWAS loci, we identified less-common and rare non-coding SNVs associated with BMD independently from GWAS common SNPs, suggesting both common and less-common variants may associate with disease risks and phenotypes in the same loci.

Students as anatomy near-peer teachers: a double-edged sword for an ancient skill.

BACKGROUND: A near-peer instructors (NPI) program was designed for 1st year medical students who successfully finished the Anatomy course, in order to develop their didactic ability and teaching skills, mostly for cadaver dissection. METHODS: Graduates of the training program were administered a voluntary survey at the end of the program, annually. Best graduates of the training program were offered a NPI position in the next academic year. They were evaluated by the first-year students, at the end of the Anatomy block. RESULTS: In a debriefing questionnaire at the end of the NPI training, on the five-point Likert scale (1 = lowest to 5 = highest), the overall rating ranged from 3.63 in 2013 to 3.71 in 2015. Learning prosection and anatomy demonstration skills scored on average from 4.30 to 4.36, respectively. The NPIs were then evaluated by first-year students at the end of the next year's Anatomy block. On the Likert scale, the average score of NPIs ranged from 4.10 in 2014 to 4.75 in 2016, on the par with the general satisfaction score for the professional preclinical teachers during the same period (which ranged from 3.80 to 4.26). CONCLUSIONS: It is suggested that students as near-peer instructors can make a valuable contribution to the teaching faculty, especially in a new medical school.

FABP4 is a leading candidate gene associated with residual feed intake in growing Holstein calves.

Ecological and economic concerns drive the need to improve feed utilization by domestic animals. Residual feed intake (RFI) is one of the most acceptable measures for feed efficiency (FE). However, phenotyping RFI-related traits is complex and expensive and requires special equipment. Advances in marker technology allow the development of various DNA-based selection tools. To assimilate these technologies for the benefit of RFI-based selection, reliable phenotypic measures are prerequisite. In the current study, we identified single nucleotide polymorphisms (SNPs) associated with RFI phenotypic consistency across different ages and diets (named RFI 1-3), using DNA samples of high or low RFI ranked Holstein calves. Using targeted sequencing of chromosomal regions associated with FE- and RFI-related traits, we identified 48 top SNPs significantly associated with at least one of three defined RFIs. Eleven of these SNPs were harbored by the fatty acid binding protein 4 (FABP4). While 10 significant SNPs found in FABP4 were common for RFI 1 and RFI 3, one SNP (FABP4_5; A

Meta-analysis of genome-wide association studies identifies two loci associated with circulating osteoprotegerin levels.

Osteoprotegerin (OPG) is involved in bone homeostasis and tumor cell survival. Circulating OPG levels are also important biomarkers of various clinical traits, such as cancers and atherosclerosis. OPG levels were measured in serum or in plasma. In a meta-analysis of genome-wide association studies in up to 10 336 individuals from European and Asian origin, we discovered that variants >100 kb upstream of the TNFRSF11B gene encoding OPG and another new locus on chromosome 17q11.2 were significantly associated with OPG variation. We also identified a suggestive locus on chromosome 14q21.2 associated with the trait. Moreover, we estimated that over half of the heritability of OPG levels could be explained by all variants examined in our study. Our findings provide further insight into the genetic regulation of circulating OPG levels.

Genetic determinants of heel bone properties: genome-wide association meta-analysis and replication in the GEFOS/GENOMOS consortium.

Quantitative ultrasound of the heel captures heel bone properties that independently predict fracture risk and, with bone mineral density (BMD) assessed by X-ray (DXA), may be convenient alternatives for evaluating osteoporosis and fracture risk. We performed a meta-analysis of genome-wide association (GWA) studies to assess the genetic determinants of heel broadband ultrasound attenuation (BUA; n = 14 260), velocity of sound (VOS; n = 15 514) and BMD (n = 4566) in 13 discovery cohorts. Independent replication involved seven cohorts with GWA data (in silico n = 11 452) and new genotyping in 15 cohorts (de novo n = 24 902). In combined random effects, meta-analysis of the discovery and replication cohorts, nine single nucleotide polymorphisms (SNPs) had genome-wide significant (P < 5 × 10(-8)) associations with heel bone properties. Alongside SNPs within or near previously identified osteoporosis susceptibility genes including ESR1 (6q25.1: rs4869739, rs3020331, rs2982552), SPTBN1 (2p16.2: rs11898505), RSPO3 (6q22.33: rs7741021), WNT16 (7q31.31: rs2908007), DKK1 (10q21.1: rs7902708) and GPATCH1 (19q13.11: rs10416265), we identified a new locus on chromosome 11q14.2 (rs597319 close to TMEM135, a gene recently linked to osteoblastogenesis and longevity) significantly associated with both BUA and VOS (P < 8.23 × 10(-14)). In meta-analyses involving 25 cohorts with up to 14 985 fracture cases, six of 10 SNPs associated with heel bone properties at P < 5 × 10(-6) also had the expected direction of association with any fracture (P < 0.05), including three SNPs with P < 0.005: 6q22.33 (rs7741021), 7q31.31 (rs2908007) and 10q21.1 (rs7902708). In conclusion, this GWA study reveals the effect of several genes common to central DXA-derived BMD and heel ultrasound/DXA measures and points to a new genetic locus with potential implications for better understanding of osteoporosis pathophysiology.

Models to explore genetics of human aging.

Genetic studies have bestowed insight into the biological mechanisms underlying inter-individual differences in susceptibility to (or resistance to) organisms' aging. Recent advances in molecular and genetic epidemiology provide tools to explore the genetic sources of the variability in biological aging in humans. To be successful, the genetic study of a complex condition such as aging requires the clear definition of essential traits that can characterize the aging process phenotypically. Phenotypes of human aging have long relied on mortality rate or exceptional longevity. Genome-wide association studies (GWAS) have been shown to present an unbiased approach to the identification of new candidate genes for human diseases. The GWAS approach can also be used for positive health phenotypes such as longevity or a delay in age-related chronic disease, as well as for other age related changes such as loss of telomere length or lens transparency. Sequencing, either in targeted regions or across the whole genome can further identify rare variation that may contribute to the biological aging mechanisms. To date, the results of the GWAS for longevity are rather disappointing, possibly in part due to the small number of individuals with GWAS data who have reached advanced old age.Human aging phenotypes are needed that can be assessed prior to death, and should be both heritable and validated as predictors of longevity. Potentially, phenotypes that focus on "successful" or "healthy" aging will be more powerful as they can be measured in large numbers of people and also are clinically relevant.We postulate that construction of an integrated phenotype of aging can be achieved capitalizing on multiple traits that may have weak correlations, but a shared underlying genetic architecture. This is based on a hypothesis that convergent results from multiple individual aging-related traits will point out the pleiotropic signals responsible for the overall rate of aging of the organism. An approach would be to investigate traits that are linked to the state of many vital functions, disability, and ultimately survival rates, to identify common biological pathways that govern aging processes in humans. New composite aging phenotypes must be validated by predicting all-cause mortality, major chronic disease and disability late in life.

A genome-wide association meta-analysis of circulating sex hormone-binding globulin reveals multiple Loci implicated in sex steroid hormone regulation.

Sex hormone-binding globulin (SHBG) is a glycoprotein responsible for the transport and biologic availability of sex steroid hormones, primarily testosterone and estradiol. SHBG has been associated with chronic diseases including type 2 diabetes (T2D) and with hormone-sensitive cancers such as breast and prostate cancer. We performed a genome-wide association study (GWAS) meta-analysis of 21,791 individuals from 10 epidemiologic studies and validated these findings in 7,046 individuals in an additional six studies. We identified twelve genomic regions (SNPs) associated with circulating SHBG concentrations. Loci near the identified SNPs included SHBG (rs12150660, 17p13.1, p = 1.8 × 10(-106)), PRMT6 (rs17496332, 1p13.3, p = 1.4 × 10(-11)), GCKR (rs780093, 2p23.3, p = 2.2 × 10(-16)), ZBTB10 (rs440837, 8q21.13, p = 3.4 × 10(-09)), JMJD1C (rs7910927, 10q21.3, p = 6.1 × 10(-35)), SLCO1B1 (rs4149056, 12p12.1, p = 1.9 × 10(-08)), NR2F2 (rs8023580, 15q26.2, p = 8.3 × 10(-12)), ZNF652 (rs2411984, 17q21.32, p = 3.5 × 10(-14)), TDGF3 (rs1573036, Xq22.3, p = 4.1 × 10(-14)), LHCGR (rs10454142, 2p16.3, p = 1.3 × 10(-07)), BAIAP2L1 (rs3779195, 7q21.3, p = 2.7 × 10(-08)), and UGT2B15 (rs293428, 4q13.2, p = 5.5 × 10(-06)). These genes encompass multiple biologic pathways, including hepatic function, lipid metabolism, carbohydrate metabolism and T2D, androgen and estrogen receptor function, epigenetic effects, and the biology of sex steroid hormone-responsive cancers including breast and prostate cancer. We found evidence of sex-differentiated genetic influences on SHBG. In a sex-specific GWAS, the loci 4q13.2-UGT2B15 was significant in men only (men p = 2.5 × 10(-08), women p = 0.66, heterogeneity p = 0.003). Additionally, three loci showed strong sex-differentiated effects: 17p13.1-SHBG and Xq22.3-TDGF3 were stronger in men, whereas 8q21.12-ZBTB10 was stronger in women. Conditional analyses identified additional signals at the SHBG gene that together almost double the proportion of variance explained at the locus. Using an independent study of 1,129 individuals, all SNPs identified in the overall or sex-differentiated or conditional analyses explained ~15.6% and ~8.4% of the genetic variation of SHBG concentrations in men and women, respectively. The evidence for sex-differentiated effects and allelic heterogeneity highlight the importance of considering these features when estimating complex trait variance.

Genetic diversity is a predictor of mortality in humans.

BACKGROUND: It has been well-established, both by population genetics theory and direct observation in many organisms, that increased genetic diversity provides a survival advantage. However, given the limitations of both sample size and genome-wide metrics, this hypothesis has not been comprehensively tested in human populations. Moreover, the presence of numerous segregating small effect alleles that influence traits that directly impact health directly raises the question as to whether global measures of genomic variation are themselves associated with human health and disease. RESULTS: We performed a meta-analysis of 17 cohorts followed prospectively, with a combined sample size of 46,716 individuals, including a total of 15,234 deaths. We find a significant association between increased heterozygosity and survival (P = 0.03). We estimate that within a single population, every standard deviation of heterozygosity an individual has over the mean decreases that person's risk of death by 1.57%. CONCLUSIONS: This effect was consistent between European and African ancestry cohorts, men and women, and major causes of death (cancer and cardiovascular disease), demonstrating the broad positive impact of genomic diversity on human survival.