Adult height variants affect birth length and growth rate in children

Characterizing differences in the muscle transcriptome between cattle with alternative LCORL-NCAPG haplotypes

BACKGROUND

The LCORL-NCAPG locus is a major quantitative trait locus (QTL) on bovine chromosome 6 (BTA6) that influences growth and carcass composition in cattle. To further understand the molecular mechanism responsible for the phenotypic changes associated with this locus, twenty-four Charolais-sired calves were selected for muscle transcriptome analysis based on alternative homozygous LCORL-NCAPG haplotypes (i.e., 12 "QQ" and 12 "qq", where "Q" is a haplotype harboring variation associated with increased growth). At 300 days of age, a biopsy of the longissimus dorsi muscle was collected from each animal for RNA sequencing.

RESULTS

Gene expression analysis identified 733 genes as differentially expressed between QQ and qq animals (q-value < 0.05). Notably, LCORL and genes known to be important regulators of growth such as IGF2 were upregulated in QQ individuals, while genes associated with adiposity such as FASN and LEP were downregulated, reflecting the increase in lean growth associated with this locus. Gene set enrichment analysis demonstrated QQ individuals had downregulation of pathways associated with adipogenesis, alongside upregulation of transcripts for cellular machinery essential for protein synthesis and energy metabolism, particularly ribosomal and mitochondrial components.

CONCLUSIONS

The differences in the muscle transcriptome between QQ and qq animals imply that muscle hypertrophy may be metabolically favored over accumulation of fat in animals with the QQ haplotype. Our findings also suggest this haplotype could be linked to a difference in LCORL expression that potentially influences the downstream transcriptional effects observed, though further research will be needed to confirm the molecular mechanisms underlying the associated changes in phenotype.

A Developmentally-Informative Genome-wide Association Study of Alcohol Use Frequency

Contemporary genome-wide association study (GWAS) methods typically do not account for variability in genetic effects throughout development. We applied genomic structural equation modeling to combine developmentally-informative phenotype data and GWAS to create polygenic scores (PGS) for alcohol use frequency that are specific to developmental stage. Longitudinal cohort studies targeted for gene-identification analyses include the Collaborative Study on the Genetics of Alcoholism (adolescence n = 1,118, early adulthood n = 2,762, adulthood n = 5,255), the National Longitudinal Study of Adolescent to Adult Health (adolescence n = 3,089, early adulthood n = 3,993, adulthood n = 5,149), and the Avon Longitudinal Study of Parents and Children (ALSPAC; adolescence n = 5,382, early adulthood n = 3,613). PGS validation analyses were conducted in the COGA sample using an alternate version of the discovery analysis with COGA removed. Results suggest that genetic liability for alcohol use frequency in adolescence may be distinct from genetic liability for alcohol use frequency later in developmental periods. The age-specific PGS predicts an increase of 4 drinking days per year per PGS standard deviation when modeled separately from the common factor PGS in adulthood. The current work was underpowered at all steps of the analysis plan. Though small sample sizes and low statistical power limit the substantive conclusions that can be drawn regarding these research questions, this work provides a foundation for future genetic studies of developmental variability in the genetic underpinnings of alcohol use behaviors and genetically-informed, age-matched phenotype prediction.

A novel approach to risk exposure and epigenetics-the use of multidimensional context to gain insights into the early origins of cardiometabolic and neurocognitive health

BACKGROUND

Each mother-child dyad represents a unique combination of genetic and environmental factors. This constellation of variables impacts the expression of countless genes. Numerous studies have uncovered changes in DNA methylation (DNAm), a form of epigenetic regulation, in offspring related to maternal risk factors. How these changes work together to link maternal-child risks to childhood cardiometabolic and neurocognitive traits remains unknown. This question is a key research priority as such traits predispose to future non-communicable diseases (NCDs). We propose viewing risk and the genome through a multidimensional lens to identify common DNAm patterns shared among diverse risk profiles.

METHODS

We identified multifactorial Maternal Risk Profiles (MRPs) generated from population-based data (n = 15,454, Avon Longitudinal Study of Parents and Children (ALSPAC)). Using cord blood HumanMethylation450 BeadChip data, we identified genome-wide patterns of DNAm that co-vary with these MRPs. We tested the prospective relation of these DNAm patterns (n = 914) to future outcomes using decision tree analysis. We then tested the reproducibility of these patterns in (1) DNAm data at age 7 and 17 years within the same cohort (n = 973 and 974, respectively) and (2) cord DNAm in an independent cohort, the Generation R Study (n = 686).

RESULTS

We identified twenty MRP-related DNAm patterns at birth in ALSPAC. Four were prospectively related to cardiometabolic and/or neurocognitive childhood outcomes. These patterns were replicated in DNAm data from blood collected at later ages. Three of these patterns were externally validated in cord DNAm data in Generation R. Compared to previous literature, DNAm patterns exhibited novel spatial distribution across the genome that intersects with chromatin functional and tissue-specific signatures.

CONCLUSIONS

To our knowledge, we are the first to leverage multifactorial population-wide data to detect patterns of variability in DNAm. This context-based approach decreases biases stemming from overreliance on specific samples or variables. We discovered molecular patterns demonstrating prospective and replicable relations to complex traits. Moreover, results suggest that patterns harbour a genome-wide organisation specific to chromatin regulation and target tissues. These preliminary findings warrant further investigation to better reflect the reality of human context in molecular studies of NCDs.

Evaluation of SNPs associated with mammographic density in European women with mammographic density in Asian women from South-East Asia

PURPOSE

Mammographic density (MD), after accounting for age and body mass index (BMI), is a strong heritable risk factor for breast cancer. Genome-wide association studies (GWAS) have identified 64 SNPs in 55 independent loci associated with MD in women of European ancestry. Their associations with MD in Asian women, however, are largely unknown.

METHOD

Using linear regression adjusting for age, BMI, and ancestry-informative principal components, we evaluated the associations of previously reported MD-associated SNPs with MD in a multi-ethnic cohort of Asian ancestry. Area and volumetric mammographic densities were determined using STRATUS (N = 2450) and Volpara™ (N = 2257). We also assessed the associations of these SNPs with breast cancer risk in an Asian population of 14,570 cases and 80,870 controls.

RESULTS

Of the 61 SNPs available in our data, 21 were associated with MD at a nominal threshold of P value < 0.05, all in consistent directions with those reported in European ancestry populations. Of the remaining 40 variants with a P-value of association > 0.05, 29 variants showed consistent directions of association as those previously reported. We found that nine of the 21 MD-associated SNPs in this study were also associated with breast cancer risk in Asian women (P < 0.05), seven of which showed a direction of associations that was consistent with that reported for MD.

CONCLUSION

Our study confirms the associations of 21 SNPs (19/55 or 34.5% out of all known MD loci identified in women of European ancestry) with area and/or volumetric densities in Asian women, and further supports the evidence of a shared genetic basis through common genetic variants for MD and breast cancer risk.

Effect of L3MBTL3/PTPN9 polymorphisms on risk to alcohol-induced ONFH in Chinese Han population

PURPOSE

Alcohol-induced osteonecrosis femoral head necrosis (ONFH) is a disease that seriously affects human health. Abnormal expression of L3MBTL3/PTPN9 gene can cause a variety of human diseases. The purpose of this study is to investigate the effect of L3MBTL3/PTPN9 gene polymorphism on the susceptibility of alcohol-induced ONFH in Chinese Han population.

METHODS

A total of 308 alcohol-induced ONFH patients and 425 healthy controls were enrolled in this case-control study. Alleles, genotypes, genetic models, haplotypes, and multifactor dimensionality reduction analyses (MDR) based on age-corrected by using odds ratio (OR) and 95% confidence interval (CI) were performed.

RESULTS

Our result revealed rs2068957 in the L3MBTL3 gene increased the risk of alcohol ONFH under the recessive model after correction. Besides, we also found that rs75393192 in the PTPN9 gene was a protective site in stratification over 40 years of age and stage. In stratified analysis of necrotic sites, we only found that rs2068957 was associated with increased susceptibility of alcohol-induced ONFH under the co-dominant model and recessive model. Haplotype "GC" in the block (rs76107647|rs10851882 in PTPN9 gene) significantly decreased the susceptibility of alcoholic ONFH.

CONCLUSIONS

Our results provide evidence that L3MBTL3/PTPN9 polymorphisms are associated with alcohol-induced ONFH risk in Chinese Han population.

Polygenic risk for depression, anxiety and neuroticism are associated with the severity and rate of change in depressive symptoms across adolescence

BACKGROUND

Adolescence marks a period where depression will commonly onset. Twin studies show that genetic influences play a role in how depression develops and changes across adolescence. Recent genome-wide association studies highlight that common genetic variants - which can be combined into polygenic risk scores (PRS) - are also implicated in depression. However, the role of PRS in adolescent depression and changes in adolescent depression is not yet understood. We aimed to examine associations between PRS for five psychiatric traits and depressive symptoms measured across adolescence using cross-sectional and growth-curve models. The five PRS were as follows: depression (DEP), major depressive disorder (MDD), anxiety (ANX), neuroticism (NEU) and schizophrenia (SCZ).

METHODS

We used data from over 6,000 participants of the Avon Longitudinal Study of Parents and Children (ALSPAC) to examine associations between the five PRS and self-reported depressive symptoms (Short Mood and Feelings Questionnaire) over 9 occasions from 10 to 24 years. The PRS were created from well-powered genome-wide association studies conducted in adult populations. We examined cross-sectional associations between the PRS at each age and then again with longitudinal trajectories of depressive symptoms in a repeated measures framework using multilevel growth-curve analysis to examine the severity and the rate of change.

RESULTS

There was strong evidence that higher PRS for DEP, MDD and NEU were associated with worse depressive symptoms throughout adolescence and into young adulthood in our cross-sectional analysis, with consistent associations observed across all nine occasions. Growth-curve analyses provided stronger associations (as measured by effect sizes) and additional insights, demonstrating that individuals with higher PRS for DEP, MDD and NEU had steeper trajectories of depressive symptoms across development, all with a greater increasing rate of change during adolescence. Evidence was less consistent for the ANX and SCZ PRS in the cross-sectional analysis, yet there was some evidence for an increasing rate of change in adolescence in the growth-curve analyses with the ANX PRS.

CONCLUSIONS

These results show that common genetic variants as indexed by varying psychiatric PRS show patterns of specificity that influence both the severity and rate of change in depressive symptoms throughout adolescence and then into young adulthood. Longitudinal data that make use of repeated measures designs have the potential to provide greater insights how genetic factors influence the onset and persistence of adolescent depression.

Novel loci and Mapuche genetic ancestry are associated with pubertal growth traits in Chilean boys

Puberty is a complex developmental process that varies considerably among individuals and populations. Genetic factors explain a large proportion of the variability of several pubertal traits. Recent genome-wide association studies (GWAS) have identified hundreds of variants involved in traits that result from body growth, like adult height. However, they do not capture many genetic loci involved in growth changes over distinct growth phases. Further, such GWAS have been mostly performed in Europeans, but it is unknown how these findings relate to other continental populations. In this study, we analyzed the genetic basis of three pubertal traits; namely, peak height velocity (PV), age at PV (APV) and height at APV (HAPV). We analyzed a cohort of 904 admixed Chilean children and adolescents with European and Mapuche Native American ancestries. Height was measured on roughly a \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$6-$$\end{document}6-month basis from childhood to adolescence between 2006 and 2019. We predict that, in average, HAPV is 4.3 cm higher in European than in Mapuche adolescents (P = 0.042), and APV is 0.73 years later in European compared with Mapuche adolescents (P = 0.023). Further, by performing a GWAS on 774, 433 single-nucleotide polymorphisms, we identified a genetic signal harboring 3 linked variants significantly associated with PV in boys (P\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$< 5 \times 10^{-8}$$\end{document}<5×10-8). This signal has never been associated with growth-related traits.

Can shorter mothers have taller children? Nutritional mobility, health equity and the intergenerational transmission of relative height

This study develops the concept of nutritional mobility, defined here as the probability that a mother ranked low in her cohort's height distribution will have a child who attains a higher rank order. We demonstrate that rank-order regression provides a robust metric of health equity, revealing differences in opportunities for each child to reach their own growth potential. We estimate four indicators of nutritional mobility and test for associations between nutritional mobility and various local economic and environmental factors. Nutritional mobility has improved over time, and the nutrition environment contributes about 2.86 times as much as a mother's height to her child's expected rank in height-for-age. Populations with the least mobility are in Latin America, and the most mobility is in more urbanized areas of Africa and Asia. Rank-order mobility is an important aspect of health equity, offering valuable insight into the role of socioecological factors in nutrition improvement across generations.

The Pathophysiological Significance of Fibulin-3

Fibulin-3 (also known as EGF-containing fibulin extracellular matrix protein 1 (EFEMP1)) is a secreted extracellular matrix glycoprotein, encoded by the EFEMP1 gene that belongs to the eight-membered fibulin protein family. It has emerged as a functionally unique member of this family, with a diverse array of pathophysiological associations predominantly centered on its role as a modulator of extracellular matrix (ECM) biology. Fibulin-3 is widely expressed in the human body, especially in elastic-fibre-rich tissues and ocular structures, and interacts with enzymatic ECM regulators, including tissue inhibitor of metalloproteinase-3 (TIMP-3). A point mutation in EFEMP1 causes an inherited early-onset form of macular degeneration called Malattia Leventinese/Doyne honeycomb retinal dystrophy (ML/DHRD). EFEMP1 genetic variants have also been associated in genome-wide association studies with numerous complex inherited phenotypes, both physiological (namely, developmental anthropometric traits) and pathological (many of which involve abnormalities of connective tissue function). Furthermore, EFEMP1 expression changes are implicated in the progression of numerous types of cancer, an area in which fibulin-3 has putative significance as a therapeutic target. Here we discuss the potential mechanistic roles of fibulin-3 in these pathologies and highlight how it may contribute to the development, structural integrity, and emergent functionality of the ECM and connective tissues across a range of anatomical locations. Its myriad of aetiological roles positions fibulin-3 as a molecule of interest across numerous research fields and may inform our future understanding and therapeutic approach to many human diseases in clinical settings.

Genetic associations of myopia severities and endophenotypes in children

OBJECTIVE

To investigate the associations of multiple single-nucleotide polymorphisms (SNPs) with the severities and endophenotypes of myopia in children.

METHODS

A total of 3300 children aged 5-10 years were recruited: 137 moderate and high myopia (SE≤-3.0D), 670 mild myopia (-3.0D<SE≤-0.5D) and 2493 controls (SE>-0.5D). 13 SNPs in 13 genes/loci were selected for genotyping in all subjects using TaqMan assays. Associations between each SNP with myopia severities and ocular traits (spherical equivalent (SE), axial length (AL) and corneal radius (CR)) were analysed.

RESULTS

When compared with controls, SNPs ZC3H11B rs4373767 (OR=1.15, p=0.038), BICC1 rs7084402 (OR=1.18, p=0.005) and GJD2 rs524952 (OR=1.14, p=0.025) showed nominal associations with overall myopia. ZC3H11B rs4373767 and BICC1 rs7084402 showed stronger associations with moderate and high myopia (rs4373767: OR=1.42, p=0.018; rs7084402: OR=1.33, p=0.025), while GJD2 rs524952 had a stronger association with mild myopia (OR=1.14, p=0.025). GJD2 rs524952 also showed a difference between emmetropia and hyperopia (p=0.018). In quantitative trait locus analysis, ZC3H11B rs4373767, KCNQ5 rs7744813 and GJD2 rs524952 were correlated with both myopic SE (β=-0.09, p=0.03; β=-0.12, p=0.007; β=-0.13, p=0.0006, respectively) and AL (β=0.07, p=0.002; β=0.09, p=0.0008; β=0.07, p=0.0003, respectively). SNTB1 rs7839488 was correlated with both AL (β=0.07, p=0.005) and CR (β=0.02, p=0.006). Moreover, ZC3H11B rs4373767-T (β=0.006; p=0.018), KCNQ5 rs7744813-A (β=0.007; p=0.015) and GJD2 rs524952-T (β=0.009; p=0.0006) were correlated with AL-CR ratio.

CONCLUSIONS AND RELEVANCE

ZC3H11B and BICC1 are genetic risk factors for moderate and high myopia, while ZC3H11B, KCNQ5, SNTB1 and GJD2 confer risk to excessive AL in children.

Genetic Polymorphisms Associated with Idiopathic Short Stature and First-Year Response to Growth Hormone Treatment

BACKGROUND/AIMS

The term idiopathic short stature (ISS) describes short stature of unknown, but likely polygenic, etiology. This study aimed to identify genetic polymorphisms associated with the ISS phenotype, and with growth response to supplemental GH.

METHODS

Using a case-control analysis we compared the prevalence of "tall" versus "short" alleles at 52 polymorphic loci (17 in growth-related candidate genes, 35 identified in prior genome-wide association studies of adult height) in 94 children with ISS followed in the Genetics and Neuroendocrinology of Short Stature International Study, versus 143 controls from the Fels Longitudinal Study.

RESULTS

Four variants were nominally associated with ISS using a genotypic model, confirmed by a simultaneous confident inference approach: compared with controls children with ISS had lower odds of "tall" alleles (odds ratio, 95% CI) for GHR (0.52, 0.29-0.96); rs2234693/ESR1 (0.50, 0.25-0.98); rs967417/BMP2 (0.39, 0.17-0.93), and rs4743034/ZNF462 (0.40, 0.18-0.89). Children with ISS also had lower odds of the "tall" allele (A) at the IGFBP3 -202 promoter polymorphism (rs2855744; 0.40, 0.20-0.80) in the simultaneous confident inference analysis. A significant association with 1st-year height SD score increase during GH treatment was observed with rs11205277, located near 4 known genes: MTMR11, SV2A, HIST2H2AA3, and SF3B4; the latter, in which heterozygous mutations occur in Nager acrofacial dysostosis, appears the most relevant gene.

CONCLUSIONS

In children with ISS we identified associations with "short" alleles at a number of height-related loci. In addition, a polymorphic variant located near SF3B4 was associated with the GH treatment response in our cohort. The findings in our small study warrant further investigation.

Associations of adult-attained height and early life energy restriction with postmenopausal breast cancer risk according to estrogen and progesterone receptor status

Adult-attained height is a marker for underlying mechanisms, such as cell growth, that may also influence postmenopausal breast cancer (BC) risk, perhaps specifically hormone-sensitive BC subtypes. Early life energy restriction may inhibit these mechanisms, resulting in shorter height and a reduced postmenopausal BC risk. Women (62,573) from the Netherlands Cohort Study completed a self-administered questionnaire in 1986 when 55-69 years old, and were followed-up for 20.3 years (case-cohort: N_subcohort = 2,438; N_cases = 3,354). Cox multivariable-adjusted hazard ratios (HR) and 95% confidence intervals (CI) were estimated for BC risk overall and by estrogen and progesterone receptor subtypes in relation to height and early life energy restriction during the Hunger Winter, War Years, and Economic Depression. Although energy restriction can only influence longitudinal growth in women exposed before and/or during the growth spurt, it may also influence BC risk when occurring after the growth spurt, possibly through different growth processes. Therefore, Cox analyses were additionally conducted according to timing of energy restriction in relation to the growth spurt. Height was associated with an increased BC risk (HR_{per 5cm} = 1.07, 95%CI:1.01-1.13), particularly hormone receptor-positive BC. Energy restriction before and/or during the growth spurt was associated with a decreased hormone receptor-positive BC risk. Energy restriction during the Hunger Winter increased the estrogen receptor-negative BC risk regardless of the timing of energy restriction. In conclusion, height and energy restriction before and/or during the growth spurt were both associated with hormone receptor-positive BC risk, in the direction as expected, indicating critical exposure windows for hormonal growth-related mechanisms.

A Large Multiethnic Genome-Wide Association Study of Adult Body Mass Index Identifies Novel Loci

Body mass index (BMI), a proxy measure for obesity, is determined by both environmental (including ethnicity, age, and sex) and genetic factors, with > 400 BMI-associated loci identified to date. However, the impact, interplay, and underlying biological mechanisms among BMI, environment, genetics, and ancestry are not completely understood. To further examine these relationships, we utilized 427,509 calendar year-averaged BMI measurements from 100,418 adults from the single large multiethnic Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort. We observed substantial independent ancestry and nationality differences, including ancestry principal component interactions and nonlinear effects. To increase the list of BMI-associated variants before assessing other differences, we conducted a genome-wide association study (GWAS) in GERA, with replication in the Genetic Investigation of Anthropomorphic Traits (GIANT) consortium combined with the UK Biobank (UKB), followed by GWAS in GERA combined with GIANT, with replication in the UKB. We discovered 30 novel independent BMI loci (P < 5.0 × 10^-8) that replicated. We then assessed the proportion of BMI variance explained by sex in the UKB using previously identified loci compared to previously and newly identified loci and found slight increases: from 3.0 to 3.3% for males and from 2.7 to 3.0% for females. Further, the variance explained by previously and newly identified variants decreased with increasing age in the GERA and UKB cohorts, echoed in the variance explained by the entire genome, which also showed gene-age interaction effects. Finally, we conducted a tissue expression QTL enrichment analysis, which revealed that GWAS BMI-associated variants were enriched in the cerebellum, consistent with prior work in humans and mice.

SOCS2 polymorphisms are not associated with clinical and biochemical phenotypes in acromegalic patients

PURPOSE

Suppressor of cytokine signaling 2 (SOCS2) is a STAT5b-regulated gene and one of its functions is to influence growth and development through negative regulatory effects on GH/IGF-1 pathway. So, we evaluate the potential influence of SOCS2 single nucleotide polymorphisms (SNPs) on clinical and laboratorial characteristics of a large cohort of Brazilian patients with acromegaly.

METHODS

Four SOCS2 SNPs (rs3782415, rs3816997, rs3825199 and rs11107116) were selected and genotyped by real-time PCR using specific Taqman probe assays. A total of 186 patients (116 women, age range 26-88 years) were evaluated.

RESULTS

No association of SOCS2 genotypes was observed with none of the following clinical and laboratorial characteristics: age, sex, body mass index, comorbidities, basal GH, oral glucose tolerance test GH nadir, IGF-I, ULNR-IGF-I.

CONCLUSION

Despite of the key role of SOCS2 in the regulation of GH receptor signaling, we did not find any significant association between SOCS2 polymorphisms and acromegaly.

Integrative genomic analysis implicates limited peripheral adipose storage capacity in the pathogenesis of human insulin resistance

Insulin resistance is a key mediator of obesity-related cardiometabolic disease, yet the mechanisms underlying this link remain obscure. Using an integrative genomic approach, we identify 53 genomic regions associated with insulin resistance phenotypes (higher fasting insulin levels adjusted for BMI, lower HDL cholesterol levels and higher triglyceride levels) and provide evidence that their link with higher cardiometabolic risk is underpinned by an association with lower adipose mass in peripheral compartments. Using these 53 loci, we show a polygenic contribution to familial partial lipodystrophy type 1, a severe form of insulin resistance, and highlight shared molecular mechanisms in common/mild and rare/severe insulin resistance. Population-level genetic analyses combined with experiments in cellular models implicate CCDC92, DNAH10 and L3MBTL3 as previously unrecognized molecules influencing adipocyte differentiation. Our findings support the notion that limited storage capacity of peripheral adipose tissue is an important etiological component in insulin-resistant cardiometabolic disease and highlight genes and mechanisms underpinning this link.

Association Study of GWAS-Derived Loci with Height in Brazilian Children: Importance of MAP3K3, MMP24 and IGF1R Polymorphisms for Height Variation

BACKGROUND/AIM

The single nucleotide polymorphisms (SNPs) rs2282978 (CDK6), rs2425019 (MMP24), rs8081612 (MAP3K3), rs2871865 (IGF1R) and rs3782415 (SOCS2) were among the SNPs most strongly associated with height in a meta-analysis of 47 genome-wide association studies (GWAS) involving 114,223 adults from six ethnic groups. The present study aimed to examine associations between these SNPs and height in Brazilian children.

METHODS

Cross-sectional heights of 1,008 healthy unrelated 4.4- to 9.7-year-old children were evaluated. All genotypes were determined by allele-specific polymerase chain reactions. Height standard deviation scores (SDS) were generated for this population and regressed on allele counts. Linear regressions were performed to estimate the effect of individual SNPs or a polygenic allelic score on height.

RESULTS

The T allele of rs8081612 (MAP3K3), the C allele of rs2871865 (IGF1R) and the G allele of rs2425019 (MMP24) were significantly associated with a 0.091-SDS greater height (95% CI 0.089-0.093, p = 0.001) by polygenic analysis. The mean height SDS difference between children with 2 'tall' alleles and children with 4 'tall' alleles was 0.24 SDS (95% CI 0.05-0.43, p = 0.01). The observed allelic effect is consistent with that found in previous GWAS.

CONCLUSIONS

Polymorphisms in MAP3K3, MMP24 and IGF1R act additively on height in children of an admixed population. These results demonstrate the importance of these loci for children's height.

Association between neonatal resuscitation and a single nucleotide polymorphism rs1835740

AIM

The aim of this work was to test whether three single nucleotide polymorphisms (SNPs) implicated in glutamate homoeostasis or signalling and cellular survival are associated with birth condition.

METHODS

This study is drawn from the Avon Longitudinal Study of Parents and Children. A total of 7611 term infants were genotyped and patient outcome data retrieved from routine medical records. Exposure measures were the presence of one or more minor alleles in one of 3 SNPs (rs2284411, rs2498804, rs1835740). The primary outcome was the need for resuscitation at birth.

RESULTS

For SNP rs1835740, infants homozygous for the minor allele compared to wild type were more likely to need resuscitation (9.2% vs. 7.0%, p = 0.041), while the odds ratio for resuscitation was associated with each increasing minor allele [OR 1.17 (1.01-1.35)]. Population attributable risk fraction was 6.5%. There was no evidence that the other two SNPs investigated were associated with birth condition.

CONCLUSION

We have tested three candidate SNPs to measure any association with birth condition. The study revealed that the rs1835740 was associated with the need for resuscitation and Apgar scores, with a substantial population impact.

Genome-wide association studies in East Asians identify new loci for waist-hip ratio and waist circumference

Sixty genetic loci associated with abdominal obesity, measured by waist circumference (WC) and waist-hip ratio (WHR), have been previously identified, primarily from studies conducted in European-ancestry populations. We conducted a meta-analysis of associations of abdominal obesity with approximately 2.5 million single nucleotide polymorphisms (SNPs) among 53,052 (for WC) and 48,312 (for WHR) individuals of Asian descent, and replicated 33 selected SNPs among 3,762 to 17,110 additional individuals. We identified four novel loci near the EFEMP1, ADAMTSL3 , CNPY2, and GNAS genes that were associated with WC after adjustment for body mass index (BMI); two loci near the NID2 and HLA-DRB5 genes associated with WHR after adjustment for BMI, and three loci near the CEP120, TSC22D2, and SLC22A2 genes associated with WC without adjustment for BMI. Functional enrichment analyses revealed enrichment of corticotropin-releasing hormone signaling, GNRH signaling, and/or CDK5 signaling pathways for those newly-identified loci. Our study provides additional insight on genetic contribution to abdominal obesity.

Assessing the Causal Relationship of Maternal Height on Birth Size and Gestational Age at Birth: A Mendelian Randomization Analysis

BACKGROUND

Observational epidemiological studies indicate that maternal height is associated with gestational age at birth and fetal growth measures (i.e., shorter mothers deliver infants at earlier gestational ages with lower birth weight and birth length). Different mechanisms have been postulated to explain these associations. This study aimed to investigate the casual relationships behind the strong association of maternal height with fetal growth measures (i.e., birth length and birth weight) and gestational age by a Mendelian randomization approach.

METHODS AND FINDINGS

We conducted a Mendelian randomization analysis using phenotype and genome-wide single nucleotide polymorphism (SNP) data of 3,485 mother/infant pairs from birth cohorts collected from three Nordic countries (Finland, Denmark, and Norway). We constructed a genetic score based on 697 SNPs known to be associated with adult height to index maternal height. To avoid confounding due to genetic sharing between mother and infant, we inferred parental transmission of the height-associated SNPs and utilized the haplotype genetic score derived from nontransmitted alleles as a valid genetic instrument for maternal height. In observational analysis, maternal height was significantly associated with birth length (p = 6.31 × 10-9), birth weight (p = 2.19 × 10-15), and gestational age (p = 1.51 × 10-7). Our parental-specific haplotype score association analysis revealed that birth length and birth weight were significantly associated with the maternal transmitted haplotype score as well as the paternal transmitted haplotype score. Their association with the maternal nontransmitted haplotype score was far less significant, indicating a major fetal genetic influence on these fetal growth measures. In contrast, gestational age was significantly associated with the nontransmitted haplotype score (p = 0.0424) and demonstrated a significant (p = 0.0234) causal effect of every 1 cm increase in maternal height resulting in ~0.4 more gestational d. Limitations of this study include potential influences in causal inference by biological pleiotropy, assortative mating, and the nonrandom sampling of study subjects.

CONCLUSIONS

Our results demonstrate that the observed association between maternal height and fetal growth measures (i.e., birth length and birth weight) is mainly defined by fetal genetics. In contrast, the association between maternal height and gestational age is more likely to be causal. In addition, our approach that utilizes the genetic score derived from the nontransmitted maternal haplotype as a genetic instrument is a novel extension to the Mendelian randomization methodology in casual inference between parental phenotype (or exposure) and outcomes in offspring.

Systems biology analysis merging phenotype, metabolomic and genomic data identifies Non-SMC Condensin I Complex, Subunit G (NCAPG) and cellular maintenance processes as major contributors to genetic variability in bovine feed efficiency

Feed efficiency is a paramount factor for livestock economy. Previous studies had indicated a substantial heritability of several feed efficiency traits. In our study, we investigated the genetic background of residual feed intake, a commonly used parameter of feed efficiency, in a cattle resource population generated from crossing dairy and beef cattle. Starting from a whole genome association analysis, we subsequently performed combined phenotype-metabolome-genome analysis taking a systems biology approach by inferring gene networks based on partial correlation and information theory approaches. Our data about biological processes enriched with genes from the feed efficiency network suggest that genetic variation in feed efficiency is driven by genetic modulation of basic processes relevant to general cellular functions. When looking at the predicted upstream regulators from the feed efficiency network, the Tumor Protein P53 (TP53) and Transforming Growth Factor beta 1 (TGFB1) genes stood out regarding significance of overlap and number of target molecules in the data set. These results further support the hypothesis that TP53 is a major upstream regulator for genetic variation of feed efficiency. Furthermore, our data revealed a significant effect of both, the Non-SMC Condensin I Complex, Subunit G (NCAPG) I442M (rs109570900) and the Growth /differentiation factor 8 (GDF8) Q204X (rs110344317) loci, on residual feed intake and feed conversion. For both loci, the growth promoting allele at the onset of puberty was associated with a negative, but favorable effect on residual feed intake. The elevated energy demand for increased growth triggered by the NCAPG 442M allele is obviously not fully compensated for by an increased efficiency in converting feed into body tissue. As a consequence, the individuals carrying the NCAPG 442M allele had an additional demand for energy uptake that is reflected by the association of the allele with increased daily energy intake as observed in our study.

A novel common variant in DCST2 is associated with length in early life and height in adulthood

Common genetic variants have been identified for adult height, but not much is known about the genetics of skeletal growth in early life. To identify common genetic variants that influence fetal skeletal growth, we meta-analyzed 22 genome-wide association studies (Stage 1; N = 28 459). We identified seven independent top single nucleotide polymorphisms (SNPs) (P < 1 × 10(-6)) for birth length, of which three were novel and four were in or near loci known to be associated with adult height (LCORL, PTCH1, GPR126 and HMGA2). The three novel SNPs were followed-up in nine replication studies (Stage 2; N = 11 995), with rs905938 in DC-STAMP domain containing 2 (DCST2) genome-wide significantly associated with birth length in a joint analysis (Stages 1 + 2; β = 0.046, SE = 0.008, P = 2.46 × 10(-8), explained variance = 0.05%). Rs905938 was also associated with infant length (N = 28 228; P = 5.54 × 10(-4)) and adult height (N = 127 513; P = 1.45 × 10(-5)). DCST2 is a DC-STAMP-like protein family member and DC-STAMP is an osteoclast cell-fusion regulator. Polygenic scores based on 180 SNPs previously associated with human adult stature explained 0.13% of variance in birth length. The same SNPs explained 2.95% of the variance of infant length. Of the 180 known adult height loci, 11 were genome-wide significantly associated with infant length (SF3B4, LCORL, SPAG17, C6orf173, PTCH1, GDF5, ZNFX1, HHIP, ACAN, HLA locus and HMGA2). This study highlights that common variation in DCST2 influences variation in early growth and adult height.

Association of ERCC1 rs3212986 & ERCC2 rs13181 polymorphisms with the risk of glioma

Objective:: Several previous studies have reported the role variant of ERCC1 rs3212986 and ERCC2 rs13181 polymorphisms in the risk of glioma, but the results of these studies are inconsistent. Therefore, we aimed to conduct a meta-analysis to investigate the role of ERCC1 rs3212986 and ERCC2 rs13181 on the risk of glioma.

Methods: A comprehensive research was conducted through the databases of Pubmed, EMBASE and the China National Knowledge Infrastructure (CNKI) platforms until June 1, 2014, including 14 eligible case-control studies.

Results: Our meta-analysis found that ERCC1 rs3212986 AA genotype was significantly associated with increased risk of glioma compared with CC genotype, and the pooled OR (95%CI) was 1.29(1.07-1.55). By subgroup analysis, ERCC1 rs3212986 AA genotype was found to be significantly correlated with increased glioma risk in Chinese population (OR=1.37, 95%CI=1.07, 1.55), Similarly, we found that ERCC2 rs13181 GT and TT genotypes were significantly associated with increased risk of glioma in Chinese population, with ORs(95%CI) of 1.47(1.17-1.85) and 1.50(1.02-2.22). But ERCC1 rs3212986 and ERCC2 rs13181 polymorphisms had no significant association with glioma risk in Caucasian populations. By begg’s funnel plot, we found that no publication bias was existed in this meta-analysis.

Conclusion: Our meta-analysis suggested that ERCC1 rs3212986 and ERCC2 rs13181 polymorphism play an important risk factor for brain tumor development in Chinese population, but no association in Caucasian populations.

Association of EFEMP1 gene polymorphisms with the risk of glioma: A hospital-based case-control study in a Chinese Han population

BACKGROUND

EGF-containing fibulin-like extracellular matrix protein1 (EFEMP1) gene was relative with the formation and development of tumors and had an anti-angiogenic function. Recently, many studies investigating the function of EFEMP1 gene, including its roles in prostate cancer and glioma, have been reported. EFEMP1 suppressed glioma growth by modulating EGFR and AKT signaling pathway or promoted growth through the regulation of Notch pathway were identified. However, the susceptibility of EFEMP1and glioma has not been well studied to date. Here, the authors were aimed to investigate whether the single nucleotide polymorphisms (SNPs) of EFEMP1 were associated with glioma susceptibility.

METHODS

The authors genotyped 14 common tagging SNPs of EFEMP1 gene via the Sequenom Mass ARRYiPLEX platform and assessed their association with glioma risk in a hospital-based case-control study in a Chinese Han population (979 cases and 1007 controls).

RESULTS

Four SNPs were significant associated with glioma risk (rs1346787, P=0.004, adjusted OR=1.49; rs3791679, P=0.014, adjusted OR=1.27; rs1346786, P=0.002, adjusted OR=1.41; rs3791675, P=0.011, adjusted OR=1.27). In further stratified analysis, all the significant SNPs except rs1346787 were associated with both low-grade gliomas and glioblastoma (GBM). In haplotype analysis, 4 haplotype blocks were identified and 2 of them were revealed significant associated with glioma, the haplotype "AA" (adjusted OR=1.44, P=0.005) in block 1 and haplotype "GG" (adjusted OR=1.65, P=0.0004) in block 2 had a 44% and 65% increased glioma risk respectively, compared with corresponding non-carriers. The results of haplotype analysis were significantly consistent with the single-locus analysis.

CONCLUSIONS

The authors' results suggested that common genetic variants in EFEMP1 gene were associated with glioma and contributed to glioma susceptibility, which might help to reveal the mechanism of gliomas and provide new insight for the diagnosis and treatment.

The correlation between reading and mathematics ability at age twelve has a substantial genetic component

Dissecting how genetic and environmental influences impact on learning is helpful for maximizing numeracy and literacy. Here we show, using twin and genome-wide analysis, that there is a substantial genetic component to children's ability in reading and mathematics, and estimate that around one half of the observed correlation in these traits is due to shared genetic effects (so-called Generalist Genes). Thus, our results highlight the potential role of the learning environment in contributing to differences in a child's cognitive abilities at age twelve.

Genome-wide association study of sexual maturation in males and females highlights a role for body mass and menarche loci in male puberty

Little is known about genes regulating male puberty. Further, while many identified pubertal timing variants associate with age at menarche, a late manifestation of puberty, and body mass, little is known about these variants' relationship to pubertal initiation or tempo. To address these questions, we performed genome-wide association meta-analysis in over 11 000 European samples with data on early pubertal traits, male genital and female breast development, measured by the Tanner scale. We report the first genome-wide significant locus for male sexual development upstream of myocardin-like 2 (MKL2) (P = 8.9 × 10(-9)), a menarche locus tagging a developmental pathway linking earlier puberty with reduced pubertal growth (P = 4.6 × 10(-5)) and short adult stature (p = 7.5 × 10(-6)) in both males and females. Furthermore, our results indicate that a proportion of menarche loci are important for pubertal initiation in both sexes. Consistent with epidemiological correlations between increased prepubertal body mass and earlier pubertal timing in girls, body mass index (BMI)-increasing alleles correlated with earlier breast development. In boys, some BMI-increasing alleles associated with earlier, and others with delayed, sexual development; these genetic results mimic the controversy in epidemiological studies, some of which show opposing correlations between prepubertal BMI and male puberty. Our results contribute to our understanding of the pubertal initiation program in both sexes and indicate that although mechanisms regulating pubertal onset in males and females may largely be shared, the relationship between body mass and pubertal timing in boys may be complex and requires further genetic studies.

Are genetic risk factors for psychosis also associated with dimension-specific psychotic experiences in adolescence?

Psychosis has been hypothesised to be a continuously distributed quantitative phenotype and disorders such as schizophrenia and bipolar disorder represent its extreme manifestations. Evidence suggests that common genetic variants play an important role in liability to both schizophrenia and bipolar disorder. Here we tested the hypothesis that these common variants would also influence psychotic experiences measured dimensionally in adolescents in the general population. Our aim was to test whether schizophrenia and bipolar disorder polygenic risk scores (PRS), as well as specific single nucleotide polymorphisms (SNPs) previously identified as risk variants for schizophrenia, were associated with adolescent dimension-specific psychotic experiences. Self-reported Paranoia, Hallucinations, Cognitive Disorganisation, Grandiosity, Anhedonia, and Parent-rated Negative Symptoms, as measured by the Specific Psychotic Experiences Questionnaire (SPEQ), were assessed in a community sample of 2,152 16-year-olds. Polygenic risk scores were calculated using estimates of the log of odds ratios from the Psychiatric Genomics Consortium GWAS stage-1 mega-analysis of schizophrenia and bipolar disorder. The polygenic risk analyses yielded no significant associations between schizophrenia and bipolar disorder PRS and the SPEQ measures. The analyses on the 28 individual SNPs previously associated with schizophrenia found that two SNPs in TCF4 returned a significant association with the SPEQ Paranoia dimension, rs17512836 (p-value = 2.57×10⁻⁴) and rs9960767 (p-value = 6.23×10⁻⁴). Replication in an independent sample of 16-year-olds (N = 3,427) assessed using the Psychotic-Like Symptoms Questionnaire (PLIKS-Q), a composite measure of multiple positive psychotic experiences, failed to yield significant results. Future research with PRS derived from larger samples, as well as larger adolescent validation samples, would improve the predictive power to test these hypotheses further. The challenges of relating adult clinical diagnostic constructs such as schizophrenia to adolescent psychotic experiences at a genetic level are discussed.

Association of IL33-IL-1 receptor-like 1 (IL1RL1) pathway polymorphisms with wheezing phenotypes and asthma in childhood

BACKGROUND

Genome-wide association studies identified IL33 and IL-1 receptor-like 1 (IL1RL1)/IL18R1 as asthma susceptibility loci. IL33 and IL1RL1 constitute a single ligand-receptor pathway.

OBJECTIVE

In 2 birth cohorts, the Prevalence and Incidence of Asthma and Mite Allergy (PIAMA) study and Avon Longitudinal Study of Parents and Children (ALSPAC), we analyzed associations of longitudinal wheezing phenotypes and asthma with single nucleotide polymorphisms (SNPs) of 8 genes encoding IL-33, IL1RL1, its coreceptor IL1RAcP, its adaptors myeloid differentiation primary response gene 88 (MyD88) and Toll-IL-11 receptor domain containing adaptor protein (TIRAP), and the downstream IL-1 receptor-associated kinase 1, IL-1 receptor-associated kinase 4, and TNF receptor-associated factor 6 (TRAF6). Furthermore, we investigated whether SNPs in this pathway show replicable evidence of gene-gene interaction.

METHODS

Ninety-four SNPs were investigated in 2007 children in the PIAMA study and 7247 children in ALSPAC. Associations with wheezing phenotypes and asthma at 8 years of age were analyzed in each cohort and subsequently meta-analyzed. Gene-gene interactions were assessed through model-based multifactor dimensionality reduction in the PIAMA study, and gene-gene interactions of 10 SNP pairs were further evaluated.

RESULTS

Intermediate-onset wheeze was associated with SNPs in several genes in the IL33-IL1RL1 pathway after applying multiple testing correction in the meta-analysis: 2 IL33 SNPs (rs4742170 and rs7037276), 1 IL-1 receptor accessory protein (IL1RAP) SNP (rs10513854), and 1 TRAF6 SNP (rs5030411). Late-onset wheeze was associated with 2 IL1RL1 SNPs (rs10208293 and rs13424006), and persistent wheeze was associated with 1 IL33 SNP (rs1342326) and 1 IL1RAP SNP (rs9290936). IL33 and IL1RL1 SNPs were nominally associated with asthma. Three SNP pairs showed interaction for asthma in the PIAMA study but not in ALSPAC.

CONCLUSIONS

IL33-IL1RL1 pathway polymorphisms are associated with asthma and specific wheezing phenotypes; that is, most SNPs are associated with intermediate-onset wheeze, a phenotype closely associated with sensitization. We speculate that IL33-IL1RL1 pathway polymorphisms affect development of wheeze and subsequent asthma through sensitization in early childhood.

Childhood height and birth weight in relation to future prostate cancer risk: a cohort study based on the copenhagen school health records register

BACKGROUND

Adult height has been positively associated with prostate cancer risk. However, the exposure window of importance is currently unknown and assessments of height during earlier growth periods are scarce. In addition, the association between birth weight and prostate cancer remains undetermined. We assessed these relationships in a cohort of the Copenhagen School Health Records Register (CSHRR).

METHODS

The CSHRR comprises 372,636 school children. For boys born between the 1930s and 1969, birth weight and annual childhood heights-measured between ages 7 and 13 years-were analyzed in relation to prostate cancer risk. Cox proportional hazards regression models were used to estimate hazard ratios (HRs) and 95% confidence intervals (CI).

RESULTS

There were 125,211 males for analysis, 2,987 of who were subsequently diagnosed with prostate cancer during 2.57 million person-years of follow-up. Height z-score was significantly associated with prostate cancer risk at all ages (HRs, 1.13 to 1.14). Height at age 13 years was more important than height change (P = 0.024) and height at age 7 years (P = 0.024), when estimates from mutually adjusted models were compared. Adjustment of birth weight did not alter the estimates. Birth weight was not associated with prostate cancer risk.

CONCLUSIONS

The association between childhood height and prostate cancer risk was driven by height at age 13 years.

IMPACT

Our findings implicate late childhood, adolescence, and adulthood growth periods as containing the exposure window(s) of interest that underlies the association between height and prostate cancer. The causal factor may not be singular given the complexity of both human growth and carcinogenesis.

Signatures of natural selection on genetic variants affecting complex human traits

It has recently been hypothesized that polygenic adaptation, resulting in modest allele frequency changes at many loci, could be a major mechanism behind the adaptation of complex phenotypes in human populations. Here we leverage the large number of variants that have been identified through genome-wide association (GWA) studies to comprehensively study signatures of natural selection on genetic variants associated with complex traits. Using population differentiation based methods, such as FST and phylogenetic branch length analyses, we systematically examined nearly 1300 SNPs associated with 38 complex phenotypes. Instead of detecting selection signatures at individual variants, we aimed to identify combined evidence of natural selection by aggregating signals across many trait associated SNPs. Our results have revealed some general features of polygenic selection on complex traits associated variants. First, natural selection acting on standing variants associated with complex traits is a common phenomenon. Second, characteristics of selection for different polygenic traits vary both temporarily and geographically. Third, some studied traits (e.g. height and urate level) could have been the primary targets of selection, as indicated by the significant correlation between the effect sizes and the estimated strength of selection in the trait associated variants; however, for most traits, the allele frequency changes in trait associated variants might have been driven by the selection on other correlated phenotypes. Fourth, the changes in allele frequencies as a result of selection can be highly stochastic, such that, polygenic adaptation may accelerate differentiation in allele frequencies among populations, but generally does not produce predictable directional changes. Fifth, multiple mechanisms (pleiotropy, hitchhiking, etc) may act together to govern the changes in allele frequencies of genetic variants associated with complex traits.

Genetic influences on trajectories of systolic blood pressure across childhood and adolescence

BACKGROUND

Blood pressure (BP) tends to increase across childhood and adolescence, but the genetic influences on rates of BP change are not known. Potentially important genetic influences could include genetic variants identified in genome-wide association studies of adults as being associated with BP, height, and body mass index. Understanding the contribution of these genetic variants to changes in BP across childhood and adolescence could yield understanding into the life course development of cardiovascular risk.

METHODS AND RESULTS

Pooling data from 2 cohorts (the Avon Longitudinal Study of Parents and Children [n=7013] and the Western Australian Pregnancy Cohort [n=1459]), we examined the associations of allelic scores of 29 single-nucleotide polymorphisms (SNPs) for adult BP, 180 height SNPs, and 32 body mass index SNPs, with trajectories of systolic BP (SBP) from 6 to 17 years of age, using linear spline multilevel models. The allelic scores of BP and body mass index SNPs were associated with SBP at 6 years of age (per-allele effect sizes, 0.097 mm Hg [SE, 0.039 mm Hg] and 0.107 mm Hg [SE, 0.037 mm Hg]); associations with age-related changes in SBP between 6 and 17 years of age were of small magnitude and imprecisely estimated. The allelic score of height SNPs was only weakly associated with SBP changes. No sex or cohort differences in genetic effects were observed.

CONCLUSIONS

Allelic scores of BP and body mass index SNPs demonstrated associations with SBP at 6 years of age with a similar magnitude but were not strongly associated with changes in SBP with age between 6 and 17 years. Further work is required to identify variants associated with changes with age in BP.

Body stature growth trajectories during childhood and the development of myopia

PURPOSE

Stature at a particular age can be considered the cumulative result of growth during a number of preceding growth trajectory periods. We investigated whether height and weight growth trajectories from birth to age 10 years were related to refractive error at ages 11 and 15 years, and eye size at age 15 years.

DESIGN

Prospective analysis in a birth cohort.

PARTICIPANTS

Children participating in the Avon Longitudinal Study of Parents and Children (ALSPAC) U.K. birth cohort (minimum N = 2676).

METHODS

Growth trajectories between birth and 10 years were modeled from a series of height and weight measurements (N = 6815). Refractive error was assessed by noncycloplegic autorefraction at ages 11 and 15 years (minimum N = 4737). Axial length (AXL) and radius of corneal curvature were measured with an IOLMaster (Carl Zeiss Meditec, Welwyn Garden City, U.K.) at age 15 years (minimum N = 2676). Growth trajectories and an allelic score for 180 genetic variants associated with adult height were tested for association with refractive error and eye size.

MAIN OUTCOME MEASURES

Noncycloplegic autorefraction at ages 11 and 15 years, and AXL and corneal curvature at age 15 years.

RESULTS

Height growth trajectory during the linear phase between 2.5 and 10 years was negatively associated with refractive error at 11 and 15 years (P<0.001), but explained <0.5% of intersubject variation. Height and weight growth trajectories, especially shortly after birth, were positively associated with AXL and corneal curvature (P<0.001), predicting 1% to 5% of trait variation. Height growth after 2.5 years was not associated with corneal curvature, whereas the association with AXL continued up to 10 years. The height allelic score was associated with corneal curvature (P = 0.03) but not with refractive error or AXL.

CONCLUSIONS

Up to the age of 10 years, shared growth mechanisms contribute to scaling of eye and body size but minimally to the development of myopia.

FINANCIAL DISCLOSURE(S)

The author(s) have no proprietary or commercial interest in any materials discussed in this article.

Development of GMDR-GPU for gene-gene interaction analysis and its application to WTCCC GWAS data for type 2 diabetes

Although genome-wide association studies (GWAS) have identified a significant number of single-nucleotide polymorphisms (SNPs) associated with many complex human traits, the susceptibility loci identified so far can explain only a small fraction of the genetic risk. Among other possible explanations, the lack of a comprehensive examination of gene–gene interaction (G×G) is often considered a source of the missing heritability. Previously, we reported a model-free Generalized Multifactor Dimensionality Reduction (GMDR) approach for detecting G×G in both dichotomous and quantitative phenotypes. However, the computational burden and less efficient implementation of the original programs make them impossible to use for GWAS. In this study, we developed a graphics processing unit (GPU)-based GMDR program (named GWAS-GPU), which is able not only to analyze GWAS data but also to run much faster than the earlier version of the GMDR program. As a demonstration of the program, we used the GMDR-GPU software to analyze a publicly available GWAS dataset on type 2 diabetes (T2D) from the Wellcome Trust Case Control Consortium. Through an exhaustive search of pair-wise interactions and a selected search of three- to five-way interactions conditioned on significant pair-wise results, we identified 24 core SNPs in six genes (FTO: rs9939973, rs9940128, rs9922047, rs1121980, rs9939609, rs9930506; TSPAN8: rs1495377; TCF7L2: rs4074720, rs7901695, rs4506565, rs4132670, rs10787472, rs11196205, rs10885409, rs11196208; L3MBTL3: rs10485400, rs4897366; CELF4: rs2852373, rs608489; RUNX1: rs445984, rs1040328, rs990074, rs2223046, rs2834970) that appear to be important for T2D. Of these core SNPs, 11 in FTO, TSPAN8, and TCF7L2 have been reported to be associated with T2D, obesity, or both, providing an independent replication of previously reported SNPs. Importantly, we identified three new susceptibility genes; i.e., L3MBTL3, CELF4, and RUNX1, for T2D, a finding that warrants further investigation with independent samples.

Describing differences in weight and length growth trajectories between white and Pakistani infants in the UK: analysis of the Born in Bradford birth cohort study using multilevel linear spline models

OBJECTIVE

To describe the growth pattern from birth to 2 years of UK-born white British and Pakistani infants.

DESIGN

Birth cohort.

SETTING

Bradford, UK.

PARTICIPANTS

314 white British boys, 383 Pakistani boys, 328 white British girls and 409 Pakistani girls.

MAIN OUTCOME MEASURES

Weight and length trajectories based on repeat measurements from birth to 2 years.

RESULTS

Linear spline multilevel models for weight and length with knot points at 4 and 9 months fitted the data well. At birth Pakistani boys were 210 g lighter (95% CI -290 to -120) and 0.5 cm shorter (-1.04 to 0.02) and Pakistani girls were 180 g lighter (-260 to -100) and 0.5 cm shorter (-0.91 to -0.03) than white British boys and girls, respectively. Pakistani infants gained length faster than white British infants between 0 and 4 months (+0.3 cm/month (0.1 to 0.5) for boys and +0.4 cm/month (0.2 to 0.6) for girls) and gained more weight per month between 9 and 24 months (+10 g/month (0 to 30) for boys and +30 g/month (20 to 40) for girls). Adjustment for maternal height attenuated ethnic differences in weight and length at birth, but not in postnatal growth. Adjustment for other confounders did not explain differences in any outcomes.

CONCLUSIONS

Pakistani infants were lighter and had shorter predicted mean length at birth than white British infants, but gained weight and length quicker in infancy. By age 2 years both ethnic groups had similar weight, but Pakistani infants were on average taller than white British infants.

Using height association studies to gain insights into human idiopathic short and syndromic stature phenotypes

Variation in adult height is not the most clinically relevant human quantitative trait, yet its study provides the foundation of many quantitative genetics theories and important statistical concepts (e.g. regression). Even today, the analysis of adult height by genome-wide association studies (GWAS) continues to significantly impact human genetics: these studies have led to the discovery of >200 loci associated with variation in adult height and have highlighted the very polygenic nature of human continuous traits. In this brief review, I discuss and provide examples on how such genetic associations, identified in individuals of normal height, could help understand the complex genetics behind such phenotypes as idiopathic short stature (ISS) or extreme/syndromic height phenotypes of unknown cause.

Coordinated genetic scaling of the human eye: shared determination of axial eye length and corneal curvature

PURPOSE

To examine the extent to which the two major determinants of refractive error, corneal curvature and axial length, are scaled relative to one another by shared genetic variants, along with their relationship to the genetic scaling of height.

METHODS

Corneal curvature, axial length, and height were measured in unrelated 14- to 17-year-old white European participants of the Avon Longitudinal Study of Parents and Children (ALSPAC; n = 1915) and in unrelated 40- to 80-year-old participants of the Singapore Chinese Eye Study (SCES; n = 1642). Univariate and bivariate heritability analyses were performed with methods that avoid confounding by common family environment, using information solely from genome-wide high-density genotypes.

RESULTS

IN ALSPAC SUBJECTS, AXIAL LENGTH, CORNEAL CURVATURE, AND HEIGHT HAD SIMILAR LOWER-BOUND HERITABILITY ESTIMATES: axial length, h(2) = 0.46 (SE = 0.16, P = 0.002); corneal curvature, h(2) = 0.42 (SE = 0.16, P = 0.004); height, h(2) = 0.48 (SE = 0.17, P = 0.002). The corresponding estimates in the SCES were 0.79 (SE = 0.18, P < 0.001), 0.35 (SE = 0.20, P = 0.036), and 0.31 (SE = 0.20, P = 0.061), respectively. The genetic correlation between corneal curvature and axial length was 0.69 (SE = 0.17, P = 0.019) for ALSPAC participants and 0.64 (SE = 0.22, P = 0.003) for SCES participants. In the subset of 1478 emmetropic ALSPAC individuals, the genetic correlation was 0.85 (SE = 0.12, P = 0.008).

CONCLUSIONS

These results imply that coordinated scaling of ocular component dimensions is largely achieved by hundreds to thousands of common genetic variants, each with a small pleiotropic effect. Furthermore, genome-wide association studies (GWAS) for either axial length or corneal curvature are likely to identify variants controlling overall eye size when using discovery cohorts dominated by emmetropes, but trait-specific variants in discovery cohorts dominated by ametropes.

A genome-wide association study for corneal curvature identifies the platelet-derived growth factor receptor α gene as a quantitative trait locus for eye size in white Europeans

PURPOSE

Corneal curvature is a key determinant of the refractive power of the eye. Variants in two genes, FKBP12-rapamycin complex-associated protein 1 (FRAP1) on chromosome 1p36.2 and platelet-derived growth factor receptor alpha (PDGFRA) on chromosome 4q12, have shown genome-wide significant association with normal variation in corneal curvature in a study of subjects of Asian origin. Variants at the PDGFRA locus have also shown genome-wide significant association with corneal astigmatism. Whether these variants influence other ocular parameters such as axial length has yet to be reported. We performed a genome-wide association study for corneal curvature in white European subjects from a population-based birth cohort, with the aim of replicating and extending the above findings.

METHODS

White European children participating in the Avon Longitudinal Study of Parents and Children (ALSPAC) birth cohort were examined at age about 15.5 years (95% confidence interval=15.45 to 15.48 years). Radius of corneal curvature and axial eye length were measured with an IOLmaster. DNA samples were genotyped with Illumina HumanHap550 arrays and untyped variants imputed using MACH, with CEU individuals from HapMap release 22, Phase II NCBI B36, Single Nucleotide Polymorphism database 126 as the reference panel. Association between corneal curvature and single nucleotide polymorphism (SNP) genotype was tested, genome-wide, using mach2qtl, with sex as a covariate (n=2023; 46.6% male).

RESULTS

The variant exhibiting the strongest evidence for association with corneal curvature (rs6554163; p=2.8×10(-6)) was located in the same linkage disequilibrium block as the previously discovered PDGFRA variants. Meta-analysis of the current and prior findings enhanced the evidence for association (rs17084051, p=4.5×10(-14)). rs6554163 genotype predicted 1.0% of variation in corneal curvature. In addition, these PDGFRA variants were associated with axial eye length, predicting 0.6% of the normal trait variation (p=5.3×10(-4)). Each copy of the minor allele of variants at the locus also increased the risk of corneal astigmatism in this white European cohort (odds ratio [OR]=1.24, 95% confidence interval=1.07-1.45; p=0.006).

CONCLUSION

As in Asians, variants at the PDGFRA locus influence corneal curvature (and corneal astigmatism). However, rather than affecting corneal curvature in isolation, this locus influences the size of the eye while maintaining its scaling.

Identification of a short region on chromosome 6 affecting direct calving ease in Piedmontese cattle breed

Calving in cattle is affected by calf morphology and by dam characteristics. It is described by two different traits: maternal calving ease, which is the ability to generate dams with good physiological predisposition to calving, and direct calving ease, which is the ability to generate calves that are easily born. The aim of this study was to identify regions of cattle genome harboring genes possibly affecting direct calving ease in the Piedmontese cattle breed. A population of 323 bulls scored for direct calving ease (EBV) was analyzed by a medium-density SNP marker panel (54,001 SNPs) to perform a genome-wide scan. The strongest signal was detected on chromosome 6 between 37.8 and 38.7 Mb where 13 SNPs associated to direct calving ease were found. Three genes are located in this region: LAP3, encoding for a leucine aminopeptidase involved in the oxytocin hydrolysis; NCAPG, encoding for a non-SMC condensin I complex, which has been associated in cattle with fetal growth and carcass size; and LCORL, which has been associated to height in humans and cattle. To further confirm the results of the genome-wide scan we genotyped additional SNPs within these genes and analyzed their association with direct calving ease. The results of this additional analysis fully confirmed the findings of the GWAS and particularly indicated LAP3 as the most probable gene involved. Linkage Disequilibrium (LD) analysis showed high correlation between SNPs located within LAP3 and LCORL indicating a possible selection signature due either to increased fitness or breeders' selection for the trait.

Pourcain B, Hofman A, Jaddoe VWV, Smith GD, Lorentzon M, Gautvik KM, Uitterlinden AG, Brommage R, Ohlsson C, Tobias JH, Rivadeneira F. Meta-analysis of genome-wide scans for total body BMD in children and adults reveals allelic heterogeneity and age-specific effects at the WNT16 locus

To identify genetic loci influencing bone accrual, we performed a genome-wide association scan for total-body bone mineral density (TB-BMD) variation in 2,660 children of different ethnicities. We discovered variants in 7q31.31 associated with BMD measurements, with the lowest P = 4.1 × 10(-11) observed for rs917727 with minor allele frequency of 0.37. We sought replication for all SNPs located ± 500 kb from rs917727 in 11,052 additional individuals from five independent studies including children and adults, together with de novo genotyping of rs3801387 (in perfect linkage disequilibrium (LD) with rs917727) in 1,014 mothers of children from the discovery cohort. The top signal mapping in the surroundings of WNT16 was replicated across studies with a meta-analysis P = 2.6 × 10(-31) and an effect size explaining between 0.6%-1.8% of TB-BMD variance. Conditional analyses on this signal revealed a secondary signal for total body BMD (P = 1.42 × 10(-10)) for rs4609139 and mapping to C7orf58. We also examined the genomic region for association with skull BMD to test if the associations were independent of skeletal loading. We identified two signals influencing skull BMD variation, including rs917727 (P = 1.9 × 10(-16)) and rs7801723 (P = 8.9 × 10(-28)), also mapping to C7orf58 (r(2) = 0.50 with rs4609139). Wnt16 knockout (KO) mice with reduced total body BMD and gene expression profiles in human bone biopsies support a role of C7orf58 and WNT16 on the BMD phenotypes observed at the human population level. In summary, we detected two independent signals influencing total body and skull BMD variation in children and adults, thus demonstrating the presence of allelic heterogeneity at the WNT16 locus. One of the skull BMD signals mapping to C7orf58 is mostly driven by children, suggesting temporal determination on peak bone mass acquisition. Our life-course approach postulates that these genetic effects influencing peak bone mass accrual may impact the risk of osteoporosis later in life.

Quantitative genomics of 30 complex phenotypes in Wagyu x Angus F₁ progeny

In the present study, a total of 91 genes involved in various pathways were investigated for their associations with six carcass traits and twenty-four fatty acid composition phenotypes in a Wagyu×Angus reference population, including 43 Wagyu bulls and their potential 791 F₁ progeny. Of the 182 SNPs evaluated, 102 SNPs that were in Hardy-Weinberg equilibrium with minor allele frequencies (MAF>0.15) were selected for parentage assignment and association studies with these quantitative traits. The parentage assignment revealed that 40 of 43 Wagyu sires produced over 96.71% of the calves in the population. Linkage disequilibrium analysis identified 75 of 102 SNPs derived from 54 genes as tagged SNPs. After Bonferroni correction, single-marker analysis revealed a total of 113 significant associations between 44 genes and 29 phenotypes (adjusted P<0.05). Multiple-marker analysis confirmed single-gene associations for 10 traits, but revealed two-gene networks for 9 traits and three-gene networks for 8 traits. Particularly, we observed that TNF (tumor necrosis factor) gene is significantly associated with both beef marbling score (P=0.0016) and palmitic acid (C16:0) (P=0.0043), RCAN1 (regulator of calcineurin 1) with rib-eye area (P=0.0103), ASB3 (ankyrin repeat and SOCS box-containing 3) with backfat (P=0.0392), ABCA1 (ATP-binding cassette A1) with both palmitic acid (C16:0) (P=0.0025) and oleic acid (C18:1n9) (P=0.0114), SLC27A1(solute carrier family 27 A1) with oleic acid (C18:1n9) (P=0.0155), CRH (corticotropin releasing hormone) with both linolenic acid (OMEGA-3) (P=0.0200) and OMEGA 6:3 RATIO (P=0.0054), SLC27A2 (solute carrier family 27 A2) with both linoleic acid (OMEGA-6) (P=0.0121) and FAT (P=0.0333), GNG3 (guanine nucleotide binding protein gamma 3 with desaturase 9 (P=0.0115), and EFEMP1 (EGF containing fibulin-like extracellular matrix protein 1), PLTP (phospholipid transfer protein) and DSEL (dermatan sulfate epimerase-like) with conjugated linoleic acid (P=0.0042-0.0044), respectively, in the Wagyu x Angus F₁ population. In addition, we observed an interesting phenomenon that crossbreeding of different breeds might change gene actions to dominant and overdominant modes, thus explaining the origin of heterosis. The present study confirmed that these important families or pathway-based genes are useful targets for improving meat quality traits and healthful beef products in cattle.