BACKGROUND: The relationship between rare variants in Ring finger protein 213 (RNF213) and intracranial atherosclerosis (ICAS) remained unelucidated. Using whole-exome sequencing (WES) and high-resolution magnetic resonance imaging (HR-MRI), this study aimed at investigating the association between rare RNF213 variants and ICAS within a Chinese community-dwelling population.
METHODS: The present study included 821 participants from Shunyi cohort. Genetic data of rare RNF213 variants were acquired by WES and were categorized by functional domains. Intracranial and extracranial atherosclerosis were assessed by brain HR-MRI and carotid ultrasound, respectively. Logistic regression and generalized linear regression were applied to evaluate the effects of rare RNF213 variants on atherosclerosis. Stratification by age were conducted with 50 years old set as the cutoff value.
RESULTS: Ninety-five participants were identified as carriers of rare RNF213 variants. Carotid plaques were observed in 367 (44.7 %) participants, while ICAS was identified in 306 (37.3 %). Rare variants of RNF213 was not associated with ECAS. Employing HR-MRI, both the presence of rare variants (β = 0.150, P = 0.025) and numerical count of variants (β = 0.182, P = 0.003) were significantly correlated with ICAS within the group of age ≤50 years. Both variant existence (β = 0.154, P = 0.014) and variant count (β = 0.188, P = 0.003) were significantly associated with plaques in middle cerebral arteries within younger subgroup, rather than basilar arteries. Furthermore, a significant association was observed between variants that located outside the N-arm domain and ICAS in the younger subgroup (OR = 2.522, P = 0.030). Statistical results remained robust after adjusted for age, gender, and cardiovascular risk factors.
CONCLUSIONS: Rare variants of RNF213 is associated with age-related ICAS in general Chinese population, highlighting the potential role of RNF213 as a genetic contributor to early-onset ICAS.

Collagen type XVII α1 (COL17A1) encodes a hemidesmosomal protein at the epidermal-dermal junction and its variants are implicated in blistering skin diseases. Recent experiments in rodents revealed that Col17a1 has critical roles in stem cells of epidermal origin and in melanoma carcinogenesis. In the present study, it was investigated whether germline variants in COL17A1 are associated with skin cancer and other cancer types using indexed consecutive autopsy cases from the Japanese Geriatric Single Nucleotide Polymorphism database (n=2,343; mean age, 80 years). The database included 12 patients with skin cancer. A total of 53 COL17A1 missense variants on an exome chip were analyzed. One variant, p.Ser1029Ala (rs118166857), which had a minor allele frequency of 1.0%, exhibited a nominal positive sign of association with skin cancer [Fisher\'s exact P=0.002, odds ratio (OR)=16.93, 95% CI: 4.44-64.64]. This variant was detected in 2/2 patients with mucosal malignant melanoma (mMM) and 1/3 patients with extramammary Paget\'s disease, and in none of the patients with non-melanoma cancer, e.g., squamous cell and basal cell carcinoma. Other cancer types were searched in the database and the p.Ser1029Ala variant was indicated to be nominally associated with breast cancer (P=0.006, OR=4.17, 95% CI: 1.72-10.11). In the two mMM cases, targeted exome sequencing of 55 cancer-predisposing genes (including tumor protein 53, BRCA1/2 and mismatch repair genes) detected no apparent pathogenic variants, but revealed variants of unknown significance in axin 2, DNA directed polymerase ζ catalytic subunit and contactin 6. Since COL17A1 provides a niche for melanocyte stem cells, it was hypothesized that the p.Ser1029Ala variant in the COL17A1 ectodomain may affect the microenvironment, e.g., the cell competition. This is a working hypothesis generated from human autopsy cases and warrants further epidemiological and molecular biological validation.

Major Depressive Disorder (MDD) is a complex neuropsychiatric disease with known genetic associations, but without known links to rare variation in the human genome. Here we aim to identify rare genetic variants associated with MDD using deep whole-genome sequencing data in an independent population.
We report the sequencing of 1,688 whole genomes in a large sample of male-male Veteran twins. Depression status was classified based on a structured diagnostic interview according to DSM-III-R diagnostic criteria. Searching only rare variants in genomic regions from recent GWAS on MDD, we used the optimised sequence kernel association test and Fisher\'s Exact test to fine map loci associated with severe depression.
Our analysis identified one gene associated with severe depression, basic helix loop helix e22 (PAdjusted = 0.03) via SKAT-O test between unrelated severely depressed cases compared to unrelated non-depressed controls. The same gene BHLHE22 had a non-silent variant rs13279074 (PAdjusted = 0.032) based on a single variant Fisher\'s Exact test between unrelated severely depressed cases compared to unrelated non-depressed controls.
The gene BHLHE22 shows compelling genetic evidence of directly impacting the severe depression phenotype. Together these results advance understanding of the genetic contribution to major depressive disorder in a new cohort and link a rare variant to severe forms of the disorder.

BACKGROUND: Bromodomain-containing protein 7 (BRD7), a member of the bromodomain-containing protein family, plays important roles in chromatin modification and transcriptional regulation. A recent model of Brd7-knockout mice presented azoospermia and male infertility, implying the potential role of BRD7 in spermatogenic failure in humans. This case-control study aimed to explore the association of the BRD7 gene with spermatogenic efficiency and the risk of spermatogenic defects in humans.
RESULTS: A total of six heterozygous variants were detected in the coding and splicing regions of the BRD7 gene in patients with azoospermia. For each of four rare variants predicted to potentially damage BRD7 function, we further identified these four variants in oligozoospermia and normozoospermia as well. However, no difference in the allele and genotype frequencies of rare variants were observed between cases with spermatogenic failure and controls with normozoospermia; the sperm products of variant carriers were similar to those of noncarriers. Moreover, similar distribution of the alleles, genotypes and haplotypes of seven tag single nucleotide polymorphisms (tagSNPs) was observed between the cases with azoospermia and oligozoospermia and controls with normozoospermia; associations of tagSNP-distinguished BRD7 alleles with sperm products were not identified.
CONCLUSIONS: The lack of an association of BRD7-linked rare and common variants with spermatogenic failure implied a limited contribution of the BRD7 gene to spermatogenic efficiency and susceptibility to male infertility in humans.
RéSUMé: CONTEXTE: Le bromodomaine contenant la protéine 7 (BRD7), un membre de la famille du bromodomaine contenant des protéines, joue des rôles importants dans la modification de la chromatine et la régulation transcriptionnelle. Un modèle récent de souris Brd7-knockout présentait une azoospermie et une infertilité mâle, ce qui implique un rôle potentiel de BRD7 dans l’altération de la spermatogenèse chez l’homme. Cette étude cas-témoins visait à explorer l’association du gène BRD7 avec l’efficacité de la spermatogenèse et le risque d’altérations spermatogéniques chez l’homme. RéSULTATS: Un total de six variants hétérozygotes ont été détectés dans les régions de codage et d’épissage du gène BRD7 chez les patients présentant une azoospermie. Pour chacun des quatre variants rares prédits pour potentiellement endommager la fonction BRD7, nous avons en outre identifié ces quatre variants dans l’oligozoospermie et la normozoospermie. Cependant, nous n’avons observé aucune différence dans les fréquences d’allèle et de génotype des variants rares entre les cas avec altérations de la spermatogenèse et les témoins avec normozoospermie ; les produits du sperme des porteurs de variants étaient semblables à ceux des non-porteurs. Par ailleurs, on a observé une distribution semblable des allèles, des génotypes et des haplotypes de sept polymorphismes simples de nucléotide de balise (tagSNPs) entre les cas avec azoospermie ou oligozoospermie et les témoins normozoospermiques; aucune association n’a pas été identifiée entre les allèles BRD7 tagSNP-distingués et des produits du sperme. CONCLUSION: L’absence d’association des variants rares liés à BRD7 et des variants communs liés à BRD7 avec les altérations de la spermatogenèse implique une contribution limitée du gène BRD7 à l’efficacité spermatogénique et à la susceptibilité à l’infertilité masculine chez l’homme.

Endometrial cancer is the fourth most commonly diagnosed cancer in women. Family history is a known risk factor for endometrial cancer. The incidence of endometrial cancer in a first-degree relative elevates the relative risk to range between 1.3 and 2.8. It is unclear to what extent or what other novel germline variants are at play in endometrial cancer. We aim to address this question by utilizing whole exome sequencing as a means to identify novel, rare variant associations between exonic regions and endometrial cancer. The MyCode community health initiative is an excellent resource for this study with germline whole exome data for 60,000 patients available in the first phase, and further 30,000 patients independently sequenced in the second phase as part of DiscovEHR study. We conducted exome-wide rare variant association using 472 cases and 4,110 controls in 60,000 patients (discovery cohort); and 261 cases and 1,531 controls from 30,000 patients (replication cohort). After binning rare germline variants into genes, case-control association tests performed using Optimal Unified Approach for Rare-Variant Association, SKAT-O. Seven genes, including RBM12, NDUFB6, ATP6V1A, RECK, SLC35E1, RFX3 (Bonferroni-corrected P < 0.05) and ATP8A1 (suggestive P < 10-5), and one long non-coding RNA, DLGAP4-AS1 (Bonferroni-corrected P < 0.05), were associated with endometrial cancer. Notably, RECK, and ATP8A1 were replicated from the replication cohort (suggestive threshold P < 0.05). Additionally, a pathway-based rare variant analysis, using pathogenic and likely pathogenic variants, identified two significant pathways, pyrimidine metabolism and protein processing in the endoplasmic reticulum (Bonferroni-corrected P < 0.05). In conclusion, our results using the single-source electronic health records (EHR) linked to genomic data highlights candidate genes and pathways associated with endometrial cancer and indicates rare variants involvement in endometrial cancer predisposition, which could help in personalized prognosis and also further our understanding of its genetic etiology.

Background: Associations of both common and rare genetic variants with fasting blood lipids have been extensively studied. However, most of the rare coding variants associated with lipids are population-specific, and exploration of genetic data from diverse population samples may enhance the identification of novel associations with rare variants. Results: We searched for novel coding genetic variants associated with fasting lipid levels in 894 samples from the Genetics of Lipid Lowering Drugs and Diet Network (GOLDN) with exome-wide sequencing-based genotype data. In single variant tests, one variant (rs11171663 in ITGA7) was associated with fasting triglyceride levels (P = 7.66E-08), explaining approximately 3.2% of the total trait variance. In gene-based tests, we found statistically significant associations between ITGA7 (P = 1.77E-07) and SLCO2A1 (P = 7.18E-07) and triglycerides, as well as between POT1 (P = 3.00E-07) and low-density lipoprotein cholesterol. In another independent replication cohort consisting of 3,183 African American samples from Hypertension Genetic Epidemiology Network (HyperGEN) and the Genetic Epidemiology Network of Arteriopathy (GENOA), the top genes achieved P-values of 0.04 (ITGA7), 0.08 (SLCO2A1), and 0.02 (POT1). In GOLDN, gene transcript levels of ITGA7 and SLCO2A1 were associated with fasting triglycerides (P = 0.07 and P = 0.02), highlighting functional relevance of our findings. Conclusion: In this study, we present preliminary evidence of novel rare variant determinants of fasting lipids, and reveal potential underlying molecular mechanisms. Moreover, these results were replicated in an independent cohort. Our findings may inform novel biomarkers of disease risk and treatment targets.

Identifying rare coding variants associated with albuminuria may open new avenues for preventing chronic kidney disease and end-stage renal disease, which are highly prevalent in individuals with diabetes. Efforts to identify genetic susceptibility variants for albuminuria have so far been limited, with the majority of studies focusing on common variants.
We performed an exome-wide association study to identify coding variants in a two-stage (discovery and replication) approach. Data from 33,985 individuals of European ancestry (15,872 with and 18,113 without diabetes) and 2605 Greenlanders were included.
We identified a rare (minor allele frequency [MAF]: 0.8%) missense (A1690V) variant in CUBN (rs141640975, β = 0.27, p = 1.3 × 10-11) associated with albuminuria as a continuous measure in the combined European meta-analysis. The presence of each rare allele of the variant was associated with a 6.4% increase in albuminuria. The rare CUBN variant had an effect that was three times stronger in individuals with type 2 diabetes compared with those without (pinteraction = 7.0 × 10-4, β with diabetes = 0.69, β without diabetes = 0.20) in the discovery meta-analysis. Gene-aggregate tests based on rare and common variants identified three additional genes associated with albuminuria (HES1, CDC73 and GRM5) after multiple testing correction (pBonferroni < 2.7 × 10-6).
The current study identifies a rare coding variant in the CUBN locus and other potential genes associated with albuminuria in individuals with and without diabetes. These genes have been implicated in renal and cardiovascular dysfunction. The findings provide new insights into the genetic architecture of albuminuria and highlight target genes and pathways for the prevention of diabetes-related kidney disease.

3q29 deletion syndrome is caused by a recurrent hemizygous 1.6 Mb deletion on the long arm of chromosome 3. The syndrome is rare (1 in 30,000 individuals) and is associated with mild to moderate intellectual disability, increased risk for autism and anxiety, and a 40-fold increased risk for schizophrenia, along with a host of physical manifestations. However, the disorder is poorly characterized, the range of manifestations is not well described, and the underlying molecular mechanism is not understood. We designed the Emory 3q29 Project to document the range of neurodevelopmental and psychiatric manifestations associated with 3q29 deletion syndrome. We will also create a biobank of samples from our 3q29 deletion carriers for mechanistic studies, which will be a publicly-available resource for qualified investigators. The ultimate goals of our study are three-fold: first, to improve management and treatment of 3q29 deletion syndrome. Second, to uncover the molecular mechanism of the disorder. Third, to enable cross-disorder comparison with other rare genetic syndromes associated with neuropsychiatric phenotypes.
We will ascertain study subjects, age 6 and older, from our existing registry ( 3q29deletion.org ). Participants and their families will travel to Atlanta, GA for phenotypic assessments, with particular emphasis on evaluation of anxiety, cognitive ability, autism symptomatology, and risk for psychosis via prodromal symptoms and syndromes. Evaluations will be performed using standardized instruments. Structural, diffusion, and resting-state functional MRI data will be collected from eligible study participants. We will also collect blood from the 3q29 deletion carrier and participating family members, to be banked at the NIMH Repository and Genomics Resource (NRGR).
The study of 3q29 deletion has the potential to transform our understanding of complex disease. Study of individuals with the deletion may provide insights into long term care and management of the disorder. Our project describes the protocol for a prospective study of the behavioral and clinical phenotype associated with 3q29 deletion syndrome. The paradigm described here could easily be adapted to study additional CNV or single gene disorders with high risk for neuropsychiatric phenotypes, and/or transferred to other study sites, providing a means for data harmonization and cross-disorder analysis.

Our understanding of genetic influences on the response of lipids to specific interventions is limited. In this study, we sought to elucidate effects of rare genetic variants on lipid response to a high-fat meal challenge and fenofibrate (FFB) therapy in the Genetics of Lipid Lowering Drugs and Diet Network (GOLDN) cohort using an exome-wide sequencing-based association study. Our results showed that the rare coding variants in ITGA7, SIPA1L2, and CEP72 are significantly associated with fasting LDL cholesterol response to FFB (P = 1.24E-07), triglyceride postprandial area under the increase (AUI) (P = 2.31E-06), and triglyceride postprandial AUI response to FFB (P = 1.88E-06), respectively. We sought to replicate the association for SIPA1L2 in the Heredity and Phenotype Intervention (HAPI) Heart Study, which included a high-fat meal challenge but not FFB treatment. The associated rare variants in GOLDN were not observed in the HAPI Heart study, and thus the gene-based result was not replicated. For functional validation, we found that gene transcript level of SIPA1L2 is associated with triglyceride postprandial AUI (P < 0.05) in GOLDN. Our study suggests unique genetic mechanisms contributing to the lipid response to the high-fat meal challenge and FFB therapy.

Next-generation sequencing-based genetic association study (GAS) is a powerful tool to identify candidate disease variants and genomic regions. Although low-coverage sequencing offers low cost but inadequacy in calling rare variants, high coverage is able to detect essentially every variant but at a high cost. Two-stage sequencing may be an economical way to conduct GAS without losing power. In two-stage sequencing, an affordable number of samples are sequenced at high coverage as the reference panel, then to impute in a larger sample is sequenced at low coverage. As unit sequencing costs continue to decrease, investigators can now conduct GAS with more flexible sequencing depths. Here, we systematically evaluate the effect of the read depth and sample size on the variant discovery power and association power for study designs using low-coverage, high-coverage, and two-stage sequencing. We consider 12 low-coverage, 12 high-coverage, and 51 two-stage design scenarios with the read depth varying from 0.5× to 80×. With state-of-the-art simulation and analysis packages and in-house scripts, we simulate the complete study process from DNA sequencing to SNP (single nucleotide polymorphism) calling and association testing. Our results show that with appropriate allocation of sequencing effort, two-stage sequencing is an effective approach for conducting GAS. We provide practical guidelines for investigators to plan the optimum sequencing-based GAS including two-stage sequencing design given their specific constraints of sequencing investment.