Chad Garner

A QTL influencing F cell production maps to a gene encoding a zinc-finger protein on chromosome 2p15

F cells measure the presence of fetal hemoglobin, a heritable quantitative trait in adults that accounts for substantial phenotypic diversity of sickle cell disease and β thalassemia. We applied a genome-wide association mapping strategy to individuals with contrasting extreme trait values and mapped a new F cell quantitative trait locus to BCL11A, which encodes a zinc-finger protein, on chromosome 2p15. The 2p15 BCL11A quantitative trait locus accounts for 15.1% of the trait variance.

Intergenic variants of HBS1L-MYB are responsible for a major quantitative trait locus on chromosome 6q23 influencing fetal hemoglobin levels in adults.

Individual variation in fetal hemoglobin (HbF, α2γ2) response underlies the remarkable diversity in phenotypic severity of sickle cell disease and β thalassemia. HbF levels and HbF-associated quantitative traits (e.g., F cell levels) are highly heritable. We have previously mapped a major quantitative trait locus (QTL) controlling F cell levels in an extended Asian-Indian kindred with β thalassemia to a 1.5-Mb interval on chromosome 6q23, but the causative gene(s) are not known. The QTL encompasses several genes including HBS1L, a member of the GTP-binding protein family that is expressed in erythroid progenitor cells. In this high-resolution association study, we have identified multiple genetic variants within and 5′ to HBS1L at 6q23 that are strongly associated with F cell levels in families of Northern European ancestry (P = 10−75). The region accounts for 17.6% of the F cell variance in northern Europeans. Although mRNA levels of HBS1L and MYB in erythroid precursors grown in vitro are positively correlated, only HBS1L expression correlates with high F cell alleles. The results support a key role for the HBS1L-related genetic variants in HbF control and illustrate the biological complexity of the mechanism of 6q QTL as a modifier of fetal hemoglobin levels in the β hemoglobinopathies.

Deep sequencing of 10,000 human genomes

Significance Large-scale initiatives toward personalized medicine are driving a massive expansion in the number of human genomes being sequenced. Therefore, there is an urgent need to define quality standards for clinical use. This includes deep coverage and sequencing accuracy of an individual’s genome. Our work represents the largest effort to date in sequencing human genomes at deep coverage with these new standards. This study identifies over 150 million human variants, a majority of them rare and unknown. Moreover, these data identify sites in the genome that are highly intolerant to variation—possibly essential for life or health. We conclude that high-coverage genome sequencing provides accurate detail on human variation for discovery and clinical applications. We report on the sequencing of 10,545 human genomes at 30×–40× coverage with an emphasis on quality metrics and novel variant and sequence discovery. We find that 84% of an individual human genome can be sequenced confidently. This high-confidence region includes 91.5% of exon sequence and 95.2% of known pathogenic variant positions. We present the distribution of over 150 million single-nucleotide variants in the coding and noncoding genome. Each newly sequenced genome contributes an average of 8,579 novel variants. In addition, each genome carries on average 0.7 Mb of sequence that is not found in the main build of the hg38 reference genome. The density of this catalog of variation allowed us to construct high-resolution profiles that define genomic sites that are highly intolerant of genetic variation. These results indicate that the data generated by deep genome sequencing is of the quality necessary for clinical use.

The Extent of Linkage Disequilibrium Caused by Selection on G6PD in Humans

The gene coding for glucose-6-phosphate dehydrogenase (G6PD) is subject to positive selection by malaria in some human populations. The G6PD A− allele, which is common in sub-Saharan Africa, is associated with deficient enzyme activity and protection from severe malaria. To delimit the impact of selection on patterns of linkage disequilibrium (LD) and nucleotide diversity, we resequenced 5.1 kb at G6PD and ∼2–3 kb at each of eight loci in a 2.5-Mb region roughly centered on G6PD in a diverse sub-Saharan African panel of 51 unrelated men (including 20 G6PD A−, 11 G6PD A+, and 20 G6PD B chromosomes). The signature of selection is evident in the absence of genetic variation at G6PD and at three neighboring loci within 0.9 Mb from G6PD among all individuals bearing G6PD A− alleles. A genomic region of ∼1.6 Mb around G6PD was characterized by long-range LD associated with the A− alleles. These patterns of nucleotide variability and LD suggest that G6PD A− is younger than previous age estimates and has increased in frequency in sub-Saharan Africa due to strong selection (0.1 < s < 0.2). These results also show that selection can lead to nonrandom associations among SNPs over great physical and genetic distances, even in African populations.

Replication of celiac disease UK genome-wide association study results in a US population

Celiac disease is a common disease with a prevalence of approximately 1%. A recent genome-wide association study (GWAS) and follow-up study identified eight loci significantly associated with celiac disease risk. We genotyped the top 1020 non-HLA single nucleotide polymorphisms (SNPs) from the GWAS study that were genotyped in the previous follow-up study. After quality control assessments, 975 SNPs were analyzed for association with 906 celiac disease cases and 3819 controls, using logistic regression. Additional genotype data were generated by imputation and analyzed across the regions showing the strongest statistical evidence for association. Twenty SNPs were associated with celiac disease with P < 0.01 in the current study as well as in the previous follow-up study, of which 16 had P < 0.001 and 11 had P < 1 x 10(-11). Five of eight regions identified in the follow-up study were strongly associated with celiac disease, including regions on 1q31, 3q25, 3q28, 4q27 and 12q24. The strongest associations were at 4q27, the region most strongly associated in the GWAS and follow-up study and containing IL2 and IL21, and at 3q28 harboring LPP. In addition, we provide new evidence for an association, not previously reported, on 2q31 harboring a strong candidate gene, ITGA4. In conclusion, in this first follow-up study of celiac cases from the USA, we provide additional evidence that five of eight previously identified regions harbor risk alleles for celiac disease, and new evidence for an association on 2q31. The underlying functional mutations responsible for these replicated associations need to be identified.

Genetic and Environmental Influences on Left Ventricular Mass A Family Study

The relations between left ventricular mass (LVM) with age, gender, body size, and blood pressure were investigated in healthy adults and children from 149 nuclear families. LVM was strongly correlated with overall weight, especially in the children. Genetic analysis indicated that the segregation of LVM was compatible with an additive polygenic model, with a heritability estimate of 34% before adjustment for weight and 28% after adjustment for weight. Genetic and/or familial environmental factors played a strong role in the correlation of LVM and weight; they accounted for all of the correlation between the 2 traits in adults and 59% of the correlation in children. Spouses exhibited a strong correlation in their weight, which suggested that common family environment may contribute to the family correlations and to the observed heritability of the trait. LVM was strongly correlated with blood pressure before adjustment for weight, but this correlation could be attributed to nonfamilial environment rather than genetic factors. After adjustment for weight, the intertrait correlations between LVM and blood pressure were nonsignificant. Thus, adjustment for weight accounts for all common determinants of LVM and blood pressure.

Genome-Wide Association Study Identifies Genetic Loci Associated with Iron Deficiency

The existence of multiple inherited disorders of iron metabolism in man, rodents and other vertebrates suggests genetic contributions to iron deficiency. To identify new genomic locations associated with iron deficiency, a genome-wide association study (GWAS) was performed using DNA collected from white men aged ≥25 y and women ≥50 y in the Hemochromatosis and Iron Overload Screening (HEIRS) Study with serum ferritin (SF) ≤ 12 µg/L (cases) and iron replete controls (SF>100 µg/L in men, SF>50 µg/L in women). Regression analysis was used to examine the association between case-control status (336 cases, 343 controls) and quantitative serum iron measures and 331,060 single nucleotide polymorphism (SNP) genotypes, with replication analyses performed in a sample of 71 cases and 161 controls from a population of white male and female veterans screened at a US Veterans Affairs (VA) medical center. Five SNPs identified in the GWAS met genome-wide statistical significance for association with at least one iron measure, rs2698530 on chr. 2p14; rs3811647 on chr. 3q22, a known SNP in the transferrin (TF) gene region; rs1800562 on chr. 6p22, the C282Y mutation in the HFE gene; rs7787204 on chr. 7p21; and rs987710 on chr. 22q11 (GWAS observed P<1.51×10−7 for all). An association between total iron binding capacity and SNP rs3811647 in the TF gene (GWAS observed P = 7.0×10−9, corrected P = 0.012) was replicated within the VA samples (observed P = 0.012). Associations with the C282Y mutation in the HFE gene also were replicated. The joint analysis of the HEIRS and VA samples revealed strong associations between rs2698530 on chr. 2p14 and iron status outcomes. These results confirm a previously-described TF polymorphism and implicate one potential new locus as a target for gene identification.

Whole-genome sequencing identifies common-to-rare variants associated with human blood metabolites

Genetic factors modifying the blood metabolome have been investigated through genome-wide association studies (GWAS) of common genetic variants and through exome sequencing. We conducted a whole-genome sequencing study of common, low-frequency and rare variants to associate genetic variations with blood metabolite levels using comprehensive metabolite profiling in 1,960 adults. We focused the analysis on 644 metabolites with consistent levels across three longitudinal data collections. Genetic sequence variations at 101 loci were associated with the levels of 246 (38%) metabolites (P ≤ 1.9 × 10−11). We identified 113 (10.7%) among 1,054 unrelated individuals in the cohort who carried heterozygous rare variants likely influencing the function of 17 genes. Thirteen of the 17 genes are associated with inborn errors of metabolism or other pediatric genetic conditions. This study extends the map of loci influencing the metabolome and highlights the importance of heterozygous rare variants in determining abnormal blood metabolic phenotypes in adults.

Lamin B‐receptor mutations in Pelger–Huët anomaly

Summary. Pelger–Huët anomaly is an inherited abnormality of neutrophils, characterized by reduced nuclear segementation and an apparently looser chromatin structure. Following linkage studies in two families, the lamin B‐receptor (LBR) was sequenced and mutations found: CCG→CTG causing proline→leucine in codon 119 of exon 3, and IVS11‐9 A→G, disrupting the splice acceptor site. The LBR gene (LBR) was also sequenced from a single English man with Pelger–Huët anomaly and a heterozygous C→G mutation was found in codon 569 of exon 14, predicted to cause a proline→arginine. Our results confirm recently published findings that LBR mutations cause Pelger–Huët.

Association Analysis of the Extended MHC Region in Celiac Disease Implicates Multiple Independent Susceptibility Loci

Celiac disease is a common autoimmune disease caused by sensitivity to the dietary protein gluten. Forty loci have been implicated in the disease. All disease loci have been characterized as low-penetrance, with the exception of the high-risk genotypes in the HLA-DQA1 and HLA-DQB1 genes, which are necessary but not sufficient to cause the disease. The very strong effects from the known HLA loci and the genetically complex nature of the major histocompatibility complex (MHC) have precluded a thorough investigation of the region. The purpose of this study was to test the hypothesis that additional celiac disease loci exist within the extended MHC (xMHC). A set of 1898 SNPs was analyzed for association across the 7.6 Mb xMHC region in 1668 confirmed celiac disease cases and 517 unaffected controls. Conditional recursive partitioning was used to create an informative indicator of the known HLA-DQA1 and HLA-DQB1 high-risk genotypes that was included in the association analysis to account for their effects. A linkage disequilibrium-based grouping procedure was utilized to estimate the number of independent celiac disease loci present in the xMHC after accounting for the known effects. There was significant statistical evidence for four new independent celiac disease loci within the classic MHC region. This study is the first comprehensive association analysis of the xMHC in celiac disease that specifically accounts for the known HLA disease genotypes and the genetic complexity of the region.