Polymorphic Allele of Human IRGM1 Is Associated with Susceptibility to Tuberculosis in African Americans

An ancestral polymorphic allele of the human autophagy-related gene IRGM1 is associated with altered gene expression and a genetic risk for Crohn's Disease (CD). We used the single nucleotide polymorphism rs10065172C/T as a marker of this polymorphic allele and genotyped 370 African American and 177 Caucasian tuberculosis (TB) cases and 180 African American and 110 Caucasian controls. Among African Americans, the TB cases were more likely to carry the CD-related T allele of rs10065172 (odds ratio of 1.54; 95% confidence interval, 1.17–2.02; P<0.01) compared to controls. Our finding suggests that this CD-related IRGM1 polymorphic allele is also associated with human susceptibility to TB disease among African Americans

Last millennium northern hemisphere summer temperatures from tree rings: Part I: The long term context

Large-scale millennial length Northern Hemisphere (NH) temperature reconstructions have been progressively improved over the last 20 years as new datasets have been developed. This paper, and its companion (Part II, Anchukaitis et al. in prep), details the latest tree-ring (TR) based NH land air temperature reconstruction from a temporal and spatial perspective. This work is the first product of a consortium called N-TREND (Northern Hemisphere Tree-Ring Network Development) which brings together dendroclimatologists to identify a collective strategy for improving large-scale summer temperature reconstructions. The new reconstruction, N-TREND2015, utilises 54 records, a significant expansion compared with previous TR studies, and yields an improved reconstruction with stronger statistical calibration metrics. N-TREND2015 is relatively insensitive to the compositing method and spatial weighting used and validation metrics indicate that the new record portrays reasonable coherence with large scale summer temperatures and is robust at all time-scales from 918 to 2004 where at least 3 TR records exist from each major continental mass. N-TREND2015 indicates a longer and warmer medieval period (∼900–1170) than portrayed by previous TR NH reconstructions and by the CMIP5 model ensemble, but with better overall agreement between records for the last 600 years. Future dendroclimatic projects should focus on developing new long records from data-sparse regions such as North America and eastern Eurasia as well as ensuring the measurement of parameters related to latewood density to complement ring-width records which can improve local based calibration substantially

Mapping Copy Number Variation by Population Scale Genome Sequencing

Genomic structural variants (SVs) are abundant in humans, differing from other forms of variation in extent, origin and functional impact. Despite progress in SV characterization, the nucleotide resolution architecture of most SVs remains unknown. We constructed a map of unbalanced SVs (that is, copy number variants) based on whole genome DNA sequencing data from 185 human genomes, integrating evidence from complementary SV discovery approaches with extensive experimental validations. Our map encompassed 22,025 deletions and 6,000 additional SVs, including insertions and tandem duplications. Most SVs (53%) were mapped to nucleotide resolution, which facilitated analysing their origin and functional impact. We examined numerous whole and partial gene deletions with a genotyping approach and observed a depletion of gene disruptions amongst high frequency deletions. Furthermore, we observed differences in the size spectra of SVs originating from distinct formation mechanisms, and constructed a map of SV hotspots formed by common mechanisms. Our analytical framework and SV map serves as a resource for sequencing-based association studies.Organismic and Evolutionary Biolog

Characterising and Predicting Haploinsufficiency in the Human Genome

Haploinsufficiency, wherein a single functional copy of a gene is insufficient to maintain normal function, is a major cause of dominant disease. Human disease studies have identified several hundred haploinsufficient (HI) genes. We have compiled a map of 1,079 haplosufficient (HS) genes by systematic identification of genes unambiguously and repeatedly compromised by copy number variation among 8,458 apparently healthy individuals and contrasted the genomic, evolutionary, functional, and network properties between these HS genes and known HI genes. We found that HI genes are typically longer and have more conserved coding sequences and promoters than HS genes. HI genes exhibit higher levels of expression during early development and greater tissue specificity. Moreover, within a probabilistic human functional interaction network HI genes have more interaction partners and greater network proximity to other known HI genes. We built a predictive model on the basis of these differences and annotated 12,443 genes with their predicted probability of being haploinsufficient. We validated these predictions of haploinsufficiency by demonstrating that genes with a high predicted probability of exhibiting haploinsufficiency are enriched among genes implicated in human dominant diseases and among genes causing abnormal phenotypes in heterozygous knockout mice. We have transformed these gene-based haploinsufficiency predictions into haploinsufficiency scores for genic deletions, which we demonstrate to better discriminate between pathogenic and benign deletions than consideration of the deletion size or numbers of genes deleted. These robust predictions of haploinsufficiency support clinical interpretation of novel loss-of-function variants and prioritization of variants and genes for follow-up studies

Genetic Association Studies of Copy-Number Variation: Should Assignment of Copy Number States Precede Testing?

Recently, structural variation in the genome has been implicated in many complex diseases. Using genomewide single nucleotide polymorphism (SNP) arrays, researchers are able to investigate the impact not only of SNP variation, but also of copy-number variants (CNVs) on the phenotype. The most common analytic approach involves estimating, at the level of the individual genome, the underlying number of copies present at each location. Once this is completed, tests are performed to determine the association between copy number state and phenotype. An alternative approach is to carry out association testing first, between phenotype and raw intensities from the SNP array at the level of the individual marker, and then aggregate neighboring test results to identify CNVs associated with the phenotype. Here, we explore the strengths and weaknesses of these two approaches using both simulations and real data from a pharmacogenomic study of the chemotherapeutic agent gemcitabine. Our results indicate that pooled marker-level testing is capable of offering a dramatic increase in power (-fold) over CNV-level testing, particularly for small CNVs. However, CNV-level testing is superior when CNVs are large and rare; understanding these tradeoffs is an important consideration in conducting association studies of structural variation

Genome-Wide Patterns of Gene Expression during Aging in the African Malaria Vector Anopheles gambiae

The primary means of reducing malaria transmission is through reduction in longevity in days of the adult female stage of the Anopheles vector. However, assessing chronological age is limited to crude physiologic methods which categorize the females binomially as either very young (nulliparous) or not very young (parous). Yet the epidemiologically relevant reduction in life span falls within the latter category. Age-grading methods that delineate chronological age, using accurate molecular surrogates based upon gene expression profiles, will allow quantification of the longevity-reducing effects of vector control tools aimed at the adult, female mosquito. In this study, microarray analyses of gene expression profiles in the African malaria vector Anopheles gambiae were conducted during natural senescence of females in laboratory conditions. Results showed that detoxification-related and stress-responsive genes were up-regulated as mosquitoes aged. A total of 276 transcripts had age-dependent expression, independently of blood feeding and egg laying events. Expression of 112 (40.6%) of these transcripts increased or decreased monotonically with increasing chronologic age. Seven candidate genes for practical age assessment were tested by quantitative gene amplification in the An. gambiae G3 strain in a laboratory experiment and the Mbita strain in field enclosures set up in western Kenya under conditions closely resembling natural ones. Results were similar between experiments, indicating that senescence is marked by changes in gene expression and that chronological age can be gauged accurately and repeatedly with this method. These results indicate that the method may be suitable for accurate gauging of the age in days of field-caught, female An. gambiae

Genome-wide Association Analysis Identifies 14 New Risk Loci for Schizophrenia

Schizophrenia is a heritable disorder with substantial public health impact. We conducted a multi-stage genome-wide association study (GWAS) for schizophrenia beginning with a Swedish national sample (5,001 cases, 6,243 controls) followed by meta-analysis with prior schizophrenia GWAS (8,832 cases, 12,067 controls) and finally by replication of SNPs in 168 genomic regions in independent samples (7,413 cases, 19,762 controls, and 581 trios). In total, 22 regions met genome-wide significance (14 novel and one previously implicated in bipolar disorder). The results strongly implicate calcium signaling in the etiology of schizophrenia, and include genome-wide significant results for CACNA1C and CACNB2 whose protein products interact. We estimate that ∼8,300 independent and predominantly common SNPs contribute to risk for schizophrenia and that these collectively account for most of its heritability. Common genetic variation plays an important role in the etiology of schizophrenia, and larger studies will allow more detailed understanding of this devastating disorder

UGT1A1 is a major locus influencing bilirubin levels in African Americans

Total serum bilirubin is associated with several clinical outcomes, including cardiovascular disease, diabetes and drug metabolism. We conducted a genome-wide association study in 619 healthy unrelated African Americans in an attempt to replicate reported findings in Europeans and Asians and to identify novel loci influencing total serum bilirubin levels. We analyzed a dense panel of over two million genotyped and imputed SNPs in additive genetic models adjusting for age, sex, and the first two significant principal components from the sample covariance matrix of genotypes. Thirty-nine SNPs spanning a 78 kb region within the UGT1A1 displayed P-values <5 × 10−8. The lowest P-value was 1.7 × 10−22 for SNP rs887829. None of SNPs in the UGT1A1 remained statistically significant in conditional association analyses that adjusted for rs887829. In addition, SNP rs10929302 located in phenobarbital response enhancer module was significantly associated with bilirubin level with a P-value of 1.37 × 10−11; this enhancer module is believed to have a critical role in phenobarbital treatment of hyperbilirubinemia. Interestingly, the lead SNP, rs887829, is in strong linkage disequilibrium (LD) (r2≥0.74) with rs10929302. Taking advantage of the lower LD and shorter haplotypes in African-ancestry populations, we identified rs887829 as a more refined proxy for the causative variant influencing bilirubin levels. Also, we replicated the reported association between variants in SEMA3C and bilirubin levels. In summary, UGT1A1 is a major locus influencing bilirubin levels and the results of this study promise to contribute to understanding of the etiology and treatment of hyperbilirubinaemia in African-ancestry populations

A Genome-Wide Association Scan on the Levels of Markers of Inflammation in Sardinians Reveals Associations That Underpin Its Complex Regulation

Identifying the genes that influence levels of pro-inflammatory molecules can help to elucidate the mechanisms underlying this process. We first conducted a two-stage genome-wide association scan (GWAS) for the key inflammatory biomarkers Interleukin-6 (IL-6), the general measure of inflammation erythrocyte sedimentation rate (ESR), monocyte chemotactic protein-1 (MCP-1), and high-sensitivity C-reactive protein (hsCRP) in a large cohort of individuals from the founder population of Sardinia. By analysing 731,213 autosomal or X chromosome SNPs and an additional ∼1.9 million imputed variants in 4,694 individuals, we identified several SNPs associated with the selected quantitative trait loci (QTLs) and replicated all the top signals in an independent sample of 1,392 individuals from the same population. Next, to increase power to detect and resolve associations, we further genotyped the whole cohort (6,145 individuals) for 293,875 variants included on the ImmunoChip and MetaboChip custom arrays. Overall, our combined approach led to the identification of 9 genome-wide significant novel independent signals—5 of which were identified only with the custom arrays—and provided confirmatory evidence for an additional 7. Novel signals include: for IL-6, in the ABO gene (rs657152, p = 2.13×10−29); for ESR, at the HBB (rs4910472, p = 2.31×10−11) and UCN119B/SPPL3 (rs11829037, p = 8.91×10−10) loci; for MCP-1, near its receptor CCR2 (rs17141006, p = 7.53×10−13) and in CADM3 (rs3026968, p = 7.63×10−13); for hsCRP, within the CRP gene (rs3093077, p = 5.73×10−21), near DARC (rs3845624, p = 1.43×10−10), UNC119B/SPPL3 (rs11829037, p = 1.50×10−14), and ICOSLG/AIRE (rs113459440, p = 1.54×10−08) loci. Confirmatory evidence was found for IL-6 in the IL-6R gene (rs4129267); for ESR at CR1 (rs12567990) and TMEM57 (rs10903129); for MCP-1 at DARC (rs12075); and for hsCRP at CRP (rs1205), HNF1A (rs225918), and APOC-I (rs4420638). Our results improve the current knowledge of genetic variants underlying inflammation and provide novel clues for the understanding of the molecular mechanisms regulating this complex process

Rare coding variants in ten genes confer substantial risk for schizophrenia

Rare coding variation has historically provided the most direct connections between gene function and disease pathogenesis. By meta-analysing the whole exomes of 24,248 schizophrenia cases and 97,322 controls, we implicate ultra-rare coding variants (URVs) in 10 genes as conferring substantial risk for schizophrenia (odds ratios of 3-50, PPeer reviewe