Better together against genetic heterogeneity: A sex-combined joint main and interaction analysis of 290 quantitative traits in the UK Biobank.

Genetic effects can be sex-specific, particularly for traits such as testosterone, a sex hormone. While sex-stratified analysis provides easily interpretable sex-specific effect size estimates, the presence of sex-differences in SNP effect implies a SNP×sex interaction. This suggests the usage of the often overlooked joint test, testing for an SNP's main and SNP×sex interaction effects simultaneously. Notably, even without individual-level data, the joint test statistic can be derived from sex-stratified summary statistics through an omnibus meta-analysis. Utilizing the available sex-stratified summary statistics of the UK Biobank, we performed such omnibus meta-analyses for 290 quantitative traits. Results revealed that this approach is robust to genetic effect heterogeneity and can outperform the traditional sex-stratified or sex-combined main effect-only tests. Therefore, we advocate using the omnibus meta-analysis that captures both the main and interaction effects. Subsequent sex-stratified analysis should be conducted for sex-specific effect size estimation and interpretation.

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

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

Comprehensive whole-genome analyses of the UK Biobank reveal significant sex differences in both genotype missingness and allele frequency on the X chromosome.

The UK Biobank is the most used dataset for genome-wide association studies (GWAS). GWAS of sex, essentially sex differences in minor allele frequencies (sdMAF), has identified autosomal SNPs with significant sdMAF, including in the UK Biobank, but the X chromosome was excluded. Our recent report identified multiple regions on the X chromosome with significant sdMAF, using short-read sequencing of other datasets. We performed a whole genome sdMAF analysis, with ~410 k white British individuals from the UK Biobank, using array genotyped, imputed or exome sequencing data. We observed marked sdMAF on the X chromosome, particularly at the boundaries between the pseudo-autosomal regions (PAR) and the non-PAR (NPR), as well as throughout the NPR, consistent with our earlier report. A small fraction of autosomal SNPs also showed significant sdMAF. Using the centrally imputed data, which relied mostly on low-coverage whole genome sequence, resulted in 2.1% of NPR SNPs with significant sdMAF. The whole exome sequencing also displays sdMAF on the X chromosome, including some NPR SNPs with heterozygous genotype calls in males. Genotyping, sequencing and imputation of X chromosomal SNPs requires further attention to ensure the integrity for downstream association analysis.

eXclusionarY: 10 years later, where are the sex chromosomes in GWASs?

10 years ago, a detailed analysis showed that only 33% of genome-wide association study (GWAS) results included the X chromosome. Multiple recommendations were made to combat such exclusion. Here, we re-surveyed the research landscape to determine whether these earlier recommendations had been translated. Unfortunately, among the genome-wide summary statistics reported in 2021 in the NHGRI-EBI GWAS Catalog, only 25% provided results for the X chromosome and 3% for the Y chromosome, suggesting that the exclusion phenomenon not only persists but has also expanded into an exclusionary problem. Normalizing by physical length of the chromosome, the average number of studies published through November 2022 with genome-wide-significant findings on the X chromosome is ∼1 study/Mb. By contrast, it ranges from ∼6 to ∼16 studies/Mb for chromosomes 4 and 19, respectively. Compared with the autosomal growth rate of ∼0.086 studies/Mb/year over the last decade, studies of the X chromosome grew at less than one-seventh that rate, only ∼0.012 studies/Mb/year. Among the studies that reported significant associations on the X chromosome, we noted extreme heterogeneities in data analysis and reporting of results, suggesting the need for clear guidelines. Unsurprisingly, among the 430 scores sampled from the PolyGenic Score Catalog, 0% contained weights for sex chromosomal SNPs. To overcome the dearth of sex chromosome analyses, we provide five sets of recommendations and future directions. Finally, until the sex chromosomes are included in a whole-genome study, instead of GWASs, we propose such studies would more properly be referred to as "AWASs," meaning "autosome-wide scans."

Major sex differences in allele frequencies for X chromosomal variants in both the 1000 Genomes Project and gnomAD.

An unexpectedly high proportion of SNPs on the X chromosome in the 1000 Genomes Project phase 3 data were identified with significant sex differences in minor allele frequencies (sdMAF). sdMAF persisted for many of these SNPs in the recently released high coverage whole genome sequence of the 1000 Genomes Project that was aligned to GRCh38, and it was consistent between the five super-populations. Among the 245,825 common (MAF>5%) biallelic X-chromosomal SNPs in the phase 3 data presumed to be of high quality, 2,039 have genome-wide significant sdMAF (p-value <5e-8). sdMAF varied by location: non-pseudo-autosomal region (NPR) = 0.83%, pseudo-autosomal regions (PAR1) = 0.29%, PAR2 = 13.1%, and X-transposed region (XTR)/PAR3 = 0.85% of SNPs had sdMAF, and they were clustered at the NPR-PAR boundaries, among others. sdMAF at the NPR-PAR boundaries are biologically expected due to sex-linkage, but have generally been ignored in association studies. For comparison, similar analyses found only 6, 1 and 0 SNPs with significant sdMAF on chromosomes 1, 7 and 22, respectively. Similar sdMAF results for the X chromosome were obtained from the high coverage whole genome sequence data from gnomAD V 3.1.2 for both the non-Finnish European and African/African American samples. Future X chromosome analyses need to take sdMAF into account.

Gut Microbiome Composition Is Associated With Future Onset of Crohn's Disease in Healthy First-Degree Relatives.

BACKGROUND & AIMS: The cause of Crohn's disease (CD) is unknown, but the current hypothesis is that microbial or environmental factors induce gut inflammation in genetically susceptible individuals, leading to chronic intestinal inflammation. Case-control studies of patients with CD have cataloged alterations in the gut microbiome composition; however, these studies fail to distinguish whether the altered gut microbiome composition is associated with initiation of CD or is the result of inflammation or drug treatment. METHODS: In this prospective cohort study, 3483 healthy first-degree relatives (FDRs) of patients with CD were recruited to identify the gut microbiome composition that precedes the onset of CD and to what extent this composition predicts the risk of developing CD. We applied a machine learning approach to the analysis of the gut microbiome composition (based on 16S ribosomal RNA sequencing) to define a microbial signature that associates with future development of CD. The performance of the model was assessed in an independent validation cohort. RESULTS: In the validation cohort, the microbiome risk score (MRS) model yielded a hazard ratio of 2.24 (95% confidence interval, 1.03-4.84; P = .04), using the median of the MRS from the discovery cohort as the threshold. The MRS demonstrated a temporal validity by capturing individuals that developed CD up to 5 years before disease onset (area under the curve > 0.65). The 5 most important taxa contributing to the MRS included Ruminococcus torques, Blautia, Colidextribacter, an uncultured genus-level group from Oscillospiraceae, and Roseburia. CONCLUSION: This study is the first to demonstrate that gut microbiome composition is associated with future onset of CD and suggests that gut microbiome is a contributor in the pathogenesis of CD.

Analyses of potential causal contributors to increased waist/hip ratio-associated cardiometabolic disease: A combined and sex-stratified Mendelian randomization study.

BACKGROUND: Increased waist/hip ratio (WHR) contributes to type 2 diabetes, fatty liver, dyslipidaemia, hypertension and coronary artery disease, with potential sex-differential effects. Postulated mediators include increased lipid flux, branched-chain amino acids, glycine and glycoprotein acetyl, but their relative contributions and sex-specific impact on WHR-associated cardiometabolic disease (CMD) are not established. METHODS: We therefore undertook combined and sex-stratified Mendelian randomization (MR) to assess the relative causal contributions of these mediators to WHR-associated CMD using summary statistics from the largest genome-wide association studies in European ancestries. RESULTS: In sex-combined MR analyses, increased WHR significantly reduces high-density lipoprotein (beta = -0.416, SE = 0.029, p = 2.87E-47), increases triglyceride (beta = 0.431, SE = 0.029, p = 1.87E-50), type 2 diabetes (odds ratio = 2.747, SE = 0.09, p = 26E-23), coronary artery disease (odds ratio = 1.478, SE = 0.045, p = 6.96E-18), alanine transaminase (beta = 0.062, SE = 0.004, p = 6.88E-67), and systolic (beta = 0.134, SE = 0.022, p = 7.81E-10) and diastolic blood pressure (beta = 0.162, SE = 0.026, p = 5.38E-10). Adjustment for the mediators attenuated WHR's effects, but the associations remained significant with concordant results in females. In males, a similar pattern was seen, except after adjusting for the effect of the ratio of monounsaturated fatty acid to total free fatty acid, the potential causal effect of WHR was no longer significant: high-density lipoprotein (beta = -0.117, SE = 0.069, p = .09) and triglyceride (beta = 0.051, SE = 0.068, p = .459). CONCLUSIONS: MR suggests WHR increases the risk of CMD independent of these mediators, with the exception of dyslipidaemia in males, which is largely driven by the monounsaturated fatty acid to total free fatty acid ratio.

Mediterranean-Like Dietary Pattern Associations With Gut Microbiome Composition and Subclinical Gastrointestinal Inflammation.

BACKGROUND & AIMS: Case-control studies have shown that patients with Crohn's disease (CD) have a microbial composition different from healthy individuals. Although the causes of CD are unknown, epidemiologic studies suggest that diet is an important contributor to CD risk, potentially via modulation of bacterial composition and gut inflammation. We hypothesized that long-term dietary clusters (DCs) are associated with gut microbiome compositions and gut inflammation. Our objectives were to identify dietary patterns and assess whether they are associated with alterations in specific gut microbial compositions and subclinical levels of gut inflammation in a cohort of healthy first-degree relatives (FDRs) of patients with CD. METHODS: As part of the Genetic, Environmental, Microbial (GEM) Project, we recruited a cohort of 2289 healthy FDRs of patients with CD. Individuals provided stool samples and answered a validated food frequency questionnaire reflecting their habitual diet during the year before sample collection. Unsupervised analysis identified 3 dietary and 3 microbial composition clusters. RESULTS: DC3, resembling the Mediterranean diet, was strongly associated with a defined microbial composition, with an increased abundance of fiber-degrading bacteria, such as Ruminococcus, as well as taxa such as Faecalibacterium. The DC3 diet was also significantly associated with lower levels of subclinical gut inflammation, defined by fecal calprotectin, compared with other dietary patterns. No significant associations were found between individual food items and fecal calprotectin, suggesting that long-term dietary patterns rather than individual food items contribute to subclinical gut inflammation. Additionally, mediation analysis demonstrated that DC3 had a direct effect on subclinical inflammation that was partially mediated by the microbiota. CONCLUSIONS: Overall, these results indicated that Mediterranean-like dietary patterns are associated with microbiome and lower intestinal inflammation. This study will help guide future dietary strategies that affect microbial composition and host gut inflammation to prevent diseases.

Genetics of osteonecrosis in children and adults with systemic lupus erythematosus.

OBJECTIVES: Genetics plays an important role in SLE risk, as well as osteonecrosis (ON), a significant and often debilitating complication of SLE. We aimed to identify genetic risk loci for ON in people with childhood-onset (cSLE) and adult-onset (aSLE) SLE. METHODS: We enrolled participants from two tertiary care centres who met classification criteria for SLE. Participants had prospectively collected clinical data and were genotyped on a multiethnic array. Un-genotyped single nucleotide polymorphisms (SNPs) were imputed, and ancestry was inferred using principal components (PCs). Our outcome was symptomatic ON confirmed by imaging. We completed time-to-ON and logistic regression of ON genome-wide association studies (GWASs) with covariates for sex, age of SLE diagnosis, five PCs for ancestry, corticosteroid use and selected SLE manifestations. We conducted separate analyses for cSLE and aSLE and meta-analysed results using inverse-variance weighting. Genome-wide significance was P < 5 × 10-8. RESULTS: The study included 940 participants with SLE, 87% female and 56% with cSLE. ON was present in 7.6% (n = 71). Median age of SLE diagnosis was 16.9 years (interquartile range [IQR]: 13.5, 29.3), with median follow-up of 8.0 years (IQR: 4.2, 15.7). Meta-GWAS of cSLE and aSLE time-to-ON of 4 431 911 SNPs identified a significant Chr.2 SNP, rs34118383 (minor allele frequency = 0.18), intronic to WIPF1 (hazard ratio = 3.2 [95% CI: 2.2, 4.8]; P = 1.0 × 10-8). CONCLUSION: We identified an intronic WIPF1 variant associated with a 3.2 times increased hazard for ON (95% CI: 2.2, 4.8; P = 1.0 × 10-8) during SLE follow-up, independent of corticosteroid exposure. The effect of the SNP on time-to-ON was similar in cSLE and aSLE. This novel discovery represents a potential ON risk locus. Our results warrant replication.

Genetics of longitudinal kidney function in children and adults with systemic lupus erythematosus.

OBJECTIVES: Genome-wide association studies (GWAS) have identified loci associated with estimated glomerular filtration rate (eGFR). Few LN risk loci have been identified to date. We tested the association of SLE and eGFR polygenic risk scores (PRS) with repeated eGFR measures from children and adults with SLE. METHODS: Patients from two tertiary care lupus clinics that met ≥4 ACR and/or SLICC criteria for SLE were genotyped on the Illumina MEGA or Omni1-Quad arrays. PRSs were calculated for SLE and eGFR, using published weighted GWA-significant alleles. eGFR was calculated using the CKD-EPI and Schwartz equations. We tested the effect of eGFR- and SLE-PRSs on eGFR mean and variance, adjusting for age at diagnosis, sex, ancestry, follow-up time, and clinical event flags. RESULTS: We included 1158 SLE patients (37% biopsy-confirmed LN) with 36 733 eGFR measures over a median of 7.6 years (IQR: 3.9-15.3). LN was associated with lower within-person mean eGFR [LN: 93.8 (s.d. 26.4) vs non-LN: 101.6 (s.d. 17.7) mL/min per 1.73 m2; P < 0.0001] and higher variance [LN median: 157.0 (IQR: 89.5, 268.9) vs non-LN median: 84.9 (IQR: 46.9, 138.2) (mL/min per 1.73 m2)2; P < 0.0001]. Increasing SLE-PRSs were associated with lower mean eGFR and greater variance, while increasing eGFR-PRS was associated with increased eGFR mean and variance. CONCLUSION: We observed significant associations between SLE and eGFR PRSs and repeated eGFR measurements, in a large cohort of children and adults with SLE. Longitudinal eGFR may serve as a powerful alternative outcome to LN categories for discovery of LN risk loci.

Linkage analysis identifies an isolated strabismus locus at 14q12 overlapping with FOXG1 syndrome region.

Strabismus is a common condition, affecting 1%-4% of individuals. Isolated strabismus has been studied in families with Mendelian inheritance patterns. Despite the identification of multiple loci via linkage analyses, no specific genes have been identified from these studies. The current study is based on a seven-generation family with isolated strabismus inherited in an autosomal dominant manner. A total of 13 individuals from a common ancestor have been included for linkage analysis. Among these, nine are affected and four are unaffected. A single linkage signal has been identified at an 8.5 Mb region of chromosome 14q12 with a multipoint LOD (logarithm of the odds) score of 4.69. Disruption of this locus is known to cause FOXG1 syndrome (or congenital Rett syndrome; OMIM #613454 and *164874), in which 84% of affected individuals present with strabismus. With the incorporation of next-generation sequencing and in-depth bioinformatic analyses, a 4 bp non-coding deletion was prioritised as the top candidate for the observed strabismus phenotype. The deletion is predicted to disrupt regulation of FOXG1, which encodes a transcription factor of the Forkhead family. Suggestive of an autoregulation effect, the disrupted sequence matches the consensus FOXG1 and Forkhead family transcription factor binding site and has been observed in previous ChIP-seq studies to be bound by Foxg1 in early mouse brain development. Future study of this specific deletion may shed light on the regulation of FOXG1 expression and may enhance our understanding of the mechanisms contributing to strabismus and FOXG1 syndrome.

Genome-wide association study identifies novel loci associated with skin autofluorescence in individuals without diabetes.

BACKGROUND: Skin autofluorescence (SAF) is a non-invasive measure reflecting accumulation of advanced glycation endproducts (AGEs) in the skin. Higher SAF levels are associated with an increased risk of developing type 2 diabetes and cardiovascular disease. An earlier genome-wide association study (GWAS) revealed a strong association between NAT2 variants and SAF. The aim of this study was to calculate SAF heritability and to identify additional genetic variants associated with SAF through genome-wide association studies (GWAS). RESULTS: In 27,534 participants without diabetes the heritability estimate of lnSAF was 33% ± 2.0% (SE) in a model adjusted for covariates. In meta-GWAS for lnSAF five SNPs, on chromosomes 8, 11, 15 and 16 were associated with lnSAF (P < 5 × 10-8): 1. rs2846707 (Chr11:102,576,358,C > T), which results in a Met30Val missense variant in MMP27 exon 1 (NM_022122.3); 2. rs2470893 (Chr15:75,019,449,C > T), in intergenic region between CYP1A1 and CYP1A2; with attenuation of the SNP-effect when coffee consumption was included as a covariate; 3. rs12931267 (Chr16:89,818,732,C > G) in intron 30 of FANCA and near MC1R; and following conditional analysis 4. rs3764257 (Chr16:89,800,887,C > G) an intronic variant in ZNF276, 17.8 kb upstream from rs12931267; finally, 30 kb downstream from NAT2 5. rs576201050 (Chr8:18,288,053,G > A). CONCLUSIONS: This large meta-GWAS revealed five SNPs at four loci associated with SAF in the non-diabetes population. Further unravelling of the genetic architecture of SAF will help in improving its utility as a tool for screening and early detection of diseases and disease complications.

Fast and Accurate Shared Segment Detection and Relatedness Estimation in Un-phased Genetic Data via TRUFFLE.

Relationship estimation and segment detection between individuals is an important aspect of disease gene mapping. Existing methods are either tailored for computational efficiency or require phasing to improve accuracy. We developed TRUFFLE, a method that integrates computational techniques and statistical principles for the identification and visualization of identity-by-descent (IBD) segments using un-phased data. By skipping the haplotype phasing step and, instead, relying on a simpler region-based approach, our method is computationally efficient while maintaining inferential accuracy. In addition, an error model corrects for segment break-ups that occur as a consequence of genotyping errors. TRUFFLE can estimate relatedness for 3.1 million pairs from the 1000 Genomes Project data in a few minutes on a typical laptop computer. Consistent with expectation, we identified only three second cousin or closer pairs across different populations, while commonly used methods identified a large number of such pairs. Similarly, within populations, we identified many fewer related pairs. Compared to methods relying on phased data, TRUFFLE has comparable accuracy but is drastically faster and has fewer broken segments. We also identified specific local genomic regions that are commonly shared within populations, suggesting selection. When applied to pedigree data, we observed 99.6% accuracy in detecting 1st to 5th degree relationships. As genomic datasets become much larger, TRUFFLE can enable disease gene mapping through implicit shared haplotypes by accurate IBD segment detection.

Anti-Microbial Antibody Response is Associated With Future Onset of Crohn's Disease Independent of Biomarkers of Altered Gut Barrier Function, Subclinical Inflammation, and Genetic Risk.

BACKGROUND AND AIMS: Altered host immune reactivity to microbial antigens is hypothesized to trigger the onset of Crohn's disease (CD). We aimed to assess whether increased serum anti-microbial antibody response in asymptomatic first-degree relatives (FDRs) of CD patients is an independent risk factor for future CD development. METHODS: We measured host serum antibody response to 6 microbial antigens at enrollment (Prometheus enzyme-linked immunosorbent assay test: anti-Saccharomyces cerevisiae antibodies immunoglobulin A/immunoglobulin G, anti-OmpC, anti-A4-Fla2, anti-FlaX, anti-CBir1) and derived the sum of positive antibodies (AS). We used samples at enrollment of prospectively followed healthy FDRs from a nested case-control cohort of the Crohn's and Colitis Canada Genetics Environment Microbial Project. Those who later developed CD (n = 77) were matched 1:4 by age, sex, follow-up duration, and geographic location with control FDRs remaining healthy (n = 307). To address our research aims, we fitted a multivariable conditional logistic regression model and performed causal mediation analysis. RESULTS: High baseline AS (≥2) (43% of cases, 11% of controls) was associated with higher risk of developing CD (adjusted odds ratio, 6.5; 95% confidence interval, 3.4-12.7; P < .001). Importantly, this association remained significant when adjusted for markers of gut barrier function, fecal calprotectin, C-reactive protein, and CD-polygenic risk score, and in subjects recruited more than 3 years before diagnosis. Causal mediation analysis showed that the effect of high AS on future CD development is partially mediated (42%) via preclinical gut inflammation. CONCLUSIONS: Our results suggest that increased anti-microbial antibody responses are associated with risk of future development of CD, independent of biomarkers of abnormal gut barrier function, subclinical inflammation, and CD-related genetic risks. This suggests that anti-microbial antibody responses are an early predisease event in the development of CD.

Enhancer-gene rewiring in the pathogenesis of Quebec platelet disorder.

Quebec platelet disorder (QPD) is an autosomal dominant bleeding disorder with a unique, platelet-dependent, gain-of-function defect in fibrinolysis, without systemic fibrinolysis. The hallmark feature of QPD is a >100-fold overexpression of PLAU, specifically in megakaryocytes. This overexpression leads to a >100-fold increase in platelet stores of urokinase plasminogen activator (PLAU/uPA); subsequent plasmin-mediated degradation of diverse α-granule proteins; and platelet-dependent, accelerated fibrinolysis. The causative mutation is a 78-kb tandem duplication of PLAU. How this duplication causes megakaryocyte-specific PLAU overexpression is unknown. To investigate the mechanism that causes QPD, we used epigenomic profiling, comparative genomics, and chromatin conformation capture approaches to study PLAU regulation in cultured megakaryocytes from participants with QPD and unaffected controls. QPD duplication led to ectopic interactions between PLAU and a conserved megakaryocyte enhancer found within the same topologically associating domain (TAD). Our results support a unique disease mechanism whereby the reorganization of sub-TAD genome architecture results in a dramatic, cell-type-specific blood disorder phenotype.

Population-based studies reveal an additive role of type IV collagen variants in hematuria and albuminuria.

Specific variants in genes that encode the α3α4α5 chains of type IV collagen cause Alport syndrome (AS), which encompass a clinical spectrum from isolated hematuria to multisystem disease affecting sight, hearing and kidney function. The commonest form is X-linked Alport syndrome (XLAS; COL4A5) with autosomal AS (COL4A3 and COL4A4) comprising a minority of cases. While historic data estimates the frequency of AS at 1:5000-10,000, recent population-based genetic studies suggest the prevalence is considerably higher. Genome-wide association studies (GWAS) have been performed in the Icelandic (deCODE) and UK (UK Biobank) populations, demonstrating an association of type IV collagen gene variants with AS relevant kidney traits. In the Icelandic population, 1 in 600 carries a 2.5-kb COL4A3 coding deletion or a COL4A3 missense variant (rs200287952[A], Gly695Arg), both of which are strongly associated with hematuria and albuminuria (P values = 1.9 × 10-5 to 2.5 × 10-20). In the UK Biobank, COL4A4 rs35138315 (Ser969X; carrier frequency 0.13%) is strongly associated with both hematuria and albuminuria (P = 1.5 × 10-73). Thus, the frequency for autosomal AS is 5-16 times higher than the historic prevalence of all forms of the disorder. Furthermore, COL4A4 rs3518315 (Ser969X) is also a reported founder mutation in families with autosomal dominant focal and segmental glomerulosclerosis and autosomal recessive forms of AS. This supports an additive mode of inheritance for specific variants, wherein a number of copies of a mutation influence disease severity in a cumulative fashion. These studies did not include the X chromosome, excluding analysis of COL4A5, which represents an area for future study.

Association of Coding Variants in Hydroxysteroid 17-beta Dehydrogenase 14 (HSD17B14) with Reduced Progression to End Stage Kidney Disease in Type 1 Diabetes.

BACKGROUND: Rare variants in gene coding regions likely have a greater impact on disease-related phenotypes than common variants through disruption of their encoded protein. We searched for rare variants associated with onset of ESKD in individuals with type 1 diabetes at advanced kidney disease stage. METHODS: Gene-based exome array analyses of 15,449 genes in five large incidence cohorts of individuals with type 1 diabetes and proteinuria were analyzed for survival time to ESKD, testing the top gene in a sixth cohort (n=2372/1115 events all cohorts) and replicating in two retrospective case-control studies (n=1072 cases, 752 controls). Deep resequencing of the top associated gene in five cohorts confirmed the findings. We performed immunohistochemistry and gene expression experiments in human control and diseased cells, and in mouse ischemia reperfusion and aristolochic acid nephropathy models. RESULTS: Protein coding variants in the hydroxysteroid 17-ß dehydrogenase 14 gene (HSD17B14), predicted to affect protein structure, had a net protective effect against development of ESKD at exome-wide significance (n=4196; P value=3.3 × 10-7). The HSD17B14 gene and encoded enzyme were robustly expressed in healthy human kidney, maximally in proximal tubular cells. Paradoxically, gene and protein expression were attenuated in human diabetic proximal tubules and in mouse kidney injury models. Expressed HSD17B14 gene and protein levels remained low without recovery after 21 days in a murine ischemic reperfusion injury model. Decreased gene expression was found in other CKD-associated renal pathologies. CONCLUSIONS: HSD17B14 gene is mechanistically involved in diabetic kidney disease. The encoded sex steroid enzyme is a druggable target, potentially opening a new avenue for therapeutic development.

Increased Intestinal Permeability Is Associated With Later Development of Crohn's Disease.

BACKGROUND & AIMS: Increased intestinal permeability has been associated with Crohn's disease (CD), but it is not clear whether it is a cause or result of the disease. We performed a prospective study to determine whether increased intestinal permeability is associated with future development of CD. METHODS: We assessed the intestinal permeability, measured by the urinary fractional excretion of lactulose-to-mannitol ratio (LMR) at recruitment in 1420 asymptomatic first-degree relatives (6-35 years old) of patients with CD (collected from 2008 through 2015). Participants were then followed up for a diagnosis of CD from 2008 to 2017, with a median follow-up time of 7.8 years. We analyzed data from 50 participants who developed CD after a median of 2.7 years during the study period, along with 1370 individuals who remained asymptomatic until October 2017. We used the Cox proportional hazards model to evaluate time-related risk of CD based on the baseline LMR. RESULTS: An abnormal LMR (>0.03) was associated with a diagnosis of CD during the follow-up period (hazard ratio, 3.03; 95% CI, 1.64-5.63; P = 3.97 × 10-4). This association remained significant even when the test was performed more than 3 years before the diagnosis of CD (hazard ratio, 1.62; 95% CI, 1.051-2.50; P = .029). CONCLUSIONS: Increased intestinal permeability is associated with later development of CD; these findings support a model in which altered intestinal barrier function contributes to pathogenesis. Abnormal gut barrier function might serve as a biomarker for risk of CD onset.

Prevalence and Clinical Features of Inflammatory Bowel Diseases Associated With Monogenic Variants, Identified by Whole-Exome Sequencing in 1000 Children at a Single Center.

BACKGROUND & AIMS: A proportion of infants and young children with inflammatory bowel diseases (IBDs) have subtypes associated with a single gene variant (monogenic IBD). We aimed to determine the prevalence of monogenic disease in a cohort of pediatric patients with IBD. METHODS: We performed whole-exome sequencing analyses of blood samples from an unselected cohort of 1005 children with IBD, aged 0-18 years (median age at diagnosis, 11.96 years) at a single center in Canada and their family members (2305 samples total). Variants believed to cause IBD were validated using Sanger sequencing. Biopsies from patients were analyzed by immunofluorescence and histochemical analyses. RESULTS: We identified 40 rare variants associated with 21 monogenic genes among 31 of the 1005 children with IBD (including 5 variants in XIAP, 3 in DOCK8, and 2 each in FOXP3, GUCY2C, and LRBA). These variants occurred in 7.8% of children younger than 6 years and 2.3% of children aged 6-18 years. Of the 17 patients with monogenic Crohn's disease, 35% had abdominal pain, 24% had nonbloody loose stool, 18% had vomiting, 18% had weight loss, and 5% had intermittent bloody loose stool. The 14 patients with monogenic ulcerative colitis or IBD-unclassified received their diagnosis at a younger age, and their most predominant feature was bloody loose stool (78%). Features associated with monogenic IBD, compared to cases of IBD not associated with a single variant, were age of onset younger than 2 years (odds ratio [OR], 6.30; P = .020), family history of autoimmune disease (OR, 5.12; P = .002), extra-intestinal manifestations (OR, 15.36; P < .0001), and surgery (OR, 3.42; P = .042). Seventeen patients had variants in genes that could be corrected with allogeneic hematopoietic stem cell transplantation. CONCLUSIONS: In whole-exome sequencing analyses of more than 1000 children with IBD at a single center, we found that 3% had rare variants in genes previously associated with pediatric IBD. These were associated with different IBD phenotypes, and 1% of the patients had variants that could be potentially corrected with allogeneic hematopoietic stem cell transplantation. Monogenic IBD is rare, but should be considered in analysis of all patients with pediatric onset of IBD.

TaxoNN: ensemble of neural networks on stratified microbiome data for disease prediction.

MOTIVATION: Research supports the potential use of microbiome as a predictor of some diseases. Motivated by the findings that microbiome data is complex in nature, and there is an inherent correlation due to hierarchical taxonomy of microbial Operational Taxonomic Units (OTUs), we propose a novel machine learning method incorporating a stratified approach to group OTUs into phylum clusters. Convolutional Neural Networks (CNNs) were used to train within each of the clusters individually. Further, through an ensemble learning approach, features obtained from each cluster were then concatenated to improve prediction accuracy. Our two-step approach comprising stratification prior to combining multiple CNNs, aided in capturing the relationships between OTUs sharing a phylum efficiently, as compared to using a single CNN ignoring OTU correlations. RESULTS: We used simulated datasets containing 168 OTUs in 200 cases and 200 controls for model testing. Thirty-two OTUs, potentially associated with risk of disease were randomly selected and interactions between three OTUs were used to introduce non-linearity. We also implemented this novel method in two human microbiome studies: (i) Cirrhosis with 118 cases, 114 controls; (ii) type 2 diabetes (T2D) with 170 cases, 174 controls; to demonstrate the model's effectiveness. Extensive experimentation and comparison against conventional machine learning techniques yielded encouraging results. We obtained mean AUC values of 0.88, 0.92, 0.75, showing a consistent increment (5%, 3%, 7%) in simulations, Cirrhosis and T2D data, respectively, against the next best performing method, Random Forest. AVAILABILITY AND IMPLEMENTATION: https://github.com/divya031090/TaxoNN_OTU. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.