Rank-based genome-wide analysis reveals the association of ryanodine receptor-2 gene variants with childhood asthma among human populations

PDE4D gene variants and haplotypes are associated with asthma and atopy in Brazilian children

PDE4D (Phosphodiesterase 4D) gene encodes a hydrolase of cyclic AMP. PDE4D genetic variants have been associated with asthma susceptibility. Therefore, this study aimed to investigate the association between PDE4D variants (and haplotypes) with asthma and atopy in a Brazilian population. The study comprised 1,246 unrelated participants from the SCAALA (Social Changes Asthma and Allergy in Latin America) program. Genotyping was performed using the Illumina 2.5 Human Omni bead chip. Multivariate logistic regression was used to investigate the association between PDE4D variants and asthma/atopy phenotypes in PLINK 1.09 software. Twenty-four SNVs in PDE4D were associated with atopy or asthma. The rs6898082 (A) variant increased asthma susceptibility (OR 2.76; CI 99% 1.26-6.03) and was also related to a greater PDE4D expression in the GTEx database. Also, the variant rs6870632 was further associated with asthma in meta-analysis with a replication cohort. In addition, the variants rs75699812 (C), rs8007656 (G), and rs958851 (T) were positively associated with atopy. Moreover, these variants formed an atopy risk haplotype (OR 1.82; CI 99% 1.15-2.88). Also, these variants were related to lower levels of IL-10. Functional in silico assessment showed that some PDE4D SNVs may have an impact on gene regulation and expression. Variants in the PDE4D are positively associated with asthma and allergy markers. It is possible that these variants lead to alteration in PDE4D expression and therefore impact immunity and pulmonary function.

Epigenome alterations in food allergy: A systematic review of candidate gene and epigenome-wide association studies

OBJECTIVE

The aim of this study was to systematically review the evidence across studies that assessed DNA methylome variations in association with food allergy (FA).

DESIGN

A systematic review of the literature and meta-analysis were carried out within several databases. However, the risk of bias in the included articles was not evaluated.

DATA SOURCES

PubMed, Cochrane Database of Systematic Reviews, and Web of Science were used to search up to July 2022.

ELIGIBILITY CRITERIA

We included targeted and epigenome-wide association studies (EWASs) that assessed DNA methylome alterations in association with FA in adult or paediatric populations.

RESULTS

Among 366 publications, only 16 were retained, which were mainly focused on FA in children. Seven candidate gene-targeted studies found associations in Th1/Th2 imbalance (IL4, IL5, IL10, INFG, IL2 and IL12B genes), regulatory T cell function (FOXP3 gene), Toll-like receptors pathway (TLR2, CD14 genes) and digestive barrier integrity (FLG gene). Nine EWAS assessed the association with peanut allergy (n = 3), cow's milk allergy (n = 2) or various food allergens (n = 4). They highlighted 11 differentially methylated loci in at least two studies (RPS6KA2, CAMTA1, CTBP2, RYR2, TRAPPC9, DOCK1, GALNTL4, HDAC4, UMODL1, ZAK and TNS3 genes). Among them, CAMTA1 and RPS6KA2, and CTBP2 are involved in regulatory T cell function and Th2 cell differentiation, respectively. Gene-functional analysis revealed two enriched gene clusters involved in immune responses and protein phosphorylation. ChIP-X Enrichment Analysis 3 showed eight significant transcription factors (RXRA, ZBTB7A, ESR1, TCF3, MYOD1, CTCF, GATA3 and CBX2). Ingenuity Pathway Analysis identified canonical pathways involved, among other, in B cell development, pathogen-induced cytokine storm signalling pathway and dendritic cell maturation.

CONCLUSION

This review highlights the involvement of epigenomic alterations of loci in Th1/Th2 and regulatory T cell differentiation in both candidate gene studies and EWAS. These alterations provide a better insight into the mechanistic aspects in FA pathogenesis and may guide the development of epigenome-based biomarkers for FA.

Seasonal characteristics of temperature variability impacts on childhood asthma hospitalization in Hefei, China: Does PM2.5 modify the association?

OBJECTIVES

Evidence of childhood asthma hospitalizations associated with temperature variability (TV) and the attributable risk are limited in China. We aim to use a comprehensive index that reflected both intra- and inter-day TV to assess the TV-childhood asthma relationship and disease burden, further to identify seasonality vulnerable populations, and to explore the effect modification of PM_2.5.

METHODS

A quasi-distributed lagged nonlinear model (DLNM) combined with a linear threshold function was applied to estimate the association between TV and childhood asthma hospitalizations during 2013-2016 in Hefei, China. Subgroup analysis was conducted by age and sex. Disease burden is reflected by the attributable fraction and attributable number. Besides, modifications of PM_2.5 were tested by introducing the cross-basis of TV and binary PM_2.5 as an interaction term.

RESULTS

The risk estimates peaked at TV_0-3 and TV_0-4 in the cool and the warm season separately, with RR of 1.051 (95%CI: 1.021-1.081) and 1.072 (95%CI: 1.008-1.125), and the effects lasted longer in the cool season. The school-age children in the warm season and all subgroups except pre-school children in the cool season were vulnerable to TV. It is estimated that the disease burden related to TV account for 6.2% (95% CI: 2.7%-9.4%) and 4% (95% CI: 0.6%-7.1%) during the cool and warm seasons in TV_0-3. In addition, the risks of TV were higher under the high PM_2.5 level compared with the low PM_2.5 level in the cool season, although no significant differences between them.

CONCLUSIONS

TV exposure significantly increases the risk and disease burden of childhood asthma hospitalizations, especially in the cool season. More medical resources should be allocated to school-age children. Giving priority to pay attention to TV in the cool season in practice could obtain the greatest public health benefits and those days with high TV and high PM_2.5 need more attention.

Unlocking the secrets of long non-coding RNAs in asthma

Asthma is a very heterozygous disease, divided in subtypes, such as eosinophilic and neutrophilic asthma. Phenotyping and endotyping of patients, especially patients with severe asthma who are refractory to standard treatment, are crucial in asthma management and are based on a combination of clinical and biological features. Nevertheless, the quest remains to find better biomarkers that distinguish asthma subtypes in a more clear and objective manner and to find new therapeutic targets to treat people with therapy-resistant asthma. In the past, research to identify asthma subtypes mainly focused on expression profiles of protein-coding genes. However, advances in RNA-sequencing technologies and the discovery of non-coding RNAs as important post-transcriptional regulators have provided an entire new field of research opportunities in asthma. This review focusses on long non-coding RNAs (lncRNAs) in asthma; these are non-coding RNAs with a length of more than 200 nucleotides. Many lncRNAs are differentially expressed in asthma, and several have been associated with asthma severity or inflammatory phenotype. Moreover, in vivo and in vitro functional studies have identified the mechanisms of action of specific lncRNAs. Although lncRNAs remain not widely studied in asthma, the current studies show the potential of lncRNAs as biomarkers and therapeutic targets as well as the need for further research.

IL13Rα1 prevents a castration resistant phenotype of prostate cancer by targeting hexokinase 2 for ubiquitin-mediated degradation

OBJECTIVE

Androgen deprivation therapy (ADT) is still the principal treatment option for prostate cancer (PCa). In addition to reactivation of androgen receptor signaling, the resistance of PCa to apoptosis during ADT also contributes to castration resistant PCa (CRPC). A previous study reported that gene transfer of IL-13Rα2 into PCa cells sensitized the cells to the IL-13R-targeted cytotoxin IL13Rα1, leading to apoptosis. Compared with IL-13Rα2, IL13Rα1 is more constitutively expressed in PCa cells, but its function in PCa remains to be established.

METHODS

We determined the role and expression of IL13Rα1 in PCa cancer cells using western blotting, flow cytometry, and cell proliferation assays. Co-immunoprecipitation and mass spectrometry were used to identify the proteins that interacted with IL13Rα1, to elucidate its function.

RESULTS

In this study, we showed that IL13Rα1 was selectively suppressed in androgen-deprived PCa cells and that its suppression tended to be associated with poor prognoses of PCa patients. IL13Rα1 overexpression promoted apoptosis and inhibited tumor growth under androgen-deprived or castrated conditions (P < 0.01). Mechanistically, IL13Rα1 recruited and facilitated ubiquitin protein ligase E3C-mediated ubiquitination and degradation of hexokinase 2 (HK2), resulting in glycolytic inhibition and eventually leading to PCa cell apoptosis. Furthermore, our data revealed that mutated ataxia-telangiectasia kinase phosphorylated and facilitated the selective ubiquitin proteasome-mediated degradation of HK2. Notably, IL13Rα1-overexpressing PCa cells were more susceptible to apoptosis and exhibited reduced tumor growth after exposure to the HK2 inhibitor, 2-deoxy-D-glucose (P < 0.01).

CONCLUSIONS

Our data identified a tumor suppressor role for IL13Rα1 in preventing the resistance of PCa cells to apoptosis during androgen deprivation by inhibiting glycolysis. IL13Rα1-mediated signaling involving HK2 may therefore provide a novel treatment target and strategy for CRPC.

Identification of ROBO2 as a Potential Locus Associated with Inhaled Corticosteroid Response in Childhood Asthma

Inhaled corticosteroids (ICS) are the most common asthma controller medication. An important contribution of genetic factors in ICS response has been evidenced. Here, we aimed to identify novel genetic markers involved in ICS response in asthma. A genome-wide association study (GWAS) of the change in lung function after 6 weeks of ICS treatment was performed in 166 asthma patients from the SLOVENIA study. Patients with an improvement in lung function ≥8% were considered as ICS responders. Suggestively associated variants (p-value ≤ 5 × 10⁻⁶) were evaluated in an independent study (n = 175). Validation of the association with asthma exacerbations despite ICS use was attempted in European (n = 2681) and admixed (n = 1347) populations. Variants previously associated with ICS response were also assessed for replication. As a result, the SNP rs1166980 from the ROBO2 gene was suggestively associated with the change in lung function (OR for G allele: 7.01, 95% CI: 3.29–14.93, p = 4.61 × 10⁻⁷), although this was not validated in CAMP. ROBO2 showed gene-level evidence of replication with asthma exacerbations despite ICS use in Europeans (minimum p-value = 1.44 × 10⁻⁵), but not in admixed individuals. The association of PDE10A-T with ICS response described by a previous study was validated. This study suggests that ROBO2 could be a potential novel locus for ICS response in Europeans.

Identification of OCA2 as a novel locus for the co-morbidity of asthma-plus-eczema

BACKGROUND

Numerous genes have been associated with the three most common allergic diseases (asthma, allergic rhinitis or eczema) but these genes explain only a part of the heritability. In the vast majority of genetic studies, complex phenotypes such as co-morbidity of two of these diseases, have not been considered. This may partly explain missing heritability.

OBJECTIVE

To identify genetic variants specifically associated with the co-morbidity of asthma-plus-eczema.

METHODS

We first conducted a meta-analysis of four GWAS (Genome-Wide Association Study) of the combined asthma-plus-eczema phenotype (total of 8807 European-ancestry subjects of whom 1208 subjects had both asthma and eczema). To assess whether the association with SNP(s) was specific to the co-morbidity, we also conducted a meta-analysis of homogeneity test of association according to disease status ("asthma-plus-eczema" vs. the presence of only one disease "asthma only or eczema only"). We then used a joint test by combining the two test statistics from the co-morbidity-SNP association and the phenotypic heterogeneity of SNP effect meta-analyses.

RESULTS

Seven SNPs were detected for specific association to the asthma-plus-eczema co-morbidity, two with significant and five with suggestive evidence using the joint test after correction for multiple testing. The two significant SNPs are located in the OCA2 gene (Oculocutaneous Albinism II), a new locus never detected for significant evidence of association with any allergic disease. This gene is a promising candidate gene, because of its link to skin and lung diseases, and to epithelial barrier and immune mechanisms.

CONCLUSION

Our study underlines the importance of studying sub-phenotypes as co-morbidities to detect new susceptibility genes.

Interactive Roles of CaMKII/Ryanodine Receptor Signaling and Inflammation in Lung Diseases

Ca²⁺/calmodulin-dependent protein kinase II (CaMKII) is a multifunctional protein kinase and has been recently recognized to play a vital role in pathological events in the pulmonary system. CaMKII has diverse downstream targets that promote vascular disease, asthma, and cancer, so improved understanding of CaMKII signaling has the potential to lead to new therapies for lung diseases. Multiple studies have demonstrated that CaMKII is involved in redox modulation of ryanodine receptors (RyRs). CaMKII can be directly activated by reactive oxygen species (ROS) which then regulates RyR activity, which is essential for Ca²⁺-dependent processes in lung diseases. Furthermore, both CaMKII and RyRs participate in the inflammation process. However, their role in the pulmonary physiology in response to ROS is still an ambiguous one. Because CaMKII and RyRs are important in pulmonary biology, cell survival, cell cycle control, and inflammation, it is possible that the relationship between ROS and CaMKII/RyRs signal complex will be necessary for understanding and treating lung diseases. Here, we review roles of CaMKII/RyRs in lung diseases to understand with how CaMKII/RyRs may act as a transduction signal to connect prooxidant conditions into specific downstream pathological effects that are relevant to rare and common forms of pulmonary disease.

The Co-existence of ADHD With Autism in Saudi Children: An Analysis Using Next-Generation DNA Sequencing

Attention-deficit/hyperactivity disorder (ADHD) is one of the most common neurodevelopmental disorders. Several studies have confirmed the co-existence of other neuropsychiatric disorders with ADHD. Out of 106 individuals suspected to have ADHD, eight Saudi Arabian pediatric patients were diagnosed with ADHD using a dual assessment procedure based on highly significant scores from the international criteria for diagnosis; (full form DMS) DSM-5. Then, these patients were examined for the co-existence of autism and ADHD using different international diagnostic protocols. Four patients with combined ADHD and autism and four ADHD patients without autism were examined for the presence of genetic variants. Six variants (chr1:98165091, chr6:32029183, chr6:32035603, chr6:32064098, chr8:2909992, chr16:84213434) were identified in 75% of the patients with ADHD and autism, indicating that these genes may have a possible role in causing autism. Five variants (The chr2:116525960, chr15:68624396, chr15:91452595, chr15:92647645, and chr16:82673047) may increase to the severity of ADHD. This study recommends screening these eleven variants in ADHD cases and their relevant controls to confirm the prevalence in the Saudi population. It is recommended that future studies examine the 11 variants in detail.

An Eight-CircRNA Assessment Model for Predicting Biochemical Recurrence in Prostate Cancer

Prostate cancer (PCa) is a high morbidity malignancy in males, and biochemical recurrence (BCR) may appear after the surgery. Our study is designed to build up a risk score model using circular RNA sequencing data for PCa. The dataset is from the GEO database, using a cohort of 144 patients in Canada. We removed the low abundance circRNAs (FPKM < 1) and obtained 546 circRNAs for the next step. BCR-related circRNAs were selected by Logistic regression using the “survival” and “survminer” R package. Least absolute shrinkage and selector operation (LASSO) regression with 10-fold cross-validation and penalty was used to construct a risk score model by “glmnet” R software package. In total, eight circRNAs (including circ_30029, circ_117300, circ_176436, circ_112897, circ_112897, circ_178252, circ_115617, circ_14736, and circ_17720) were involved in our risk score model. Further, we employed differentially expressed mRNAs between high and low risk score groups. The following Gene Ontology (GO) analysis were visualized by Omicshare Online tools. As per the GO analysis results, tumor immune microenvironment related pathways are significantly enriched. “CIBERSORT” and “ESTIMATE” R package were used to detect tumor-infiltrating immune cells and compare the level of microenvironment scores between high and low risk score groups. What’s more, we verified two of eight circRNA’s (circ_14736 and circ_17720) circular characteristics and tested their biological function with qPCR and CCK8 in vitro. circ_14736 and circ_17720 were detected in exosomes of PCa patients’ plasma. This is the first bioinformatics study to establish a prognosis model for prostate cancer using circRNA. These circRNAs were associated with CD8⁺ T cell activities and may serve as a circRNA-based liquid biopsy panel for disease prognosis.

IGSF3 mutation identified in patient with severe COPD alters cell function and motility

Cigarette smoking (CS) and genetic susceptibility determine the risk for development, progression, and severity of chronic obstructive pulmonary diseases (COPD). We posited that an incidental balanced reciprocal chromosomal translocation was linked to a patient's risk of severe COPD. We determined that 46,XX,t(1;4)(p13.1;q34.3) caused a breakpoint in the immunoglobulin superfamily member 3 (IGSF3) gene, with markedly decreased expression. Examination of COPDGene cohort identified 14 IGSF3 SNPs, of which rs1414272 and rs12066192 were directly and rs6703791 inversely associated with COPD severity, including COPD exacerbations. We confirmed that IGSF3 is a tetraspanin-interacting protein that colocalized with CD9 and integrin B1 in tetraspanin-enriched domains. IGSF3-deficient patient-derived lymphoblastoids exhibited multiple alterations in gene expression, especially in the unfolded protein response and ceramide pathways. IGSF3-deficient lymphoblastoids had high ceramide and sphingosine-1 phosphate but low glycosphingolipids and ganglioside levels, and they were less apoptotic and more adherent, with marked changes in multiple TNFRSF molecules. Similarly, IGSF3 knockdown increased ceramide in lung structural cells, rendering them more adherent, with impaired wound repair and weakened barrier function. These findings suggest that, by maintaining sphingolipid and membrane receptor homeostasis, IGSF3 is required for cell mobility-mediated lung injury repair. IGSF3 deficiency may increase susceptibility to CS-induced lung injury in COPD.

A Bioinformatics Crash Course for Interpreting Genomics Data

Reductions in genotyping costs and improvements in computational power have made conducting genome-wide association studies (GWAS) standard practice for many complex diseases. GWAS is the assessment of genetic variants across the genome of many individuals to determine which, if any, genetic variants are associated with a specific trait. As with any analysis, there are evolving best practices that should be followed to ensure scientific rigor and reliability in the conclusions. This article presents a brief summary for many of the key bioinformatics considerations when either planning or evaluating GWAS. This review is meant to serve as a guide to those without deep expertise in bioinformatics and GWAS and give them tools to critically evaluate this popular approach to investigating complex diseases. In addition, a checklist is provided that can be used by investigators to evaluate whether a GWAS has appropriately accounted for the many potential sources of bias and generally followed current best practices.

Genetic profiling for disease stratification in chronic obstructive pulmonary disease and asthma

PURPOSE OF REVIEW

In asthma and chronic obstructive pulmonary disease (COPD), the movement towards genetic profiling with a push towards 'personalized medicine' has been hindered by complex environment--gene interactions and lack of tools to identify clear causal genetic traits. In this review, we will discuss the need for genetic profiling in asthma and COPD, what methods are currently used in the clinics and the recent finding using new sequencing methods.

RECENT FINDINGS

Over the past 10-15 years, genome-wide association studies analysis of common variants has provide little in the way of new genetic profiling markers for asthma and COPD. Whole exome/genome sequencing has provided a new method to identify lowly abundant alleles, which might have a much higher impact. Although, low population numbers due to high costs has hindered early studies, recent studies have reached genome wide significance.

SUMMARY

The use of genetic profiling of COPD in the clinic is current limited to the identification of Alpha-1 antitrypsin deficiency, while being absent in asthma. Advances in sequencing technology provide new avenues to identify disease causes or therapy response altering variants that in the short-term will allow for the development of screening procedures for disease to identify patients at risk of developing asthma or COPD.

Admixture mapping of asthma in southwestern Europeans with North African ancestry influences

The prevalence of asthma symptoms in Canary Islanders, a southwestern European population from Spain, is almost three times higher than the country average. Because the genetic risks identified so far explain <5% of asthma heritability, here we aimed to discover new asthma loci by completing the first admixture mapping study in Canary Islanders leveraging their distinctive genetic makeup, where significant northwest African influences coexist in the European genetic diversity landscape. A 2-stage study was conducted in 1,491 unrelated individuals self-declaring having a Canary Islands origin for the 4 grandparents. Local ancestry estimates were obtained for the shared positions with reference data from putative ancestral populations from Europe, North Africa, and sub-Saharan Africa. Case-control deviations in local ancestry were tested for each ancestry separately using logistic regressions adjusted for principal components, followed by fine-mapping analyses based on imputed genotypes and analyses of the likely deleterious exonic variants. The admixture mapping analysis of asthma detected that local North African ancestry in a locus spanning 365 kb of chromosome 16q23.3 was associated with asthma risk at study-wide significance [lowest P = 1.12 × 10^-4; odds ratio (OR) = 2.05; 95% confidence interval (CI) = 1.41-3.00]. Fine-mapping studies identified a variant associated with asthma, and results were replicated in independent samples (rs3852738, OR = 1.34; 95% CI = 1.13-1.59, P = 7.58 × 10^-4). Whole exome sequencing data from a subset of individuals revealed an enrichment of likely deleterious variants among asthma cases in 16q23.3, particularly in the phospholipase Cγ2 (PLCG2) gene (P = 3.67 × 10^-4). By completing the first mapping study of asthma in admixed populations from Europe, our results revealed a new plausible asthma locus.

A GWAS approach identifies Dapp1 as a determinant of air pollution-induced airway hyperreactivity

Asthma is a chronic inflammatory disease of the airways with contributions from genes, environmental exposures, and their interactions. While genome-wide association studies (GWAS) in humans have identified ~200 susceptibility loci, the genetic factors that modulate risk of asthma through gene-environment (GxE) interactions remain poorly understood. Using the Hybrid Mouse Diversity Panel (HMDP), we sought to identify the genetic determinants of airway hyperreactivity (AHR) in response to diesel exhaust particles (DEP), a model traffic-related air pollutant. As measured by invasive plethysmography, AHR under control and DEP-exposed conditions varied 3-4-fold in over 100 inbred strains from the HMDP. A GWAS with linear mixed models mapped two loci significantly associated with lung resistance under control exposure to chromosomes 2 (p = 3.0x10^-6) and 19 (p = 5.6x10^-7). The chromosome 19 locus harbors Il33 and is syntenic to asthma association signals observed at the IL33 locus in humans. A GxE GWAS for post-DEP exposure lung resistance identified a significantly associated locus on chromosome 3 (p = 2.5x10^-6). Among the genes at this locus is Dapp1, an adaptor molecule expressed in immune-related and mucosal tissues, including the lung. Dapp1-deficient mice exhibited significantly lower AHR than control mice but only after DEP exposure, thus functionally validating Dapp1 as one of the genes underlying the GxE association at this locus. In summary, our results indicate that some of the genetic determinants for asthma-related phenotypes may be shared between mice and humans, as well as the existence of GxE interactions in mice that modulate lung function in response to air pollution exposures relevant to humans.

The genetic factors that modulate risk of asthma through gene-environment (GxE) interactions are poorly understood, due in large part to the inherent difficulties in carrying out such studies in humans. To address these challenges, we used the Hybrid Mouse Diversity Panel to elucidate the genetic architecture of asthma-related phenotypes in mice and identify loci that are associated with airway hyperreactivity (AHR) under control exposure conditions and in response to diesel exhaust particles (DEP), as a model traffic-related air pollutant. In the absence of exposure, we identified two loci on chromosomes 2 and 19 for AHR. The locus on chromosome 19 harbors Il33 and is syntenic to association signals observed for asthma at the IL33 locus in humans. In response to DEP exposure, we mapped AHR to a region on chromosome 3 and used a genetically modified mouse model to functionally demonstrate that Dapp1 is one of the genes underlying the GxE association at this locus. Collectively, our results support the concept that some of the genetic determinants for asthma-related phenotypes may be shared between mice and humans as well as the existence of GxE interactions in mice that modulate lung function in response to air pollution exposures relevant to humans.

BAL Cell Gene Expression in Severe Asthma Reveals Mechanisms of Severe Disease and Influences of Medications

Rationale: Gene expression of BAL cells, which samples the cellular milieu within the lower respiratory tract, has not been well studied in severe asthma.Objectives: To identify new biomolecular mechanisms underlying severe asthma by an unbiased, detailed interrogation of global gene expression.Methods: BAL cell expression was profiled in 154 asthma and control subjects. Of these participants, 100 had accompanying airway epithelial cell gene expression. BAL cell expression profiles were related to participant (age, sex, race, and medication) and sample traits (cell proportions), and then severity-related gene expression determined by correlating transcripts and coexpression networks to lung function, emergency department visits or hospitalizations in the last year, medication use, and quality-of-life scores.Measurements and Main Results: Age, sex, race, cell proportions, and medications strongly influenced BAL cell gene expression, but leading severity-related genes could be determined by carefully identifying and accounting for these influences. A BAL cell expression network enriched for cAMP signaling components most differentiated subjects with severe asthma from other subjects. Subsequently, an in vitro cellular model showed this phenomenon was likely caused by a robust upregulation in cAMP-related expression in nonsevere and β-agonist-naive subjects given a β-agonist before cell collection. Interestingly, ELISAs performed on BAL lysates showed protein levels may partly disagree with expression changes.Conclusions: Gene expression in BAL cells is influenced by factors seldomly considered. Notably, β-agonist exposure likely had a strong and immediate impact on cellular gene expression, which may not translate to important disease mechanisms or necessarily match protein levels. Leading severity-related genes were discovered in an unbiased, system-wide analysis, revealing new targets that map to asthma susceptibility loci.

Molecular Subtyping and Prognostic Assessment Based on Tumor Mutation Burden in Patients with Lung Adenocarcinomas

The distinct molecular subtypes of lung cancer are defined by monogenic biomarkers, such as EGFR, KRAS, and ALK rearrangement. Tumor mutation burden (TMB) is a potential biomarker for response to immunotherapy, which is one of the measures for genomic instability. The molecular subtyping based on TMB has not been well characterized in lung adenocarcinomas in the Chinese population. Here we performed molecular subtyping based on TMB with the published whole exome sequencing data of 101 lung adenocarcinomas and compared the different features of the classified subtypes, including clinical features, somatic driver genes, and mutational signatures. We found that patients with lower TMB have a longer disease-free survival, and higher TMB is associated with smoking and aging. Analysis of somatic driver genes and mutational signatures demonstrates a significant association between somatic RYR2 mutations and the subtype with higher TMB. Molecular subtyping based on TMB is a potential prognostic marker for lung adenocarcinoma. Signature 4 and the mutation of RYR2 are highlighted in the TMB-High group. The mutation of RYR2 is a significant biomarker associated with high TMB in lung adenocarcinoma.

A genome-wide association study of tramadol metabolism from post-mortem samples

Phase I tramadol metabolism requires cytochrome p450 family 2, subfamily D, polypeptide 6 (CYP2D6) to form O-desmethyltramadol (M1). CYP2D6 genetic variants may infer metabolizer phenotype; however, drug ADME (absorption, distribution, metabolism, and excretion) and response depend on protein pathway(s), not CYP2D6 alone. There is a paucity of data regarding the contribution of trans-acting proteins to idiosyncratic phenotypes following drug exposure. A genome-wide association study identified five markers (rs79983226/kgp11274252, rs9384825, rs62435418/kgp10370907, rs72732317/kgp3743668, and rs184199168/exm1592932) associated with the conversion of tramadol to M1 (M1:T). These SNPs reside within five genes previously implicated with adverse reactions. Analysis of accompanying toxicological meta-data revealed a significant positive linear relationship between M1:T and degree of sample polypharmacy. Taken together, these data identify candidate loci for potential clinical inferences of phenotype following exposure to tramadol and highlight sample polypharmacy as a possible diagnostic covariate in post-mortem genetic studies.

African ancestry is associated with cluster-based childhood asthma subphenotypes

BACKGROUND

Childhood asthma is a syndrome composed of heterogeneous phenotypes; furthermore, intrinsic biologic variation among racial/ethnic populations suggests possible genetic ancestry variation in childhood asthma. The objective of the study is to identify clinically homogeneous asthma subphenotypes in a diverse sample of asthmatic children and to assess subphenotype-specific genetic ancestry in African-American asthmatic children.

METHODS

A total of 1211 asthmatic children including 813 in the Childhood Asthma Management Program and 398 in the Childhood Asthma Research and Education program were studied. Unsupervised cluster analysis on clinical phenotypes was conducted to identify homogeneous subphenotypes. Subphenotype-specific genetic ancestry was estimated for 167 African-American asthmatic children. Genetic ancestry association with subphenotypes/clinical phenotypes were determined.

RESULTS

Three distinct subphenotypes were identified: a moderate atopic dermatitis (AD) group with negative skin prick test (SPT) and preserved lung function; a high AD group with positive SPT and airway hyperresponsiveness; and a low AD group with positive SPT and lower lung function. African ancestry at asthma genome-wide association study (GWAS) SNPs differed between subphenotypes (64, 89, and 94% for the three subphenotypes, respectively) and was inversely correlated with AD; each additional 10% increase in African ancestry was associated with 1.5 fold higher in IgE and 6.3 higher odds of positive SPT (all p-values < 0.0001).

CONCLUSIONS

By conducting phenotype-based cluster analysis and assessing subphenotype-specific genetic ancestry, we were able to identify homogeneous subphenotypes for childhood asthma that showed significant variation in genetic ancestry of African-American asthmatic children. This finding demonstrates the utility of these complementary approaches to understand and refine childhood asthma subphenotypes and enable more targeted therapy.

MKLN1 splicing defect in dogs with lethal acrodermatitis

Lethal acrodermatitis (LAD) is a genodermatosis with monogenic autosomal recessive inheritance in Bull Terriers and Miniature Bull Terriers. The LAD phenotype is characterized by poor growth, immune deficiency, and skin lesions, especially at the paws. Utilizing a combination of genome wide association study and haplotype analysis, we mapped the LAD locus to a critical interval of ~1.11 Mb on chromosome 14. Whole genome sequencing of an LAD affected dog revealed a splice region variant in the MKLN1 gene that was not present in 191 control genomes (chr14:5,731,405T>G or MKLN1:c.400+3A>C). This variant showed perfect association in a larger combined Bull Terrier/Miniature Bull Terrier cohort of 46 cases and 294 controls. The variant was absent from 462 genetically diverse control dogs of 62 other dog breeds. RT-PCR analysis of skin RNA from an affected and a control dog demonstrated skipping of exon 4 in the MKLN1 transcripts of the LAD affected dog, which leads to a shift in the MKLN1 reading frame. MKLN1 encodes the widely expressed intracellular protein muskelin 1, for which diverse functions in cell adhesion, morphology, spreading, and intracellular transport processes are discussed. While the pathogenesis of LAD remains unclear, our data facilitate genetic testing of Bull Terriers and Miniature Bull Terriers to prevent the unintentional production of LAD affected dogs. This study may provide a starting point to further clarify the elusive physiological role of muskelin 1 in vivo.

Lethal acrodermatitis (LAD) is an autosomal recessive hereditary disease in dogs. It is characterized by poor growth, immune deficiency and characteristic skin lesions of the paws and of the face. We mapped the LAD locus to a ~1.11 Mb segment on canine chromosome 14. Whole genome sequence data of an LAD affected dog and 191 controls revealed a candidate causative variant in the MKLN1 gene, encoding muskelin 1. The identified variant, a single nucleotide substitution, MKLN1:c.400+3A>C, altered the 5’-splice site at the beginning of intron 4. We experimentally confirmed that this variant leads to complete skipping of exon 4 in the MKLN1 mRNA in skin. Various cellular functions have been postulated for muskelin 1 including roles in intracellular transport processes, cell morphology, cell spreading, and cell adhesion. Our data from dogs reveal a novel in vivo role for muskelin 1 that is related to the immune system and skin. MKLN1 thus represents a novel candidate gene for human patients with unsolved acrodermatitis and/or immune deficiency phenotypes. LAD affected dogs may serve as models to gain more insights into the function of muskelin 1.

Genome-wide association analysis identifies new candidate risk loci for familial intracranial aneurysm in the French-Canadian population

Intracranial Aneurysm (IA) is a common disease with a worldwide prevalence of 1–3%. In the French-Canadian (FC) population, where there is an important founder effect, the incidence of IA is higher and is frequently seen in families. In this study, we genotyped a cohort of 257 mostly familial FC IA patients and 1,992 FC controls using the Illumina NeuroX SNP-chip. The most strongly associated loci were tested in 34 Inuit IA families and in 32 FC IA patients and 106 FC controls that had been exome sequenced (WES). After imputation, one locus at 3p14.2 (FHIT, rs1554600, p = 4.66 × 10^–9) reached a genome-wide significant level of association and a subsequent validation in Nunavik Inuit cohort further confirmed the significance of the FHIT variant association (rs780365, FBAT-O, p = 0.002839). Additionally, among the other promising loci (p < 5 × 10⁻⁶), the one at 3q13.2 (rs78125721, p = 4.77 × 10⁻⁷), which encompasses CCDC80, also showed an increased mutation burden in the WES data (CCDC80, SKAT-O, p = 0.0005). In this study, we identified two new potential IA loci in the FC population: FHIT, which is significantly associated with hypertensive IA, and CCDC80, which has potential genetic and functional relevance to IA pathogenesis, providing evidence on the additional risk loci for familial IA. We also replicated the previous IA GWAS risk locus 18q11.2, and suggested a potential locus at 8p23.1 that warrants further study.

The association between ambient temperature and childhood asthma: a systematic review

The objectives of this study are to review available information on the association between ambient temperature and childhood asthma, and to elucidate the possible underlying mechanisms of this relationship. A systematic review was conducted based on the papers retrieved from four databases, including PubMed, ProQuest, ScienceDirect, and Scopus. Papers examining the association of absolute temperature or temperature variation with childhood asthma published from 1 January 2000 to 31 December 2016 were included. Thirteen papers have quantified the effect of absolute temperature on childhood asthma, and six papers have examined the effect of intra- or inter-day temperature variation on childhood asthma. All studies were conducted in urban areas. Aeroallergen sensitizations were only considered in the analyses of one study. Discrepancy existed in the significance of the relationship between absolute temperature and childhood asthma, and also in the shape of this relationship (i.e. linear or non-linear) and whether temperature effects were lagged. Increasing evidence is suggesting non-linear relationship between absolute temperature and childhood asthma. Future research should investigate the burden of childhood asthma specifically attributable to extreme temperatures and temperature variation using advanced statistical approach, particularly in rural areas, after properly considering aeroallergens and air pollution. Projecting future burden of childhood asthma under climate change scenarios is also warranted.

Multiancestry association study identifies new asthma risk loci that colocalize with immune-cell enhancer marks

We examined common variation in asthma risk by conducting a meta-analysis of worldwide asthma genome-wide association studies (23,948 cases, 118,538 controls) from ethnically-diverse populations. We identified five new asthma loci, uncovered two additional novel associations at two known asthma loci, established asthma associations at two loci implicated previously in comorbidity of asthma plus hay fever, and confirmed nine known loci. Investigation of pleiotropy showed large overlaps in genetic variants with autoimmune and inflammatory diseases. Enrichment of asthma risk loci in enhancer marks, especially in immune cells, suggests a major role of these loci in the regulation of immune-related mechanisms.

Differential Expression of Serum MicroRNAs Supports CD4⁺ T Cell Differentiation into Th2/Th17 Cells in Severe Equine Asthma

MicroRNAs (miRNAs) regulate post-transcriptional gene expression and may be exported from cells via exosomes or in partnership with RNA-binding proteins. MiRNAs in body fluids can act in a hormone-like manner and play important roles in disease initiation and progression. Hence, miRNAs are promising candidates as biomarkers. To identify serum miRNA biomarkers in the equine model of asthma we investigated small RNA derived from the serum of 34 control and 37 asthmatic horses. These samples were used for next generation sequencing, novel miRNA identification and differential miRNA expression analysis. We identified 11 significantly differentially expressed miRNAs between case and control horses: eca-miR-128, eca-miR-744, eca-miR-197, eca-miR-103, eca-miR-107a, eca-miR-30d, eca-miR-140-3p, eca-miR-7, eca-miR-361-3p, eca-miR-148b-3p and eca-miR-215. Pathway enrichment using experimentally validated target genes of the human homologous miRNAs showed a significant enrichment in the regulation of epithelial-to-mesenchymal transition (key player in airway remodeling in asthma) and the phosphatidylinositol (3,4,5)-triphosphate (PIP3) signaling pathway (modulator of CD4⁺ T cell maturation and function). Downregulated miR-128 and miR-744 supports a Th2/Th17 type immune response in severe equine asthma.

Multiple Transcriptome Data Analysis Reveals Biologically Relevant Atopic Dermatitis Signature Genes and Pathways

Several studies have identified genes that are differentially expressed in atopic dermatitis (AD) compared to normal skin. However, there is also considerable variation in the list of differentially expressed genes (DEGs) reported by different groups and the exact cause of AD is still not fully understood. Using a rank-based approach, we analyzed gene expression data from five different microarray studies, comprising a total of 127 samples and more than 250,000 transcripts. A total of 89 AD gene expression signatures '89ADGES', including FLG gene, were identified to show dysregulation consistently across these studies. Using a Support Vector Machine, we showed that the '89ADGES' discriminates AD from normal skin with 98% predictive accuracy. Functional annotation of these genes implicated their roles in immune responses (e.g., betadefensin, microseminoprotein), keratinocyte differentiation/epidermal development (e.g., FLG, CORIN, AQP, LOR, KRT16), inflammation (e.g., IL37, IL27RA, CCL18) and lipid metabolism (e.g., AKR1B10, FAD7, FAR2). Subsequently, we validated a subset of signature genes using quantitative PCR in a mouse model. Using a bioinformatic approach, we identified keratinocyte pathway over-represented (P = <0.0006) among the 89 signature genes. Keratinocytes are known to play a major role in barrier function due to their location in the epidermis. Our result suggests that besides immune- mediated pathway, skin barrier pathways such as the keratinocyte differentiation pathway play a key role in AD pathogenesis. A better understanding of the role of keratinocytes in AD will be important for developing novel "barrier therapy" for this disease.

Genome-Wide Association Study of Staphylococcus aureus Carriage in a Community-Based Sample of Mexican-Americans in Starr County, Texas

Staphylococcus aureus is the number one cause of hospital-acquired infections. Understanding host pathogen interactions is paramount to the development of more effective treatment and prevention strategies. Therefore, whole exome sequence and chip-based genotype data were used to conduct rare variant and genome-wide association analyses in a Mexican-American cohort from Starr County, Texas to identify genes and variants associated with S. aureus nasal carriage. Unlike most studies of S. aureus that are based on hospitalized populations, this study used a representative community sample. Two nasal swabs were collected from participants (n = 858) 11–17 days apart between October 2009 and December 2013, screened for the presence of S. aureus, and then classified as either persistent, intermittent, or non-carriers. The chip-based and exome sequence-based single variant association analyses identified 1 genome-wide significant region (KAT2B) for intermittent and 11 regions suggestively associated with persistent or intermittent S. aureus carriage. We also report top findings from gene-based burden analyses of rare functional variation. Notably, we observed marked differences between signals associated with persistent and intermittent carriage. In single variant analyses of persistent carriage, 7 of 9 genes in suggestively associated regions and all 5 top gene-based findings are associated with cell growth or tight junction integrity or are structural constituents of the cytoskeleton, suggesting that variation in genes associated with persistent carriage impact cellular integrity and morphology.

Mapping asthma-associated variants in admixed populations

Admixed populations arise when two or more previously isolated populations interbreed. Mapping asthma susceptibility loci in an admixed population using admixture mapping (AM) involves screening the genome of individuals of mixed ancestry for chromosomal regions that have a higher frequency of alleles from a parental population with higher asthma risk as compared with parental population with lower asthma risk. AM takes advantage of the admixture created in populations of mixed ancestry to identify genomic regions where an association exists between genetic ancestry and asthma (in contrast to between the genotype of the marker and asthma). The theory behind AM is that chromosomal segments of affected individuals contain a significantly higher-than-average proportion of alleles from the high-risk parental population and thus are more likely to harbor disease-associated loci. Criteria to evaluate the applicability of AM as a gene mapping approach include: (1) the prevalence of the disease differences in ancestral populations from which the admixed population was formed; (2) a measurable difference in disease-causing alleles between the parental populations; (3) reduced linkage disequilibrium (LD) between unlinked loci across chromosomes and strong LD between neighboring loci; (4) a set of markers with noticeable allele-frequency differences between parental populations that contributes to the admixed population (single nucleotide polymorphisms (SNPs) are the markers of choice because they are abundant, stable, relatively cheap to genotype, and informative with regard to the LD structure of chromosomal segments); and (5) there is an understanding of the extent of segmental chromosomal admixtures and their interactions with environmental factors. Although genome-wide association studies have contributed greatly to our understanding of the genetic components of asthma, the large and increasing degree of admixture in populations across the world create many challenges for further efforts to map disease-causing genes. This review, summarizes the historical context of admixed populations and AM, and considers current opportunities to use AM to map asthma genes. In addition, we provide an overview of the potential limitations and future directions of AM in biomedical research, including joint admixture and association mapping for asthma and asthma-related disorders.

Genes associated with RSV lower respiratory tract infection and asthma: the application of genetic epidemiological methods to understand causality

Infants with respiratory syncytial virus (RSV) lower respiratory tract infections (LRIs) are at increased risk for childhood asthma. The objectives of this article are to review the genes associated with both RSV LRI and asthma, review analytic approaches to assessing shared genetic risk and propose a future perspective on how these approaches can help us to understand the role of infant RSV infection as both an important risk factor for asthma and marker of shared genetic etiology between the two conditions. The review of shared genes and thus pathways associated with severity of response to RSV infection and asthma risk can help us to understand mechanisms of disease and ultimately propose new and novel targets for primary prevention of both diseases.

Genome-Wide Association Study Identifies Novel Loci Associated With Diisocyanate-Induced Occupational Asthma

Diisocyanates, reactive chemicals used to produce polyurethane products, are the most common causes of occupational asthma. The aim of this study is to identify susceptibility gene variants that could contribute to the pathogenesis of diisocyanate asthma (DA) using a Genome-Wide Association Study (GWAS) approach. Genome-wide single nucleotide polymorphism (SNP) genotyping was performed in 74 diisocyanate-exposed workers with DA and 824 healthy controls using Omni-2.5 and Omni-5 SNP microarrays. We identified 11 SNPs that exceeded genome-wide significance; the strongest association was for the rs12913832 SNP located on chromosome 15, which has been mapped to the HERC2 gene (p = 6.94 × 10(-14)). Strong associations were also found for SNPs near the ODZ3 and CDH17 genes on chromosomes 4 and 8 (rs908084, p = 8.59 × 10(-9) and rs2514805, p = 1.22 × 10(-8), respectively). We also prioritized 38 SNPs with suggestive genome-wide significance (p < 1 × 10(-6)). Among them, 17 SNPs map to the PITPNC1, ACMSD, ZBTB16, ODZ3, and CDH17 gene loci. Functional genomics data indicate that 2 of the suggestive SNPs (rs2446823 and rs2446824) are located within putative binding sites for the CCAAT/Enhancer Binding Protein (CEBP) and Hepatocyte Nuclear Factor 4, Alpha transcription factors (TFs), respectively. This study identified SNPs mapping to the HERC2, CDH17, and ODZ3 genes as potential susceptibility loci for DA. Pathway analysis indicated that these genes are associated with antigen processing and presentation, and other immune pathways. Overlap of 2 suggestive SNPs with likely TF binding sites suggests possible roles in disruption of gene regulation. These results provide new insights into the genetic architecture of DA and serve as a basis for future functional and mechanistic studies.

Contributing factors to the development of childhood asthma: working toward risk minimization

Asthma is the most common chronic disease in childhood, and considerable research has been undertaken to find ways to prevent its development and reduce its prevalence. For such interventions to be successful, risk factors for asthma emergence should be identified and clearly defined. Data are robust for some of them, including atopy, viral infections and exposure to airborne irritants, whereas it is less conclusive for others, such as aeroallergen exposure and bacterial infections. Several interventions for asthma prevention, including avoidance and pharmacotherapy, have been attempted. However, most of them have furnished equivocal results. Various issues hinder the establishment of risk factors for asthma development and reduce the effectiveness of interventions, including the complexity of the disease and the fluidity of the developing systems in childhood. In this review, we revisit the evidence on pediatric asthma risk factors and prevention and discuss issues that perplex this field.

Childhood asthma is associated with polymorphic markers of PROC on 2q14 in addition to 17q21 locus

BACKGROUND

Childhood asthma is caused by both genetic and environmental factors. The first genomewide association study (GWAS) for asthma revealed putative candidates on nine chromosomal regions in Caucasians, with 17q21 locus being the most widely replicated one. However, there was no replication study for the other loci. This study investigated genetic associations between childhood asthma and autosomal single nucleotide polymorphisms (SNPs) on eight loci reported in the first GWAS among Hong Kong Chinese.

METHODS

510 asthmatic children and 510 non-allergic controls were recruited. 110 tagging SNPs selected based on r(2 ) ≥ 0.80 and minor allele frequency ≥0.05 for Han Chinese among all SNPs located 50-kb upstream and downstream of significant autosomal SNPs were genotyped by TaqMan allelic discrimination assays. Transcription factor binding of SNPs was determined by electrophoretic mobility shift assay (EMSA).

RESULTS

Asthma was significantly associated with SNPs on 17q21 and 2q14 loci. Twelve SNPs on 17q21 were associated with asthma, with rs6503527 being the most significant SNP. Five SNPs of protein C gene (PROC) on 2q14 were associated with asthma, with rs6755028 being the most significant SNP. Plasma protein C concentrations were higher in asthmatic patients than controls, and five PROC SNPs were associated with plasma protein C concentrations. EMSA showed specific differential binding of rs878461 to nuclear extracts from bronchial epithelial and hepatocarcinoma cell lines.

CONCLUSIONS

Our findings identify PROC on 2q14 as a novel candidate for childhood asthma and replicate the genetic association for 17q21 locus. Rs878461 of PROC may increase asthma susceptibility by altering transcription factor binding.

Genome-wide association studies in Africans and African Americans: expanding the framework of the genomics of human traits and disease

Genomic research is one of the tools for elucidating the pathogenesis of diseases of global health relevance and paving the research dimension to clinical and public health translation. Recent advances in genomic research and technologies have increased our understanding of human diseases, genes associated with these disorders, and the relevant mechanisms. Genome-wide association studies (GWAS) have proliferated since the first studies were published several years ago and have become an important tool in helping researchers comprehend human variation and the role genetic variants play in disease. However, the need to expand the diversity of populations in GWAS has become increasingly apparent as new knowledge is gained about genetic variation. Inclusion of diverse populations in genomic studies is critical to a more complete understanding of human variation and elucidation of the underpinnings of complex diseases. In this review, we summarize the available data on GWAS in recent African ancestry populations within the western hemisphere (i.e. African Americans and peoples of the Caribbean) and continental African populations. Furthermore, we highlight ways in which genomic studies in populations of recent African ancestry have led to advances in the areas of malaria, HIV, prostate cancer, and other diseases. Finally, we discuss the advantages of conducting GWAS in recent African ancestry populations in the context of addressing existing and emerging global health conditions.

Genetics of allergic diseases

Genome-wide association studies (GWAS) have been employed in the field of allergic disease, and significant associations have been published for nearly 100 asthma genes/loci. An outcome of GWAS in allergic disease has been the formation of national and international collaborations leading to consortia meta-analyses, and an appreciation for the specificity of genetic associations to sub-phenotypes of allergic disease. Molecular genetics has undergone a technological revolution, leading to next-generation sequencing strategies that are increasingly employed to hone in on the causal variants associated with allergic diseases. Unmet needs include the inclusion of diverse cohorts and strategies for managing big data.

Immune mechanisms and development of childhood asthma

Early life influences are crucial for the development of distinct childhood asthma phenotypes, which are currently included under the term asthma syndrome. Improved characterisation of different childhood asthma phenotypes will help to elucidate specific underlying immune mechanisms--namely, endotypes. Besides genetics, epigenetics and environmental factors have an effect on innate and adaptive immune regulatory networks. Crucial determining factors for complex immune regulation and barrier function include family history of atopy, respiratory infections, microbiome, and nutrition. Recent diagnostic approaches, including biomarkers, might offer a unique opportunity to improve definitions of asthma sub-phenotypes, prediction of outcome, and treatment options, by referring to the underlying pathophysiology. For prevention and patient-individualised medicine, a multifactorial approach incorporating deep phenotyping and mathematical models for analysis to extend our present knowledge is needed.

Finding novel molecular connections between developmental processes and disease

Identifying molecular connections between developmental processes and disease can lead to new hypotheses about health risks at all stages of life. Here we introduce a new approach to identifying significant connections between gene sets and disease genes, and apply it to several gene sets related to human development. To overcome the limits of incomplete and imperfect information linking genes to disease, we pool genes within disease subtrees in the MeSH taxonomy, and we demonstrate that such pooling improves the power and accuracy of our approach. Significance is assessed through permutation. We created a web-based visualization tool to facilitate multi-scale exploration of this large collection of significant connections (http://gda.cs.tufts.edu/development). High-level analysis of the results reveals expected connections between tissue-specific developmental processes and diseases linked to those tissues, and widespread connections to developmental disorders and cancers. Yet interesting new hypotheses may be derived from examining the unexpected connections. We highlight and discuss the implications of three such connections, linking dementia with bone development, polycystic ovary syndrome with cardiovascular development, and retinopathy of prematurity with lung development. Our results provide additional evidence that plays a key role in the early pathogenesis of polycystic ovary syndrome. Our evidence also suggests that the VEGF pathway and downstream NFKB signaling may explain the complex relationship between bronchopulmonary dysplasia and retinopathy of prematurity, and may form a bridge between two currently-competing hypotheses about the molecular origins of bronchopulmonary dysplasia. Further data exploration and similar queries about other gene sets may generate a variety of new information about the molecular relationships between additional diseases.

Understanding the roles that genes involved in normal human development can play in disease processes is an important part of predicting disease risk and designing novel treatment approaches. In this study, we have identified classes of disease that are associated with a surprisingly large number of genes involved in any of several tissue-specific developmental processes. To do so, we developed a novel approach whose strength comes from pooling genetic information across related diseases, overcoming problems ordinarily posed by limited information about individual gene-disease relationships. We demonstrate the method's efficacy both by examining its ability to highlight connections between gene sets and disease classes that are known to be related, and by demonstrating that the approach recovers expected broad classes of connections, such as those between heart development and cardiovascular disorders. However, by examining unexpected connections in this data set, we are able to develop new understanding of some surprising disease relationships, such as the one between dementia and osteoporosis. Such connections may lead to a better overall understanding of the role of development in lifelong health, as well as to the design of new methods to treat a range of diseases.