On the effect heterogeneity of established disease susceptibility loci for Alzheimer's disease across different genetic ancestries

INTRODUCTION

Genome-wide association studies have identified numerous disease susceptibility loci (DSLs) for Alzheimer's disease (AD). However, only a limited number of studies have investigated the dependence of the genetic effect size of established DSLs on genetic ancestry.

METHODS

We utilized the whole genome sequencing data from the Alzheimer's Disease Sequencing Project (ADSP) including 35,569 participants. A total of 25,459 subjects in four distinct populations (African ancestry, non-Hispanic White, admixed Hispanic, and Asian) were analyzed.

RESULTS

We found that nine DSLs showed significant heterogeneity across populations. Single nucleotide polymorphism (SNP) rs2075650 in translocase of outer mitochondrial membrane 40 (TOMM40) showed the largest heterogeneity (Cochran's Q = 0.00, I2 = 90.08), followed by other SNPs in apolipoprotein C1 (APOC1) and apolipoprotein E (APOE). Two additional loci, signal-induced proliferation-associated 1 like 2 (SIPA1L2) and solute carrier 24 member 4 (SLC24A4), showed significant heterogeneity across populations.

DISCUSSION

We observed substantial heterogeneity for the APOE-harboring 19q13.32 region with TOMM40/APOE/APOC1 genes. The largest risk effect was seen among African Americans, while Asians showed a surprisingly small risk effect.

Gasdermin B, an asthma-susceptibility gene, promotes MAVS-TBK1 signaling and airway inflammation

RATIONALE

Respiratory virus-induced inflammation is the leading cause of asthma exacerbation, frequently accompanied by induction of IFN-stimulated genes (ISGs). How asthma genetic susceptible genes modulate cellular response upon viral infection through fine-tuning ISGs induction and subsequent airway inflammation in genetically susceptible asthmatics remains largely unknown.

OBJECTIVES

To decipher the functions of GSDMB in respiratory virus-induced lung inflammation.

METHODS

In two independent cohorts, we analyzed expression correlation between GSDMB and ISGs. In human bronchial epithelial cell line or primary cells, we generated GSDMB-overexpressing and -deficient cells. A series of qPCR, ELISA and co-immunoprecipitation assays were performed to determine the function and mechanism of GSDMB for ISGs induction. We also generated a novel transgenic mouse line with inducible expression of human unique GSDMB gene in airway epithelial cells and applied respiratory syncytial virus (RSV) infection to determine the role of GSDMB on RSV-induced lung inflammation in vivo.

MEASUREMENTS AND MAIN RESULTS

Gasdermin B encoded by GSDMB, one of the most significant asthma-susceptible genes at 17q21, acts as a novel RNA sensor, promoting MAVS-TBK1 signaling and subsequent inflammation. In airway epithelium, GSDMB is induced by respiratory viral infections. Expression of GSDMB and ISGs significantly correlated in respiratory epithelium from two independent asthma cohorts. Notably, inducible expression of human GSDMB gene in mouse airway epithelium leads to enhanced ISGs induction, increased airway inflammation with mucus hyper-secretion upon RSV infection.

CONCLUSIONS

GSDMB promotes ISGs expression and airway inflammation upon respiratory virus infection, thereby conferring asthma risk in risk allele carriers.

Fast computation of the eigensystem of genomic similarity matrices

The computation of a similarity measure for genomic data is a standard tool in computational genetics. The principal components of such matrices are routinely used to correct for biases due to confounding by population stratification, for instance in linear regressions. However, the calculation of both a similarity matrix and its singular value decomposition (SVD) are computationally intensive. The contribution of this article is threefold. First, we demonstrate that the calculation of three matrices (called the covariance matrix, the weighted Jaccard matrix, and the genomic relationship matrix) can be reformulated in a unified way which allows for the application of a randomized SVD algorithm, which is faster than the traditional computation. The fast SVD algorithm we present is adapted from an existing randomized SVD algorithm and ensures that all computations are carried out in sparse matrix algebra. The algorithm only assumes that row-wise and column-wise subtraction and multiplication of a vector with a sparse matrix is available, an operation that is efficiently implemented in common sparse matrix packages. An exception is the so-called Jaccard matrix, which does not have a structure applicable for the fast SVD algorithm. Second, an approximate Jaccard matrix is introduced to which the fast SVD computation is applicable. Third, we establish guaranteed theoretical bounds on the accuracy (in [Formula: see text] norm and angle) between the principal components of the Jaccard matrix and the ones of our proposed approximation, thus putting the proposed Jaccard approximation on a solid mathematical foundation, and derive the theoretical runtime of our algorithm. We illustrate that the approximation error is low in practice and empirically verify the theoretical runtime scalings on both simulated data and data of the 1000 Genome Project.

A consistent pattern of slide effects in Illumina DNA methylation BeadChip array data

Background: Recent studies have identified thousands of associations between DNA methylation CpGs and complex diseases/traits, emphasizing the critical role of epigenetics in understanding disease aetiology and identifying biomarkers. However, association analyses based on methylation array data are susceptible to batch/slide effects, which can lead to inflated false positive rates or reduced statistical powerResults: We use multiple DNA methylation datasets based on the popular Illumina Infinium MethylationEPIC BeadChip array to describe consistent patterns and the joint distribution of slide effects across CpGs, confirming and extending previous results. The susceptible CpGs overlap with the Illumina Infinium HumanMethylation450 BeadChip array content.Conclusions: Our findings reveal systematic patterns in slide effects. The observations provide further insights into the characteristics of these effects and can improve existing adjustment approaches.

A Polygenic Risk Score for Idiopathic Pulmonary Fibrosis and Interstitial Lung Abnormalities

Rationale: In addition to rare genetic variants and the MUC5B locus, common genetic variants contribute to idiopathic pulmonary fibrosis (IPF) risk. The predictive power of common variants outside the MUC5B locus for IPF and interstitial lung abnormalities (ILAs) is unknown. Objectives: We tested the predictive value of IPF polygenic risk scores (PRSs) with and without the MUC5B region on IPF, ILA, and ILA progression. Methods: We developed PRSs that included (PRS-M5B) and excluded (PRS-NO-M5B) the MUC5B region (500-kb window around rs35705950-T) using an IPF genome-wide association study. We assessed PRS associations with area under the receiver operating characteristic curve (AUC) metrics for IPF, ILA, and ILA progression. Measurements and Main Results: We included 14,650 participants (1,970 IPF; 1,068 ILA) from six multi-ancestry population-based and case-control cohorts. In cases excluded from genome-wide association study, the PRS-M5B (odds ratio [OR] per SD of the score, 3.1; P = 7.1 × 10-95) and PRS-NO-M5B (OR per SD, 2.8; P = 2.5 × 10-87) were associated with IPF. Participants in the top PRS-NO-M5B quintile had ∼sevenfold odds for IPF compared with those in the first quintile. A clinical model predicted IPF (AUC, 0.61); rs35705950-T and PRS-NO-M5B demonstrated higher AUCs (0.73 and 0.7, respectively), and adding both genetic predictors to a clinical model yielded the highest performance (AUC, 0.81). The PRS-NO-M5B was associated with ILA (OR, 1.25) and ILA progression (OR, 1.16) in European ancestry participants. Conclusions: A common genetic variant risk score complements the MUC5B variant to identify individuals at high risk of interstitial lung abnormalities and pulmonary fibrosis.

Benchmarking omics-based prediction of asthma development in children

BACKGROUND

Asthma is a heterogeneous disease with high morbidity. Advancement in high-throughput multi-omics approaches has enabled the collection of molecular assessments at different layers, providing a complementary perspective of complex diseases. Numerous computational methods have been developed for the omics-based patient classification or disease outcome prediction. Yet, a systematic benchmarking of those methods using various combinations of omics data for the prediction of asthma development is still lacking.

OBJECTIVE

We aimed to investigate the computational methods in disease status prediction using multi-omics data.

METHOD

We systematically benchmarked 18 computational methods using all the 63 combinations of six omics data (GWAS, miRNA, mRNA, microbiome, metabolome, DNA methylation) collected in The Vitamin D Antenatal Asthma Reduction Trial (VDAART) cohort. We evaluated each method using standard performance metrics for each of the 63 omics combinations.

RESULTS

Our results indicate that overall Logistic Regression, Multi-Layer Perceptron, and MOGONET display superior performance, and the combination of transcriptional, genomic and microbiome data achieves the best prediction. Moreover, we find that including the clinical data can further improve the prediction performance for some but not all the omics combinations.

CONCLUSIONS

Specific omics combinations can reach the optimal prediction of asthma development in children. And certain computational methods showed superior performance than other methods.

FGF20 and PGM2 variants are associated with childhood asthma in family-based whole-genome sequencing studies

BACKGROUND

Asthma is a heterogeneous common respiratory disease that remains poorly understood. The established genetic associations fail to explain the high estimated heritability, and the prevalence of asthma differs between populations and geographic regions. Robust association analyses incorporating different genetic ancestries and whole-genome sequencing data may identify novel genetic associations.

METHODS

We performed family-based genome-wide association analyses of childhood-onset asthma based on whole-genome sequencing (WGS) data for the 'The Genetic Epidemiology of Asthma in Costa Rica' study (GACRS) and the Childhood Asthma Management Program (CAMP). Based on parent-child trios with children diagnosed with asthma, we performed a single variant analysis using an additive and a recessive genetic model and a region-based association analysis of low-frequency and rare variants.

RESULTS

Based on 1180 asthmatic trios (894 GACRS trios and 286 CAMP trios, a total of 3540 samples with WGS data), we identified three novel genetic loci associated with childhood-onset asthma: rs4832738 on 4p14 ($P=1.72\ast{10}^{-9}$, recessive model), rs1581479 on 8p22 ($P=1.47\ast{10}^{-8}$, additive model) and rs73367537 on 10q26 ($P=1.21\ast{10}^{-8}$, additive model in GACRS only). Integrative analyses suggested potential novel candidate genes underlying these associations: PGM2 on 4p14 and FGF20 on 8p22.

CONCLUSION

Our family-based whole-genome sequencing analysis identified three novel genetic loci for childhood-onset asthma. Gene expression data and integrative analyses point to PGM2 on 4p14 and FGF20 on 8p22 as linked genes. Furthermore, region-based analyses suggest independent potential low-frequency/rare variant associations on 8p22. Follow-up analyses are needed to understand the functional mechanisms and generalizability of these associations.

A comparison between similarity matrices for principal component analysis to assess population stratification in sequenced genetic data sets

Genetic similarity matrices are commonly used to assess population substructure (PS) in genetic studies. Through simulation studies and by the application to whole-genome sequencing (WGS) data, we evaluate the performance of three genetic similarity matrices: the unweighted and weighted Jaccard similarity matrices and the genetic relationship matrix. We describe different scenarios that can create numerical pitfalls and lead to incorrect conclusions in some instances. We consider scenarios in which PS is assessed based on loci that are located across the genome ('globally') and based on loci from a specific genomic region ('locally'). We also compare scenarios in which PS is evaluated based on loci from different minor allele frequency bins: common (>5%), low-frequency (5-0.5%) and rare (<0.5%) single-nucleotide variations (SNVs). Overall, we observe that all approaches provide the best clustering performance when computed based on rare SNVs. The performance of the similarity matrices is very similar for common and low-frequency variants, but for rare variants, the unweighted Jaccard matrix provides preferable clustering features. Based on visual inspection and in terms of standard clustering metrics, its clusters are the densest and the best separated in the principal component analysis of variants with rare SNVs compared with the other methods and different allele frequency cutoffs. In an application, we assessed the role of rare variants on local and global PS, using WGS data from multiethnic Alzheimer's disease data sets and European or East Asian populations from the 1000 Genome Project.

Zinc finger protein 33B demonstrates sex interaction with atopy-related markers in childhood asthma

BACKGROUND

Sex differences related to immune responses can influence atopic manifestations in childhood asthma. While genome-wide association studies have investigated a sex-specific genetic architecture of the immune response, gene-by-sex interactions have not been extensively analysed for atopy-related markers including allergy skin tests, IgE and eosinophils in asthmatic children.

METHODS

We performed a genome-wide gene-by-sex interaction analysis for atopy-related markers using whole-genome sequencing data based on 889 trios from the Genetic Epidemiology of Asthma in Costa Rica Study (GACRS) and 284 trios from the Childhood Asthma Management Program (CAMP). We also tested the findings in UK Biobank participants with self-reported childhood asthma. Furthermore, downstream analyses in GACRS integrated gene expression to disentangle observed associations.

RESULTS

Single nucleotide polymorphism (SNP) rs1255383 at 10q11.21 demonstrated a genome-wide significant gene-by-sex interaction (pinteraction=9.08×10-10) for atopy (positive skin test) with opposite direction of effects between females and males. In the UK Biobank participants with a history of childhood asthma, the signal was consistently observed with the same sex-specific effect directions for high eosinophil count (pinteraction=0.0058). Gene expression of ZNF33B (zinc finger protein 33B), located at 10q11.21, was moderately associated with atopy in girls, but not in boys.

CONCLUSIONS

We report SNPs in/near a zinc finger gene as novel sex-differential loci for atopy-related markers with opposite effect directions in females and males. A potential role for ZNF33B should be studied further as an important driver of sex-divergent features of atopy in childhood asthma.

Unsupervised outlier detection applied to SARS-CoV-2 nucleotide sequences can identify sequences of common variants and other variants of interest

As of June 2022, the GISAID database contains more than 11 million SARS-CoV-2 genomes, including several thousand nucleotide sequences for the most common variants such as delta or omicron. These SARS-CoV-2 strains have been collected from patients around the world since the beginning of the pandemic. We start by assessing the similarity of all pairs of nucleotide sequences using the Jaccard index and principal component analysis. As shown previously in the literature, an unsupervised cluster analysis applied to the SARS-CoV-2 genomes results in clusters of sequences according to certain characteristics such as their strain or their clade. Importantly, we observe that nucleotide sequences of common variants are often outliers in clusters of sequences stemming from variants identified earlier on during the pandemic. Motivated by this finding, we are interested in applying outlier detection to nucleotide sequences. We demonstrate that nucleotide sequences of common variants (such as alpha, delta, or omicron) can be identified solely based on a statistical outlier criterion. We argue that outlier detection might be a useful surveillance tool to identify emerging variants in real time as the pandemic progresses.

Assessing the contribution of rare genetic variants to phenotypes of chronic obstructive pulmonary disease using whole-genome sequence data

RATIONALE

Genetic variation has a substantial contribution to chronic obstructive pulmonary disease (COPD) and lung function measurements. Heritability estimates using genome-wide genotyping data can be biased if analyses do not appropriately account for the nonuniform distribution of genetic effects across the allele frequency and linkage disequilibrium (LD) spectrum. In addition, the contribution of rare variants has been unclear.

OBJECTIVES

We sought to assess the heritability of COPD and lung function using whole-genome sequence data from the Trans-Omics for Precision Medicine program.

METHODS

Using the genome-based restricted maximum likelihood method, we partitioned the genome into bins based on minor allele frequency and LD scores and estimated heritability of COPD, FEV1% predicted and FEV1/FVC ratio in 11 051 European ancestry and 5853 African-American participants.

MEASUREMENTS AND MAIN RESULTS

In European ancestry participants, the estimated heritability of COPD, FEV1% predicted and FEV1/FVC ratio were 35.5%, 55.6% and 32.5%, of which 18.8%, 19.7%, 17.8% were from common variants, and 16.6%, 35.8%, and 14.6% were from rare variants. These estimates had wide confidence intervals, with common variants and some sets of rare variants showing a statistically significant contribution (P-value < 0.05). In African-Americans, common variant heritability was similar to European ancestry participants, but lower sample size precluded calculation of rare variant heritability.

CONCLUSIONS

Our study provides updated and unbiased estimates of heritability for COPD and lung function, and suggests an important contribution of rare variants. Larger studies of more diverse ancestry will improve accuracy of these estimates.

A robust and adaptive framework for interaction testing in quantitative traits between multiple genetic loci and exposure variables

The identification and understanding of gene-environment interactions can provide insights into the pathways and mechanisms underlying complex diseases. However, testing for gene-environment interaction remains a challenge since a.) statistical power is often limited and b.) modeling of environmental effects is nontrivial and such model misspecifications can lead to false positive interaction findings. To address the lack of statistical power, recent methods aim to identify interactions on an aggregated level using, for example, polygenic risk scores. While this strategy can increase the power to detect interactions, identifying contributing genes and pathways is difficult based on these relatively global results. Here, we propose RITSS (Robust Interaction Testing using Sample Splitting), a gene-environment interaction testing framework for quantitative traits that is based on sample splitting and robust test statistics. RITSS can incorporate sets of genetic variants and/or multiple environmental factors. Based on the user's choice of statistical/machine learning approaches, a screening step selects and combines potential interactions into scores with improved interpretability. In the testing step, the application of robust statistics minimizes the susceptibility to main effect misspecifications. Using extensive simulation studies, we demonstrate that RITSS controls the type 1 error rate in a wide range of scenarios, and we show how the screening strategy influences statistical power. In an application to lung function phenotypes and human height in the UK Biobank, RITSS identified highly significant interactions based on subcomponents of genetic risk scores. While the contributing single variant interaction signals are weak, our results indicate interaction patterns that result in strong aggregated effects, providing potential insights into underlying gene-environment interaction mechanisms.

Region-based analysis of rare genomic variants in whole-genome sequencing datasets reveal two novel Alzheimer's disease-associated genes: DTNB and DLG2

Alzheimer's disease (AD) is a genetically complex disease for which nearly 40 loci have now been identified via genome-wide association studies (GWAS). We attempted to identify groups of rare variants (alternate allele frequency <0.01) associated with AD in a region-based, whole-genome sequencing (WGS) association study (rvGWAS) of two independent AD family datasets (NIMH/NIA; 2247 individuals; 605 families). Employing a sliding window approach across the genome, we identified several regions that achieved association p values <10^-6, using the burden test or the SKAT statistic. The genomic region around the dystobrevin beta (DTNB) gene was identified with the burden and SKAT test and replicated in case/control samples from the ADSP study reaching genome-wide significance after meta-analysis (p_meta = 4.74 × 10^-8). SKAT analysis also revealed region-based association around the Discs large homolog 2 (DLG2) gene and replicated in case/control samples from the ADSP study (p_meta = 1 × 10^-6). In conclusion, in a region-based rvGWAS of AD we identified two novel AD genes, DLG2 and DTNB, based on association with rare variants.

Novel recessive locus for body mass index in childhood asthma

Most genome-wide association studies of obesity and body mass index (BMI) have so far assumed an additive mode of inheritance in their analysis, although association testing supports a recessive effect for some of the established loci, for example, rs1421085 in FTO In two whole-genome sequencing (WGS) studies of children with asthma and their parents (892 Costa Rican trios and 286 North American trios), we discovered an association between a locus (rs9292139) in LOC102724122 and BMI that reaches genome-wide significance under a recessive model in the combined analysis. As the association does not achieve significance under an additive model, our finding illustrates the benefits of the recessive model in WGS analyses.

Genome-wide association analysis of COVID-19 mortality risk in SARS-CoV-2 genomes identifies mutation in the SARS-CoV-2 spike protein that colocalizes with P.1 of the Brazilian strain

SARS‐CoV‐2 mortality has been extensively studied in relation to host susceptibility. How sequence variations in the SARS‐CoV‐2 genome affect pathogenicity is poorly understood. Starting in October 2020, using the methodology of genome‐wide association studies (GWAS), we looked at the association between whole‐genome sequencing (WGS) data of the virus and COVID‐19 mortality as a potential method of early identification of highly pathogenic strains to target for containment. Although continuously updating our analysis, in December 2020, we analyzed 7548 single‐stranded SARS‐CoV‐2 genomes of COVID‐19 patients in the GISAID database and associated variants with mortality using a logistic regression. In total, evaluating 29,891 sequenced loci of the viral genome for association with patient/host mortality, two loci, at 12,053 and 25,088 bp, achieved genome‐wide significance (p values of 4.09e−09 and 4.41e−23, respectively), though only 25,088 bp remained significant in follow‐up analyses. Our association findings were exclusively driven by the samples that were submitted from Brazil (p value of 4.90e−13 for 25,088 bp). The mutation frequency of 25,088 bp in the Brazilian samples on GISAID has rapidly increased from about 0.4 in October/December 2020 to 0.77 in March 2021. Although GWAS methodology is suitable for samples in which mutation frequencies varies between geographical regions, it cannot account for mutation frequencies that change rapidly overtime, rendering a GWAS follow‐up analysis of the GISAID samples that have been submitted after December 2020 as invalid. The locus at 25,088 bp is located in the P.1 strain, which later (April 2021) became one of the distinguishing loci (precisely, substitution V1176F) of the Brazilian strain as defined by the Centers for Disease Control. Specifically, the mutations at 25,088 bp occur in the S2 subunit of the SARS‐CoV‐2 spike protein, which plays a key role in viral entry of target host cells. Since the mutations alter amino acid coding sequences, they potentially imposing structural changes that could enhance viral infectivity and symptom severity. Our analysis suggests that GWAS methodology can provide suitable analysis tools for the real‐time detection of new more transmissible and pathogenic viral strains in databases such as GISAID, though new approaches are needed to accommodate rapidly changing mutation frequencies over time, in the presence of simultaneously changing case/control ratios. Improvements of the associated metadata/patient information in terms of quality and availability will also be important to fully utilize the potential of GWAS methodology in this field.

A unifying framework for rare variant association testing in family-based designs, including higher criticism approaches, SKATs, and burden tests

MOTIVATION

Analysis of rare variants in family-based studies remains a challenge. Transmission-based approaches provide robustness against population stratification, but the evaluation of the significance of test statistics based on asymptotic theory can be imprecise. Also, power will depend heavily on the choice of the test statistic and on the underlying genetic architecture of the locus, which will be generally unknown.

RESULTS

In our proposed framework, we utilize the FBAT haplotype algorithm to obtain the conditional offspring genotype distribution under the null hypothesis given the sufficient statistic. Based on this conditional offspring genotype distribution, the significance of virtually any association test statistic can be evaluated based on simulations or exact computations, without the need for asymptotic approximations. Besides standard linear burden-type statistics, this enables our approach to also evaluate other test statistics such as variance components statistics, higher criticism approaches, and maximum-single-variant-statistics, where asymptotic theory might be involved or does not provide accurate approximations for rare variant data. Based on these P-values, combined test statistics such as the aggregated Cauchy association test (ACAT) can also be utilized. In simulation studies, we show that our framework outperforms existing approaches for family-based studies in several scenarios. We also applied our methodology to a TOPMed whole-genome sequencing dataset with 897 asthmatic trios from Costa Rica.

AVAILABILITY AND IMPLEMENTATION

FBAT software is available at https://sites.google.com/view/fbatwebpage. Simulation code is available at https://github.com/julianhecker/FBAT_rare_variant_test_simulations. Whole-genome sequencing data for 'NHLBI TOPMed: The Genetic Epidemiology of Asthma in Costa Rica' is available at https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000988.v4.p1.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A novel locus for exertional dyspnoea in childhood asthma

Most children diagnosed with asthma have respiratory symptoms such as cough, dyspnoea and wheezing, which are also important markers of overall respiratory function. A decade of genome-wide association studies (GWAS) have investigated genetic susceptibility to asthma itself, but few have focused on important respiratory symptoms that characterise childhood asthma.Using whole-genome sequencing (WGS) data for 894 asthmatic trios from a Costa Rican cohort, we performed family-based association tests (FBATs) to assess the association between genetic variants and multiple asthma-relevant respiratory phenotypes: cough, phlegm, wheezing, exertional dyspnoea and exertional chest tightness. We tested whether genome-wide significant associations were replicated in two additional studies: 1) 286 asthmatic trios from the Childhood Asthma Management Program (CAMP), and 2) 2691 African American current or former smokers from the COPDGene study.In the 894 Costa Rican trios, we identified a genome-wide significant association (p=2.16×10^-9) between exertional dyspnoea and the single nucleotide polymorphism (SNP) rs10165869, located on chromosome 2q37.3, that was replicated in the CAMP cohort (p=0.023) with the same direction of association (combined p=3.28×10^-10). This association was not found in the African American participants from COPDGene. We also found suggestive evidence for an association between SNP rs10165869 and the atypical chemokine receptor 3 (ACKR3).Our finding encourages the secondary association analysis of a wider range of phenotypes that characterise respiratory symptoms in other airway diseases/studies.

locStra: Fast analysis of regional/global stratification in whole-genome sequencing studies

locStra is an R -package for the analysis of regional and global population stratification in whole-genome sequencing (WGS) studies, where regional stratification refers to the substructure defined by the loci in a particular region on the genome. Population substructure can be assessed based on the genetic covariance matrix, the genomic relationship matrix, and the unweighted/weighted genetic Jaccard similarity matrix. Using a sliding window approach, the regional similarity matrices are compared with the global ones, based on user-defined window sizes and metrics, for example, the correlation between regional and global eigenvectors. An algorithm for the specification of the window size is provided. As the implementation fully exploits sparse matrix algebra and is written in C++, the analysis is highly efficient. Even on single cores, for realistic study sizes (several thousand subjects, several million rare variants per subject), the runtime for the genome-wide computation of all regional similarity matrices does typically not exceed one hour, enabling an unprecedented investigation of regional stratification across the entire genome. The package is applied to three WGS studies, illustrating the varying patterns of regional substructure across the genome and its beneficial effects on association testing.

The effects of Hurricane Irma on seagrass meadows in previously eutrophic estuaries in Southwest Florida (USA)

In six contiguous estuaries in Southwest Florida (USA) focused management actions over the past several decades have reduced watershed nutrient loads, resulting in an additional 11,672 ha of seagrass meadows between 1999 and 2016, an improvement of 32%. However, in September of 2017, Hurricane Irma made landfall in the state of Florida, affecting the open water and watersheds of each of these six estuaries. In response, seagrass coverage declined by 1203 ha between 2016 and 2018, a system-wide decrease of 3%. The range of decreases associated with Hurricane Irma varied from less than a 1% loss of seagrass coverage in St. Joseph Sound to declines of 7 and 11% in Clearwater Harbor and Lemon Bay, respectively. Areas with the largest losses between 2016 and 2018 were those systems where seagrass coverage had declined in prior years, indicating the effects of Hurricane Irma might have been intensified by prior impacts.

Identification of Novel Alzheimer's Disease Loci Using Sex-Specific Family-Based Association Analysis of Whole-Genome Sequence Data

With the advent of whole genome-sequencing (WGS) studies, family-based designs enable sex-specific analysis approaches that can be applied to only affected individuals; tests using family-based designs are attractive because they are completely robust against the effects of population substructure. These advantages make family-based association tests (FBATs) that use siblings as well as parents especially suited for the analysis of late-onset diseases such as Alzheimer's Disease (AD). However, the application of FBATs to assess sex-specific effects can require additional filtering steps, as sensitivity to sequencing errors is amplified in this type of analysis. Here, we illustrate the implementation of robust analysis approaches and additional filtering steps that can minimize the chances of false positive-findings due to sex-specific sequencing errors. We apply this approach to two family-based AD datasets and identify four novel loci (GRID1, RIOK3, MCPH1, ZBTB7C) showing sex-specific association with AD risk. Following stringent quality control filtering, the strongest candidate is ZBTB7C (Pinter = 1.83 × 10-7), in which the minor allele of rs1944572 confers increased risk for AD in females and protection in males. ZBTB7C encodes the Zinc Finger and BTB Domain Containing 7C, a transcriptional repressor of membrane metalloproteases (MMP). Members of this MMP family were implicated in AD neuropathology.

Genetic Advances in Chronic Obstructive Pulmonary Disease. Insights from COPDGene

Chronic obstructive pulmonary disease (COPD) is a common and progressive disease that is influenced by both genetic and environmental factors. For many years, knowledge of the genetic basis of COPD was limited to Mendelian syndromes, such as alpha-1 antitrypsin deficiency and cutis laxa, caused by rare genetic variants. Over the past decade, the proliferation of genome-wide association studies, the accessibility of whole-genome sequencing, and the development of novel methods for analyzing genetic variation data have led to a substantial increase in the understanding of genetic variants that play a role in COPD susceptibility and COPD-related phenotypes. COPDGene (Genetic Epidemiology of COPD), a multicenter, longitudinal study of over 10,000 current and former cigarette smokers, has been pivotal to these breakthroughs in understanding the genetic basis of COPD. To date, over 20 genetic loci have been convincingly associated with COPD affection status, with additional loci demonstrating association with COPD-related phenotypes such as emphysema, chronic bronchitis, and hypoxemia. In this review, we discuss the contributions of the COPDGene study to the discovery of these genetic associations as well as the ongoing genetic investigations of COPD subtypes, protein biomarkers, and post-genome-wide association study analysis.

Whole-exome sequencing of 81 individuals from 27 multiply affected bipolar disorder families

Bipolar disorder (BD) is a highly heritable neuropsychiatric disease characterized by recurrent episodes of depression and mania. Research suggests that the cumulative impact of common alleles explains 25-38% of phenotypic variance, and that rare variants may contribute to BD susceptibility. To identify rare, high-penetrance susceptibility variants for BD, whole-exome sequencing (WES) was performed in three affected individuals from each of 27 multiply affected families from Spain and Germany. WES identified 378 rare, non-synonymous, and potentially functional variants. These spanned 368 genes, and were carried by all three affected members in at least one family. Eight of the 368 genes harbored rare variants that were implicated in at least two independent families. In an extended segregation analysis involving additional family members, five of these eight genes harbored variants showing full or nearly full cosegregation with BD. These included the brain-expressed genes RGS12 and NCKAP5, which were considered the most promising BD candidates on the basis of independent evidence. Gene enrichment analysis for all 368 genes revealed significant enrichment for four pathways, including genes reported in de novo studies of autism (p_adj < 0.006) and schizophrenia (p_adj = 0.015). These results suggest a possible genetic overlap with BD for autism and schizophrenia at the rare-sequence-variant level. The present study implicates novel candidate genes for BD development, and may contribute to an improved understanding of the biological basis of this common and often devastating disease.

Whole Genome Sequencing Identifies CRISPLD2 as a Lung Function Gene in Children With Asthma

BACKGROUND

Asthma is a common respiratory disorder with a highly heterogeneous nature that remains poorly understood. The objective was to use whole genome sequencing (WGS) data to identify regions of common genetic variation contributing to lung function in individuals with a diagnosis of asthma.

METHODS

WGS data were generated for 1,053 individuals from trios and extended pedigrees participating in the family-based Genetic Epidemiology of Asthma in Costa Rica study. Asthma affection status was defined through a physician's diagnosis of asthma, and most participants with asthma also had airway hyperresponsiveness (AHR) to methacholine. Family-based association tests for single variants were performed to assess the associations with lung function phenotypes.

RESULTS

A genome-wide significant association was identified between baseline FEV1/FVC ratio and a single-nucleotide polymorphism in the top hit cysteine-rich secretory protein LCCL domain-containing 2 (CRISPLD2) (rs12051168; P = 3.6 × 10-8 in the unadjusted model) that retained suggestive significance in the covariate-adjusted model (P = 5.6 × 10-6). Rs12051168 was also nominally associated with other related phenotypes: baseline FEV1 (P = 3.3 × 10-3), postbronchodilator (PB) FEV1 (7.3 × 10-3), and PB FEV1/FVC ratio (P = 2.7 × 10-3). The identified baseline FEV1/FVC ratio and rs12051168 association was meta-analyzed and replicated in three independent cohorts in which most participants with asthma also had confirmed AHR (combined weighted z-score P = .015) but not in cohorts without information about AHR.

CONCLUSIONS

These findings suggest that using specific asthma characteristics, such as AHR, can help identify more genetically homogeneous asthma subgroups with genotype-phenotype associations that may not be observed in all children with asthma. CRISPLD2 also may be important for baseline lung function in individuals with asthma who also may have AHR.

Shared Genetic Etiology of Obesity-Related Traits and Barrett's Esophagus/Adenocarcinoma: Insights from Genome-Wide Association Studies

BACKGROUND

Obesity is a major risk factor for esophageal adenocarcinoma (EA) and its precursor Barrett's esophagus (BE). Research suggests that individuals with high genetic risk to obesity have a higher BE/EA risk. To facilitate understanding of biological factors that lead to progression from BE to EA, the present study investigated the shared genetic background of BE/EA and obesity-related traits.

METHODS

Cross-trait linkage disequilibrium score regression was applied to summary statistics from genome-wide association meta-analyses on BE/EA and on obesity traits. Body mass index (BMI) was used as a proxy for general obesity, and waist-to-hip ratio (WHR) for abdominal obesity. For single marker analyses, all genome-wide significant risk alleles for BMI and WHR were compared with summary statistics of the BE/EA meta-analyses.

RESULTS

Sex-combined analyses revealed a significant genetic correlation between BMI and BE/EA (r_g = 0.13, P = 2 × 10^-04) and a r_g of 0.12 between WHR and BE/EA (P = 1 × 10^-02). Sex-specific analyses revealed a pronounced genetic correlation between BMI and EA in females (r_g = 0.17, P = 1.2 × 10^-03), and WHR and EA in males (r_g = 0.18, P = 1.51 × 10^-02). On the single marker level, significant enrichment of concordant effects was observed for BMI and BE/EA risk variants (P = 8.45 × 10^-03) and WHR and BE/EA risk variants (P = 2 × 10^-02).

CONCLUSIONS

Our study provides evidence for sex-specific genetic correlations that might reflect specific biological mecha-nisms. The data demonstrate that shared genetic factors are particularly relevant in progression from BE to EA.

IMPACT

Our study quantifies the genetic correlation between BE/EA and obesity. Further research is now warranted to elucidate these effects and to understand the shared pathophysiology.

A flexible and nearly optimal sequential testing approach to randomized testing: QUICK-STOP

In the analysis of current life science datasets, we often encounter scenarios in which the application of asymptotic theory to hypothesis testing can be problematic. Besides improved asymptotic results, permutation/simulation-based tests are a general approach to address this issue. However, these randomized tests can impose a massive computational burden, for example, in scenarios in which large numbers of statistical tests are computed, and the specified significance level is very small. Stopping rules aim to assess significance with the smallest possible number of draws while controlling the probabilities of errors due to statistical uncertainty. In this communication, we derive a general stopping rule, QUICK-STOP, based on the sequential testing theory that is easy to implement, controls the error probabilities rigorously, and is nearly optimal in terms of expected draws. In a simulation study, we show that our approach outperforms current stopping approaches for general randomized tests by factor 10 and does not impose an additional computational burden. We illustrate our approach by applying our stopping rule to a single-variant analysis of a whole-genome sequencing study for lung function.

Evaluation of single-view contrast-enhanced mammography as novel reading strategy: a non-inferiority feasibility study

BACKGROUND

Guidelines recommend screening of high-risk women using breast magnetic resonance imaging (MRI). Contrast-enhanced mammography (CEM) has matured, providing excellent diagnostic accuracy. To lower total radiation dose, evaluation of single-view (1 V) CEM exams might be considered instead of double-view (2 V) readings as an alternative reading strategy in women who cannot undergo MRI.

METHODS

This retrospective non-inferiority feasibility study evaluates whether the use of 1 V results in an acceptable sensitivity for detecting breast cancer (non-inferiority margin, - 10%). CEM images from May 2013 to December 2017 were included. 1 V readings were performed by consensus opinion of three radiologists, followed by 2 V readings being performed after 6 weeks. Cases were considered "malignant" if the final BI-RADS score was ≥ 4, enabling calculation of sensitivity, specificity, and area under the receiver operating characteristic curve (AUC). Histopathological results or follow-up served as a gold standard.

RESULTS

A total of 368 cases were evaluated. Mean follow-up for benign or negative cases was 20.9 months. Sensitivity decreased by 9.6% from 92.9 to 83.3% when only 1 V was used for evaluation (p < 0.001). The lower limit of the 90% confidence interval around the difference in sensitivity between 1 V and 2 V readings was - 15% and lies below the predefined non-inferiority margin of - 10%. Hence, non-inferiority of 1 V to 2 V reading cannot be concluded. AUC for 1 V was significantly lower, 0.861 versus 0.899 for 2 V (p = 0.0174).

CONCLUSION

Non-inferiority of 1 V evaluations as an alternative reading strategy to standard 2 V evaluations could not be concluded. 1 V evaluations had lower diagnostic performance compared with 2 V evaluations.

KEY POINTS

• To lower radiation exposure used in contrast-enhanced mammography, we studied a hypothetical alternative strategy: single-view readings (1 V) versus (standard) double-view readings (2 V). • Based on our predefined margin of - 10%, non-inferiority of 1 V could not be concluded. • 1 V evaluation is not recommended as an alternative reading strategy to lower CEM-related radiation exposure.

Effect of population stratification on SNP-by-environment interaction

Proportions of false-positive rates in genome-wide association analysis are affected by population stratification, and if it is not correctly adjusted, the statistical analysis can produce the large false-negative finding. Therefore various approaches have been proposed to adjust such problems in genome-wide association studies. However, in spite of its importance, a few studies have been conducted in genome-wide single nucleotide polymorphism (SNP)-by-environment interaction studies. In this report, we illustrate in which scenarios can lead to the false-positive rates in association mapping and approach to maintaining the overall type-1 error rate.

PolyGEE: a generalized estimating equation approach to the efficient and robust estimation of polygenic effects in large-scale association studies

To quantify polygenic effects, i.e. undetected genetic effects, in large-scale association studies, we propose a generalized estimating equation (GEE) based estimation framework. We develop a marginal model for single-variant association test statistics of complex diseases that generalizes existing approaches such as LD Score regression and that is applicable to population-based designs, to family-based designs or to arbitrary combinations of both. We extend the standard GEE approach so that the parameters of the proposed marginal model can be estimated based on working-correlation/linkage-disequilibrium (LD) matrices from external reference panels. Our method achieves substantial efficiency gains over standard approaches, while it is robust against misspecification of the LD structure, i.e. the LD structure of the reference panel can differ substantially from the true LD structure in the study population. In simulation studies and in applications to population-based and family-based studies, we illustrate the features of the proposed GEE framework. Our results suggest that our approach can be up to 100% more efficient than existing methodology.

A comparison of popular TDT-generalizations for family-based association analysis

The transmission disequilibrium test (TDT) is the gold standard for testing the association between a genetic variant and disease in samples consisting of affected individuals and their parents. In practice, more complex pedigree structures, that is siblings with no parents, or three-generational pedigrees with possibly missing genotypes, are common. There are several generalizations of the TDT that are suitable for use with arbitrary pedigree structures. We consider three such frequently used generalizations, family-based association test, pedigree disequilibrium test, and generalized disequilibrium test, that have accompanying software and compare them regarding validity and power in the single variant setting. We use simulations to study the effects of population admixture, populations whose genotypes are not in Hardy-Weinberg equilibrium (HWE), different pedigree structures, and the presence of linkage. Whereas our results show that some TDT generalizations can have a substantially increased Type 1 error, these tests are often used in substantive research without caveats about the validity of their Type 1 error. For the association analysis of rare variants in sequencing studies, region-based extensions of the TDT generalizations, that rely on the postulated robustness of the single variant tests, have been proposed. We discuss the implications of our results for these region-based extensions.

Exome sequencing in large, multiplex bipolar disorder families from Cuba

Bipolar disorder (BD) is a major psychiatric illness affecting around 1% of the global population. BD is characterized by recurrent manic and depressive episodes, and has an estimated heritability of around 70%. Research has identified the first BD susceptibility genes. However, the underlying pathways and regulatory networks remain largely unknown. Research suggests that the cumulative impact of common alleles with small effects explains only around 25-38% of the phenotypic variance for BD. A plausible hypothesis therefore is that rare, high penetrance variants may contribute to BD risk. The present study investigated the role of rare, nonsynonymous, and potentially functional variants via whole exome sequencing in 15 BD cases from two large, multiply affected families from Cuba. The high prevalence of BD in these pedigrees renders them promising in terms of the identification of genetic risk variants with large effect sizes. In addition, SNP array data were used to calculate polygenic risk scores for affected and unaffected family members. After correction for multiple testing, no significant increase in polygenic risk scores for common, BD-associated genetic variants was found in BD cases compared to healthy relatives. Exome sequencing identified a total of 17 rare and potentially damaging variants in 17 genes. The identified variants were shared by all investigated BD cases in the respective pedigree. The most promising variant was located in the gene SERPING1 (p.L349F), which has been reported previously as a genome-wide significant risk gene for schizophrenia. The present data suggest novel candidate genes for BD susceptibility, and may facilitate the discovery of disease-relevant pathways and regulatory networks.

Widespread recovery of seagrass coverage in Southwest Florida (USA): Temporal and spatial trends and management actions responsible for success

In Southwest Florida, a variety of human impacts had caused widespread losses of seagrass coverage from historical conditions. St. Joseph Sound and Clearwater Harbor lost approximately 24 and 51%, respectively, of their seagrass coverage between 1950 and 1999, while Tampa Bay and Sarasota Bay had lost 46% and 15%, respectively, of their seagrass coverage between 1950 and the 1980s. However, over the period of 1999 to 2016, the largest of the six estuaries, Tampa Bay, added 408 ha of seagrass per year, while the remaining five estuaries examined in this paper added approximately 269 ha per year. In total, seagrass coverage in these six estuaries increased 12,171 ha between the 1980s and 2016. Focused resource management plans have held the line on nitrogen loads from non-point sources, allowing seagrass resources to expand in response to reductions in point source loads that have been implemented over the past few decades.

Analysis of the Influence of microRNAs in Lithium Response in Bipolar Disorder

Bipolar disorder (BD) is a common, highly heritable neuropsychiatric disease characterized by recurrent episodes of mania and depression. Lithium is the best-established long-term treatment for BD, even though individual response is highly variable. Evidence suggests that some of this variability has a genetic basis. This is supported by the largest genome-wide association study (GWAS) of lithium response to date conducted by the International Consortium on Lithium Genetics (ConLiGen). Recently, we performed the first genome-wide analysis of the involvement of miRNAs in BD and identified nine BD-associated miRNAs. However, it is unknown whether these miRNAs are also associated with lithium response in BD. In the present study, we therefore tested whether common variants at these nine candidate miRNAs contribute to the variance in lithium response in BD. Furthermore, we systematically analyzed whether any other miRNA in the genome is implicated in the response to lithium. For this purpose, we performed gene-based tests for all known miRNA coding genes in the ConLiGen GWAS dataset (n = 2,563 patients) using a set-based testing approach adapted from the versatile gene-based test for GWAS (VEGAS2). In the candidate approach, miR-499a showed a nominally significant association with lithium response, providing some evidence for involvement in both development and treatment of BD. In the genome-wide miRNA analysis, 71 miRNAs showed nominally significant associations with the dichotomous phenotype and 106 with the continuous trait for treatment response. A total of 15 miRNAs revealed nominal significance in both phenotypes with miR-633 showing the strongest association with the continuous trait (p = 9.80E-04) and miR-607 with the dichotomous phenotype (p = 5.79E-04). No association between miRNAs and treatment response to lithium in BD in either of the tested conditions withstood multiple testing correction. Given the limited power of our study, the investigation of miRNAs in larger GWAS samples of BD and lithium response is warranted.

Family-based tests for associating haplotypes with general phenotype data: Improving the FBAT-haplotype algorithm

For family-based association studies, Horvath et al. proposed an algorithm for the association analysis between haplotypes and arbitrary phenotypes when the phase of the haplotypes is unknown, that is, genotype data is given. Their approach to haplotype analysis maintains the original features of the TDT/FBAT-approach, that is, complete robustness against genetic confounding and misspecification of the phenotype. The algorithm has been implemented in the FBAT and PBAT software package and has been used in numerous substantive manuscripts. Here, we propose a simplification of the original algorithm that maintains the original approach but reduces the computational burden of the approach substantially and gives valuable insights regarding the conditional distribution. With the modified algorithm, the application to whole-genome sequencing (WGS) studies becomes feasible; for example, in sliding window approaches or spatial-clustering approaches. The reduction of the computational burden that our modification provides is especially dramatic when both parental genotypes are missing. For example, for eight variants and 441 nuclear families with mostly offspring-only families, in a WGS study at the APOE locus, the running time decreased from approximately 21 hr for the original algorithm to 0.11 sec after our modification.

On the association analysis of genome-sequencing data: A spatial clustering approach for partitioning the entire genome into nonoverlapping windows

For the association analysis of whole-genome sequencing (WGS) studies, we propose an efficient and fast spatial-clustering algorithm. Compared to existing analysis approaches for WGS data, that define the tested regions either by sliding or consecutive windows of fixed sizes along variants, a meaningful grouping of nearby variants into consecutive regions has the advantage that, compared to sliding window approaches, the number of tested regions is likely to be smaller. In comparison to consecutive, fixed-window approaches, our approach is likely to group nearby variants together. Given existing biological evidence that disease-associated mutations tend to physically cluster in specific regions along the chromosome, the identification of meaningful groups of nearby located variants could thus lead to a potential power gain for association analysis. Our algorithm defines consecutive genomic regions based on the physical positions of the variants, assuming an inhomogeneous Poisson process and groups together nearby variants. As parameters are estimated locally, the algorithm takes the differing variant density along the chromosome into account and provides locally optimal partitioning of variants into consecutive regions. An R-implementation of the algorithm is provided. We discuss the theoretical advances of our algorithm compared to existing, window-based approaches and show the performance and advantage of our introduced algorithm in a simulation study and by an application to Alzheimer's disease WGS data. Our analysis identifies a region in the ITGB3 gene that potentially harbors disease susceptibility loci for Alzheimer's disease. The region-based association signal of ITGB3 replicates in an independent data set and achieves formally genome-wide significance. Software Implementation: An implementation of the algorithm in R is available at: https://github.com/heidefier/cluster_wgs_data.

Identification of shared risk loci and pathways for bipolar disorder and schizophrenia

Bipolar disorder (BD) is a highly heritable neuropsychiatric disease characterized by recurrent episodes of mania and depression. BD shows substantial clinical and genetic overlap with other psychiatric disorders, in particular schizophrenia (SCZ). The genes underlying this etiological overlap remain largely unknown. A recent SCZ genome wide association study (GWAS) by the Psychiatric Genomics Consortium identified 128 independent genome-wide significant single nucleotide polymorphisms (SNPs). The present study investigated whether these SCZ-associated SNPs also contribute to BD development through the performance of association testing in a large BD GWAS dataset (9747 patients, 14278 controls). After re-imputation and correction for sample overlap, 22 of 107 investigated SCZ SNPs showed nominal association with BD. The number of shared SCZ-BD SNPs was significantly higher than expected (p = 1.46x10-8). This provides further evidence that SCZ-associated loci contribute to the development of BD. Two SNPs remained significant after Bonferroni correction. The most strongly associated SNP was located near TRANK1, which is a reported genome-wide significant risk gene for BD. Pathway analyses for all shared SCZ-BD SNPs revealed 25 nominally enriched gene-sets, which showed partial overlap in terms of the underlying genes. The enriched gene-sets included calcium- and glutamate signaling, neuropathic pain signaling in dorsal horn neurons, and calmodulin binding. The present data provide further insights into shared risk loci and disease-associated pathways for BD and SCZ. This may suggest new research directions for the treatment and prevention of these two major psychiatric disorders.

Analysis of Rare Variants in the Alcohol Dependence Candidate Gene GATA4

Background

Common variants in the gene GATA binding protein 4 (GATA4) show association with alcohol dependence (AD). The aim of this study was to identify rare variants in GATA4 in order to elucidate the role of this gene in AD susceptibility. Identification of rare variants may provide a more complete picture of the allelic architecture at this risk locus.

Methods

Sanger sequencing of all 6 coding exons of GATA4 was performed in 528 patients and 517 controls. Four in silico prediction tools were used to determine the effect of a DNA variant on the amino acid sequence and protein function. Five variants were included in the replication step. Of these, 4 were successfully genotyped in our replication cohort of 655 patients and 1,501 controls. All patients fulfilled DSM‐IV criteria for AD, and all individuals were of German descent.

Results

In the discovery step, 19 different heterozygous variants were identified. Four patient‐specific and potentially functionally relevant variants were followed up. Only the variant S379S (c.1137C>T) remained patient specific (1/1,166 patients vs. 0/1,997 controls). None of the variants showed a statistically significant association with AD.

Conclusions

The present study elucidated the role of GATA4 in AD susceptibility by identifying rare variants via Sanger sequencing and subsequent replication. Although novel patient‐specific rare variants of GATA4 were identified, none received support in the independent replication step. However, given previous robust findings of association with common variants, GATA4 remains a promising candidate gene for AD.

Augmentation of Transient Donor Cell Chimerism and Alloantigen-Specific Regulation of Lung Transplants in Miniature Swine

Donor alloantigen infusion induces T cell regulation and transplant tolerance in small animals. Here, we study donor splenocyte infusion in a large animal model of pulmonary transplantation. Major histocompatibility complex-mismatched single lung transplantation was performed in 28 minipigs followed by a 28-day course of methylprednisolone and tacrolimus. Some animals received a perioperative donor or third party splenocyte infusion, with or without low-dose irradiation (IRR) before surgery. Graft survival was significantly prolonged in animals receiving both donor splenocytes and IRR compared with controls with either donor splenocytes or IRR only. In animals with donor splenocytes and IRR, increased donor cell chimerism and CD4(+) CD25(high+) T cell frequencies were detected in peripheral blood associated with decreased interferon-γ production of leukocytes. Secondary third-party kidney transplants more than 2 years after pulmonary transplantation were acutely rejected despite maintained tolerance of the lung allografts. As a cellular control, additional animals received third-party splenocytes or donor splenocyte protein extracts. While animals treated with third-party splenocytes showed significant graft survival prolongation, the subcellular antigen infusion showed no such effect. In conclusion, minipigs conditioned with preoperative IRR and donor, or third-party, splenocyte infusions may develop long-term donor-specific pulmonary allograft survival in the presence of high levels of circulating regulatory T cells.

Utilizing the Jaccard index to reveal population stratification in sequencing data: a simulation study and an application to the 1000 Genomes Project

MOTIVATION

Population stratification is one of the major sources of confounding in genetic association studies, potentially causing false-positive and false-negative results. Here, we present a novel approach for the identification of population substructure in high-density genotyping data/next generation sequencing data. The approach exploits the co-appearances of rare genetic variants in individuals. The method can be applied to all available genetic loci and is computationally fast. Using sequencing data from the 1000 Genomes Project, the features of the approach are illustrated and compared to existing methodology (i.e. EIGENSTRAT). We examine the effects of different cutoffs for the minor allele frequency on the performance of the approach. We find that our approach works particularly well for genetic loci with very small minor allele frequencies. The results suggest that the inclusion of rare-variant data/sequencing data in our approach provides a much higher resolution picture of population substructure than it can be obtained with existing methodology. Furthermore, in simulation studies, we find scenarios where our method was able to control the type 1 error more precisely and showed higher power.

AVAILABILITY AND IMPLEMENTATION

CONTACT

dmitry.prokopenko@uni-bonn.de

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

On the Recombination Rate Estimation in the Presence of Population Substructure

As recombination events are not uniformly distributed along the human genome, the estimation of fine-scale recombination maps, e.g. HapMap Project, has been one of the major research endeavors over the last couple of years. For simulation studies, these estimates provide realistic reference scenarios to design future study and to develop novel methodology. To achieve a feasible framework for the estimation of such recombination maps, existing methodology uses sample probabilities for a two-locus model with recombination, with recent advances allowing for computationally fast implementations. In this work, we extend the existing theoretical framework for the recombination rate estimation to the presence of population substructure. We show under which assumptions the existing methodology can still be applied. We illustrate our extension of the methodology by an extensive simulation study.

Using Network Methodology to Infer Population Substructure

One of the main caveats of association studies is the possible affection by bias due to population stratification. Existing methods rely on model-based approaches like structure and ADMIXTURE or on principal component analysis like EIGENSTRAT. Here we provide a novel visualization technique and describe the problem of population substructure from a graph-theoretical point of view. We group the sequenced individuals into triads, which depict the relational structure, on the basis of a predefined pairwise similarity measure. We then merge the triads into a network and apply community detection algorithms in order to identify homogeneous subgroups or communities, which can further be incorporated as covariates into logistic regression. We apply our method to populations from different continents in the 1000 Genomes Project and evaluate the type 1 error based on the empirical p-values. The application to 1000 Genomes data suggests that the network approach provides a very fine resolution of the underlying ancestral population structure. Besides we show in simulations, that in the presence of discrete population structures, our developed approach maintains the type 1 error more precisely than existing approaches.

An IgE epitope of Bet v 1 and fagales PR10 proteins as defined by a human monoclonal IgE

BACKGROUND

Analyses of the molecular basis underlying allergenicity and allergen cross-reactivity, as well as improvement of allergy diagnostics and therapeutics, are hampered by the lack of human monoclonal IgE antibodies and knowledge about their epitopes. Here, we addressed the consecutive generation and epitope delineation of a human monoclonal IgE against the prototypic allergen Bet v 1.

METHODS

Phage-display scFv hybrid libraries of allergic donor-derived VH epsilon and synthetic VL were established from 107 mononuclear cells. An obtained scFv was converted into human immunoglobulin formats including IgE. Using variants of Bet v 1, the epitope of the antibody was mapped and extrapolated to other PR10 proteins.

RESULTS

The obtained antibodies exhibited pronounced reactivity with Bet v 1, but were not reactive with the homologous PR10 protein Mal d 1. The epitope as defined by the IgE paratope and a set of chimeric Bet v 1 fusion proteins and fragments could be assigned to a C-terminal helix-structured motif comprised by amino acid residues 132-154, including the critical residue E149. Grafting this motif re-established the reactivity of the per se nonreactive Mal d 1 framework. Cross-reactivities predicted by primary structure analyses of different isoforms and PR10 proteins were verified by allergen chip-based analyses.

CONCLUSIONS

The obtained results demonstrate that hybrid IgE repertoires represent a source for human antibodies with genuine paratopes. The IgE-derived information about the IgE epitope nature of Bet v 1 and homologues allows for detailed insights into molecular aspects of allergenicity and cross-reactivity within the PR10 protein family.

Ageing of intravenous glucose solutions

The pH and acid content of 63 glucose-based solutions was measured up to eight years after their expiry. The pH ranged from 2.27 to 4.33. Newly manufactured solutions with 5 per cent or less of glucose which were packaged in plastic bags contained in the order of 0.01 mEq/L of acid per gram of glucose. Solutions with higher glucose concentrations had lower pH and greater acid. As solutions in bags became older, pH decreased and acid increased. In contrast, 25 per cent and 50 per cent glucose solutions in glass containers were much less acid and showed only slight changes with age in pH or acid. The greater rate of acid increase in glucose solutions in plastic bags may be due to the bags being permeable to oxygen and thus allowing slow oxidation of the glucose.

Economic evaluation of donepezil in moderate to severe Alzheimer disease

OBJECTIVE

To investigate the costs to society of Alzheimer disease (AD) care in a multinational, randomized, placebo-controlled trial of donepezil in patients with moderate to severe AD.

METHODS

A total of 290 patients with AD (screening standardized Mini-Mental State Examination score 5 to 17) were randomized to receive either donepezil (n = 144; 5 mg/day for 28 days, followed by 10 mg/day as per clinician's judgment) or placebo (n = 146) for 24 weeks. The authors collected data on patient and caregiver health resource utilization prospectively using the Canadian Utilization of Services Tracking questionnaire. Costs were calculated for patients and caregivers in each group based on resource utilization multiplied by the unit prices for each resource. A cost (the average Ontario minimum wage for 1998 [Can 6.85 dollars per hour]) was assigned to unpaid time that caregivers spent assisting the patient with activities of daily living (ADL).

RESULTS

Patient and caregiver demographics at baseline were similar across the two groups. After adjusting for baseline total cost per patient, the mean total societal cost per patient for the 24-week period was donepezil, Can 9,904 dollars (US 6,686 dollars) and placebo, Can 10,236 dollars (US 6,910 dollars). This net cost saving of Can 332 dollars (US 224 dollars) included the average 24-week cost of donepezil treatment. Most of the cost-saving with donepezil treatment was due to less use of residential care by patients, and caregivers spending less time assisting patients with ADL.

CONCLUSION

This cost-consequence analysis reveals economic benefits of treatment of moderate to severe AD with donepezil.