The Association Study between Twenty One Polymorphisms in Seven Candidate Genes and Coronary Heart Diseases in Chinese Han Population

Previous genome-wide association studies (GWAS) in multiple populations identified several genetic loci for coronary heart diseases (CHD). Here we utilized a 2-stage candidate gene association strategy in Chinese Han population to shed light on the putative association between several metabolic-related candidate genes and CHD. At the 1^st stage, 190 patients with CHD and 190 controls were genotyped through the MassARRAY platform. At the 2^nd stage, a larger sample including 400 patients and 392 controls was genotyped by the High Resolution Melt (HRM) method to confirm or rule out the associations with CHD. MLXIP expression level was quantified by the real time PCR in 65 peripheral blood samples. From the 21 studied single nucleotide polymorphisms (SNPs) of seven candidate genes: MLXIPL, MLXIP, MLX, ADIPOR1, VDR, SREBF1 and NR1H3, only one tag SNP rs4758685 (T→C) was found to be statistically associated with CHD (P-value = 0.02, Odds ratio (OR) of 0.83). After adjustment for the age, sex, lipid levels and diabetes, the association remained significant (P-value = 0.03). After adjustment for the hypertension, P-value became 0.20 although there was a significant difference in the allele distribution between the CHD patients with hypertension and the controls (P-value = 0.04, 406 vs 582). In conclusion, among the 21 tested SNPs, we identified a novel association between rs4758685 of MLXIP gene and CHD. The C allele of common variant rs4758685 interacted with hypertension, and was found to be protective against CHD in both allelic and genotypic models in Chinese Han population.     

Exome sequencing in multiplex families with left-sided cardiac defects has high yield for disease gene discovery

Congenital heart disease (CHD) is a common group of birth defects with a strong genetic contribution to their etiology, but historically the diagnostic yield from exome studies of isolated CHD has been low. Pleiotropy, variable expressivity, and the difficulty of accurately phenotyping newborns contribute to this problem. We hypothesized that performing exome sequencing on selected individuals in families with multiple members affected by left-sided CHD, then filtering variants by population frequency, in silico predictive algorithms, and phenotypic annotations from publicly available databases would increase this yield and generate a list of candidate disease-causing variants that would show a high validation rate. In eight of the nineteen families in our study (42%), we established a well-known gene/phenotype link for a candidate variant or performed confirmation of a candidate variant’s effect on protein function, including variants in genes not previously described or firmly established as disease genes in the body of CHD literature: BMP10, CASZ1, ROCK1 and SMYD1. Two plausible variants in different genes were found to segregate in the same family in two instances suggesting oligogenic inheritance. These results highlight the need for functional validation and demonstrate that in the era of next-generation sequencing, multiplex families with isolated CHD can still bring high yield to the discovery of novel disease genes.     

Haploinsufficiency of TAB2 Causes Congenital Heart Defects in Humans

Congenital heart defects (CHDs) are the most common major developmental anomalies and the most frequent cause for perinatal mortality, but their etiology remains often obscure. We identified a locus for CHDs on 6q24-q25. Genotype-phenotype correlations in 12 patients carrying a chromosomal deletion on 6q delineated a critical 850 kb region on 6q25.1 harboring five genes. Bioinformatics prioritization of candidate genes in this locus for a role in CHDs identified the TGF-β-activated kinase 1/MAP3K7 binding protein 2 gene (TAB2) as the top-ranking candidate gene. A role for this candidate gene in cardiac development was further supported by its conserved expression in the developing human and zebrafish heart. Moreover, a critical, dosage-sensitive role during development was demonstrated by the cardiac defects observed upon titrated knockdown of tab2 expression in zebrafish embryos. To definitively confirm the role of this candidate gene in CHDs, we performed mutation analysis of TAB2 in 402 patients with a CHD, which revealed two evolutionarily conserved missense mutations. Finally, a balanced translocation was identified, cosegregating with familial CHD. Mapping of the breakpoints demonstrated that this translocation disrupts TAB2. Taken together, these data clearly demonstrate a role for TAB2 in human cardiac development.     

The contribution of de novo and rare inherited copy number changes to congenital heart disease in an unselected sample of children with conotruncal defects or hypoplastic left heart disease

Congenital heart disease (CHD) is the most common congenital malformation, with evidence of a strong genetic component. We analyzed data from 223 consecutively ascertained families, each consisting of at least one child affected by a conotruncal defect (CNT) or hypoplastic left heart disease (HLHS) and both parents. The NimbleGen HD2-2.1 comparative genomic hybridization platform was used to identify de novo and rare inherited copy number variants (CNVs). Excluding 10 cases with 22q11.2 DiGeorge deletions, we validated de novo CNVs in 8 % of 148 probands with CNTs, 12.7 % of 71 probands with HLHS and none in 4 probands with both. Only 2 % of control families showed a de novo CNV. We also identified a group of ultra-rare inherited CNVs that occurred de novo in our sample, contained a candidate gene for CHD, recurred in our sample or were present in an affected sibling. We confirmed the contribution to CHD of copy number changes in genes such as GATA4 and NODAL and identified several genes in novel recurrent CNVs that may point to novel CHD candidate loci. We also found CNVs previously associated with highly variable phenotypes and reduced penetrance, such as dup 1q21.1, dup 16p13.11, dup 15q11.2-13, dup 22q11.2, and del 2q23.1. We found that the presence of extra-cardiac anomalies was not related to the frequency of CNVs, and that there was no significant difference in CNV frequency or specificity between the probands with CNT and HLHS. In agreement with other series, we identified likely causal CNVs in 5.6 % of our total sample, half of which were de novo.     

Identification of embryonic lethal genes in humans by autozygosity mapping and exome sequencing in consanguineous families

BackgroundIdentifying genetic variants that lead to discernible phenotypes is the core of Mendelian genetics. An approach that considers embryonic lethality as a bona fide Mendelian phenotype has the potential to reveal novel genetic causes, which will further our understanding of early human development at a molecular level. Consanguineous families in which embryonic lethality segregates as a recessive Mendelian phenotype offer a unique opportunity for high throughput novel gene discovery as has been established for other recessive postnatal phenotypes.ResultsWe have studied 24 eligible families using autozygosity mapping and whole-exome sequencing. In addition to revealing mutations in genes previously linked to embryonic lethality in severe cases, our approach revealed seven novel candidate genes (THSD1, PIGC, UBN1, MYOM1, DNAH14, GALNT14, and FZD6). A founder mutation in one of these genes, THSD1, which has been linked to vascular permeability, accounted for embryonic lethality in three of the study families. Unlike the other six candidate genes, we were able to identify a second mutation in THSD1 in a family with a less severe phenotype consisting of hydrops fetalis and persistent postnatal edema, which provides further support for the proposed link between this gene and embryonic lethality.ConclusionsOur study represents an important step towards the systematic analysis of “embryonic lethal genes” in humans.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0681-6) contains supplementary material, which is available to authorized users.     

An Integrative Analysis to Identify Driver Genes in Esophageal Squamous Cell Carcinoma

BackgroundFew driver genes have been well established in esophageal squamous cell carcinoma (ESCC). Identification of the genomic aberrations that contribute to changes in gene expression profiles can be used to predict driver genes.MethodsWe searched for driver genes in ESCC by integrative analysis of gene expression microarray profiles and copy number data. To narrow down candidate genes, we performed survival analysis on expression data and tested the genetic vulnerability of each genes using public RNAi screening data. We confirmed the results by performing RNAi experiments and evaluating the clinical relevance of candidate genes in an independent ESCC cohort.ResultsWe found 10 significantly recurrent copy number alterations accompanying gene expression changes, including loci 11q13.2, 7p11.2, 3q26.33, and 17q12, which harbored CCND1, EGFR, SOX2, and ERBB2, respectively. Analysis of survival data and RNAi screening data suggested that GRB7, located on 17q12, was a driver gene in ESCC. In ESCC cell lines harboring 17q12 amplification, knockdown of GRB7 reduced the proliferation, migration, and invasion capacities of cells. Moreover, siRNA targeting GRB7 had a synergistic inhibitory effect when combined with trastuzumab, an anti-ERBB2 antibody. Survival analysis of the independent cohort also showed that high GRB7 expression was associated with poor prognosis in ESCC.ConclusionOur integrative analysis provided important insights into ESCC pathogenesis. We identified GRB7 as a novel ESCC driver gene and potential new therapeutic target.     

Temporal and tissue specific gene expression patterns of the zebrafish kinesin-1 heavy chain family,kif5s,during development

Homo- and heterodimers of Kif5 proteins form the motor domain of Kinesin-1, a major plus-end directed microtubule motor. Kif5s have been implicated in the intracellular transport of organelles, vesicles, proteins, and RNAs in many cell types. There are three mammalian KIF5s. KIF5A and KIF5C proteins are strictly neural in mouse whereas, KIF5B is ubiquitously expressed. Mouse knockouts indicate crucial roles for KIF5 in development and human mutations in KIF5A lead to the neurodegenerative disease Hereditary Spastic Paraplegia. However, the developmental functions and the extent to which individual kif5 functions overlap have not been elucidated. Zebrafish possess five kif5 genes: kif5Aa, kif5Ab, kif5Ba, kif5Bb, and kif5C. Here we report their tissue specific expression patterns in embryonic and larval stages. Specifically, we find that kif5As are strictly zygotic and exhibit neural-specific expression. In contrast, kif5Bs exhibit strong maternal contribution and are ubiquitously expressed. Lastly, kif5C exhibits weak maternal expression followed by enrichment in neural populations. In addition, kif5s show distinct expression domains in the larval retina.     

The Landscape of Candidate Driver Genes Differs between Male and Female Breast Cancer

The rapidly growing collection of diverse genome-scale data from multiple tumor types sheds light on various aspects of the underlying tumor biology. With the objective to identify genes of importance for breast tumorigenesis in men and to enable comparisons with genes important for breast cancer development in women, we applied the computational framework COpy Number and EXpression In Cancer (CONEXIC) to detect candidate driver genes among all altered passenger genes. Unique to this approach is that each driver gene is associated with several gene modules that are believed to be altered by the driver. Thirty candidate drivers were found in the male breast cancers and 67 in the female breast cancers. We identified many known drivers of breast cancer and other types of cancer, in the female dataset (e.g. GATA3, CCNE1, GRB7, CDK4). In contrast, only three known cancer genes were found among male breast cancers; MAP2K4, LHP, and ZNF217. Many of the candidate drivers identified are known to be involved in processes associated with tumorigenesis, including proliferation, invasion and differentiation. One of the modules identified in male breast cancer was regulated by THY1, a gene involved in invasion and related to epithelial-mesenchymal transition. Furthermore, men with THY1 positive breast cancers had significantly inferior survival. THY1 may thus be a promising novel prognostic marker for male breast cancer. Another module identified among male breast cancers, regulated by SPAG5, was closely associated with proliferation. Our data indicate that male and female breast cancers display highly different landscapes of candidate driver genes, as only a few genes were found in common between the two. Consequently, the pathobiology of male breast cancer may differ from that of female breast cancer and can be associated with differences in prognosis; men diagnosed with breast cancer may consequently require different management and treatment strategies than women.     

Identification of rare genetic variants in novel loci associated with Paget’s disease of bone

In genome-wide association studies, single nucleotide polymorphisms located in five novel loci were associated with PDB. We aimed at identifying rare genetic variants of candidate genes located in these loci and search for genetic association with PDB in the French-Canadian population. Exons, promoter and exon–intron junctions from patients with familial PDB and healthy individuals were sequenced in candidate genes, located within novel loci associated with PDB in our population. Rare variant was defined by a minor allele frequency <0.05 or absent from dbSNP (NCBI). We sequenced seven genes in 1p13 locus, three genes in 7q33, three genes in 8q22, and five genes in 15q24 locus. We identified 126 rare variants in at least one patient with PDB of whom 55 were located in 1p13 locus, 32 in 7q33, 10 in 8q22 and 29 in 15q24 locus. We located 71 of these 126 rare variants in an intron, 30 in an exon and 9 in an untranslated region. 60 % of these variants were located in functionally relevant gene regions. Among the 12 missense rare variants in PDB, two (rs62620995 in TM7SF4; rs62641691 in CD276) were predicted to be damaging by in silico analysis tools. Rs62620995, which altered a conserved amino acid (p.Leu397Phe) in the TM7SF4 gene, encoding the DC-STAMP protein involved in osteoclastogenesis through RANK signaling pathway, was found to have a marginal association with PDB (p = 0.09). Rs35500845, located in the CTHRC1 gene, which encodes a regulator of collagen matrix deposition, was also associated with PDB in the French-Canadian population (p = 0.046).     

Chromosome microarray testing for patients with congenital heart defects reveals novel disease causing loci and high diagnostic yield

Background

Congenital heart defects (CHD), as the most common congenital anomaly, have been reported to be frequently associated with pathogenic copy number variants (CNVs). Currently, patients with CHD are routinely offered chromosomal microarray (CMA) testing, but the diagnostic yield of CMA on CHD patients has not been extensively evaluated based on a large patient cohort. In this study, we retrospectively assessed the detected CNVs in a total of 514 CHD cases (a 422-case clinical cohort from Boston Children''s Hospital (BCH) and a 92-case research cohort from Shanghai Children’s Medical Center (SCMC)) and conducted a genotype-phenotype analysis. Furthermore, genes encompassed in pathogenic/likely pathogenic CNVs were prioritized by integrating several tools and public data sources for novel CHD candidate gene identification.

Results

Based on the BCH cohort, the overall diagnostic yield of CMA testing for CHD patients was 12.8(pathogenic CNVs)-18.5% (pathogenic and likely pathogenic CNVs). The diagnostic yield of CMA for syndromic CHD was 14.1-20.6% (excluding aneuploidy cases), whereas the diagnostic yield for isolated CHD was 4.3-9.3%. Four recurrent genomic loci (4q terminal region, 15q11.2, 16p12.2 and Yp11.2) were more significantly enriched in cases than in controls. These regions are considered as novel CHD loci. We further identified 20 genes as the most likely novel CHD candidate genes through gene prioritization analysis.

Conclusion

The high clinical diagnostic yield of CMA in this study provides supportive evidence for CMA as the first-line genetic diagnostic tool for CHD patients. The CNVs detected in our study suggest a number of CHD candidate genes that warrant further investigation.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1127) contains supplementary material, which is available to authorized users.     

Molecular basis of albinism in India: Evaluation of seven potential candidate genes and some new findings

Albinism represents a group of genetic disorders with a broad spectrum of hypopigmentary phenotypes dependent on the genetic background of the patients. Oculocutaneous albinism (OCA) patients have little or no pigment in their eyes, skin and hair, whereas ocular albinism (OA) primarily presents the ocular symptoms, and the skin and hair color may vary from near normal to very fair. Mutations in genes directly or indirectly regulating melanin production are responsible for different forms of albinism with overlapping clinical features. In this study, 27 albinistic individuals from 24 families were screened for causal variants by a PCR-sequencing based approach. TYR, OCA2, TYRP1, SLC45A2, SLC24A5, TYRP2 and SILV were selected as candidate genes. We identified 5 TYR and 3 OCA2 mutations, majority in homozygous state, in 8 unrelated patients including a case of autosomal recessive ocular albinism (AROA). A homozygous 4-nucleotide novel insertion in SLC24A5 was detected in a person showing with extreme cutaneous hypopigmentation. A potential causal variant was identified in the TYRP2 gene in a single patient. Haplotype analyses in the patients carrying homozygous mutations in the classical OCA genes suggested founder effect. This is the first report of an Indian AROA patient harboring a mutation in OCA2. Our results also reveal for the first time that mutations in SLC24A5 could contribute to extreme hypopigmentation in humans.     

Eleven Candidate Susceptibility Genes for Common Familial Colorectal Cancer

Hereditary factors are presumed to play a role in one third of colorectal cancer (CRC) cases. However, in the majority of familial CRC cases the genetic basis of predisposition remains unexplained. This is particularly true for families with few affected individuals. To identify susceptibility genes for this common phenotype, we examined familial cases derived from a consecutive series of 1514 Finnish CRC patients. Ninety-six familial CRC patients with no previous diagnosis of a hereditary CRC syndrome were included in the analysis. Eighty-six patients had one affected first-degree relative, and ten patients had two or more. Exome sequencing was utilized to search for genes harboring putative loss-of-function variants, because such alterations are likely candidates for disease-causing mutations. Eleven genes with rare truncating variants in two or three familial CRC cases were identified: UACA, SFXN4, TWSG1, PSPH, NUDT7, ZNF490, PRSS37, CCDC18, PRADC1, MRPL3, and AKR1C4. Loss of heterozygosity was examined in all respective cancer samples, and was detected in seven occasions involving four of the candidate genes. In all seven occasions the wild-type allele was lost (P = 0.0078) providing additional evidence that these eleven genes are likely to include true culprits. The study provides a set of candidate predisposition genes which may explain a subset of common familial CRC. Additional genetic validation in other populations is required to provide firm evidence for causality, as well as to characterize the natural history of the respective phenotypes.     

Genetic Variation in the NOC Gene Is Associated with Body Mass Index in Chinese Subjects

Circadian clock genes are critical regulators of energy homeostasis and metabolism. However, whether variation in the circadian genes is associated with metabolic phenotypes in humans remains to be explored. In this study, we systemically genotyped 20 tag single nucleotide polymorphisms (SNPs) in 8 candidate genes involved in circadian clock, including CLOCK, BMAL1(ARNTL), PER1, PER2, CRY1, CRY2, CSNK1E,, and NOC(CCRN4L) in 1,510 non-diabetic Chinese subjects in Taipei and Yunlin populations in Taiwan. Their associations with metabolic phenotypes were analyzed. We found that genetic variation in the NOC gene, rs9684900 was associated with body mass index (BMI) (P = 0.0016, Bonferroni corrected P = 0.032). Another variant, rs135764 in the CSNK1E gene was associated with fasting glucose (P = 0.0023, Bonferroni corrected P = 0.046). These associations were consistent in both Taipei and Yunlin populations. Significant epistatic and joint effects between SNPs on BMI and related phenotypes were observed. Furthermore, NOC mRNA levels in human abdominal adipose tissue were significantly increased in obese subjects compared to non-obese controls.

Conclusion

Genetic variation in the NOC gene is associated with BMI in Chinese subjects.     

A Systems Genetics Approach Identifies CXCL14, ITGAX,and LPCAT2 as Novel Aggressive Prostate Cancer Susceptibility Genes

Although prostate cancer typically runs an indolent course, a subset of men develop aggressive, fatal forms of this disease. We hypothesize that germline variation modulates susceptibility to aggressive prostate cancer. The goal of this work is to identify susceptibility genes using the C57BL/6-Tg(TRAMP)8247Ng/J (TRAMP) mouse model of neuroendocrine prostate cancer. Quantitative trait locus (QTL) mapping was performed in transgene-positive (TRAMPxNOD/ShiLtJ) F2 intercross males (n = 228), which facilitated identification of 11 loci associated with aggressive disease development. Microarray data derived from 126 (TRAMPxNOD/ShiLtJ) F2 primary tumors were used to prioritize candidate genes within QTLs, with candidate genes deemed as being high priority when possessing both high levels of expression-trait correlation and a proximal expression QTL. This process enabled the identification of 35 aggressive prostate tumorigenesis candidate genes. The role of these genes in aggressive forms of human prostate cancer was investigated using two concurrent approaches. First, logistic regression analysis in two human prostate gene expression datasets revealed that expression levels of five genes (CXCL14, ITGAX, LPCAT2, RNASEH2A, and ZNF322) were positively correlated with aggressive prostate cancer and two genes (CCL19 and HIST1H1A) were protective for aggressive prostate cancer. Higher than average levels of expression of the five genes that were positively correlated with aggressive disease were consistently associated with patient outcome in both human prostate cancer tumor gene expression datasets. Second, three of these five genes (CXCL14, ITGAX, and LPCAT2) harbored polymorphisms associated with aggressive disease development in a human GWAS cohort consisting of 1,172 prostate cancer patients. This study is the first example of using a systems genetics approach to successfully identify novel susceptibility genes for aggressive prostate cancer. Such approaches will facilitate the identification of novel germline factors driving aggressive disease susceptibility and allow for new insights into these deadly forms of prostate cancer.     

Screening of indigenous goats for prolificacy associated DNA markers of sheep

The present study was undertaken to explore the genetic basis of caprine prolificacy and to screen indigenous goats for prolificacy associated markers of sheep in BMPR1B, GDF9 and BMP15 genes. To detect the associated mutations and identify novel allelic variants in the candidate genes, representative samples were collected from the breeding tract of indigenous goat breeds varying in prolificacy and geographic distribution. DNA was extracted and PCR amplification was done using primers designed or available in literature for the coding DNA sequence of candidate genes. Direct sequencing was done to identify the genetic variations. Mutations in the candidate genes associated with fecundity in sheep were not detected in Indian goats. Three non-synonymous SNPs (C818T, A959C and G1189A) were identified in exon 2 of GDF9 gene out of which mutation A959C has been associated with prolificacy in exotic goats. Two novel SNPs (G735A and C808G) were observed in exon 2 of BMP15 gene.     

12 Survival-related differentially expressed genes based on the TARGET-osteosarcoma database

The Therapeutically Applicable Research to Generate Effective Treatments (TARGET) project aims to determine molecular changes that drive childhood cancers, including osteosarcoma. The main purpose of the program is to use the open-source database to develop novel, effective, and less toxic therapies. We downloaded TARGET-OS RNA-Sequencing data through R studio and merged the mRNA expression of genes with clinical information (vital status, survival time and gender). Further, we analyzed differential gene expressions between dead and alive patients based on TARGET-OS project. By this study, we found 5758 differentially expressed genes between deceased and alive patients with a false discovery rate below 0.05; 4469 genes were upregulated in deceased patients compared to alive, whereas 1289 genes were downregulated. The survival-related genes were obtained using Kaplan–Meier survival analysis and Cox univariate regression (KM < 0.05 and Cox P-value < 0.05). Out of 5758 differentially expressed genes, only 217 have been associated with overall survival. Eight survival-related downregulated genes (ERCC4, CLUAP1, CTNNBIP1, GCA, RAB40C, SIRPA, USP11, and TCN2) and four survival-related upregulated genes (MUC1, COL13A1, JAG2 and KAZALD1) were selected for further analysis as potential independent prognostic candidate genes. This study may help to discover novel prognostic markers and potential therapeutic targets for osteosarcoma.     

Selection of biomarkers for HIV-1 latency by integrated analysis

A major obstacle in the treatment of human immunodeficiency virus type 1 (HIV-1) is its ability to establish latent infection. To find novel biomarkers associated with the mechanism of HIV-1 latent infection, we identified 70 candidate genes in HIV-1 latently infected cells through the integrated analysis in a previous study. It is important to select more effective biomarkers among 70 candidates and to verify the possibility of selected biomarkers for HIV-1 latency. We identified the 24 and 25 genes from 70 candidate genes in significantly enriched categories selected by Database for Annotation, Visualization and Integrated Discovery (DAVID) software and Gene Set Enrichment Analysis (GSEA) software, respectively. Also, we investigated genes regulated in both HIV-1 latently infected cell lines and PBMCs from HIV-1 infected patients and found the genes with a common pattern of expression levels in both cell lines and PBMCs. Consequently, we identified nine genes, APBB2, GMPR, IGF2BP3, LRP1, MAD2L2, MX1, OXR1, PTK2B, and TNFSF13B, via integrated analysis. Especially, APBB2 and MAD2L2 were identified in both DAVID and GSEA software. Our findings suggest that nine genes were identified via integrated analysis as potential biomarkers and in particular, APBB2 and MAD2L2 may be considered as more significant biomarkers for HIV-1 latency.