Survival bias may unduly impact genetic association with complex diseases; gene-specific survival effects may further complicate such investigations. Coronary artery disease (CAD) is a complex phenotype for which little is understood about gene-specific survival effects; yet, such information can offer insight into refining genetic associations, improving replications, and can provide candidate genes for both mortality risk and improved survivorship in CAD. Building on our previous work, the purpose of this current study was to: evaluate LSAMP SNP-specific hazards for all-cause mortality post-catheterization in a larger cohort of our CAD cases; and, perform additional replication in an independent dataset. We examined two LSAMP SNPs-rs1462845 and rs6788787-using CAD case-only Cox proportional hazards regression for additive genetic effects, censored on time-to-all-cause mortality or last follow-up among Caucasian subjects from the Catheterization Genetics Study (CATHGEN; n = 2,224) and the Intermountain Heart Collaborative Study (IMHC; n = 3,008). Only after controlling for age, sex, body mass index, histories of smoking, type 2 diabetes, hyperlipidemia and hypertension (HR = 1.11, 95%CI = 1.01-1.22, p = 0.032), rs1462845 conferred significantly increased hazards of all-cause mortality among CAD cases. Even after controlling for multiple covariates, but in only the primary cohort, rs6788787 conferred significantly improved survival (HR = 0.80, 95% CI = 0.69-0.92, p = 0.002). Post-hoc analyses further stratifying by sex and disease severity revealed replicated effects for rs1462845: even after adjusting for aforementioned covariates and coronary interventional procedures, males with severe burden of CAD had significantly amplified hazards of death with the minor variant of rs1462845 in both cohorts (HR = 1.29, 95% CI = 1.08-1.55, p = 0.00456; replication HR = 1.25, 95% CI = 1.05-1.49, p = 0.013). Kaplan-Meier curves revealed unique cohort-specific genotype effects on survival. Additional analyses demonstrated that the homozygous risk genotype ('A/A') fully explained the increased hazard in both cohorts. None of the post-hoc analyses in control subjects were significant for any model. This suggests that genetic effects of rs1462845 on survival are unique to CAD presence. This represents formal, replicated evidence of genetic contribution of rs1462845 to increased risk for all-cause mortality; the contribution is unique to CAD case status and specific to males with severe burden of CAD.

The balance of myeloid populations and lymphoid populations must be well controlled. Here we found that osteopontin (OPN) skewed this balance during pathogenic conditions such as infection and autoimmunity. Notably, two isoforms of OPN exerted distinct effects in shifting this balance through cell-type-specific regulation of apoptosis. Intracellular OPN (iOPN) diminished the population size of myeloid progenitor cells and myeloid cells, and secreted OPN (sOPN) increase the population size of lymphoid cells. The total effect of OPN on skewing the leukocyte population balance was observed as host sensitivity to early systemic infection with Candida albicans and T cell-mediated colitis. Our study suggests previously unknown detrimental roles for two OPN isoforms in causing the imbalance of leukocyte populations.

Given the importance of cardiovascular disease (CVD) to public health and the demonstrated heritability of both disease status and its related risk factors, identifying the genetic variation underlying these susceptibilities is a critical step in understanding the pathogenesis of CVD and informing prevention and treatment strategies. Although one can look for genetic variation underlying susceptibility to CVD per se, it can be difficult to define the disease phenotype for such a qualitative analysis and CVD itself represents a convergence of diverse etiologic pathways. Alternatively, one can study the genetics of intermediate traits that are known risk factors for CVD, which can be measured quantitatively. Using the latter strategy, we have measured 21 cardiovascular-related biomarkers in an extended multigenerational pedigree, the CARRIAGE family (Carolinas Region Interaction of Aging, Genes, and Environment). These biomarkers belong to inflammatory and immune, connective tissue, lipid, and hemostasis pathways. Of these, 18 met our quality control standards. Using the pedigree and biomarker data, we have estimated the broad sense heritability (H2) of each biomarker (ranging from 0.09-0.56). A genome-wide panel of 6,015 SNPs was used subsequently to map these biomarkers as quantitative traits. Four showed noteworthy evidence for linkage in multipoint analysis (LOD score ≥ 2.6): paraoxonase (chromosome 8p11, 21), the chemokine RANTES (22q13.33), matrix metalloproteinase 3 (MMP3, 17p13.3), and granulocyte colony stimulating factor (GCSF, 8q22.1). Identifying the causal variation underlying each linkage score will help to unravel the genetic architecture of these quantitative traits and, by extension, the genetic architecture of cardiovascular risk.

Coronary artery disease (CAD), and one of its intermediate risk factors, dyslipidemia, possess a demonstrable genetic component, although the genetic architecture is incompletely defined. We previously reported a linkage peak on chromosome 5q31-33 for early-onset CAD where the strength of evidence for linkage was increased in families with higher mean low density lipoprotein-cholesterol (LDL-C). Therefore, we sought to fine-map the peak using association mapping of LDL-C as an intermediate disease-related trait to further define the etiology of this linkage peak. The study populations consisted of 1908 individuals from the CATHGEN biorepository of patients undergoing cardiac catheterization; 254 families (N = 827 individuals) from the GENECARD familial study of early-onset CAD; and 162 aorta samples harvested from deceased donors. Linkage disequilibrium-tagged SNPs were selected with an average of one SNP per 20 kb for 126.6-160.2 MB (region of highest linkage) and less dense spacing (one SNP per 50 kb) for the flanking regions (117.7-126.6 and 160.2-167.5 MB) and genotyped on all samples using a custom Illumina array. Association analysis of each SNP with LDL-C was performed using multivariable linear regression (CATHGEN) and the quantitative trait transmission disequilibrium test (QTDT; GENECARD). SNPs associated with the intermediate quantitative trait, LDL-C, were then assessed for association with CAD (i.e., a qualitative phenotype) using linkage and association in the presence of linkage (APL; GENECARD) and logistic regression (CATHGEN and aortas).

Amplification of the q21-q23 region on chromosome 1 is frequently found in sarcomas and a variety of other solid tumours. Previous analyses of sarcomas have indicated the presence of at least two separate amplicons within this region, one located in 1q21 and one located near the apolipoprotein A-II (APOA2) gene in 1q23. In this study we have mapped and characterized the amplicon in 1q23 in more detail.

Fuchs corneal dystrophy (FCD) is a degenerative genetic disorder of the corneal endothelium that represents one of the most common causes of corneal transplantation in the United States. Despite its high prevalence (4% over the age of 40), the underlying genetic basis of FCD is largely unknown. Here we report missense mutations in TCF8, a transcription factor whose haploinsufficiency causes posterior polymorphous corneal dystrophy (PPCD), in a cohort of late-onset FCD patients. In contrast to PPCD-causing mutations, all of which are null, FCD-associated mutations encode rare missense changes suggested to cause loss of function by an in vivo complementation assay. Importantly, segregation of a recurring p.Q840P mutation in a large, multigenerational FCD pedigree showed this allele to be sufficient but not necessary for pathogenesis. Execution of a genome-wide scan conditioned for the presence of the 840P allele identified an additional late-onset FCD locus on chromosome 9p, whereas haplotype analysis indicated that the presence of the TCF8 allele and the disease haplotype on 9p leads to a severe FCD manifestation with poor prognosis. Our data suggest that PPCD and FCD are allelic variants of the same disease continuum and that genetic interaction between genes that cause corneal dystrophies can modulate the expressivity of the phenotype.

Transcription factors are key mediators of human complex disease processes. Identifying the target genes of transcription factors will increase our understanding of the biological network leading to disease risk. The prediction of transcription factor binding sites (TFBSs) is one method to identify these target genes; however, current prediction methods need improvement. We chose the transcription factor upstream stimulatory factor 1 ( USF1 ) to evaluate the performance of our novel TFBS prediction method because of its known genetic association with coronary artery disease (CAD) and the recent availability of USF1 chromatin immunoprecipitation microarray (ChIP-chip) results. The specific goals of our study were to develop a novel and accurate genome-scale method for predicting USF1 binding sites and associated target genes to aid in the study of CAD. Previously published USF1 ChIP-chip data for 1 per cent of the genome were used to develop and evaluate several kernel logistic regression prediction models. A combination of genomic features (phylogenetic conservation, regulatory potential, presence of a CpG island and DNaseI hypersensitivity), as well as position weight matrix (PWM) scores, were used as variables for these models. Our most accurate predictor achieved an area under the receiver operator characteristic curve of 0.827 during cross-validation experiments, significantly outperforming standard PWM-based prediction methods. When applied to the whole human genome, we predicted 24,010 USF1 binding sites within 5 kilobases upstream of the transcription start site of 9,721 genes. These predictions included 16 of 20 genes with strong evidence of USF1 regulation. Finally, in the spirit of genomic convergence, we integrated independent experimental CAD data with these USF1 binding site prediction results to develop a prioritised set of candidate genes for future CAD studies. We have shown that our novel prediction method, which employs genomic features related to the presence of regulatory elements, enables more accurate and efficient prediction of USF1 binding sites. This method can be extended to other transcription factors identified in human disease studies to help further our understanding of the biology of complex disease.

Multiple sclerosis (MS) is an autoimmune disorder where T cells attack neurons in the central nervous system (CNS) leading to demyelination and neurological deficits. A driver of increased MS risk is the soluble form of the interleukin-7 receptor alpha chain gene (sIL7R) produced by alternative splicing of IL7R exon 6. Here, we identified the RNA helicase DDX39B as a potent activator of this exon and consequently a repressor of sIL7R, and we found strong genetic association of DDX39B with MS risk. Indeed, we showed that a genetic variant in the 5' UTR of DDX39B reduces translation of DDX39B mRNAs and increases MS risk. Importantly, this DDX39B variant showed strong genetic and functional epistasis with allelic variants in IL7R exon 6. This study establishes the occurrence of biological epistasis in humans and provides mechanistic insight into the regulation of IL7R exon 6 splicing and its impact on MS risk.

During pandemics, individuals exhibit differences in risk and clinical outcomes. Here, we developed single-cell high-throughput human in vitro susceptibility testing (scHi-HOST), a method for rapidly identifying genetic variants that confer resistance and susceptibility. We applied this method to influenza A virus (IAV), the cause of four pandemics since the start of the 20th century. scHi-HOST leverages single-cell RNA sequencing (scRNA-seq) to simultaneously assign genetic identity to cells in mixed infections of cell lines of European, African, and Asian origin, reveal associated genetic variants for viral burden, and identify expression quantitative trait loci. Integration of scHi-HOST with human challenge and experimental validation demonstrated that a missense variant in endoplasmic reticulum aminopeptidase 1 (ERAP1; rs27895) increased IAV burden in cells and human volunteers. rs27895 exhibits population differentiation, likely contributing to greater permissivity of cells from African populations to IAV. scHi-HOST is a broadly applicable method and resource for decoding infectious-disease genetics.

Lymphoblastoid cell lines (LCLs) are generated by transforming primary B cells with Epstein-Barr virus (EBV) and are used extensively as model systems in viral oncology, immunology, and human genetics research. In this study, we characterized single-cell transcriptomic profiles of five LCLs and present a simple discrete-time simulation to explore the influence of stochasticity on LCL clonal evolution. Single-cell RNA sequencing (scRNA-seq) revealed substantial phenotypic heterogeneity within and across LCLs with respect to immunoglobulin isotype; virus-modulated host pathways involved in survival, activation, and differentiation; viral replication state; and oxidative stress. This heterogeneity is likely attributable to intrinsic variance in primary B cells and host-pathogen dynamics. Stochastic simulations demonstrate that initial primary cell heterogeneity, random sampling, time in culture, and even mild differences in phenotype-specific fitness can contribute substantially to dynamic diversity in populations of nominally clonal cells.

Background DNA methylation is implicated in many chronic diseases and may contribute to mortality. Therefore, we conducted an epigenome-wide association study (EWAS) for all-cause mortality with whole-transcriptome data in a cardiovascular cohort (CATHGEN [Catheterization Genetics]). Methods and Results Cases were participants with mortality≥7 days postcatheterization whereas controls were alive with≥2 years of follow-up. The Illumina Human Methylation 450K and EPIC arrays (Illumina, San Diego, CA) were used for the discovery and validation sets, respectively. A linear model approach with empirical Bayes estimators adjusted for confounders was used to assess difference in methylation (Δβ). In the discovery set (55 cases, 49 controls), 25 629 (6.5%) probes were differently methylated (P<0.05). In the validation set (108 cases, 108 controls), 3 probes were differentially methylated with a false discovery rate-adjusted P<0.10: cg08215811 (SLC4A9; log2 fold change=-0.14); cg17845532 (MATK; fold change=-0.26); and cg17944110 (castor zinc finger 1 [CASZ1]; FC=0.26; P<0.0001; false discovery rate-adjusted P=0.046-0.080). Meta-analysis identified 6 probes (false discovery rate-adjusted P<0.05): the 3 above, cg20428720 (intergenic), cg17647904 (NCOR2), and cg23198793 (CAPN3). Messenger RNA expression of 2 MATK isoforms was lower in cases (fold change=-0.24 [P=0.007] and fold change=-0.61 [P=0.009]). The CASZ1, NCOR2, and CAPN3 transcripts did not show differential expression (P>0.05); the SLC4A9 transcript did not pass quality control. The cg17944110 probe is located within a potential regulatory element; expression of predicted targets (using GeneHancer) of the regulatory element, UBIAD1 (P=0.01) and CLSTN1 (P=0.03), were lower in cases. Conclusions We identified 6 novel methylation sites associated with all-cause mortality. Methylation in CASZ1 may serve as a regulatory element associated with mortality in cardiovascular patients. Larger studies are necessary to confirm these observations.

White matter injuries (WMIs) are the leading cause of neurologic impairment in infants born premature. There are no treatment options available. The most common forms of WMIs in infants occur prior to the onset of normal myelination, making its pathophysiology distinctive, thus requiring a tailored approach to treatment. Neonates present a unique opportunity to repair WMIs due to a transient abundance of neural stem/progenitor cells (NSPCs) present in the germinal matrix with oligodendrogenic potential. We identified an endogenous oxysterol, 20-αHydroxycholesterol (20HC), in human maternal breast milk that induces oligodendrogenesis through a sonic hedgehog (shh), Gli-dependent mechanism. Following WMI in neonatal mice, injection of 20HC induced subventricular zone-derived oligodendrogenesis and improved myelination in the periventricular white matter, resulting in improved motor outcomes. Targeting the oligodendrogenic potential of postnatal NSPCs in neonates with WMIs may be further developed into a novel approach to mitigate this devastating complication of preterm birth.

Chiari Type I Malformation (CMI) is characterized by displacement of the cerebellar tonsils below the base of the skull, resulting in significant neurologic morbidity. Although multiple lines of evidence support a genetic contribution to disease, no genes have been identified. We therefore conducted the largest whole genome linkage screen to date using 367 individuals from 66 families with at least two individuals presenting with nonsyndromic CMI with or without syringomyelia. Initial findings across all 66 families showed minimal evidence for linkage due to suspected genetic heterogeneity. In order to improve power to localize susceptibility genes, stratified linkage analyses were performed using clinical criteria to differentiate families based on etiologic factors. Families were stratified on the presence or absence of clinical features associated with connective tissue disorders (CTDs) since CMI and CTDs frequently co-occur and it has been proposed that CMI patients with CTDs represent a distinct class of patients with a different underlying disease mechanism. Stratified linkage analyses resulted in a marked increase in evidence of linkage to multiple genomic regions consistent with reduced genetic heterogeneity. Of particular interest were two regions (Chr8, Max LOD = 3.04; Chr12, Max LOD = 2.09) identified within the subset of "CTD-negative" families, both of which harbor growth differentiation factors (GDF6, GDF3) implicated in the development of Klippel-Feil syndrome (KFS). Interestingly, roughly 3-5% of CMI patients are diagnosed with KFS. In order to investigate the possibility that CMI and KFS are allelic, GDF3 and GDF6 were sequenced leading to the identification of a previously known KFS missense mutation and potential regulatory variants in GDF6. This study has demonstrated the value of reducing genetic heterogeneity by clinical stratification implicating several convincing biological candidates and further supporting the hypothesis that multiple, distinct mechanisms are responsible for CMI.

Fuchs endothelial corneal dystrophy (FECD) is a genetically heterogeneous disorder that has been primarily studied in patients of European or Asian ancestry. Given the sparse literature on African Americans with FECD, we sought to characterize the genetic variation in three known FECD candidate genes in African American patients with FECD.

Fuchs endothelial corneal dystrophy (FECD) is a common, late-onset disorder of the corneal endothelium. Although progress has been made in understanding the genetic basis of FECD by studying large families in which the phenotype is transmitted in an autosomal dominant fashion, a recently reported genome-wide association study identified common alleles at a locus on chromosome 18 near TCF4 which confer susceptibility to FECD. Here, we report the findings of our independent validation study for TCF4 using the largest FECD dataset to date (450 FECD cases and 340 normal controls). Logistic regression with sex as a covariate was performed for three genetic models: dominant (DOM), additive (ADD), and recessive (REC). We found significant association with rs613872, the target marker reported by Baratz et al.(2010), for all three genetic models (DOM: P = 9.33×10(-35); ADD: P = 7.48×10(-30); REC: P = 5.27×10(-6)). To strengthen the association study, we also conducted a genome-wide linkage scan on 64 multiplex families, composed primarily of affected sibling pairs (ASPs), using both parametric and non-parametric two-point and multipoint analyses. The most significant linkage region localizes to chromosome 18 from 69.94cM to 85.29cM, with a peak multipoint HLOD = 2.5 at rs1145315 (75.58cM) under the DOM model, mapping 1.5 Mb proximal to rs613872. In summary, our study presents evidence to support the role of the intronic TCF4 single nucleotide polymorphism rs613872 in late-onset FECD through both association and linkage studies.

To provide a definitive linkage map for multiple sclerosis, we have genotyped the Illumina BeadArray linkage mapping panel (version 4) in a data set of 730 multiplex families of Northern European descent. After the application of stringent quality thresholds, data from 4,506 markers in 2,692 individuals were included in the analysis. Multipoint nonparametric linkage analysis revealed highly significant linkage in the major histocompatibility complex (MHC) on chromosome 6p21 (maximum LOD score [MLS] 11.66) and suggestive linkage on chromosomes 17q23 (MLS 2.45) and 5q33 (MLS 2.18). This set of markers achieved a mean information extraction of 79.3% across the genome, with a Mendelian inconsistency rate of only 0.002%. Stratification based on carriage of the multiple sclerosis-associated DRB1*1501 allele failed to identify any other region of linkage with genomewide significance. However, ordered-subset analysis suggested that there may be an additional locus on chromosome 19p13 that acts independent of the main MHC locus. These data illustrate the substantial increase in power that can be achieved with use of the latest tools emerging from the Human Genome Project and indicate that future attempts to systematically identify susceptibility genes for multiple sclerosis will have to involve large sample sizes and an association-based methodology.

Variation in the major histocompatibility complex (MHC) on chromosome 6p21 is known to influence susceptibility to multiple sclerosis with the strongest effect originating from the HLA-DRB1 gene in the class II region. The possibility that other genes in the MHC independently influence susceptibility to multiple sclerosis has been suggested but remains unconfirmed.

The transcription factor GATA2 plays an essential role in the establishment and maintenance of adult hematopoiesis. It is expressed in hematopoietic stem cells, as well as the cells that make up the aortic vasculature, namely aortic endothelial cells and smooth muscle cells. We have shown that GATA2 expression is predictive of location within the thoracic aorta; location is suggested to be a surrogate for disease susceptibility. The GATA2 gene maps beneath the Chromosome 3q linkage peak from our family-based sample set (GENECARD) study of early-onset coronary artery disease. Given these observations, we investigated the relationship of several known and novel polymorphisms within GATA2 to coronary artery disease. We identified five single nucleotide polymorphisms that were significantly associated with early-onset coronary artery disease in GENECARD. These results were validated by identifying significant association of two of these single nucleotide polymorphisms in an independent case-control sample set that was phenotypically similar to the GENECARD families. These observations identify GATA2 as a novel susceptibility gene for coronary artery disease and suggest that the study of this transcription factor and its downstream targets may uncover a regulatory network important for coronary artery disease inheritance.

Androgen receptor splice variant 7 (AR-V7) results in a truncated receptor, which leads to ligand-independent constitutive activation that is not inhibited by anti-androgen therapies, including abiraterone or enzalutamide. Given that previous reports suggested that circulating tumor cell (CTC) AR-V7 detection is a poor prognostic indicator for the clinical efficacy of secondary hormone therapies, we conducted a prospective multicenter validation study.

Air pollution is a worldwide contributor to cardiovascular disease mortality and morbidity. Traffic-related air pollution is a widespread environmental exposure and is associated with multiple cardiovascular outcomes such as coronary atherosclerosis, peripheral arterial disease, and myocardial infarction. Despite the recognition of the importance of both genetic and environmental exposures to the pathogenesis of cardiovascular disease, studies of how these two contributors operate jointly are rare. We performed a genome-wide interaction study (GWIS) to examine gene-traffic exposure interactions associated with coronary atherosclerosis. Using race-stratified cohorts of 538 African-Americans (AA) and 1562 European-Americans (EA) from a cardiac catheterization cohort (CATHGEN), we identify gene-by-traffic exposure interactions associated with the number of significantly diseased coronary vessels as a measure of chronic atherosclerosis. We found five suggestive (P<1x10-5) interactions in the AA GWIS, of which two (rs1856746 and rs2791713) replicated in the EA cohort (P < 0.05). Both SNPs are in the PIGR-FCAMR locus and are eQTLs in lymphocytes. The protein products of both PIGR and FCAMR are implicated in inflammatory processes. In the EA GWIS, there were three suggestive interactions; none of these replicated in the AA GWIS. All three were intergenic; the most significant interaction was in a regulatory region associated with SAMSN1, a gene previously associated with atherosclerosis and B cell activation. In conclusion, we have uncovered several novel genes associated with coronary atherosclerosis in individuals chronically exposed to increased ambient concentrations of traffic air pollution. These genes point towards inflammatory pathways that may modify the effects of air pollution on cardiovascular disease risk.