Power of QTL detection using association tests with family controls

The power of testing for a population-wide association between a biallelic quantitative trait locus and a linked biallelic marker locus is predicted both empirically and deterministically for several tests. The tests were based on the analysis of variance (ANOVA) and on a number of transmission disequilibrium tests (TDT). Deterministic power predictions made use of family information, and were functions of population parameters including linkage disequilibrium, allele frequencies, and recombination rate. Deterministic power predictions were very close to the empirical power from simulations in all scenarios considered in this study. The different TDTs had very similar power, intermediate between one-way and nested ANOVAs. One-way ANOVA was the only test that was not robust against spurious disequilibrium. Our general framework for predicting power deterministically can be used to predict power in other association tests. Deterministic power calculations are a powerful tool for researchers to plan and evaluate experiments and obviate the need for elaborate simulation studies.

A description of the Mei2-like protein family; structure, phylogenetic distribution and biological context

The Schizosaccharomyces pombe Mei2 gene encodes an RNA recognition motif (RRM) protein that stimulates meiosis upon binding a specific non-coding RNA and subsequent accumulation in a “mei2-dot” in the nucleus. We present here the first systematic characterization of the family of proteins with characteristic Mei2-like amino acid sequences. Mei2-like proteins are an ancient eukaryotic protein family with three identifiable RRMs. The C-terminal RRM (RRM3) is unique to Mei2-like proteins and is the most highly conserved of the three RRMs. RRM3 also contains conserved sequence elements at its C-terminus not found in other RRM domains. Single copy Mei2-like genes are present in some fungi, in alveolates such as Paramecium and in the early branching eukaryote Entamoeba histolytica, while plants contain small families of Mei2-like genes. While the C-terminal RRM is highly conserved between plants and fungi, indicating conservation of molecular mechanisms, plant Mei2-like genes have changed biological context to regulate various aspects of developmental pattern formation.

An investigation of the C77G and C772T variations within the human protein tyrosine phosphatase receptor type C gene for association with multiple sclerosis in an Australian population

Multiple sclerosis (MS) is a common cause of neurological disability in young adults. The disease generally manifests in early to middle adulthood and causes various neurological deficits. Autoreactive T lymphocytes and their associated antigens have long been presumed important features of MS pathogenesis. The Protein tyrosine phosphatase receptor type C gene (PTPRC) encodes the T-cell receptor CD45. Variations within PTPRC have been previously associated with diseases of autoimmune origin such as type 1 diabetes mellitus and Graves' disease. We set out to investigate two variants within the PTPRC gene, C77G and C772T in subjects with MS and matched healthy controls to determine whether significant differences exist in these markers in an Australian population. We employed high resolution melt analysis (HRM) and restriction length polymorphism (RFLP) techniques to determine genotypic and allelic frequencies. Our study found no significant difference between frequencies for PTPRC C77G by either genotype (Χ2 = 0.65, P = 0.72) or allele (Χ2 = 0.48, P = 0.49). Similarly, we did not find evidence to suggest an association between PTPRC C772T by genotype (Χ2 = 1.06, P = 0.59) or allele (Χ2 = 0.20, P = 0.66). Linkage disequilibrium (LD) analysis showed strong linkage disequilibrium between the two tested markers (D' = 0.9970, SD = 0.0385). This study reveals no evidence to suggest that these markers are associated with MS in the tested Australian Caucasian population. Although the PTPRC gene has a significant role in regulating CD4+ and CD8+ autoreactive T-cells, interferon-beta responsiveness, and potentially other important processes, our study does not support a role for the two tested variants of this gene in MS susceptibility in the Australian population.

Proteolytic activation of Chlamydia trachomatis HTRA is mediated by PDZ1 domain interactions with protease domain loops L3 and LC and beta strand β5

Chlamydia trachomatis is a bacterial pathogen responsible for one of the most prevalent sexually transmitted infections worldwide. Its unique development cycle has limited our understanding of its pathogenic mechanisms. However, CtHtrA has recently been identified as a potential C. trachomatis virulence factor. CtHtrA is a tightly regulated quality control protein with a monomeric structural unit comprised of a chymotrypsin-like protease domain and two PDZ domains. Activation of proteolytic activity relies on the C-terminus of the substrate allosterically binding to the PDZ1 domain, which triggers subsequent conformational change and oligomerization of the protein into 24-mers enabling proteolysis. This activation is mediated by a cascade of precise structural arrangements, but the specific CtHtrA residues and structural elements required to facilitate activation are unknown. Using in vitro analysis guided by homology modeling, we show that the mutation of residues Arg362 and Arg224, predicted to disrupt the interaction between the CtHtrA PDZ1 domain and loop L3, and between loop L3 and loop LD, respectively, are critical for the activation of proteolytic activity. We also demonstrate that mutation to residues Arg299 and Lys160, predicted to disrupt PDZ1 domain interactions with protease loop LC and strand β5, are also able to influence proteolysis, implying their involvement in the CtHtrA mechanism of activation. This is the first investigation of protease loop LC and strand β5 with respect to their potential interactions with the PDZ1 domain. Given their high level of conservation in bacterial HtrA, these structural elements may be equally significant in the activation mechanism of DegP and other HtrA family members.

Association of the NuMA region on chromosome 11q13 with breast cancer susceptibility

The development of breast cancer is a complex process that involves multiple genes at many stages, from initial cell cycle dysregulation to disease progression. To identify genetic variations that influence this process, we conducted a large-scale association study using a collection of German cases and controls and >25,000 SNPs located within 16,000 genes. One of the loci identified was located on chromosome 11q13 [odds ratio (OR)=1.85, P=0.017]. The initial association was subsequently tested in two independent breast cancer collections. In both sample sets, the frequency of the susceptibility allele was increased in the cases (OR=1.6, P=0.01). The susceptibility allele was also associated with an increase in cancer family history (P=0.1). Fine mapping showed that the region of association extends approximately 300 kb and spans several genes, including the gene encoding the nuclear mitotic apparatus protein (NuMA). A nonsynonymous SNP (A794G) in NuMA was identified that showed a stronger association with breast cancer risk than the initial marker SNP (OR=2.8, P=0.005 initial sample; OR=2.1, P=0.002 combined). NuMA is a cell cycle-related protein essential for normal mitosis that is degraded in early apoptosis. NuMA-retinoic acid receptor alpha fusion proteins have been described in acute promyelocytic leukemia. Although the potential functional relevance of the A794G variation requires further biological validation, we conclude that variations in the NuMA gene are likely responsible for the observed increased breast cancer risk.

G-protein β3 subunit gene (GNB3) variant in causation of essential hypertension

Essential hypertensives display enhanced signal transduction through pertussis toxin-sensitive G proteins. The T allele of a C825T variant in exon 10 of the G protein β3 subunit gene (GNB3) induces formation of a splice variant (Gβ3-s) with enhanced activity. The T allele of GNB3 was shown recently to be associated with hypertension in unselected German patients (frequency=0.31 versus 0.25 in control). To confirm and extend this finding in a different setting, we performed an association study in Australian white hypertensives. This involved an extensively examined cohort of 110 hypertensives, each of whom were the offspring of 2 hypertensive parents, and 189 normotensives whose parents were both normotensive beyond age 50 years. Genotyping was performed by polymerase chain reaction and digestion with BseDI, which either cut (C allele) or did not cut (T allele) the 268-bp polymerase chain reaction product. T allele frequency in the hypertensive group was 0.43 compared with 0.25 in the normotensive group (χ2=22; P=0.00002; odds ratio=2.3; 95% CI=1.7 to 3.3). The T allele tracked with higher pretreatment blood pressure: diastolic=105±7, 109±16, and 128±28 mm Hg (mean±SD) for CC, CT, and 7T, respectively (P=0.001 by 1-way ANOVA). Blood pressures were higher in female hypertensives with a T allele (P=0.006 for systolic and 0.0003 for diastolic by ANOVA) than they were in male hypertensives. In conclusion, the present study of a group with strong family history supports a role for a genetically determined, physiologically active splice variant of the G protein β3 subunit gene in the causation of essential hypertension.

Ran is a potential therapeutic target for cancer cells with molecular changes associated with activation of the PI3K/Akt/mTORC1 and Ras/MEK/ERK pathways

Purpose Cancer cells have been shown to be more susceptible to Ran knockdown than normal cells. We now investigate whether Ran is a potential therapeutic target of cancers with frequently found mutations that lead to higher Ras/MEK/ERK [mitogen-activated protein/extracellular signal-regulated kinase (ERK; MEK)] and phosphoinositide 3-kinase (PI3K)/Akt/mTORC1 activities. Experimental Design Apoptosis was measured by flow cytometry [propidium iodide (PI) and Annexin V staining] and MTT assay in cancer cells grown under different conditions after knockdown of Ran. The correlations between Ran expression and patient survival were examined in breast and lung cancers. Results Cancer cells with their PI3K/Akt/mTORC1 and Ras/MEK/ERK pathways inhibited are less susceptible to Ran silencing-induced apoptosis. K-Ras-mutated, c-Met-amplified, and Pten-deleted cancer cells are also more susceptible to Ran silencing-induced apoptosis than their wild-type counterparts and this effect is reduced by inhibitors of the PI3K/Akt/mTORC1 and MEK/ERK pathways. Overexpression of Ran in clinical specimens is significantly associated with poor patient outcome in both breast and lung cancers. This association is dramatically enhanced in cancers with increased c-Met or osteopontin expression, or with oncogenic mutations of K-Ras or PIK3CA, all of which are mutations that potentially correlate with activation of the PI3K/Akt/mTORC1 and/or Ras/MEK/ERK pathways. Silencing Ran also results in dysregulation of nucleocytoplasmic transport of transcription factors and downregulation of Mcl-1 expression, at the transcriptional level, which are reversed by inhibitors of the PI3K/Akt/mTORC1 and MEK/ERK pathways. Conclusion Ran is a potential therapeutic target for treatment of cancers with mutations/changes of expression in protooncogenes that lead to activation of the PI3K/Akt/mTORC1 and Ras/MEK/ERK pathways. ©2011 AACR.

DRB2 is required for microRNA biogenesis in Arabidopsis thaliana

Background The Arabidopsis thaliana (Arabidopsis) DOUBLE-STRANDED RNA BINDING (DRB) protein family consists of five members, DRB1 to DRB5. The biogenesis of two developmentally important small RNA (sRNA) species, the microRNAs (miRNAs) and trans-acting small interfering RNAs (tasiRNAs) by DICER-LIKE (DCL) endonucleases requires the assistance of DRB1 and DRB4 respectively. The importance of miRNA-directed target gene expression in plant development is exemplified by the phenotypic consequence of loss of DRB1 activity (drb1 plants). Principal Findings Here we report that the developmental phenotype of the drb235 triple mutant plant is the result of deregulated miRNA biogenesis in the shoot apical meristem (SAM) region. The expression of DRB2, DRB3 and DRB5 in wild-type seedlings is restricted to the SAM region. Small RNA sequencing of the corresponding tissue of drb235 plants revealed altered miRNA accumulation. Approximately half of the miRNAs detected remained at levels equivalent to those of wild-type plants. However, the accumulation of the remaining miRNAs was either elevated or reduced in the triple mutant. Examination of different single and multiple drb mutants revealed a clear association between the loss of DRB2 activity and altered accumulation for both the elevated and reduced miRNA classes. Furthermore, we show that the constitutive over-expression of DRB2 outside of its wild-type expression domain can compensate for the loss of DRB1 activity in drb1 plants. Conclusions/Significance Our results suggest that in the SAM region, DRB2 is both antagonistic and synergistic to the role of DRB1 in miRNA biogenesis, adding an additional layer of gene regulatory complexity in this developmentally important tissue.

Gamma tocotrienol targets tyrosine phosphatase SHP2 in mammospheres resulting in cell death through RAS/ERK pathway

Background There is increasing evidence supporting the concept of cancer stem cells (CSCs), which are responsible for the initiation, growth and metastasis of tumors. CSCs are thus considered the target for future cancer therapies. To achieve this goal, identifying potential therapeutic targets for CSCs is essential. Methods We used a natural product of vitamin E, gamma tocotrienol (gamma-T3), to treat mammospheres and spheres from colon and cervical cancers. Western blotting and real-time RT-PCR were employed to identify the gene and protein targets of gamma-T3 in mammospheres. Results We found that mammosphere growth was inhibited in a dose dependent manner, with total inhibition at high doses. Gamma-T3 also inhibited sphere growth in two other human epithelial cancers, colon and cervix. Our results suggested that both Src homology 2 domain-containing phosphatase 1 (SHP1) and 2 (SHP2) were affected by gamma-T3 which was accompanied by a decrease in K- and H-Ras gene expression and phosphorylated ERK protein levels in a dose dependent way. In contrast, expression of self-renewal genes TGF-beta and LIF, as well as ESR signal pathways were not affected by the treatment. These results suggest that gamma-T3 specifically targets SHP2 and the RAS/ERK signaling pathway. Conclusions SHP1 and SHP2 are potential therapeutic targets for breast CSCs and gamma-T3 is a promising natural drug for future breast cancer therapy.

Association of STAT3 and TNFRSF1A with ankylosing spondylitis in Han Chinese

Objectives: Recent association studies by the Australo-Anglo-American Spondyloarthritis Consortium (TASC) in Caucasian European populations from Australia, North America and the UK have identified a number of genes as being associated with ankylosing spondylitis (AS). A candidate gene study in a Han Chinese population was performed based on these findings to identify associated genes in this population. Methods: A case-control study was performed in a Han Chinese population of patients with AS (n=775) and controls (n=1587) from Shanghai and Nanjing. All patients met the modified New York criteria for AS. The cases and controls were genotyped for 115 single nucleotide polymorphisms (SNPs) tagging IL23R, ERAP1, STAT3, JAK2, TNFRSF1A and TRADD, as well as other confirmation SNPs from the TASC study, using the Sequenom iPlex and the ABI OpenArray platforms. Statistical analysis of genotyped SNPs was performed using the Cochran - Armitage test for trend and meta-analysis was performed using METAL. SNPs in AS-associated genes in this study were then imputed using MaCH, and association with AS tested by logistic regression. Results: SNPs in TNFRSF1A (rs4149577, p=8.2×10-4), STAT3 (rs2293152, p=0.0015; rs1053005, p=0.017) and ERAP1 (rs27038, p=0.0091; rs27037, p=0.0092) were significantly associated with AS in Han Chinese. Association was also observed between AS and the intergenic region 2p15 (rs10865331, p=0.023). The lack of association between AS and IL23R in Han Chinese was confirmed (all SNPs p>0.1). Conclusions: The study results demonstrate for the first time that genetic polymorphisms in STAT3, TNFRSF1A and 2p15 are associated with AS in Han Chinese, suggesting common pathogenic mechanisms for the disease in Chinese and Caucasian European populations. Furthermore, previous findings demonstrating that ERAP1, but not IL23R, is associated with AS in Chinese patients were confirmed.

The chromosome 16q region associated with ankylosing spondylitis includes the candidate gene tumour necrosis factor receptor type 1-associated death domain (TRADD)

Objective: To replicate and refine the reported association of ankylosing spondylitis (AS) with two nonsynonymous single nucleotide polymorphisms (nsSNPs) on chromosome 16q22.1. Methods: Firstly, 730 independent UK patients with AS were genotyped for rs9939768 and rs6979 and allele frequencies were compared with 2879 previously typed historic disease controls. Secondly, the two data sets were combined in meta-analyses. Finally, 5 tagging SNPs, located between rs9939768 and rs6979, were analysed in 1604 cases and 1020 controls. Results: The association of rs6979 with AS was replicated, p=0.03, OR=1.14 (95% CI 1.01 to 1.28), and a trend for association with rs9939768 detected, p=0.06, OR=1.25 (95% CI 0.99 to 1.57). Meta-analyses revealed association of both SNPs with AS, p=0.0008, OR=1.31 (95% CI 1.12 to 1.54) and p=0.0009, OR=1.15 (95% CI 1.06 to 1.23) for rs9939768 and rs6979, respectively. New associations with rs9033 and rs868213 (p=0.00002, OR=1.23 (95% CI 1.12 to 1.36) and p=0.00002 OR=1.45 (95% CI 1.22 to 1.72), respectively, were identified. Conclusions: The region on chromosome 16 that has been replicated in the present work is interesting as the highly plausible candidate gene, tumour necrosis factor receptor type 1 (TNFR1)-associated death domain (TRADD), is located between rs9033 and rs868213. It will require additional work to identify the primary genetic association(s) with AS.

Association of sporadic chondrocalcinosis with a -4-basepair G-to-A transition in the 5′-untranslated region of ANKH that promotes enhanced expression of ANKH protein and excess generation of extracellular inorganic pyrophosphate

Objective Certain mutations in ANKH, which encodes a multiple-pass transmembrane protein that regulates inorganic pyrophosphate (PPi) transport, are linked to autosomal-dominant familial chondrocalcinosis. This study investigated the potential for ANKH sequence variants to promote sporadic chondrocalcinosis. Methods ANKH variants identified by genomic sequencing were screened for association with chondrocalcinosis in 128 patients with severe sporadic chondrocalcinosis or pseudogout and in ethnically matched healthy controls. The effects of specific variants on expression of common markers were evaluated by in vitro transcription/translation. The function of these variants was studied in transfected human immortalized CH-8 articular chondrocytes. Results Sporadic chondrocalcinosis was associated with a G-to-A transition in the ANKH 5′-untranslated region (5′-UTR) at 4 bp upstream of the start codon (in homozygotes of the minor allele, genotype relative risk 6.0, P = 0.0006; overall genotype association P = 0.02). This -4-bp transition, as well as 2 mutations previously linked with familial and sporadic chondrocalcinosis (+14 bp C-to-T and C-terminal GAG deletion, respectively), but not the French familial chondrocalcinosis kindred 143-bp T-to-C mutation, increased reticulocyte ANKH transcription/ANKH translation in vitro. Transfection of complementary DNA for both the wild-type ANKH and the -4-bp ANKH protein variant promoted increased extracellular PPi in CH-8 cells, but unexpectedly, these ANKH mutants had divergent effects on the expression of extracellular PPi and the chondrocyte hypertrophy marker, type X collagen. Conclusion A subset of sporadic chondrocalcinosis appears to be heritable via a -4-bp G-to-A ANKH 5′-UTR transition that up-regulates expression of ANKH and extracellular PPi in chondrocyte cells. Distinct ANKH mutations associated with heritable chondrocalcinosis may promote disease by divergent effects on extracellular PPi and chondrocyte hypertrophy, which is likely to mediate differences in the clinical phenotypes and severity of the disease.

Autosomal dominant spondylocostal dysostosis is caused by mutation in TBX6

In humans, congenital spinal defects occur with an incidence of 0.5-1 per 1000 live births. One of the most severe syndromes with such defects is spondylocostal dysostosis (SCD). Over the past decade, the genetic basis of several forms of autosomal recessive SCD cases has been solved with the identification of four causative genes (DLL3, MESP2, LFNG and HES7). Autosomal dominant forms of SCD have also been reported, but to date no genetic etiology has been described for these. Here, we have used exome capture and next-generation sequencing to identify a stoploss mutation in TBX6 that segregates with disease in two generations of one family. We show that this mutation has a deleterious effect on the transcriptional activation activity of the TBX6 protein, likely due to haploinsufficiency. In mouse, Tbx6 is essential for the patterning of the vertebral precursor tissues, somites; thus, mutation of TBX6 is likely to be causative of SCD in this family. This is the first identification of the genetic cause of an autosomal dominant form of SCD, and also demonstrates the potential of exome sequencing to identify genetic causes of dominant diseases even in small families with few affected individuals.

Breakthroughs in genetic studies of ankylosing spondylitis

Ankylosing spondylitis (AS), the prototypic seronegative arthropathy, is known to be highly heritable, with >90% of the risk of developing the disease determined genetically. As with most common heritable diseases, progress in identifying the genes involved using family-based or candidate gene approaches has been slow. The recent development of the genome-wide association study approach has revolutionized genetic studies of such diseases. Early studies in ankylosing spondylitis have produced two major breakthroughs in the identification of genes contributing roughly one third of the population attributable risk of the disease, and pointing directly to a potential therapy. These exciting findings highlight the potential of future more comprehensive genetic studies of determinants of disease risk and clinical manifestations, and are the biggest advance in our understanding of the causation of the disease since the discovery of the association with HLA-B27.