Resolving candidate genes of mouse skeletal muscle QTL via RNA-Seq and expression network analyses

BACKGROUND:

We have recently identified a number of Quantitative Trait Loci (QTL) contributing to the 2-fold muscle weight difference between the LG/J and SM/J mouse strains and refined their confidence intervals. To facilitate nomination of the candidate genes responsible for these differences we examined the transcriptome of the tibialis anterior (TA) muscle of each strain by RNA-Seq.

13,726 genes were expressed in mouse skeletal muscle. Intersection of a set of 1061 differentially expressed transcripts with a mouse muscle Bayesian Network identified a coherent set of differentially expressed genes that we term the LG/J and SM/J Regulatory Network (LSRN). The integration of the QTL, transcriptome and the network analyses identified eight key drivers of the LSRN (Kdr, Plbd1, Mgp, Fah, Prss23, 2310014F06Rik, Grtp1, Stk10) residing within five QTL regions, which were either polymorphic or differentially expressed between the two strains and are strong candidates for quantitative trait genes (QTGs) underlying muscle mass. The insight gained from network analysis including the ability to make testable predictions is illustrated by annotating the LSRN with knowledge-based signatures and showing that the SM/J state of the network corresponds to a more oxidative state. We validated this prediction by NADH tetrazolium reductase staining in the TA muscle revealing higher oxidative potential of the SM/J compared to the LG/J strain (p<0.03).

Thus, integration of fine resolution QTL mapping, RNA-Seq transcriptome information and mouse muscle Bayesian Network analysis provides a novel and unbiased strategy for nomination of muscle QTGs.

Association Study of 167 Candidate Genes for Schizophrenia Selected by a Multi-Domain Evidence-Based Prioritization Algorithm and Neurodevelopmental Hypothesis

Integrating evidence from multiple domains is useful in prioritizing disease candidate genes for subsequent testing. We ranked all known human genes (n = 3819) under linkage peaks in the Irish Study of High-Density Schizophrenia Families using three different evidence domains: 1) a meta-analysis of microarray gene expression results using the Stanley Brain collection, 2) a schizophrenia protein-protein interaction network, and 3) a systematic literature search. Each gene was assigned a domain-specific p-value and ranked after evaluating the evidence within each domain. For comparison to this
ranking process, a large-scale candidate gene hypothesis was also tested by including genes with Gene Ontology terms related to neurodevelopment. Subsequently, genotypes of 3725 SNPs in 167 genes from a custom Illumina iSelect array were used to evaluate the top ranked vs. hypothesis selected genes. Seventy-three genes were both highly ranked and involved in neurodevelopment (category 1) while 42 and 52 genes were exclusive to neurodevelopment (category 2) or highly ranked (category 3), respectively. The most significant associations were observed in genes PRKG1, PRKCE, and CNTN4 but no individual SNPs were significant after correction for multiple testing. Comparison of the approaches showed an excess of significant tests using the hypothesis-driven neurodevelopment category. Random selection of similar sized genes from two independent genome-wide association studies (GWAS) of schizophrenia showed the excess was unlikely by chance. In a further meta-analysis of three GWAS datasets, four candidate SNPs reached nominal significance. Although gene ranking using integrated sources of prior information did not enrich for significant results in the current experiment, gene selection using an a priori hypothesis (neurodevelopment) was superior to random selection. As such, further development of gene ranking strategies using more carefully selected sources of information is warranted.

Identification of potentially pathogenic variants in candidate breast and ovarian cancer susceptibility genes in BRCAX patients

Mutations within the BRCA1 and BRCA2 genes account for approximately 20% of hereditary breast cancers, with a further 10%–15% being attributable to rare mutations in moderate-risk genes and common variants in low-risk genes. The genes harbouring mutations in the remaining ∼65% of hereditary breast cancers are unknown. The identification of mutation carriers in hereditary breast and ovarian cancer (hboc) families is critical for determining who is most at risk of developing the disease and therefore who should be offered risk-reducing procedures or more intensive screening, or both.

Many of the high- and moderate-risk genes for hereditary breast cancers encode proteins that work in concert to maintain genomic stability and in dna damage signalling and repair. A novel BRCA1 protein complex identified within the research group whose target genes are involved in dna repair provided novel candidates for hboc susceptibility genes. These 12 candidate genes were sequenced in a cohort of 675 affected individuals from the Kathleen Cunningham Foundation Consortium for Research into Familial Breast Cancer (kConFab) with hereditary breast or ovarian cancer, but with no mutations in known susceptibility genes (BRCAx patients). This analysis identified 20 individuals (each from a different BRCAx family) with different potentially pathogenic variants across 6 of the candidate hboc susceptibility genes. The family members of each BRCAx index case were tested for the presence of the specific mutation identified in the proband to examine segregation with disease. To further expand on the potential role of the novel candidate hboc susceptibility genes identified in this study, the genetic variation of a second cohort of 520 Northern Irish BRCAx patients is being characterized using a 61-gene panel.

Can Certain Genotypes Predispose to Poor Asthma Control in Children? A Pharmacogenetic Study of 9 Candiate Genes in Children with Difficult Asthma

Objective: We tested the hypothesis that patients with difficult asthma have an increased frequency of certain genotypes that predispose them to asthma exacerbations and poor asthma control.

Methods: A total of 180 Caucasian children with confirmed asthma diagnosis were selected from two phenotypic groups; difficult (n = 112) versus mild/moderate asthma (n = 68) groups. All patients were screened for 19 polymorphisms in 9 candidate genes to evaluate their association with difficult asthma.

Key Results: The results indicated that LTA4H A-9188.G, TNFa G-308.A and IL-4Ra A1727.G polymorphisms were significantly associated with the development of difficult asthma in paediatric patients (p,0.001, p = 0.019 and p = 0.037, respectively). Haplotype analysis also revealed two haplotypes (ATA haplotype of IL-4Ra A1199.C, IL-4Ra T1570.C and IL- 4Ra A1727.G and CA haplotype of TNFa C-863.A and TNFa G-308.A polymorphisms) which were significantly associated with difficult asthma in children (p = 0.04 and p = 0.018, respectively).

Conclusions and Clinical Relevance: The study revealed multiple SNPs and haplotypes in LTA4H, TNFa and IL4-Ra genes which constitute risk factors for the development of difficult asthma in children. Of particular interest is the LTA4H A- 9188.G polymorphism which has been reported, for the first time, to have strong association with severe asthma in children. Our results suggest that screening for patients with this genetic marker could help characterise the heterogeneity of responses to leukotriene-modifying medications and, hence, facilitate targeting these therapies to the subset of patients who are most likely to gain benefit. ©2013 Almomani et al.

Integration of global SNP-based mapping and expression arrays reveals key regions, mechanisms, and genes important in the pathogenesis of multiple myeloma.

Multiple myeloma is characterized by genomic alterations frequently involving gains and losses of chromosomes. Single nucleotide polymorphism (SNP)-based mapping arrays allow the identification of copy number changes at the sub-megabase level and the identification of loss of heterozygosity (LOH) due to monosomy and uniparental disomy (UPD). We have found that SNP-based mapping array data and fluorescence in situ hybridization (FISH) copy number data correlated well, making the technique robust as a tool to investigate myeloma genomics. The most frequently identified alterations are located at 1p, 1q, 6q, 8p, 13, and 16q. LOH is found in these large regions and also in smaller regions throughout the genome with a median size of 1 Mb. We have identified that UPD is prevalent in myeloma and occurs through a number of mechanisms including mitotic nondisjunction and mitotic recombination. For the first time in myeloma, integration of mapping and expression data has allowed us to reduce the complexity of standard gene expression data and identify candidate genes important in both the transition from normal to monoclonal gammopathy of unknown significance (MGUS) to myeloma and in different subgroups within myeloma. We have documented these genes, providing a focus for further studies to identify and characterize those that are key in the pathogenesis of myeloma.

Sequence variation in the CHAT locus shows no association with late-onset Alzheimer's disease

There is substantial evidence for a susceptibility gene for late-onset Alzheimer's disease (AD) on chromosome 10. One of the characteristic features of AD is the degeneration and dysfunction of the cholinergic system. The genes encoding choline acetyltransferase (ChAT) and its vesicular transporter (VAChT), CHAT and SLC18A3 respectively, map to the linked region of chromosome 10 and are therefore both positional and obvious functional candidate genes for late-onset AD. We have screened both genes for sequence variants and investigated each for association with late-onset AD in up to 500 late-onset AD cases and 500 control DNAs collected in the UK. We detected a total of 17 sequence variants. Of these, 14 were in CHAT, comprising three non-synonymous variants (D7N in the S exon, A120T in exon 5 and L243F in exon 8), one synonymous change (H547H), nine single-nucleotide polymorphisms in intronic, untranslated or promoter regions, and a variable number of tandem repeats in intron 7. Three non-coding SNPs were detected in SLC18A3. None demonstrated any reproducible association with late-onset AD in our samples. Levels of linkage disequilibrium were generally low across the CHAT locus but two of the coding variants, D7N and A120T, proved to be in complete linkage disequilibrium.

CRYBB1 mutation associated with congenital cataract and microcornea

The molecular characterization of a UK family with an autosomal dominant congenital cataract associated with microcornea is reported. METHODS: Family history and clinical data were recorded. This phenotype was linked to a 7.6 cM region of chromosome 22q11.2-q12.2, spanning the beta-crystallin gene cluster (ZMax of 3.91 for marker D22S1114 at theta=0). Candidate genes were PCR amplified and screened for mutations on both strands using direct sequencing. RESULTS: Sequencing of the coding regions and flanking intronic sequences of CRYBB2 and CRYBB1 showed the presence of a novel, heterozygous X253R change in exon 6 of CRYBB1. SSCP analysis confirmed that this sequence change segregated with the disease phenotype in all available family members and was not found in 109 ethnically matched controls. CONCLUSIONS: X253R is predicted to elongate the COOH-terminal extension of the protein and would be expected to disrupt beta-crystallin interactions. This is the first documented involvement of CRYBB1 in ocular development beyond cataractogenesis.

Genetic disorders and spermatogenesis.

Male infertility affects one man in twenty and a genetic basis seems likely in at least 30% of those men. Genetic regulation of fertility involves the inter-related processes of testicular development, spermatogenesis (involving germ cell mitosis, meiosis and spermatid maturation), and their endocrine and paracrine regulation. In regard to spermatogenesis, particular attention has been given to the Yq11 region, where some spermatogenesis genes ('azoospermia factors') appear to be located. Several candidate genes have been identified but have not been shown to have a defined or essential role in spermatogenesis. Microdeletions of Yq11 are found in approximately 15% of azoospermic or severely oligospermic men. The complexity of the genetic control of male fertility is demonstrated by the evidence for genes involved in spermatogenesis and sexual differentiation on the X chromosome and autosomes. Better understanding of the genetic regulation of normal spermatogenesis will provide new probes for clinical studies; however, at present the majority of spermatogenic failure remains without an identified genetic linkage. The advent of intracytoplasmic sperm injection permits fertility in many previously sterile men and presents the possibility of their transmission of infertility; appropriate counselling is required.

Investigation of the association of BMP gene variants with nephropathy in Type 1 diabetes mellitus

Aims Diabetic nephropathy is a leading cause of end-stage renal disease. The transforming growth factor beta-bone morphogenic protein (BMP) pathway is implicated in the pathogenesis of diabetic nephropathy. The BMP2, BMP4 and BMP7 genes are located near linkage peaks for renal dysfunction, and we hypothesize that genetic polymorphisms in these biological and positional candidate genes may be risk factors for diabetic kidney disease.

Targeted genome-wide investigation identifies novel SNP's associated with diabetic nephropathy

Loci contributing to complex disease have been identified by focusing on genome-wide scans utilising non-synonymous single nucleotide polymorphisms (nsSNPs). We employed Illumina’s HNS12 BeadChip (13,917 high-value SNPs) which was specifically designed to capture nsSNPs and ideally complements more dense genome-wide association studies that fail to consider many of these putatively functional variants. The HNS12 panel also includes 870 tag SNPs covering the major histocompatibility region. All individuals genotyped in this study were Caucasians with (cases) and without (controls) diabetic nephropathy. About 449 individuals with type 2 diabetes (203 cases, 246 controls) were genotyped in the initial study. 1,467 individuals with type 1 diabetes (718 cases, 749 controls) were genotyped in the follow up study. 11,152 SNPs were successfully analysed and ranked for association with diabetic nephropathy based on significance (P) values. The top ranked 32 SNPs were subsequently genotyped using MassARRAY iPLEX™ and TaqMan technologies to investigate association of these polymorphisms with nephropathy in individuals with type 1 diabetes. The top ranked nsSNP, rs1543547 (P = 10-5), is located in RAET1L, a major histocompatibility class I-related gene at 6q25.1. Of particular interest, multiple nsSNPs within the top ranked (0.2%) SNPs are within several plausible candidate genes for nephropathy on 3q21.3 and 6p21.3.