Multiple SCN5A variant enhancers modulate its cardiac gene expression and the QT interval

The rationale for genome-wide association study (GWAS) results is sequence variation in cis-regulatory elements (CREs) modulating a target gene's expression as the major cause of trait variation. To understand the complete molecular landscape of one of these GWAS loci, we performed in vitro reporter screens in cardiomyocyte cell lines for CREs overlapping nearly all common variants associated with any of five independent QT interval (QTi)-associated GWAS hits at the SCN5A-SCN10A locus. We identified 13 causal CRE variants using allelic reporter activity, cardiomyocyte nuclear extract-based binding assays, overlap with human cardiac tissue DNaseI hypersensitive regions, and predicted impact of sequence variants on DNaseI sensitivity. Our analyses identified at least one high-confidence causal CRE variant for each of the five sentinel hits that could collectively predict SCN5A cardiac gene expression and QTi association. Although all 13 variants could explain SCN5A gene expression, the highest statistical significance was obtained with seven variants (inclusive of the five above). Thus, multiple, causal, mutually associated CRE variants can underlie GWAS signals.

Association of breast cancer risk with genetic variants showing differential allelic expression: Identification of a novel breast cancer susceptibility locus at 4q21

There are significant inter-individual differences in the levels of gene expression. Through modulation of gene expression, cis-acting variants represent an important source of phenotypic variation. Consequently, cis-regulatory SNPs associated with differential allelic expression are functional candidates for further investigation as disease-causing variants. To investigate whether common variants associated with differential allelic expression were involved in breast cancer susceptibility, a list of genes was established on the basis of their involvement in cancer related pathways and/or mechanisms. Thereafter, using data from a genome-wide map of allelic expression associated SNPs, 313 genetic variants were selected and their association with breast cancer risk was then evaluated in 46,451 breast cancer cases and 42,599 controls of European ancestry ascertained from 41 studies participating in the Breast Cancer Association Consortium. The associations were evaluated with overall breast cancer risk and with estrogen receptor negative and positive disease. One novel breast cancer susceptibility locus on 4q21 (rs11099601) was identified (OR = 1.05, P = 5.6x10-6). rs11099601 lies in a 135 kb linkage disequilibrium block containing several genes, including, HELQ, encoding the protein HEL308 a DNA dependant ATPase and DNA Helicase involved in DNA repair, MRPS18C encoding the Mitochondrial Ribosomal Protein S18C and FAM175A (ABRAXAS), encoding a BRCA1 BRCT domain-interacting protein involved in DNA damage response and double-strand break (DSB) repair. Expression QTL analysis in breast cancer tissue showed rs11099601 to be associated with HELQ (P = 8.28x10-14), MRPS18C (P = 1.94x10-27) and FAM175A (P = 3.83x10-3), explaining about 20%, 14% and 1%, respectively of the variance inexpression of these genes in breast carcinomas.

Rare double sex and mab-3-related transcription factor 1 regulatory variants in severe spermatogenic failure

The double sex and mab-3-related transcription factor 1 (DMRT1) gene has long been linked to sex-determining pathways across vertebrates and is known to play an essential role in gonadal development and maintenance of spermatogenesis in mice. In humans, the genomic region harboring the DMRT gene cluster has been implicated in disorders of sex development and recently DMRT1 deletions were shown to be associated with non-obstructive azoospermia (NOA). In this work, we have employed different methods to screen a cohort of Portuguese NOA patients for DMRT1 exonic insertions and deletions [by multiplex ligation probe assay (MLPA); n = 68] and point mutations (by Sanger sequencing; n = 155). We have found three novel patient-specific non-coding variants in heterozygosity that were absent from 357 geographically matched controls. One of these is a complex variant with a putative regulatory role (c.-223_-219CGAAA>T), located in the promoter region within a conserved sequence involved in Dmrt1 repression. Moreover, while DMRT1 domains are highly conserved across vertebrates and show reduced levels of diversity in human populations, two rare synonymous substitutions (rs376518776 and rs34946058) and two rare non-coding variants that potentially affect DMRT1 expression and splicing (rs144122237 and rs200423545) were overrepresented in patients when compared with 376 Portuguese controls (301 fertile and 75 normozoospermic). Overall our previous and present results suggest a role of changes in DMRT1 dosage in NOA potentially also through a process of gene misregulation, even though DMRT1 deleterious variants seem to be rare.