Sandra Lai

Genome-wide association study shows BCL11A associated with persistent fetal hemoglobin and amelioration of the phenotype of -thalassemia

β-Thalassemia and sickle cell disease both display a great deal of phenotypic heterogeneity, despite being generally thought of as simple Mendelian diseases. The reasons for this are not well understood, although the level of fetal hemoglobin (HbF) is one well characterized ameliorating factor in both of these conditions. To better understand the genetic basis of this heterogeneity, we carried out genome-wide scans with 362,129 common SNPs on 4,305 Sardinians to look for genetic linkage and association with HbF levels, as well as other red blood cell-related traits. Among major variants affecting HbF levels, SNP rs11886868 in the BCL11A gene was strongly associated with this trait (P < 10−35). The C allele frequency was significantly higher in Sardinian individuals with elevated HbF levels, detected by screening for β-thalassemia, and patients with attenuated forms of β-thalassemia vs. those with thalassemia major. We also show that the same BCL11A variant is strongly associated with HbF levels in a large cohort of sickle cell patients. These results indicate that BCL11A variants, by modulating HbF levels, act as an important ameliorating factor of the β-thalassemia phenotype, and it is likely they could help ameliorate other hemoglobin disorders. We expect our findings will help to characterize the molecular mechanisms of fetal globin regulation and could eventually contribute to the development of new therapeutic approaches for β-thalassemia and sickle cell anemia.

Common variants in the GDF5-UQCC region are associated with variation in human height

Identifying genetic variants that influence human height will advance our understanding of skeletal growth and development. Several rare genetic variants have been convincingly and reproducibly associated with height in mendelian syndromes, and common variants in the transcription factor gene HMGA2 are associated with variation in height in the general population. Here we report genome-wide association analyses, using genotyped and imputed markers, of 6,669 individuals from Finland and Sardinia, and follow-up analyses in an additional 28,801 individuals. We show that common variants in the osteoarthritis-associated locus GDF5-UQCC contribute to variation in height with an estimated additive effect of 0.44 cm (overall P < 10−15). Our results indicate that there may be a link between the genetic basis of height and osteoarthritis, potentially mediated through alterations in bone growth and development.

The GLUT9 gene is associated with serum uric acid levels in Sardinia and Chianti cohorts.

High serum uric acid levels elevate pro-inflammatory–state gout crystal arthropathy and place individuals at high risk for cardiovascular morbidity and mortality. Genome-wide scans in the genetically isolated Sardinian population identified variants associated with serum uric acid levels as a quantitative trait. They mapped within GLUT9, a Chromosome 4 glucose transporter gene predominantly expressed in liver and kidney. SNP rs6855911 showed the strongest association (p = 1.84 × 10−16), along with eight others (p = 7.75 × 10−16 to 6.05 × 10−11). Individuals homozygous for the rare allele of rs6855911 (minor allele frequency = 0.26) had 0.6 mg/dl less uric acid than those homozygous for the common allele; the results were replicated in an unrelated cohort from Tuscany. Our results suggest that polymorphisms in GLUT9 could affect glucose metabolism and uric acid synthesis and/or renal reabsorption, influencing serum uric acid levels over a wide range of values.

Genome-wide Association Study of Vitamin B6, Vitamin B12, Folate, and Homocysteine Blood Concentrations

The B vitamins are components of one-carbon metabolism (OCM) that contribute to DNA synthesis and methylation. Homocysteine, a by-product of OCM, has been associated with coronary heart disease, stroke and neurological disease. To investigate genetic factors that affect circulating vitamin B6, vitamin B12, folate and homocysteine, a genome-wide association analysis was conducted in the InCHIANTI (N = 1175), SardiNIA (N = 1115), and BLSA (N = 640) studies. The top loci were replicated in an independent sample of 687 participants in the Progetto Nutrizione study. Polymorphisms in the ALPL gene (rs4654748, p = 8.30 x 10(-18)) were associated with vitamin B6 and FUT2 (rs602662, [corrected] p = 2.83 x 10(-20)) with vitamin B12 serum levels. The association of MTHFR, a gene consistently associated with homocysteine, was confirmed in this meta-analysis. The ALPL gene likely influences the catabolism of vitamin B6 while FUT2 interferes with absorption of vitamin B12. These findings highlight mechanisms that affect vitamin B6, vitamin B12 and homocysteine serum levels.

Variants within the immunoregulatory CBLB gene are associated with Multiple Sclerosis

A genome-wide association scan of ∼6.6 million genotyped or imputed variants in 882 Sardinian individuals with multiple sclerosis (cases) and 872 controls suggested association of CBLB gene variants with disease, which was confirmed in 1,775 cases and 2,005 controls (rs9657904, overall P = 1.60 × 10−10, OR = 1.40). CBLB encodes a negative regulator of adaptive immune responses, and mice lacking the ortholog are prone to experimental autoimmune encephalomyelitis, the animal model of multiple sclerosis.

Common variants in the SLCO1B3 locus are associated with bilirubin levels and unconjugated hyperbilirubinemia

Bilirubin, resulting largely from the turnover of hemoglobin, is found in the plasma in two main forms: unconjugated or conjugated with glucuronic acid. Unconjugated bilirubin is transported into hepatocytes. There, it is glucuronidated by UGT1A1 and secreted into the bile canaliculi. We report a genome wide association scan in 4300 Sardinian individuals for total serum bilirubin levels. In addition to the two known loci previously involved in the regulation of bilirubin levels, UGT1A1 (P = 6.2 x 10(-62)) and G6PD (P = 2.5 x 10(-8)), we observed a strong association on chromosome 12 within the SLCO1B3 gene (P = 3.9 x 10(-9)). Our findings were replicated in an independent sample of 1860 Sardinians and in 832 subjects from the Old Order Amish (combined P < 5 x 10(-14)). We also show that SLC01B3 variants contribute to idiopathic mild unconjugated hyperbilirubinemia. Thus, SLC01B3 appears to be involved in the regulation of serum bilirubin levels in healthy individuals and in some bilirubin-related disorders that are only partially explained by other known gene variants.

Phosphodiesterase 8B Gene Variants Are Associated with Serum TSH Levels and Thyroid Function

Thyroid-stimulating hormone (TSH) controls thyroid growth and hormone secretion through binding to its G protein-coupled receptor (TSHR) and production of cyclic AMP (cAMP). Serum TSH is a sensitive indicator of thyroid function, and overt abnormalities in thyroid function lead to common endocrine disorders affecting approximately 10% of individuals over a life span. By genotyping 362,129 SNPs in 4,300 Sardinians, we identified a strong association (p = 1.3 x 10(-11)) between alleles of rs4704397 and circulating TSH levels; each additional copy of the minor A allele was associated with an increase of 0.13 muIU/ml in TSH. The single-nucleotide polymorphism (SNP) is located in intron 1 of PDE8B, encoding a high-affinity cAMP-specific phosphodiesterase. The association was replicated in 4,158 individuals, including additional Sardinians and two genetically distant cohorts from Tuscany and the Old Order Amish (overall p value = 1.9 x 10(-20)). In addition to association of TSH levels with SNPs in PDE8B, our genome scan provided evidence for association with PDE10A and several biologically interesting candidates in a focused analysis of 24 genes. In particular, we found evidence for association of TSH levels with SNPs in the THRB (rs1505287, p = 7.3 x 10(-5)), GNAQ (rs10512065, p = 2.0 x 10(-4)), TG (rs2252696, p = 2.2 x 10(-3)), POU1F1 (rs1976324, p = 3.9 x 10(-3)), PDE4D (rs27178, p = 8.3 x 10(-3)), and TSHR (rs4903957, p = 8.6 x 10(-3)) loci. Overall, the results suggest a primary effect of PDE8B variants on cAMP levels in the thyroid. This would affect production of T4 and T3 and feedback to alter TSH release by the pituitary. PDE8B may thus provide a candidate target for the treatment of thyroid dysfunction.

Genome-wide association scan of trait depression

BACKGROUND Independent of temporal circumstances, some individuals have greater susceptibility to depressive affects, such as feelings of guilt, sadness, hopelessness, and loneliness. Identifying the genetic variants that contribute to these individual differences can point to biological pathways etiologically involved in psychiatric disorders. METHODS Genome-wide association scans for the depression scale of the Revised NEO Personality Inventory in community-based samples from a genetically homogeneous area of Sardinia, Italy (n = 3972) and from the Baltimore Longitudinal Study of Aging in the United States (n = 839). RESULTS Meta-analytic results for genotyped or imputed single nucleotide polymorphisms indicate that the strongest association signals for trait depression were found in RORA (rs12912233; p = 6 × 10⁻⁷), a gene involved in circadian rhythm. A plausible biological association was also found with single nucleotide polymorphisms within GRM8 (rs17864092; p = 5 × 10⁻⁶), a metabotropic receptor for glutamate, a major excitatory neurotransmitter in the central nervous system. CONCLUSIONS These findings suggest shared genetic basis underlying the continuum from personality traits to psychopathology.

Low-Pass DNA Sequencing of 1200 Sardinians Reconstructs European Y-Chromosome Phylogeny

Examining Y The evolution of human populations has long been studied with unique sequences from the nonrecombining, male-specific Y chromosome (see the Perspective by Cann). Poznik et al. (p. 562) examined 9.9 Mb of the Y chromosome from 69 men from nine globally divergent populations—identifying population and individual specific sequence variants that elucidate the evolution of the Y chromosome. Sequencing of maternally inherited mitochondrial DNA allowed comparison between the relative rates of evolution, which suggested that the coalescence, or origin, of the human Y chromosome and mitochondria both occurred approximately 120 thousand years ago. Francalacci et al. (p. 565) investigated the sequence divergence of 1204 Y chromosomes that were sampled within the isolated and genetically informative Sardinian population. The sequence analyses, along with archaeological records, were used to calibrate and increase the resolution of the human phylogenetic tree. Local human demographic history is inferred from in-depth DNA sequence analysis of Sardinian mens Y chromosomes. [Also see Perspective by Cann] Genetic variation within the male-specific portion of the Y chromosome (MSY) can clarify the origins of contemporary populations, but previous studies were hampered by partial genetic information. Population sequencing of 1204 Sardinian males identified 11,763 MSY single-nucleotide polymorphisms, 6751 of which have not previously been observed. We constructed a MSY phylogenetic tree containing all main haplogroups found in Europe, along with many Sardinian-specific lineage clusters within each haplogroup. The tree was calibrated with archaeological data from the initial expansion of the Sardinian population ~7700 years ago. The ages of nodes highlight different genetic strata in Sardinia and reveal the presumptive timing of coalescence with other human populations. We calculate a putative age for coalescence of ~180,000 to 200,000 years ago, which is consistent with previous mitochondrial DNA–based estimates.

COL4A1 is associated with arterial stiffness by genome-wide association scan.

Background—Pulse wave velocity (PWV), a noninvasive index of central arterial stiffness, is a potent predictor of cardiovascular mortality and morbidity. Heritability and linkage studies have pointed toward a genetic component affecting PWV. We conducted a genome-wide association study to identify single-nucleotide polymorphisms (SNPs) associated with PWV. Methods and Results—The study cohort included participants from the SardiNIA study for whom PWV measures were available. Genotyping was performed in 4221 individuals, using either the Affymetrix 500K or the Affymetrix 10K mapping array sets (with imputation of the missing genotypes). Associations with PWV were evaluated using an additive genetic model that included age, age2, and sex as covariates. The findings were tested for replication in an independent internal Sardinian cohort of 1828 individuals, using a custom chip designed to include the top 43 nonredundant SNPs associated with PWV. Of the loci that were tested for association with PWV, the nonsynonymous SNP rs3742207 in the COL4A1 gene on chromosome 13 and SNP rs1495448 in the MAGI1 gene on chromosome 3 were successfully replicated (P=7.08×10−7 and P=1.06×10−5, respectively, for the combined analyses). The association between rs3742207 and PWV was also successfully replicated (P=0.02) in an independent population, the Old-Order Amish, leading to an overall P=5.16×10−8. Conclusions—A genome-wide association study identified a SNP in the COL4A1 gene that was significantly associated with PWV in 2 populations. Collagen type 4 is the major structural component of basement membranes, suggesting that previously unrecognized cell-matrix interactions may exert an important role in regulating arterial stiffness.