Sandra Lai

Background: Growth differentiation factor 15 (GDF‐15) levels are emerging as a candidate biomarker of aging. The present study aimed to: (1) characterize the association of GDF‐15 with the continuum of arterial stiffening, assessed as carotid‐femoral pulse wave velocity, as age increases; (2) determine the predictive role of serum GDF‐15 levels on mortality; and (3) identify genetic determinants of serum GDF‐15 levels. Methods and Results: Serum levels of GDF‐15 and established cardiovascular risk factors, including pulse wave velocity, were assessed in a large (4736 individual) Sardinian population. Serum levels of GDF‐15, which can be reliably measured repeatedly over time, increase with age; are associated with a stiffer aorta; “mediate” a large proportion of the age‐associated increase in arterial stiffness; pose risks because of their association with greater mortality; and are significantly associated with the variant rs11549407, which causes thalassemia major in homozygosity. Conclusions: Because of its consistent ability to predict functional and clinical outcomes, including all‐cause mortality, we conclude that GDF‐15 serum levels serve as a robust biomarker for the continuum from health to the emergence of clinical disease during aging and, subsequently, to the likelihood of mortality.

Fresh blood immunophenotyping by flow cytometry, based on the reliable simultaneous detection of several markers in a cell, is the method of choice to study the circulating human immune system. Especially in large and multicenter studies, high sample quality is difficult to achieve, and adequate collection and storage of samples with fine-tuned whole blood cryopreservation is mandatory. Here, we compared the quality of immunophenotypic data obtained from fresh blood with those obtained after five cryopreservation methods by quantifying the levels of 41 immune cell populations. They comprised B and T lymphocyte subsets and their maturation stages, as well as monocytes and granulocytes. Three methods used fixative solutions and two other methods used dimethyl sulfoxide solutions to preserve cell viability. The fixative methods prevented detection of markers critical for identification of B and T cell subsets, including CD27, CXCR3, and CCR6. The other two methods permitted reliable discrimination of most immune-cell populations in thawed samples, though some cell frequencies varied compared to the corresponding fresh sample. Of those two methods, the one preserving blood in media containing dimethyl sulfoxide produced results that were most similar to those with fresh samples.

To investigate the genetic regulation of platelet (PLT) levels we carried out a whole-genome association analysis in 6,528 Sardinians from the general population of the Lanusei valley. We found 6 variants significantly influencing PLT levels, including a novel rare missense mutation (p.Pro27Ser) in the GP1BB protein that is associated with PLT reduction (P=1.17×10 ⁻¹⁶ ). This mutation is rare in the SardiNIA population cohort (frequency of 0.45%), even rarer in the rest of the Sardinian island (frequency of 0.16%), and not reported elsewhere. Notably, GP1BB is involved in Bernard-Soulier syndrome (BSS), a rare autosomal recessive bleeding disorder caused by a defect in the platelet GPIb-IX-V protein complex. Consistently, the 57 identified individuals heterozygous for the p.P27S mutation showed mild thrombocytopenia, morphologically enlarged platelets (P=2.13×10 ⁻¹⁰ ), and reduced expression of two GPIb-IX-V-complex components: GPIbα (−26.51%, P=3.66×10 ⁻⁸ ) and GPIX (−24.69%, P=2.66×10 ⁻⁶ ). Molecular modeling infers a corresponding reduction in the stability of GP1BB. These observations predict that in homozygosity as well as in individuals carrying specific compound heterozygous configurations, this variant likely causes BSS.

Genomewide association studies of autoimmune diseases have mapped hundreds of susceptibility regions in the genome. However, only for a few association signals has the causal gene been identified, and for even fewer have the causal variant and underlying mechanism been defined. Coincident associations of DNA variants affecting both the risk of autoimmune disease and quantitative immune variables provide an informative route to explore disease mechanisms and drug-targetable pathways. Using case–control samples from Sardinia, Italy, we performed a genomewide association study in multiple sclerosis followed by TNFSF13B locus–specific association testing in systemic lupus erythematosus (SLE). Extensive phenotyping of quantitative immune variables, sequence-based fine mapping, cross-population and cross-phenotype analyses, and gene-expression studies were used to identify the causal variant and elucidate its mechanism of action. Signatures of positive selection were also investigated. A variant in TNFSF13B, encoding the cytokine and drug target B-cell activating factor (BAFF), was associated with multiple sclerosis as well as SLE. The disease-risk allele was also associated with up-regulated humoral immunity through increased levels of soluble BAFF, B lymphocytes, and immunoglobulins. The causal variant was identified: an insertion–deletion variant, GCTGT→A (in which A is the risk allele), yielded a shorter transcript that escaped microRNA inhibition and increased production of soluble BAFF, which in turn up-regulated humoral immunity. Population genetic signatures indicated that this autoimmunity variant has been evolutionarily advantageous, most likely by augmenting resistance to malaria. A TNFSF13B variant was associated with multiple sclerosis and SLE, and its effects were clarified at the population, cellular, and molecular levels. (Funded by the Italian Foundation for Multiple Sclerosis and others.)

Melinda Mills, Nicola Barban, Harold Snieder, Marcel den Hoed and colleagues perform a meta-analysis of data from over 300,000 individuals for age at first birth and number of children ever born. They identify 12 significant loci that associate with these traits, providing insights into the genetic basis of human reproductive behavior. The genetic architecture of human reproductive behavior—age at first birth (AFB) and number of children ever born (NEB)—has a strong relationship with fitness, human development, infertility and risk of neuropsychiatric disorders. However, very few genetic loci have been identified, and the underlying mechanisms of AFB and NEB are poorly understood. We report a large genome-wide association study of both sexes including 251,151 individuals for AFB and 343,072 individuals for NEB. We identified 12 independent loci that are significantly associated with AFB and/or NEB in a SNP-based genome-wide association study and 4 additional loci associated in a gene-based effort. These loci harbor genes that are likely to have a role, either directly or by affecting non-local gene expression, in human reproduction and infertility, thereby increasing understanding of these complex traits.

Despite progress in identifying genes associated with breast cancer, many more risk loci exist. Genome-wide association analyses in genetically-homogeneous populations, such as that of Sardinia (Italy), could represent an additional approach to detect low penetrance alleles. We performed a genome-wide association study comparing 1431 Sardinian patients with non-familial, BRCA1/2-mutation-negative breast cancer to 2171 healthy Sardinian blood donors. DNA was genotyped using GeneChip Human Mapping 500 K Arrays or Genome-Wide Human SNP Arrays 6.0. To increase genomic coverage, genotypes of additional SNPs were imputed using data from HapMap Phase II. After quality control filtering of genotype data, 1367 cases (9 men) and 1658 controls (1156 men) were analyzed on a total of 2,067,645 SNPs. Overall, 33 genomic regions (67 candidate SNPs) were associated with breast cancer risk at the p < 10⁻⁶ level. Twenty of these regions contained defined genes, including one already associated with breast cancer risk: TOX3. With a lower threshold for preliminary significance to p < 10⁻⁵, we identified 11 additional SNPs in FGFR2, a well-established breast cancer-associated gene. Ten candidate SNPs were selected, excluding those already associated with breast cancer, for technical validation as well as replication in 1668 samples from the same population. Only SNP rs345299, located in intron 1 of VAV3, remained suggestively associated (p-value, 1.16x10⁻⁵), but it did not associate with breast cancer risk in pooled data from two large, mixed-population cohorts. This study indicated the role of TOX3 and FGFR2 as breast cancer susceptibility genes in BRCA1/2-wild-type breast cancer patients from Sardinian population.

The utility of genotype imputation in genome-wide association studies is increasing as progressively larger reference panels are improved and expanded through whole-genome sequencing. Developing general guidelines for optimally cost-effective imputation, however, requires evaluation of performance issues that include the relative utility of study-specific compared with general/multipopulation reference panels; genotyping with various array scaffolds; effects of different ethnic backgrounds; and assessment of ranges of allele frequencies. Here we compared the effectiveness of study-specific reference panels to the commonly used 1000 Genomes Project (1000G) reference panels in the isolated Sardinian population and in cohorts of European ancestry including samples from Minnesota (USA). We also examined different combinations of genome-wide and custom arrays for baseline genotypes. In Sardinians, the study-specific reference panel provided better coverage and genotype imputation accuracy than the 1000G panels and other large European panels. In fact, even gene-centered custom arrays (interrogating ~200 000 variants) provided highly informative content across the entire genome. Gain in accuracy was also observed for Minnesotans using the study-specific reference panel, although the increase was smaller than in Sardinians, especially for rare variants. Notably, a combined panel including both study-specific and 1000G reference panels improved imputation accuracy only in the Minnesota sample, and only at rare sites. Finally, we found that when imputation is performed with a study-specific reference panel, cutoffs different from the standard thresholds of MACH-Rsq and IMPUTE-INFO metrics should be used to efficiently filter badly imputed rare variants. This study thus provides general guidelines for researchers planning large-scale genetic studies.

Here 106 genomic loci associated with age at menarche, a marker of puberty timing in females, are identified; these loci show enrichment for genes involved in nuclear hormone receptor function, body mass index, and rare disorders of puberty, and for genes located in imprinted regions, with parent-of-origin specific effects at several loci. The age at which females first experience menstruation, called menarche, is a heritable trait associated with risks for obesity, type 2 diabetes, cardiovascular disease, breast cancer and general mortality. This large-scale genome-wide association study identifies 123 signals at 106 genomic loci associated with age at menarche. New findings include parent-of-origin-specific allelic associations (both maternally and paternally driven) at three imprinted loci and the implication of retinoic acid and GABAB receptor II signalling and lysine-specific histone demethylation. These data bring new insights into the genetic architecture of puberty timing and suggest a model involving thousands of genetic variants. Age at menarche is a marker of timing of puberty in females. It varies widely between individuals, is a heritable trait and is associated with risks for obesity, type 2 diabetes, cardiovascular disease, breast cancer and all-cause mortality1. Studies of rare human disorders of puberty and animal models point to a complex hypothalamic-pituitary-hormonal regulation2,3, but the mechanisms that determine pubertal timing and underlie its links to disease risk remain unclear. Here, using genome-wide and custom-genotyping arrays in up to 182,416 women of European descent from 57 studies, we found robust evidence (P < 5 × 10−8) for 123 signals at 106 genomic loci associated with age at menarche. Many loci were associated with other pubertal traits in both sexes, and there was substantial overlap with genes implicated in body mass index and various diseases, including rare disorders of puberty. Menarche signals were enriched in imprinted regions, with three loci (DLK1-WDR25, MKRN3-MAGEL2 and KCNK9) demonstrating parent-of-origin-specific associations concordant with known parental expression patterns. Pathway analyses implicated nuclear hormone receptors, particularly retinoic acid and γ-aminobutyric acid-B2 receptor signalling, among novel mechanisms that regulate pubertal timing in humans. Our findings suggest a genetic architecture involving at least hundreds of common variants in the coordinated timing of the pubertal transition.