Vectors and methods for hairpin RNA and artificial microRNA-mediated gene silencing in plants

In plant cells, DICER-LIKE4 processes perfectly double-stranded RNA (dsRNA) into short interfering (si) RNAs, and DICER-LIKE1 generates micro (mi) RNAs from primary miRNA transcripts (pri-miRNA) that form fold-back structures of imperfectly dsRNA. Both si and miRNAs direct the endogenous endonuclease, ARGONAUTE1 to cleave complementary target single-stranded RNAs and either small RNA (sRNA)-directed pathway can be harnessed to silence genes in plants. A routine way of inducing and directing RNA silencing by siRNAs is to express self-complementary single-stranded hairpin RNA (hpRNA), in which the duplexed region has the same sequence as part of the target gene's mRNA. Artificial miRNA (amiRNA)-mediated silencing uses an endogenous pri-miRNA, in which the original miRNA/miRNA* sequence has been replaced with a sequence complementary to the new target gene. In this chapter, we describe the plasmid vector systems routinely used by our research group for the generation of either hpRNA-derived siRNAs or amiRNAs.

Efficient silencing of endogenous microRNAs using artificial microRNAs in Arabidopsis thaliana

We report here that the expression of endogenous microRNAs (miRNAs) can be efficiently silenced in Arabidopsis thaliana (Arabidopsis) using artificial miRNA (amiRNA) technology. We demonstrate that an amiRNA designed to target a mature miRNA directs silencing against all miRNA family members, whereas an amiRNA designed to target the stem-loop region of a miRNA precursor transcript directs silencing against only the individual family member targeted. Furthermore, our results indicate that amiRNAs targeting both the mature miRNA and stem-loop sequence direct RNA silencing through cleavage of the miRNA precursor transcript, which presumably occurs in the nucleus of a plant cell during the initial stages of miRNA biogenesis. This suggests that small RNA (sRNA)-guided RNA cleavage in plants occurs not only in the cytoplasm, but also in the nucleus. Many plant miRNA gene families have been identified via sequencing and bioinformatic analysis, but, to date, only a small tranche of these have been functionally characterized due to a lack of effective forward or reverse genetic tools. Our findings therefore provide a new and powerful reverse-genetic tool for the analysis of miRNA function in plants. © The Author 2010. Published by the Molecular Plant Shanghai Editorial Office in association with Oxford University Press on behalf of CSPP and IPPE, SIBS, CAS.

Rice ragged stunt oryzavirus genome segments S7 and S10 encode non-structural proteins of M(r) 68 025 (Pns7) and M(r) 32364 (Pns10) : brief report

The nucleotide sequences of genome segments S7 and S10 of a Thai-isolate of rice ragged stunt virus (RRSV) were determined. The 1938 bp S7 sequence contains a single large open reading frame (ORF) spanning nucleotides 20 to 1 843 that is predicted to encode a protein of M(r) 68 025. The 1 162 bp S10 sequence has a major ORF spanning nucleotides 142 to 1 032 that is predicted to encode a protein of M(r) 32364. This S10 ORF is preceded by a small ORF (nt 20-55) which is probably a minicistron. Coupled in vitro transcription-translation from the two major ORFs gave protein products of the expected sizes. However, no protein was visualised from S10 when the small ORF sequence was included. Proteins were expressed in Escherichia coli from the full length ORF of S7 (P7) and from a segment of the S10 ORF (P10) fused to the ORF of glutathione S-transferase (GST). Neither fusion protein was recognised by polyclonal antibodies raised against RRSV particles. Furthermore, polyclonal antibodies raised against GST-P7 fusion protein did not recognise any virion structural polypeptides. These data strongly suggest that the proteins P7 and P10 do not form part of RRSV particle. This is further supported by observed sequence homology (though very weak) of predicted.

The M(r) 43K major capsid protein of rice ragged stunt oryzavirus is a post-translationally processed product of a M(r) 67 348 polypeptide encoded by genome segment 8

The nucleotide sequence of DNA complementary to rice ragged stunt oryzavirus (RRSV) genome segment 8 (S8) of an isolate from Thailand was determined. RRSV S8 is 1 914 bp in size and contains a single large open reading frame (ORF) spanning nucleotides 23 to 1 810 which is capable of encoding a protein of M(r) 67 348. The N-terminal amino acid sequence of a ~43K virion polypeptide matched to that inferred for an internal region of the S8 coding sequence. These data suggest that the 43K protein is encoded by S8 and is derived by a proteolytic cleavage. Predicted polypeptide sizes from this possible cleavage of S8 protein are 26K and 42K. Polyclonal antibodies raised against a maltose binding protein (MBP)-S8 fusion polypeptide (expressed in Escherichia coli) recognised four RRSV particle associated polypeptides of M(r) 67K, 46K, 43K and 26K and all except the 26K polypeptide were also highly immunoreactive to polyclonal antibodies raised against purified RRSV particles. Cleavage of the MBP-S8 fusion polypeptide with protease Factor X produced the expected 40K MBP and two polypeptides of apparent M(r) 46K and 26K. Antibodies to purified RRSV particles reacted strongly with the intact fusion protein and the 46K cleavage product but weakly to the 26K product. Furthermore, in vitro transcription and translation of the S8 coding region revealed a post-translational self cleavage of the 67K polypeptide to 46K and 26K products. These data indicate that S8 encodes a structural polypeptide, the majority of which is auto- catalytically cleaved to 26K and 46K proteins. The data also suggest that the 26K protein is the self cleaving protease and that the 46K product is further processed or undergoes stable conformational changes to a ~43K major capsid protein.

Sequence and identification of the barley yellow dwarf virus coat protein gene

The nucleotide sequence of the coat protein gene of barley yellow dwarf virus (BYDV, PAV serotype) was determined, and the amino acid sequence was deduced. The open reading frame, encoding a protein of relative molecular mass (Mr) 22,047, was confirmed as the coat protein gene by comparison with amino acid sequences of tryptic peptides derived from dissociated virions. In addition, a fragment of this gene expressed in Escherichia coli produced a product which was recognized by antibodies prepared against purified BYDV virions. An overlapping reading frame encoding an Mr 17,147 protein is contained completely within the coat protein gene. © 1988.

Functional role for senataxin, defective in ataxia oculomotor apraxia type 2, in transcriptional regulation

Ataxia oculomotor apraxia type 2 (AOA2) is an autosomal recessive neurodegenerative disorder characterized by cerebellar ataxia and oculomotor apraxia. The gene mutated in AOA2, SETX, encodes senataxin, a putative DNA/RNA helicase which shares high homology to the yeast Sen1p protein and has been shown to play a role in the response to oxidative stress. To investigate further the function of senataxin, we identified novel senataxin-interacting proteins, the majority of which are involved in transcription and RNA processing, including RNA polymerase II. Binding of RNA polymerase II to candidate genes was significantly reduced in senataxin deficient cells and this was accompanied by decreased transcription of these genes, suggesting a role for senataxin in the regulation/modulation of transcription. RNA polymerase II-dependent transcription termination was defective in cells depleted of senataxin in keeping with the observed interaction of senataxin with poly(A) binding proteins 1 and 2. Splicing efficiency of specific mRNAs and alternate splice-site selection of both endogenous genes and artificial minigenes were altered in senataxin depleted cells. These data suggest that senataxin, similar to its yeast homolog Sen1p, plays a role in coordinating transcriptional events, in addition to its role in DNA repair.

Disrupted sensory gating in pathological gambling

Background Some neurochemical evidence as well as recent studies on molecular genetics suggest that pathologic gambling may be related to dysregulated dopamine neurotransmission. Methods The current study examined sensory (motor) gating in pathologic gamblers as a putative measure of endogenous brain dopamine activity with prepulse inhibition of the acoustic startle eye-blink response and the auditory P300 event-related potential. Seventeen pathologic gamblers and 21 age- and gender-matched healthy control subjects were assessed. Both prepulse inhibition measures were recorded under passive listening and two-tone prepulse discrimination conditions. Results Compared to the control group, pathologic gamblers exhibited disrupted sensory (motor) gating on all measures of prepulse inhibition. Sensory motor gating deficits of eye-blink responses were most profound at 120-millisecond prepulse lead intervals in the passive listening task and at 240-millisecond prepulse lead intervals in the two-tone prepulse discrimination task. Sensory gating of P300 was also impaired in pathologic gamblers, particularly at 500-millisecond lead intervals, when performing the discrimination task on the prepulse. Conclusions In the context of preclinical studies on the disruptive effects of dopamine agonists on prepulse inhibition, our findings suggest increased endogenous brain dopamine activity in pathologic gambling in line with previous neurobiological findings.

Fine-mapping of the HNF1B multicancer locus identifies candidate variants that mediate endometrial cancer risk

Common variants in the hepatocyte nuclear factor 1 homeobox B (HNF1B) gene are associated with the risk of Type II diabetes and multiple cancers. Evidence to date indicates that cancer risk may be mediated via genetic or epigenetic effects on HNF1B gene expression. We previously found single-nucleotide polymorphisms (SNPs) at the HNF1B locus to be associated with endometrial cancer, and now report extensive fine-mapping and in silico and laboratory analyses of this locus. Analysis of 1184 genotyped and imputed SNPs in 6608 Caucasian cases and 37 925 controls, and 895 Asian cases and 1968 controls, revealed the best signal of association for SNP rs11263763 (P = 8.4 × 10−14, odds ratio = 0.86, 95% confidence interval = 0.82–0.89), located within HNF1B intron 1. Haplotype analysis and conditional analyses provide no evidence of further independent endometrial cancer risk variants at this locus. SNP rs11263763 genotype was associated with HNF1B mRNA expression but not with HNF1B methylation in endometrial tumor samples from The Cancer Genome Atlas. Genetic analyses prioritized rs11263763 and four other SNPs in high-to-moderate linkage disequilibrium as the most likely causal SNPs. Three of these SNPs map to the extended HNF1B promoter based on chromatin marks extending from the minimal promoter region. Reporter assays demonstrated that this extended region reduces activity in combination with the minimal HNF1B promoter, and that the minor alleles of rs11263763 or rs8064454 are associated with decreased HNF1B promoter activity. Our findings provide evidence for a single signal associated with endometrial cancer risk at the HNF1B locus, and that risk is likely mediated via altered HNF1B gene expression.

Autosomal dominant familial calcium pyrophosphate dihydrate deposition disease is caused by mutation in the transmembrane protein ANKH

Familial autosomal dominant calcium pyrophosphate dihydrate (CPPD) chondrocalcinosis has previously been mapped to chromosome 5pl5. We have identified a mutation in the ANKH gene that segregates with the disease in a family with this condition. ANKH encodes a putative transmembrane inorganic pyrophosphate (PPi) transport channel. We postulate that loss of function of ANKH causes elevated extracellular PPi levels, predisposing to CPPD crystal deposition.

The dominant 55kDa allergen of the subtropical Bahia grass (Paspalum notatum) pollen is a group 13 pollen allergen, Pas n 13

Bahia grass, Paspalum notatum, is an important pollen allergen source with a long season of pollination and wide distribution in subtropical and temperate regions. We aimed to characterize the 55. kDa allergen of Bahia grass pollen (BaGP) and ascertain its clinical importance. BaGP extract was separated by 2D-PAGE and immunoblotted with serum IgE of a grass pollen-allergic patient. The amino-terminal protein sequence of the predominant allergen isoform at 55. kDa had similarity with the group 13 allergens of Timothy grass and maize pollen, Phl p 13 and Zea m 13. Four sequences obtained by rapid amplification of the allergen cDNA ends represented multiple isoforms of Pas n 13. The predicted full length cDNA for Pas n 13 encoded a 423 amino acid glycoprotein including a signal peptide of 28 residues and with a predicted pI of 7.0. Tandem mass spectrometry of tryptic peptides of 2D gel spots identified peptides specific to the deduced amino acid sequence for each of the four Pas n 13 cDNA, representing 47% of the predicted mature protein sequence of Pas n 13. There was 80.6% and 72.6% amino acid identity with Zea m 13 and Phl p 13, respectively. Reactivity with a Phl p 13-specific monoclonal antibody AF6 supported designation of this allergen as Pas n 13. The allergen was purified from BaGP extract by ammonium sulphate precipitation, hydrophobic interaction and size exclusion chromatography. Purified Pas n 13 reacted with serum IgE of 34 of 71 (48%) grass pollen-allergic patients and specifically inhibited IgE reactivity with the 55. kDa band of BaGP for two grass pollen-allergic donors. Four isoforms of Pas n 13 from pI 6.3-7.8 had IgE-reactivity with grass pollen allergic sera. The allergenic activity of purified Pas n 13 was demonstrated by activation of basophils from whole blood of three grass pollen-allergic donors tested but not control donors. Pas n 13 is thus a clinically relevant pollen allergen of the subtropical Bahia grass likely to be important in eliciting seasonal allergic rhinitis and asthma in grass pollen-allergic patients.

Genetic determinants of heel bone properties: Genome-wide association meta-analysis and replication in the GEFOS/GENOMOS consortium

Quantitative ultrasound of the heel captures heel bone properties that independently predict fracture risk and, with bone mineral density (BMD) assessed by X-ray (DXA), may be convenient alternatives for evaluating osteoporosis and fracture risk. We performed a meta-analysis of genome-wide association (GWA) studies to assess the genetic determinants of heel broadband ultrasound attenuation (BUA; n 5 14 260), velocity of sound (VOS; n 5 15 514) and BMD (n 5 4566) in 13 discovery cohorts. Independent replication involved seven cohorts with GWA data (in silico n 5 11 452) and new genotyping in 15 cohorts (de novo n 5 24 902). In combined random effects, meta-analysis of the discovery and replication cohorts, nine single nucleotide polymorphisms (SNPs) had genome-wide significant (P < 5 3 108) associations with heel bone properties. Alongside SNPs within or near previously identified osteoporosis susceptibility genes including ESR1 (6q25.1: rs4869739, rs3020331, rs2982552), SPTBN1 (2p16.2: rs11898505), RSPO3 (6q22.33: rs7741021), WNT16 (7q31.31: rs2908007), DKK1 (10q21.1: rs7902708) and GPATCH1 (19q13.11: rs10416265), we identified a new locus on chromosome 11q14.2 (rs597319 close to TMEM135, a gene recently linked to osteoblastogenesis and longevity) significantly associated with both BUA and VOS (P < 8.23 3 1014). In meta-analyses involving 25 cohorts with up to 14 985 fracture cases, six of 10 SNPs associated with heel bone properties at P < 5 3 106 also had the expected direction of association with any fracture (P < 0.05), including threeSNPswithP < 0.005: 6q22.33 (rs7741021), 7q31.31 (rs2908007) and 10q21.1 (rs7902708). In conclusion, thisGWAstudy reveals the effect of several genescommon to central DXA-derivedBMDand heel ultrasound/DXAmeasures and points to anewgenetic locus with potential implications for better understanding of osteoporosis pathophysiology.

Founder p.Arg 446* mutation in the PDHX gene explains over half of cases with congenital lactic acidosis in Roma children

Investigation of 31 of Roma patients with congenital lactic acidosis (CLA) from Bulgaria identified homozygosity for the R446* mutation in the PDHX gene as the most common cause of the disorder in this ethnic group. It accounted for around 60% of patients in the study and over 25% of all CLA cases referred to the National Genetic Laboratory in Bulgaria. The detection of a homozygous patient from Hungary and carriers among population controls from Romania and Slovakia suggests a wide spread of the mutation in the European Roma population. The clinical phenotype of the twenty R446* homozygotes was relatively homogeneous, with lactic acidosis crisis in the first days or months of life as the most common initial presentation (15/20 patients) and delayed psychomotor development and/or seizures in infancy as the leading manifestations in a smaller group (5/20 patients). The subsequent clinical picture was dominated by impaired physical growth and a very consistent pattern of static cerebral palsy-like encephalopathy with spasticity and severe to profound mental retardation seen in over 80% of cases. Most patients had a positive family history. We propose testing for the R446* mutation in PDHX as a rapid first screening in Roma infants with metabolic acidosis. It will facilitate and accelerate diagnosis in a large proportion of cases, allow early rehabilitation to alleviate the chronic clinical course, and prevent further affected births in high-risk families.

Genome-wide association study of intraocular pressure identifies the GLCCI1/ICA1 region as a glaucoma susceptibility locus

To discover quantitative trait loci for intraocular pressure, amajor risk factor for glaucoma and the only modifiable one,weperformed agenome-wide association studyonadiscoverycohort of2175individualsfromSydney, Australia. We found a novel association between intraocular pressure and a common variant at 7p21 near to GLCCI1 and ICA1. The findings in this region were confirmed through two UK replication cohorts totalling 4866 individuals (rs59072263, Pcombined = 1.10 × 10-8). A copy of the G allele at this SNP is associated with an increase in mean IOP of 0.45 mmHg (95%CI = 0.30-0.61 mmHg). These results lend support to the implication of vesicle trafficking and glucocorticoid inducibility pathways in the determination of intraocular pressure and in the pathogenesis of primary open-angle glaucoma.

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.

Using genome-Wwide complex trait analysis to quantify 'missing heritability' in Parkinson's disease

Genome-wide association studies (GWASs) have been successful at identifying single-nucleotide polymorphisms (SNPs) highly associated with common traits; however, a great deal of the heritable variation associated with common traits remains unaccounted for within the genome. Genome-wide complex trait analysis (GCTA) is a statistical method that applies a linear mixed model to estimate phenotypic variance of complex traits explained by genome-wide SNPs, including those not associated with the trait in a GWAS. We applied GCTA to 8 cohorts containing 7096 case and 19 455 control individuals of European ancestry in order to examine the missing heritability present in Parkinson's disease (PD). We meta-analyzed our initial results to produce robust heritability estimates for PD types across cohorts. Our results identify 27% (95% CI 17-38, P = 8.08E - 08) phenotypic variance associated with all types of PD, 15% (95% CI -0.2 to 33, P = 0.09) phenotypic variance associated with early-onset PD and 31% (95% CI 17-44, P = 1.34E - 05) phenotypic variance associated with late-onset PD. This is a substantial increase from the genetic variance identified by top GWAS hits alone (between 3 and 5%) and indicates there are substantially more risk loci to be identified. Our results suggest that although GWASs are a useful tool in identifying the most common variants associated with complex disease, a great deal of common variants of small effect remain to be discovered. © Published by Oxford University Press 2012.