Diffusion imaging genomics provides novel insight into early mechanisms of cerebral small vessel disease.

Cerebral small vessel disease (cSVD) is a leading cause of stroke and dementia. Genetic risk loci for white matter hyperintensities (WMH), the most common MRI-marker of cSVD in older age, were recently shown to be significantly associated with white matter (WM) microstructure on diffusion tensor imaging (signal-based) in young adults. To provide new insights into these early changes in WM microstructure and their relation with cSVD, we sought to explore the genetic underpinnings of cutting-edge tissue-based diffusion imaging markers across the adult lifespan. We conducted a genome-wide association study of neurite orientation dispersion and density imaging (NODDI) markers in young adults (i-Share study: N = 1 758, (mean[range]) 22.1[18-35] years), with follow-up in young middle-aged (Rhineland Study: N = 714, 35.2[30-40] years) and late middle-aged to older individuals (UK Biobank: N = 33 224, 64.3[45-82] years). We identified 21 loci associated with NODDI markers across brain regions in young adults. The most robust association, replicated in both follow-up cohorts, was with Neurite Density Index (NDI) at chr5q14.3, a known WMH locus in VCAN. Two additional loci were replicated in UK Biobank, at chr17q21.2 with NDI, and chr19q13.12 with Orientation Dispersion Index (ODI). Transcriptome-wide association studies showed associations of STAT3 expression in arterial and adipose tissue (chr17q21.2) with NDI, and of several genes at chr19q13.12 with ODI. Genetic susceptibility to larger WMH volume, but not to vascular risk factors, was significantly associated with decreased NDI in young adults, especially in regions known to harbor WMH in older age. Individually, seven of 25 known WMH risk loci were associated with NDI in young adults. In conclusion, we identified multiple novel genetic risk loci associated with NODDI markers, particularly NDI, in early adulthood. These point to possible early-life mechanisms underlying cSVD and to processes involving remyelination, neurodevelopment and neurodegeneration, with a potential for novel approaches to prevention.

Contribution of genetic variants in the development of familial premature coronary artery disease in a cohort of cardiac patients.

Coronary artery disease (CAD), the most prevalent cardiovascular disease, is the leading cause of death worldwide. Heritable factors play a significant role in the pathogenesis of CAD. It has been proposed that approximately one-third of patients with CAD have a positive family history, and individuals with such history are at ~1.5-fold increased risk of CAD in their lifespans. Accordingly, the long-recognized familial clustering of CAD is a strong risk factor for this disease. Our study aimed to identify candidate genetic variants contributing to CAD by studying a cohort of 60 large Iranian families with at least two members in different generations afflicted with premature CAD (PCAD), defined as established disease at ≤45 years in men and ≤55 years in women. Exome sequencing was performed for a subset of the affected individuals, followed by prioritization and Sanger sequencing of candidate variants in all available family members. Subsequently, apparently healthy carriers of potential risk variants underwent coronary computed tomography angiography (CCTA), followed by co-segregation analysis of the combined data. Putative causal variants were identified in seven genes, ABCG8, CD36, CYP27A1, PIK3C2G, RASSF9, RYR2, and ZFYVE21, co-segregating with familial PCAD in seven unrelated families. Among these, PIK3C2G, RASSF9, and ZFYVE21 are novel candidate CAD susceptibility genes. Our findings indicate that rare variants in genes identified in this study are involved in CAD development.

Identification of circulating proteins associated with general cognitive function among middle-aged and older adults.

Identifying circulating proteins associated with cognitive function may point to biomarkers and molecular process of cognitive impairment. Few studies have investigated the association between circulating proteins and cognitive function. We identify 246 protein measures quantified by the SomaScan assay as associated with cognitive function (p < 4.9E-5, n up to 7289). Of these, 45 were replicated using SomaScan data, and three were replicated using Olink data at Bonferroni-corrected significance. Enrichment analysis linked the proteins associated with general cognitive function to cell signaling pathways and synapse architecture. Mendelian randomization analysis implicated higher levels of NECTIN2, a protein mediating viral entry into neuronal cells, with higher Alzheimer's disease (AD) risk (p = 2.5E-26). Levels of 14 other protein measures were implicated as consequences of AD susceptibility (p < 2.0E-4). Proteins implicated as causes or consequences of AD susceptibility may provide new insight into the potential relationship between immunity and AD susceptibility as well as potential therapeutic targets.

Genetic Insights on the Relation of Vascular Risk Factors and Cervical Artery Dissection.

BACKGROUND: The association between vascular risk factors and cervical artery dissections (CeADs), a leading cause of ischemic stroke (IS) in the young, remains controversial. OBJECTIVES: This study aimed to explore the causal relation of vascular risk factors with CeAD risk and recurrence and compare it to their relation with non-CeAD IS. METHODS: This study used 2-sample Mendelian randomization analyses to explore the association of blood pressure (BP), lipid levels, type 2 diabetes, waist-to-hip ratio, smoking, and body mass index with CeAD and non-CeAD IS. To simulate effects of the most frequently used BP-lowering drugs, this study constructed genetic proxies and tested their association with CeAD and non-CeAD IS. In analyses among patients with CeAD, the investigators studied the association between weighted genetic risk scores of vascular risk factors and the risk of multiple or early recurrent dissections. RESULTS: Genetically determined higher systolic BP (OR: 1.51; 95% CI: 1.32-1.72) and diastolic BP (OR: 2.40; 95% CI: 1.92-3.00) increased the risk of CeAD (P < 0.0001). Genetically determined higher body mass index was inconsistently associated with a lower risk of CeAD. Genetic proxies for ß-blocker effects were associated with a lower risk of CeAD (OR: 0.65; 95% CI: 0.50-0.85), whereas calcium-channel blockers were associated with a lower risk of non-CeAD IS (OR: 0.75; 95% CI: 0.63-0.90). Weighted genetic risk scores for systolic BP and diastolic BP were associated with an increased risk of multiple or early recurrent CeAD. CONCLUSIONS: These results are supportive of a causal association between higher BP and increased CeAD risk and recurrence and provide genetic evidence for lower CeAD risk under ß-blockers. This may inform secondary prevention strategies and trial design for CeAD.

Dyslipidemia, inflammation, calcification, and adiposity in aortic stenosis: a genome-wide study.

AIMS: Although highly heritable, the genetic etiology of calcific aortic stenosis (AS) remains incompletely understood. The aim of this study was to discover novel genetic contributors to AS and to integrate functional, expression, and cross-phenotype data to identify mechanisms of AS. METHODS AND RESULTS: A genome-wide meta-analysis of 11.6 million variants in 10 cohorts involving 653 867 European ancestry participants (13 765 cases) was performed. Seventeen loci were associated with AS at P ≤ 5 × 10-8, of which 15 replicated in an independent cohort of 90 828 participants (7111 cases), including CELSR2-SORT1, NLRP6, and SMC2. A genetic risk score comprised of the index variants was associated with AS [odds ratio (OR) per standard deviation, 1.31; 95% confidence interval (CI), 1.26-1.35; P = 2.7 × 10-51] and aortic valve calcium (OR per standard deviation, 1.22; 95% CI, 1.08-1.37; P = 1.4 × 10-3), after adjustment for known risk factors. A phenome-wide association study indicated multiple associations with coronary artery disease, apolipoprotein B, and triglycerides. Mendelian randomization supported a causal role for apolipoprotein B-containing lipoprotein particles in AS (OR per g/L of apolipoprotein B, 3.85; 95% CI, 2.90-5.12; P = 2.1 × 10-20) and replicated previous findings of causality for lipoprotein(a) (OR per natural logarithm, 1.20; 95% CI, 1.17-1.23; P = 4.8 × 10-73) and body mass index (OR per kg/m2, 1.07; 95% CI, 1.05-1.9; P = 1.9 × 10-12). Colocalization analyses using the GTEx database identified a role for differential expression of the genes LPA, SORT1, ACTR2, NOTCH4, IL6R, and FADS. CONCLUSION: Dyslipidemia, inflammation, calcification, and adiposity play important roles in the etiology of AS, implicating novel treatments and prevention strategies.

Application of Genomic Sequencing to Refine Patient Stratification for Adjuvant Therapy in Renal Cell Carcinoma.

PURPOSE: Patients with resected localized clear-cell renal cell carcinoma (ccRCC) remain at variable risk of recurrence. Incorporation of biomarkers may refine risk prediction and inform adjuvant treatment decisions. We explored the role of tumor genomics in this setting, leveraging the largest cohort to date of localized ccRCC tissues subjected to targeted gene sequencing. EXPERIMENTAL DESIGN: The somatic mutation status of 12 genes was determined in 943 ccRCC cases from a multinational cohort of patients, and associations to outcomes were examined in a Discovery (n = 469) and Validation (n = 474) framework. RESULTS: Tumors containing a von-Hippel Lindau (VHL) mutation alone were associated with significantly improved outcomes in comparison with tumors containing a VHL plus additional mutations. Within the Discovery cohort, those with VHL+0, VHL+1, VHL+2, and VHL+≥3 tumors had disease-free survival (DFS) rates of 90.8%, 80.1%, 68.2%, and 50.7% respectively, at 5 years. This trend was replicated in the Validation cohort. Notably, these genomically defined groups were independent of tumor mutational burden. Amongst patients eligible for adjuvant therapy, those with a VHL+0 tumor (29%) had a 5-year DFS rate of 79.3% and could, therefore, potentially be spared further treatment. Conversely, patients with VHL+2 and VHL+≥3 tumors (32%) had equivalent DFS rates of 45.6% and 35.3%, respectively, and should be prioritized for adjuvant therapy. CONCLUSIONS: Genomic characterization of ccRCC identified biologically distinct groups of patients with divergent relapse rates. These groups account for the ∼80% of cases with VHL mutations and could be used to personalize adjuvant treatment discussions with patients as well as inform future adjuvant trial design.

Gene-mapping study of extremes of cerebral small vessel disease reveals TRIM47 as a strong candidate.

Cerebral small vessel disease is a leading cause of stroke and a major contributor to cognitive decline and dementia, but our understanding of specific genes underlying the cause of sporadic cerebral small vessel disease is limited. We report a genome-wide association study and a whole-exome association study on a composite extreme phenotype of cerebral small vessel disease derived from its most common MRI features: white matter hyperintensities and lacunes. Seventeen population-based cohorts of older persons with MRI measurements and genome-wide genotyping (n = 41 326), whole-exome sequencing (n = 15 965), or exome chip (n = 5249) data contributed 13 776 and 7079 extreme small vessel disease samples for the genome-wide association study and whole-exome association study, respectively. The genome-wide association study identified significant association of common variants in 11 loci with extreme small vessel disease, of which the chr12q24.11 locus was not previously reported to be associated with any MRI marker of cerebral small vessel disease. The whole-exome association study identified significant associations of extreme small vessel disease with common variants in the 5' UTR region of EFEMP1 (chr2p16.1) and one probably damaging common missense variant in TRIM47 (chr17q25.1). Mendelian randomization supports the causal association of extensive small vessel disease severity with increased risk of stroke and Alzheimer's disease. Combined evidence from summary-based Mendelian randomization studies and profiling of human loss-of-function allele carriers showed an inverse relation between TRIM47 expression in the brain and blood vessels and extensive small vessel disease severity. We observed significant enrichment of Trim47 in isolated brain vessel preparations compared to total brain fraction in mice, in line with the literature showing Trim47 enrichment in brain endothelial cells at single cell level. Functional evaluation of TRIM47 by small interfering RNAs-mediated knockdown in human brain endothelial cells showed increased endothelial permeability, an important hallmark of cerebral small vessel disease pathology. Overall, our comprehensive gene-mapping study and preliminary functional evaluation suggests a putative role of TRIM47 in the pathophysiology of cerebral small vessel disease, making it an important candidate for extensive in vivo explorations and future translational work.

Differentially methylated CpGs in response to growth hormone administration in children with idiopathic short stature.

BACKGROUND: Recombinant human growth hormone (rhGH) has shown a great growth-promoting potential in children with idiopathic short stature (ISS). However, the response to rhGH differs across individuals, largely due to genetic and epigenetic heterogeneity. Since epigenetic marks on the methylome can be dynamically influenced by GH, we performed a comprehensive pharmacoepigenomics analysis of DNA methylation changes associated with long-term rhGH administration in children with ISS. RESULTS: We measured DNA methylation profiles before and after GH treatment (with a duration of ~ 18 months in average) on 47 healthy children using customized methylC-seq capture sequencing. Their changes were compared and associated with changes in plasma IGF1 by adjusting sex, age, treatment duration and estimated blood proportions. We observed a considerable inter-individual heterogeneity of DNA methylation changes responding to GH treatment. We identified 267 response-associated differentially methylated cytosines (DMCs) that were enriched in promoter regions, CpG islands and blood cell-type-specific regulatory elements. Furthermore, the genes associated with these DMCs were enriched in the biology process of "cell development," "neuron differentiation" and "developmental growth," and in the TGF-beta signaling pathway, PPAR Alpha pathway, endoderm differentiation pathway, adipocytokine signaling pathway as well as PI3K-Akt signaling pathway, and cAMP signaling pathway. CONCLUSION: Our study provides a first insight in DNA methylation changes associated with rhGH administration, which may help understand mechanisms of epigenetic regulation on GH-responsive genes.

Cohort profile: genomic data for 26 622 individuals from the Canadian Longitudinal Study on Aging (CLSA).

PURPOSE: The Canadian Longitudinal Study on Aging (CLSA) Comprehensive cohort was established to provide unique opportunities to study the genetic and environmental contributions to human disease as well as ageing process. The aim of this report was to describe the genomic data included in CLSA. PARTICIPANTS: A total of 26 622 individuals from the CLSA Comprehensive cohort of men and women aged 45-85 recruited between 2010 and 2015 underwent genome-wide genotyping of DNA samples collected from blood. Comprehensive quality control metrics were measured for genetic markers and samples, respectively. The genotypes were imputed to the TOPMed reference panel. Sex chromosome abnormalities were identified by copy number profiling. Classical human leukocyte antigen gene haplotypes were imputed at two-field (four-digit). FINDINGS TO DATE: Of the 26 622 genotyped participants, 24 655 (92.6%) were identified as having European ancestry. These genomic data were linked to physical, lifestyle, medical, economic, environmental and psychosocial factors collected longitudinally in CLSA. The combined analysis, including CLSA genomic data, uncovered over 100 novel loci associated with key parameters to define glaucoma. The CLSA genomic dataset validated the contribution of a polygenic risk score to screen individuals with high fracture risk. It is also a valuable resource to directly identify common genetic variations associated with conditions related to complex traits. Taking advantage of the comprehensive interview and physical information collected in CLSA, this genomic dataset has been linked to psychosocial factors to investigate both the independent and interactive effects on cardiovascular disease. FUTURE PLANS: The CLSA overall is ongoing. Follow-up data will continue to be collected from participants in the current genomic subcohort, including the DNA methylation and metabolomic data. Ongoing studies focus on elucidating the role of genetic factors in cognitive decline and cardiovascular diseases. This genomic data resource is available on request through the CLSA data access application process.

Genomic Studies Across the Lifespan Point to Early Mechanisms Determining Subcortical Volumes.

BACKGROUND: Subcortical brain structures play a key role in pathological processes of age-related neurodegenerative disorders. Mounting evidence also suggests that early-life factors may have an impact on the development of common late-life neurological diseases, including genetic factors that can influence both brain maturation and neurodegeneration. METHODS: Using large population-based brain imaging datasets across the lifespan (N ≤ 40,628), we aimed to 1) estimate the heritability of subcortical volumes in young (18-35 years), middle (35-65 years), and older (65+ years) age, and their genetic correlation across age groups; 2) identify whether genetic loci associated with subcortical volumes in older persons also show associations in early adulthood, and explore underlying genes using transcriptome-wide association studies; and 3) explore their association with neurological phenotypes. RESULTS: Heritability of subcortical volumes consistently decreased with increasing age. Genetic risk scores for smaller caudate nucleus, putamen, and hippocampus volume in older adults were associated with smaller volumes in young adults. Individually, 10 loci associated with subcortical volumes in older adults also showed associations in young adults. Within these loci, transcriptome-wide association studies showed that expression of several genes in brain tissues (especially MYLK2 and TUFM) was associated with subcortical volumes in both age groups. One risk variant for smaller caudate nucleus volume (TUFM locus) was associated with lower cognitive performance. Genetically predicted Alzheimer's disease was associated with smaller subcortical volumes in middle and older age. CONCLUSIONS: Our findings provide novel insights into the genetic determinants of subcortical volumes across the lifespan. More studies are needed to decipher the underlying biology and clinical impact.

Common Genetic Variants Contribute to Risk of Transposition of the Great Arteries.

RATIONALE: Dextro-transposition of the great arteries (D-TGA) is a severe congenital heart defect which affects approximately 1 in 4,000 live births. While there are several reports of D-TGA patients with rare variants in individual genes, the majority of D-TGA cases remain genetically elusive. Familial recurrence patterns and the observation that most cases with D-TGA are sporadic suggest a polygenic inheritance for the disorder, yet this remains unexplored. OBJECTIVE: We sought to study the role of common single nucleotide polymorphisms (SNPs) in risk for D-TGA. METHODS AND RESULTS: We conducted a genome-wide association study in an international set of 1,237 patients with D-TGA and identified a genome-wide significant susceptibility locus on chromosome 3p14.3, which was subsequently replicated in an independent case-control set (rs56219800, meta-analysis P=8.6x10-10, OR=0.69 per C allele). SNP-based heritability analysis showed that 25% of variance in susceptibility to D-TGA may be explained by common variants. A genome-wide polygenic risk score derived from the discovery set was significantly associated to D-TGA in the replication set (P=4x10-5). The genome-wide significant locus (3p14.3) co-localizes with a putative regulatory element that interacts with the promoter of WNT5A, which encodes the Wnt Family Member 5A protein known for its role in cardiac development in mice. We show that this element drives reporter gene activity in the developing heart of mice and zebrafish and is bound by the developmental transcription factor TBX20. We further demonstrate that TBX20 attenuates Wnt5a expression levels in the developing mouse heart. CONCLUSIONS: This work provides support for a polygenic architecture in D-TGA and identifies a susceptibility locus on chromosome 3p14.3 near WNT5A. Genomic and functional data support a causal role of WNT5A at the locus.

An innate contribution of human nicotinic receptor polymorphisms to COPD-like lesions.

Chronic Obstructive Pulmonary Disease is a generally smoking-linked major cause of morbidity and mortality. Genome-wide Association Studies identified a locus including a non-synonymous single nucleotide polymorphism in CHRNA5, rs16969968, encoding the nicotinic acetylcholine receptor α5 subunit, predisposing to both smoking and Chronic Obstructive Pulmonary Disease. Here we report that nasal polyps from rs16969968 non-smoking carriers exhibit airway epithelium remodeling and inflammation. These hallmarks of Chronic Obstructive Pulmonary Disease occur spontaneously in mice expressing human rs16969968. They are significantly amplified after exposure to porcine pancreatic elastase, an emphysema model, and to oxidative stress with a polymorphism-dependent alteration of lung function. Targeted rs16969968 expression in epithelial cells leads to airway remodeling in vivo, increased proliferation and production of pro-inflammatory cytokines through decreased calcium entry and increased adenylyl-cyclase activity. We show that rs16969968 directly contributes to Chronic Obstructive Pulmonary Disease-like lesions, sensitizing the lung to the action of oxidative stress and injury, and represents a therapeutic target.

A small number of early introductions seeded widespread transmission of SARS-CoV-2 in Québec, Canada.

BACKGROUND: Québec was the Canadian province most impacted by COVID-19, with 401,462 cases as of September 24th, 2021, and 11,347 deaths due mostly to a very severe first pandemic wave. In April 2020, we assembled the Coronavirus Sequencing in Québec (CoVSeQ) consortium to sequence SARS-CoV-2 genomes in Québec to track viral introduction events and transmission within the province. METHODS: Using genomic epidemiology, we investigated the arrival of SARS-CoV-2 to Québec. We report 2921 high-quality SARS-CoV-2 genomes in the context of > 12,000 publicly available genomes sampled globally over the first pandemic wave (up to June 1st, 2020). By combining phylogenetic and phylodynamic analyses with epidemiological data, we quantify the number of introduction events into Québec, identify their origins, and characterize the spatiotemporal spread of the virus. RESULTS: Conservatively, we estimated approximately 600 independent introduction events, the majority of which happened from spring break until 2 weeks after the Canadian border closed for non-essential travel. Subsequent mass repatriations did not generate large transmission lineages (> 50 sequenced cases), likely due to mandatory quarantine measures in place at the time. Consistent with common spring break and "snowbird" destinations, most of the introductions were inferred to have originated from Europe via the Americas. Once introduced into Québec, viral lineage sizes were overdispersed, with a few lineages giving rise to most infections. Consistent with founder effects, the earliest lineages to arrive tended to spread most successfully. Fewer than 100 viral introductions arrived during spring break, of which 7-12 led to the largest transmission lineages of the first wave (accounting for 52-75% of all sequenced infections). These successful transmission lineages dispersed widely across the province. Transmission lineage size was greatly reduced after March 11th, when a quarantine order for returning travellers was enacted. While this suggests the effectiveness of early public health measures, the biggest transmission lineages had already been ignited prior to this order. CONCLUSIONS: Combined, our results reinforce how, in the absence of tight travel restrictions or quarantine measures, fewer than 100 viral introductions in a week can ensure the establishment of extended transmission chains.

Integrated genomics point to immune vulnerabilities in pleural mesothelioma.

Pleural mesothelioma is an aggressive malignancy with limited effective therapies. In order to identify therapeutic targets, we integrated SNP genotyping, sequencing and transcriptomics from tumours and low-passage patient-derived cells. Previously unrecognised deletions of SUFU locus (10q24.32), observed in 21% of 118 tumours, resulted in disordered expression of transcripts from Hedgehog pathways and the T-cell synapse including VISTA. Co-deletion of Interferon Type I genes and CDKN2A was present in half of tumours and was a predictor of poor survival. We also found previously unrecognised deletions in RB1 in 26% of cases and show sub-micromolar responses to downstream PLK1, CHEK1 and Aurora Kinase inhibitors in primary mesothelioma cells. Defects in Hippo pathways that included RASSF7 amplification and NF2 or LATS1/2 mutations were present in 50% of tumours and were accompanied by micromolar responses to the YAP1 inhibitor Verteporfin. Our results suggest new therapeutic avenues in mesothelioma and indicate targets and biomarkers for immunotherapy.

Y disruption, autosomal hypomethylation and poor male lung cancer survival.

Lung cancer is the most frequent cause of cancer death worldwide. It affects more men than women, and men generally have worse survival outcomes. We compared gene co-expression networks in affected and unaffected lung tissue from 126 consecutive patients with Stage IA-IV lung cancer undergoing surgery with curative intent. We observed marked degradation of a sex-associated transcription network in tumour tissue. This disturbance, detected in 27.7% of male tumours in the discovery dataset and 27.3% of male tumours in a further 123-sample replication dataset, was coincident with partial losses of the Y chromosome and extensive autosomal DNA hypomethylation. Central to this network was the epigenetic modifier and regulator of sexually dimorphic gene expression, KDM5D. After accounting for prognostic and epidemiological covariates including stage and histology, male patients with tumour KDM5D deficiency showed a significantly increased risk of death (Hazard Ratio [HR] 3.80, 95% CI 1.40-10.3, P = 0.009). KDM5D deficiency was confirmed as a negative prognostic indicator in a further 1100 male lung tumours (HR 1.67, 95% CI 1.4-2.0, P = 1.2 × 10-10). Our findings identify tumour deficiency of KDM5D as a prognostic marker and credible mechanism underlying sex disparity in lung cancer.

Rare loss-of-function variants in type I IFN immunity genes are not associated with severe COVID-19.

A recent report found that rare predicted loss-of-function (pLOF) variants across 13 candidate genes in TLR3- and IRF7-dependent type I IFN pathways explain up to 3.5% of severe COVID-19 cases. We performed whole-exome or whole-genome sequencing of 1,864 COVID-19 cases (713 with severe and 1,151 with mild disease) and 15,033 ancestry-matched population controls across 4 independent COVID-19 biobanks. We tested whether rare pLOF variants in these 13 genes were associated with severe COVID-19. We identified only 1 rare pLOF mutation across these genes among 713 cases with severe COVID-19 and observed no enrichment of pLOFs in severe cases compared to population controls or mild COVID-19 cases. We found no evidence of association of rare LOF variants in the 13 candidate genes with severe COVID-19 outcomes.

Morphological findings in frozen non-neoplastic kidney tissues of patients with kidney cancer from large-scale multicentric studies on renal cancer.

There are unexplained geographical variations in the incidence of kidney cancer with the high rates reported in Baltic countries, as well as eastern and central Europe. Having access to a large and well-annotated collection of "tumor/non-tumor" pairs of kidney cancer patients from the Czech Republic, Romania, Serbia, UK, and Russia, we aimed to analyze the morphology of non-neoplastic renal tissue in nephrectomy specimens. By applying digital pathology, we performed a microscopic examination of 1012 frozen non-neoplastic kidney tissues from patients with renal cell carcinoma. Four components of renal parenchyma were evaluated and scored for the intensity of interstitial inflammation and fibrosis, tubular atrophy, glomerulosclerosis, and arterial wall thickening, globally called chronic renal parenchymal changes. Moderate or severe changes were observed in 54 (5.3%) of patients with predominance of occurrence in Romania (OR = 2.67, CI 1.07-6.67) and Serbia (OR = 4.37, CI 1.20-15.96) in reference to those from Russia. Further adjustment for comorbidities, tumor characteristics, and stage did not change risk estimates. In multinomial regression model, relative probability of non-glomerular changes was 5.22 times higher for Romania and Serbia compared to Russia. Our findings show that the frequency of chronic renal parenchymal changes, with the predominance of chronic interstitial nephritis pattern, in kidney cancer patients varies by country, significantly more frequent in countries located in central and southeastern Europe where the incidence of kidney cancer has been reported to be moderate to high. The observed association between these pathological features and living in certain geographic areas requires a larger population-based study to confirm this association on a large scale.

A Distributed Whole Genome Sequencing Benchmark Study.

Population sequencing often requires collaboration across a distributed network of sequencing centers for the timely processing of thousands of samples. In such massive efforts, it is important that participating scientists can be confident that the accuracy of the sequence data produced is not affected by which center generates the data. A study was conducted across three established sequencing centers, located in Montreal, Toronto, and Vancouver, constituting Canada's Genomics Enterprise (www.cgen.ca). Whole genome sequencing was performed at each center, on three genomic DNA replicates from three well-characterized cell lines. Secondary analysis pipelines employed by each site were applied to sequence data from each of the sites, resulting in three datasets for each of four variables (cell line, replicate, sequencing center, and analysis pipeline), for a total of 81 datasets. These datasets were each assessed according to multiple quality metrics including concordance with benchmark variant truth sets to assess consistent quality across all three conditions for each variable. Three-way concordance analysis of variants across conditions for each variable was performed. Our results showed that the variant concordance between datasets differing only by sequencing center was similar to the concordance for datasets differing only by replicate, using the same analysis pipeline. We also showed that the statistically significant differences between datasets result from the analysis pipeline used, which can be unified and updated as new approaches become available. We conclude that genome sequencing projects can rely on the quality and reproducibility of aggregate data generated across a network of distributed sites.

Shared genomic segment analysis with equivalence testing.

An important aspect of disease gene mapping is replication, that is, a putative finding in one group of individuals is confirmed in another set of individuals. As it can happen by chance that individuals share an estimated disease position, we developed a statistical approach to determine the p-value for multiple individuals or families to share a possibly small number of candidate susceptibility variants. Here, we focus on candidate variants for dominant traits that have been obtained by our previously developed heterozygosity analysis, and we are testing the sharing of candidate variants obtained for different individuals. Our approach allows for multiple pathogenic variants in a gene to contribute to disease, and for estimated disease variant positions to be imprecise. Statistically, the method developed here falls into the category of equivalence testing, where the classical null and alternative hypotheses of homogeneity and heterogeneity are reversed. The null hypothesis situation is created by permuting genomic locations of variants for one individual after another. We applied our methodology to the ALSPAC data set of 1,927 whole-genome sequenced individuals, where some individuals carry a pathogenic variant for the BRCA1 gene, but no two individuals carry the same variant. Our shared genomic segment analysis found significant evidence for BRCA1 pathogenic variants within ±5 kb of a given DNA variant.