Canadian COVID-19 host genetics cohort replicates known severity associations.

The HostSeq initiative recruited 10,059 Canadians infected with SARS-CoV-2 between March 2020 and March 2023, obtained clinical information on their disease experience and whole genome sequenced (WGS) their DNA. We analyzed the WGS data for genetic contributors to severe COVID-19 (considering 3,499 hospitalized cases and 4,975 non-hospitalized after quality control). We investigated the evidence for replication of loci reported by the International Host Genetics Initiative (HGI); analyzed the X chromosome; conducted rare variant gene-based analysis and polygenic risk score testing. Population stratification was adjusted for using meta-analysis across ancestry groups. We replicated two loci identified by the HGI for COVID-19 severity: the LZTFL1/SLC6A20 locus on chromosome 3 and the FOXP4 locus on chromosome 6 (the latter with a variant significant at P < 5E-8). We found novel significant associations with MRAS and WDR89 in gene-based analyses, and constructed a polygenic risk score that explained 1.01% of the variance in severe COVID-19. This study provides independent evidence confirming the robustness of previously identified COVID-19 severity loci by the HGI and identifies novel genes for further investigation.

Characterization of direct and/or indirect genetic associations for multiple traits in longitudinal studies of disease progression.

When quantitative longitudinal traits are risk factors for disease progression and subject to random biological variation, joint model analysis of time-to-event and longitudinal traits can effectively identify direct and/or indirect genetic association of single nucleotide polymorphisms (SNPs) with time-to-event. We present a joint model that integrates: (1) a multivariate linear mixed model describing trajectories of multiple longitudinal traits as a function of time, SNP effects, and subject-specific random effects and (2) a frailty Cox survival model that depends on SNPs, longitudinal trajectory effects, and subject-specific frailty accounting for dependence among multiple time-to-event traits. Motivated by complex genetic architecture of type 1 diabetes complications (T1DC) observed in the Diabetes Control and Complications Trial (DCCT), we implement a 2-stage approach to inference with bootstrap joint covariance estimation and develop a hypothesis testing procedure to classify direct and/or indirect SNP association with each time-to-event trait. By realistic simulation study, we show that joint modeling of 2 time-to-T1DC (retinopathy and nephropathy) and 2 longitudinal risk factors (HbA1c and systolic blood pressure) reduces estimation bias in genetic effects and improves classification accuracy of direct and/or indirect SNP associations, compared to methods that ignore within-subject risk factor variability and dependence among longitudinal and time-to-event traits. Through DCCT data analysis, we demonstrate feasibility for candidate SNP modeling and quantify effects of sample size and Winner's curse bias on classification for 2 SNPs identified as having indirect associations with time-to-T1DC traits. Joint analysis of multiple longitudinal and multiple time-to-event traits provides insight into complex traits architecture.

Penalized maximum likelihood inference under the mixture cure model in sparse data.

INTRODUCTION: When a study sample includes a large proportion of long-term survivors, mixture cure (MC) models that separately assess biomarker associations with long-term recurrence-free survival and time to disease recurrence are preferred to proportional-hazards models. However, in samples with few recurrences, standard maximum likelihood can be biased. OBJECTIVE AND METHODS: We extend Firth-type penalized likelihood (FT-PL) developed for bias reduction in the exponential family to the Weibull-logistic MC, using the Jeffreys invariant prior. Via simulation studies based on a motivating cohort study, we compare parameter estimates of the FT-PL method to those by ML, as well as type 1 error (T1E) and power obtained using likelihood ratio statistics. RESULTS: In samples with relatively few events, the Firth-type penalized likelihood estimates (FT-PLEs) have mean bias closer to zero and smaller mean squared error than maximum likelihood estimates (MLEs), and can be obtained in samples where the MLEs are infinite. Under similar T1E rates, FT-PL consistently exhibits higher statistical power than ML in samples with few events. In addition, we compare FT-PL estimation with two other penalization methods (a log-F prior method and a modified Firth-type method) based on the same simulations. DISCUSSION: Consistent with findings for logistic and Cox regressions, FT-PL under MC regression yields finite estimates under stringent conditions, and better bias-and-variance balance than the other two penalizations. The practicality and strength of FT-PL for MC analysis is illustrated in a cohort study of breast cancer prognosis with long-term follow-up for recurrence-free survival.

Two-phase sample selection strategies for design and analysis in post-genome-wide association fine-mapping studies.

Post-GWAS analysis, in many cases, focuses on fine-mapping targeted genetic regions discovered at GWAS-stage; that is, the aim is to pinpoint potential causal variants and susceptibility genes for complex traits and disease outcomes using next-generation sequencing (NGS) technologies. Large-scale GWAS cohorts are necessary to identify target regions given the typically modest genetic effect sizes. In this context, two-phase sampling design and analysis is a cost-reduction technique that utilizes data collected during phase 1 GWAS to select an informative subsample for phase 2 sequencing. The main goal is to make inference for genetic variants measured via NGS by efficiently combining data from phases 1 and 2. We propose two approaches for selecting a phase 2 design under a budget constraint. The first method identifies sampling fractions that select a phase 2 design yielding an asymptotic variance covariance matrix with certain optimal characteristics, for example, smallest trace, via Lagrange multipliers (LM). The second relies on a genetic algorithm (GA) with a defined fitness function to identify exactly a phase 2 subsample. We perform comprehensive simulation studies to evaluate the empirical properties of the proposed designs for a genetic association study of a quantitative trait. We compare our methods against two ranked designs: residual-dependent sampling and a recently identified optimal design. Our findings demonstrate that the proposed designs, GA in particular, can render competitive power in combined phase 1 and 2 analysis compared with alternative designs while preserving type 1 error control. These results are especially evident under the more practical scenario where design values need to be defined a priori and are subject to misspecification. We illustrate the proposed methods in a study of triglyceride levels in the North Finland Birth Cohort of 1966. R code to reproduce our results is available at github.com/egosv/TwoPhase_postGWAS.

Non-invasive Electroarthrography Measures Load-Induced Cartilage Streaming Potentials via Electrodes Placed on Skin Surrounding an Articular Joint.

OBJECTIVE: We aimed to demonstrate that electroarthrography (EAG) measures streaming potentials originating in the cartilage extracellular matrix during load bearing through electrodes adhered to skin surrounding an articular joint. DESIGN: Equine metacarpophalangeal joints were subjected to simulated physiological loads while (1) replacing synovial fluid with immersion buffers of different electrolyte concentrations and (2) directly degrading cartilage with trypsin. RESULTS: An inverse relationship between ionic strength and EAG coefficient was detected. Compared to native synovial fluid, EAG coefficients increased (P < 0.05) for 5 of 6 electrodes immersed in 0.1X phosphate-buffered saline (PBS) (0.014 M NaCl), decreased (P < 0.05) for 4 of 6 electrodes in 1X PBS (0.14 M NaCl), and decreased (P < 0.05) for all 6 electrodes in 10X PBS (1.4 M NaCl). This relationship corresponds to similar studies where streaming potentials were directly measured on cartilage. EAG coefficients, obtained after trypsin degradation, were reduced (P < 0.05) in 6 of 8, and 7 of 8 electrodes, during simulated standing and walking, respectively. Trypsin degradation was confirmed by direct cartilage assessments. Streaming potentials, measured by directly contacting cartilage, indicated lower cartilage stiffness (P < 10-5). Unconfined compression data revealed reduced Em, representing proteoglycan matrix stiffness (P = 0.005), no change in Ef, representing collagen network stiffness (P = 0.15), and no change in permeability (P = 0.24). Trypsin depleted proteoglycan as observed by both dimethylmethylene blue assay (P = 0.0005) and safranin-O stained histological sections. CONCLUSION: These data show that non-invasive EAG detects streaming potentials produced by cartilage during joint compression and has potential to become a diagnostic tool capable of detecting early cartilage degeneration.

Statistical power in COVID-19 case-control host genomic study design.

The identification of genetic variation that directly impacts infection susceptibility to SARS-CoV-2 and disease severity of COVID-19 is an important step towards risk stratification, personalized treatment plans, therapeutic, and vaccine development and deployment. Given the importance of study design in infectious disease genetic epidemiology, we use simulation and draw on current estimates of exposure, infectivity, and test accuracy of COVID-19 to demonstrate the feasibility of detecting host genetic factors associated with susceptibility and severity in published COVID-19 study designs. We demonstrate that limited phenotypic data and exposure/infection information in the early stages of the pandemic significantly impact the ability to detect most genetic variants with moderate effect sizes, especially when studying susceptibility to SARS-CoV-2 infection. Our insights can aid in the interpretation of genetic findings emerging in the literature and guide the design of future host genetic studies.

Gene-based and pathway-based testing for rare-variant association in affected sib pairs.

Next generation sequencing technologies have made it possible to investigate the role of rare variants (RVs) in disease etiology. Because RVs associated with disease susceptibility tend to be enriched in families with affected individuals, study designs based on affected sib pairs (ASP) can be more powerful than case-control studies. We construct tests of RV-set association in ASPs for single genomic regions as well as for multiple regions. Single-region tests can efficiently detect a gene region harboring susceptibility variants, while multiple-region extensions are meant to capture signals dispersed across a biological pathway, potentially as a result of locus heterogeneity. Within ascertained ASPs, the test statistics contrast the frequencies of duplicate rare alleles (usually appearing on a shared haplotype) against frequencies of a single rare allele copy (appearing on a nonshared haplotype); we call these allelic parity tests. Incorporation of minor allele frequency estimates from reference populations can markedly improve test efficiency. Under various genetic penetrance models, application of the tests in simulated ASP data sets demonstrates good type I error properties as well as power gains over approaches that regress ASP rare allele counts on sharing state, especially in small samples. We discuss robustness of the allelic parity methods to the presence of genetic linkage, misspecification of reference population allele frequencies, sequencing error and de novo mutations, and population stratification. As proof of principle, we apply single- and multiple-region tests in a motivating study data set consisting of whole exome sequencing of sisters ascertained with early onset breast cancer.

DNA methylation age calculators reveal association with diabetic neuropathy in type 1 diabetes.

BACKGROUND: Many CpGs become hyper or hypo-methylated with age. Multiple methods have been developed by Horvath et al. to estimate DNA methylation (DNAm) age including Pan-tissue, Skin & Blood, PhenoAge, and GrimAge. Pan-tissue and Skin & Blood try to estimate chronological age in the normal population whereas PhenoAge and GrimAge use surrogate markers associated with mortality to estimate biological age and its departure from chronological age. Here, we applied Horvath's four methods to calculate and compare DNAm age in 499 subjects with type 1 diabetes (T1D) from the Diabetes Control and Complications Trial/Epidemiology of Diabetes Interventions and Complications (DCCT/EDIC) study using DNAm data measured by Illumina EPIC array in the whole blood. Association of the four DNAm ages with development of diabetic complications including cardiovascular diseases (CVD), nephropathy, retinopathy, and neuropathy, and their risk factors were investigated. RESULTS: Pan-tissue and GrimAge were higher whereas Skin & Blood and PhenoAge were lower than chronological age (p < 0.0001). DNAm age was not associated with the risk of CVD or retinopathy over 18-20 years after DNAm measurement. However, higher PhenoAge (ß = 0.023, p = 0.007) and GrimAge (ß = 0.029, p = 0.002) were associated with higher albumin excretion rate (AER), an indicator of diabetic renal disease, measured over time. GrimAge was also associated with development of both diabetic peripheral neuropathy (OR = 1.07, p = 9.24E-3) and cardiovascular autonomic neuropathy (OR = 1.06, p = 0.011). Both HbA1c (ß = 0.38, p = 0.026) and T1D duration (ß = 0.01, p = 0.043) were associated with higher PhenoAge. Employment (ß = - 1.99, p = 0.045) and leisure time (ß = - 0.81, p = 0.022) physical activity were associated with lower Pan-tissue and Skin & Blood, respectively. BMI (ß = 0.09, p = 0.048) and current smoking (ß = 7.13, p = 9.03E-50) were positively associated with Skin & Blood and GrimAge, respectively. Blood pressure, lipid levels, pulse rate, and alcohol consumption were not associated with DNAm age regardless of the method used. CONCLUSIONS: Various methods of measuring DNAm age are sub-optimal in detecting people at higher risk of developing diabetic complications although some work better than the others.

Osteosarcoma and soft-tissue sarcomas with an immune infiltrate express PD-L1: relation to clinical outcome and Th1 pathway activation.

Immune checkpoint proteins, such as PD-L1 and PD-1, are important in several cancers; however, their role in osteosarcoma (OSA) and soft tissue sarcoma (STS) remains unclear. Our aims were to determine whether subsets of OSA/STS harbor tumor-infiltrating lymphocytes (TILs) and express PD-L1, and how PD-L1 expression is related to clinical outcome. Tissue sections of 25 cases each of untreated undifferentiated pleomorphic sarcoma (UPS), myxofibrosarcoma (MFS), liposarcoma (LPS) and 24 of leiomyosarcoma (LMS) were subjected to immunohistochemistry (IHC) for immune cells, PD-L1 and PD-1. RT-qPCR was utilized to quantify levels of PD-L1 mRNA from 33 UPS, 57 MFS and 79 OSA primary-untreated specimens. PD-L1 mRNA levels were tested for their correlation with overall survival in patients presenting without metastases. Transcriptome analysis evaluated biological pathway differences between high and low PD-L1 expressers. A subset of UPS and MFS contained TILs and expressed PD-L1 and PD-1; LMS and LPS did not. PD-L1 levels by IHC and RT-qPCR were positively correlated. PD-L1 over-expression was associated with better survival for UPS and OSA, but not MFS. The Th1 pathway was significantly activated in UPS with high levels of PD-L1 and improved survival. Some sarcoma subtypes harbor TILs and express PD-L1. Patients with UPS and OSA with high levels of PD-L1 had better overall survival than those with low expression levels. Important biological pathways distinguish PD-L1 high and low groups. The stratification of patients with OSA/STS with respect to potential immune therapies may be improved through investigation of the expression of immune cells and checkpoint proteins.

Genome-Wide Association Study of Diabetic Kidney Disease Highlights Biology Involved in Glomerular Basement Membrane Collagen.

BACKGROUND: Although diabetic kidney disease demonstrates both familial clustering and single nucleotide polymorphism heritability, the specific genetic factors influencing risk remain largely unknown. METHODS: To identify genetic variants predisposing to diabetic kidney disease, we performed genome-wide association study (GWAS) analyses. Through collaboration with the Diabetes Nephropathy Collaborative Research Initiative, we assembled a large collection of type 1 diabetes cohorts with harmonized diabetic kidney disease phenotypes. We used a spectrum of ten diabetic kidney disease definitions based on albuminuria and renal function. RESULTS: Our GWAS meta-analysis included association results for up to 19,406 individuals of European descent with type 1 diabetes. We identified 16 genome-wide significant risk loci. The variant with the strongest association (rs55703767) is a common missense mutation in the collagen type IV alpha 3 chain (COL4A3) gene, which encodes a major structural component of the glomerular basement membrane (GBM). Mutations in COL4A3 are implicated in heritable nephropathies, including the progressive inherited nephropathy Alport syndrome. The rs55703767 minor allele (Asp326Tyr) is protective against several definitions of diabetic kidney disease, including albuminuria and ESKD, and demonstrated a significant association with GBM width; protective allele carriers had thinner GBM before any signs of kidney disease, and its effect was dependent on glycemia. Three other loci are in or near genes with known or suggestive involvement in this condition (BMP7) or renal biology (COLEC11 and DDR1). CONCLUSIONS: The 16 diabetic kidney disease-associated loci may provide novel insights into the pathogenesis of this condition and help identify potential biologic targets for prevention and treatment.

gpart: human genome partitioning and visualization of high-density SNP data by identifying haplotype blocks.

SUMMARY: For the analysis of high-throughput genomic data produced by next-generation sequencing (NGS) technologies, researchers need to identify linkage disequilibrium (LD) structure in the genome. In this work, we developed an R package gpart which provides clustering algorithms to define LD blocks or analysis units consisting of SNPs. The visualization tool in gpart can display the LD structure and gene positions for up to 20 000 SNPs in one image. The gpart functions facilitate construction of LD blocks and SNP partitions for vast amounts of genome sequencing data within reasonable time and memory limits in personal computing environments. AVAILABILITY AND IMPLEMENTATION: The R package is available at https://bioconductor.org/packages/gpart. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Fresh Cut Versus Stored Cut Paraffin-embedded Tissue: Effect on Immunohistochemical Staining for Common Breast Cancer Markers.

The proper handling of unstained paraffin slides for immunohistochemistry has been a matter of debate, with several studies demonstrating loss of antigenicity with prolonged storage at room temperature, 4°C and -20°C. The purpose of this study was to determine whether long-term storage of unstained slides at -80°C would impact the staining intensity and expression distribution of markers used to molecularly subtype breast cancer specimens [estrogen receptor (ER), human epidermal growth factor receptor 2 (HER2), cytokeratin 5 (CK5), epidermal growth factor receptor (EGFR), and Ki67]. The staining pattern of previously unstained breast tumor slides (n=39 to 64) stored at -80°C for a minimum of 9.93 years (avg., 12.8 y) was compared with the staining pattern of fresh cut slides from the same tumors. The Allred scoring method was used to score ER (0 to 2, negative; 3 to 8, positive), CK5 (≥4, positive), and EGFR (≥4, positive). ASCO/CAP guidelines were used to assess HER2 (0/1+, 2+, or 3+). Ki67 scores were determined based on the proportion of cells stained of any intensity, with 20% staining used as a cut-off. Agreement was assessed using concordance rates and chance-corrected agreement statistics. The chance-corrected agreements were as follows: 0.94 (38/39) for ER, 0.92 (53/55) for CK5, 0.87 (61/64) for EGFR, 0.86 (37/39) for HER2, and 0.67 (46/54) for Ki67. Long-term storage of cut unstained slides at -80°C does not significantly impact the scoring interpretation of ER, CK5, EGFR, and HER2.

Integrating epigenetic, genetic, and phenotypic data to uncover gene-region associations with triglycerides in the GOLDN study.

BACKGROUND: There has been significant interest in investigating genome-wide and epigenome-wide associations with lipids. Testing at the gene or region level may improve power in such studies. METHODS: We analyze chromosome 11 cytosine-phosphate-guanine (CpG) methylation levels and single-nucleotide polymorphism (SNP) genotypes from the original Genetics of Lipid Lowering Drugs and Diet Network (GOLDN) study, aiming to explore the association between triglyceride levels and genetic/epigenetic factors. We apply region-based tests of association to methylation and genotype data, in turn, which seek to increase power by reducing the dimension of the gene-region variables. We also investigate whether integrating 2 omics data sets (methylation and genotype) into the triglyceride association analysis helps or hinders detection of candidate gene regions. RESULTS: Gene-region testing identified 1 CpG region that had been previously reported in the GOLDN study data and another 2 gene regions that are also associated with triglyceride levels. Testing on the combined genetic and epigenetic data detected the same genes as using epigenetic or genetic data alone. CONCLUSIONS: Region-based testing can uncover additional association signals beyond those detected using single-variant testing.

Tumoral BRD4 expression in lymph node-negative breast cancer: association with T-bet+ tumor-infiltrating lymphocytes and disease-free survival.

BACKGROUND: We previously observed that T-bet+ tumor-infiltrating T lymphocytes (T-bet+ TILs) in primary breast tumors were associated with adverse clinicopathological features, yet favorable clinical outcome. We identified BRD4 (Bromodomain-Containing Protein 4), a member of the  Bromodomain and Extra Terminal domain (BET) family, as a gene that distinguished T-bet+/high and T-bet-/low tumors. In clinical studies, BET inhibitors have been shown to suppress inflammation in various cancers, suggesting a potential link between BRD4 and immune infiltration in cancer. Hence, we examined the BRD4 expression and clinicopathological features of breast cancer. METHODS: The cohort consisted of a prospectively ascertained consecutive series of women with axillary node-negative breast cancer with long follow-up. Gene expression microarray data were used to detect mRNAs differentially expressed between T-bet+/high (n = 6) and T-bet-/low (n = 41) tumors. Tissue microarrays (TMAs) constructed from tumors of 612 women were used to quantify expression of BRD4 by immunohistochemistry, which was analyzed for its association with T-bet+ TILs, Jagged1, clinicopathological features, and disease-free survival. RESULTS: Microarray analysis indicated that BRD4 mRNA expression was up to 44-fold higher in T-bet+/high tumors compared to T-bet-/low tumors (p = 5.38E-05). Immunohistochemical expression of BRD4 in cancer cells was also shown to be associated with T-bet+ TILs (p = 0.0415) as well as with Jagged1 mRNA and protein expression (p = 0.0171, 0.0010 respectively). BRD4 expression correlated with larger tumor size (p = 0.0049), pre-menopausal status (p = 0.0018), and high Ki-67 proliferative index (p = 0.0009). Women with high tumoral BRD4 expression in the absence of T-bet+ TILs exhibited a significantly poorer outcome (log rank test p = 0.0165) relative to other subgroups. CONCLUSIONS: The association of BRD4 expression with T-bet+ TILs, and T-bet+ TIL-dependent disease-free survival suggests a potential link between BRD4-mediated tumor development and tumor immune surveillance, possibly through BRD4's regulation of Jagged1 signaling pathways. Further understanding BRD4's role in different immune contexts may help to identify an appropriate subset of breast cancer patients who may benefit from BET inhibitors without the risk of diminishing the anti-tumoral immune activity.

Meta-genome-wide association studies identify a locus on chromosome 1 and multiple variants in the MHC region for serum C-peptide in type 1 diabetes.

AIMS/HYPOTHESIS: The aim of this study was to identify genetic variants associated with beta cell function in type 1 diabetes, as measured by serum C-peptide levels, through meta-genome-wide association studies (meta-GWAS). METHODS: We performed a meta-GWAS to combine the results from five studies in type 1 diabetes with cross-sectionally measured stimulated, fasting or random C-peptide levels, including 3479 European participants. The p values across studies were combined, taking into account sample size and direction of effect. We also performed separate meta-GWAS for stimulated (n = 1303), fasting (n = 2019) and random (n = 1497) C-peptide levels. RESULTS: In the meta-GWAS for stimulated/fasting/random C-peptide levels, a SNP on chromosome 1, rs559047 (Chr1:238753916, T>A, minor allele frequency [MAF] 0.24-0.26), was associated with C-peptide (p = 4.13 × 10-8), meeting the genome-wide significance threshold (p < 5 × 10-8). In the same meta-GWAS, a locus in the MHC region (rs9260151) was close to the genome-wide significance threshold (Chr6:29911030, C>T, MAF 0.07-0.10, p = 8.43 × 10-8). In the stimulated C-peptide meta-GWAS, rs61211515 (Chr6:30100975, T/-, MAF 0.17-0.19) in the MHC region was associated with stimulated C-peptide (ß [SE] = - 0.39 [0.07], p = 9.72 × 10-8). rs61211515 was also associated with the rate of stimulated C-peptide decline over time in a subset of individuals (n = 258) with annual repeated measures for up to 6 years (p = 0.02). In the meta-GWAS of random C-peptide, another MHC region, SNP rs3135002 (Chr6:32668439, C>A, MAF 0.02-0.06), was associated with C-peptide (p = 3.49 × 10-8). Conditional analyses suggested that the three identified variants in the MHC region were independent of each other. rs9260151 and rs3135002 have been associated with type 1 diabetes, whereas rs559047 and rs61211515 have not been associated with a risk of developing type 1 diabetes. CONCLUSIONS/INTERPRETATION: We identified a locus on chromosome 1 and multiple variants in the MHC region, at least some of which were distinct from type 1 diabetes risk loci, that were associated with C-peptide, suggesting partly non-overlapping mechanisms for the development and progression of type 1 diabetes. These associations need to be validated in independent populations. Further investigations could provide insights into mechanisms of beta cell loss and opportunities to preserve beta cell function.

Two-phase designs for joint quantitative-trait-dependent and genotype-dependent sampling in post-GWAS regional sequencing.

We evaluate two-phase designs to follow-up findings from genome-wide association study (GWAS) when the cost of regional sequencing in the entire cohort is prohibitive. We develop novel expectation-maximization-based inference under a semiparametric maximum likelihood formulation tailored for post-GWAS inference. A GWAS-SNP (where SNP is single nucleotide polymorphism) serves as a surrogate covariate in inferring association between a sequence variant and a normally distributed quantitative trait (QT). We assess test validity and quantify efficiency and power of joint QT-SNP-dependent sampling and analysis under alternative sample allocations by simulations. Joint allocation balanced on SNP genotype and extreme-QT strata yields significant power improvements compared to marginal QT- or SNP-based allocations. We illustrate the proposed method and evaluate the sensitivity of sample allocation to sampling variation using data from a sequencing study of systolic blood pressure.

A new haplotype block detection method for dense genome sequencing data based on interval graph modeling of clusters of highly correlated SNPs.

Motivation: Linkage disequilibrium (LD) block construction is required for research in population genetics and genetic epidemiology, including specification of sets of single nucleotide polymorphisms (SNPs) for analysis of multi-SNP based association and identification of haplotype blocks in high density sequencing data. Existing methods based on a narrow sense definition do not allow intermediate regions of low LD between strongly associated SNP pairs and tend to split high density SNP data into small blocks having high between-block correlation. Results: We present Big-LD, a block partition method based on interval graph modeling of LD bins which are clusters of strong pairwise LD SNPs, not necessarily physically consecutive. Big-LD uses an agglomerative approach that starts by identifying small communities of SNPs, i.e. the SNPs in each LD bin region, and proceeds by merging these communities. We determine the number of blocks using a method to find maximum-weight independent set. Big-LD produces larger LD blocks compared to existing methods such as MATILDE, Haploview, MIG ++, or S-MIG ++ and the LD blocks better agree with recombination hotspot locations determined by sperm-typing experiments. The observed average runtime of Big-LD for 13 288 240 non-monomorphic SNPs from 1000 Genomes Project autosome data (286 East Asians) is about 5.83 h, which is a significant improvement over the existing methods. Availability and implementation: Source code and documentation are available for download at http://github.com/sunnyeesl/BigLD. Contact: yyoo@snu.ac.kr. Supplementary information: Supplementary data are available at Bioinformatics online.

Model-Free Linkage Analysis of a Binary Trait.

Genetic linkage analysis aims to detect chromosomal regions containing genetic variants that influence risk of specific inherited diseases. The presence of linkage is indicated when a disease or trait cosegregates through the families with genetic markers at a particular region of the genome. Two main types of genetic linkage analysis are in common use, namely model-based linkage analysis and model-free linkage analysis. In this chapter, we focus solely on the latter type and specifically on binary traits or phenotypes, such as the presence or absence of a specific disease. Model-free linkage analysis is based on allele-sharing, where patterns of genetic similarity among affected relatives are compared to chance expectations. Because the model-free methods do not require the specification of the inheritance parameters of a genetic model, they are preferred by many researchers at early stages in the study of a complex disease. We introduce the history of model-free linkage analysis in Subheading 1. Table 1 describes a standard model-free linkage analysis workflow. We describe three popular model-free linkage analysis methods, the nonparametric linkage (NPL) statistic, the affected sib-pair (ASP) likelihood ratio test, and a likelihood approach for pedigrees. The theory behind each linkage test is described in this section together with a simple example of the relevant calculations. Table 4 provides a summary of popular genetic analysis software packages that implement model-free linkage models. In Subheading 2, we work through the methods on a rich example providing sample software code and output. Subheading 3 contains notes with additional details on various topics that may need further consideration during analysis.

Multiple linear combination (MLC) regression tests for common variants adapted to linkage disequilibrium structure.

By jointly analyzing multiple variants within a gene, instead of one at a time, gene-based multiple regression can improve power, robustness, and interpretation in genetic association analysis. We investigate multiple linear combination (MLC) test statistics for analysis of common variants under realistic trait models with linkage disequilibrium (LD) based on HapMap Asian haplotypes. MLC is a directional test that exploits LD structure in a gene to construct clusters of closely correlated variants recoded such that the majority of pairwise correlations are positive. It combines variant effects within the same cluster linearly, and aggregates cluster-specific effects in a quadratic sum of squares and cross-products, producing a test statistic with reduced degrees of freedom (df) equal to the number of clusters. By simulation studies of 1000 genes from across the genome, we demonstrate that MLC is a well-powered and robust choice among existing methods across a broad range of gene structures. Compared to minimum P-value, variance-component, and principal-component methods, the mean power of MLC is never much lower than that of other methods, and can be higher, particularly with multiple causal variants. Moreover, the variation in gene-specific MLC test size and power across 1000 genes is less than that of other methods, suggesting it is a complementary approach for discovery in genome-wide analysis. The cluster construction of the MLC test statistics helps reveal within-gene LD structure, allowing interpretation of clustered variants as haplotypic effects, while multiple regression helps to distinguish direct and indirect associations.

Identification of low frequency and rare variants for hypertension using sparse-data methods.

Availability of genomic sequence data provides opportunities to study the role of low-frequency and rare variants in the etiology of complex disease. In this study, we conduct association analyses of hypertension status in the cohort of 1943 unrelated Mexican Americans provided by Genetic Analysis Workshop 19, focusing on exonic variants in MAP4 on chromosome 3. Our primary interest is to compare the performance of standard and sparse-data approaches for single-variant tests and variant-collapsing tests for sets of rare and low-frequency variants. We analyze both the real and the simulated phenotypes.