miRPlant : an integrated tool for identification of plant miRNA from RNA sequencing data

Background Small RNA sequencing is commonly used to identify novel miRNAs and to determine their expression levels in plants. There are several miRNA identification tools for animals such as miRDeep, miRDeep2 and miRDeep*. miRDeep-P was developed to identify plant miRNA using miRDeep’s probabilistic model of miRNA biogenesis, but it depends on several third party tools and lacks a user-friendly interface. The objective of our miRPlant program is to predict novel plant miRNA, while providing a user-friendly interface with improved accuracy of prediction. Result We have developed a user-friendly plant miRNA prediction tool called miRPlant. We show using 16 plant miRNA datasets from four different plant species that miRPlant has at least a 10% improvement in accuracy compared to miRDeep-P, which is the most popular plant miRNA prediction tool. Furthermore, miRPlant uses a Graphical User Interface for data input and output, and identified miRNA are shown with all RNAseq reads in a hairpin diagram. Conclusions We have developed miRPlant which extends miRDeep* to various plant species by adopting suitable strategies to identify hairpin excision regions and hairpin structure filtering for plants. miRPlant does not require any third party tools such as mapping or RNA secondary structure prediction tools. miRPlant is also the first plant miRNA prediction tool that dynamically plots miRNA hairpin structure with small reads for identified novel miRNAs. This feature will enable biologists to visualize novel pre-miRNA structure and the location of small RNA reads relative to the hairpin. Moreover, miRPlant can be easily used by biologists with limited bioinformatics skills.

Global dissemination of a multidrug resistant Escherichia coli clone

Escherichia coli sequence type 131 (ST131) is a globally disseminated, multidrug resistant (MDR) clone responsible for a high proportion of urinary tract and bloodstream infections. The rapid emergence and successful spread of E. coli ST131 is strongly associated with several factors, including resistance to fluoroquinolones, high virulence gene content, the possession of the type 1 fimbriae FimH30 allele, and the production of the CTX-M-15 extended spectrum β-lactamase (ESBL). Here, we used genome sequencing to examine the molecular epidemiology of a collection of E. coli ST131 strains isolated from six distinct geographical locations across the world spanning 2000–2011. The global phylogeny of E. coli ST131, determined from whole-genome sequence data, revealed a single lineage of E. coli ST131 distinct from other extraintestinal E. coli strains within the B2 phylogroup. Three closely related E. coli ST131 sublineages were identified, with little association to geographic origin. The majority of single-nucleotide variants associated with each of the sublineages were due to recombination in regions adjacent to mobile genetic elements (MGEs). The most prevalent sublineage of ST131 strains was characterized by fluoroquinolone resistance, and a distinct virulence factor and MGE profile. Four different variants of the CTX-M ESBL–resistance gene were identified in our ST131 strains, with acquisition of CTX-M-15 representing a defining feature of a discrete but geographically dispersed ST131 sublineage. This study confirms the global dispersal of a single E. coli ST131 clone and demonstrates the role of MGEs and recombination in the evolution of this important MDR pathogen.

Association of CYP1B1 codon 432 mutant allele in head and neck squamous cell cancer is reflected by somatic mutations of p53 in tumor tissue

Tobacco use is causally associated with head and neck squamous cell cancer (HNSCC). Here, we present the results of a case-control study that investigated the effects that the genetic variants of the cytochrome (CYP)1A1, CYP1B1, glutathione-S-transferase (GST)M1, GSTT1, and GSTP1 genes have on modifying the risk of smoking-related HNSCC. Allelisms of the CYP1A1, GSTT1, GSTM1, and GSTT1 genes alone were not associated with an increased risk. CYP1B1 codon 432 polymorphism was found to be a putative susceptibility factor in smoking-related HNSCC. The frequency of CYP1B1 polymorphism was significantly higher (P < 0.001) in the group of smoking cases when compared with smoking controls. Additionally, an odds ratio (OR) of 4.53 (2.62-7.98) was discovered when investigating smoking and nonsmoking cases for the susceptible genotype CYP1B1*2/*2, when compared with the presence of the genotype wild type. In combination with polymorphic variants of the GST genes, a synergistic-effect OR was observed. The calculated OR for the combined genotype CYP1B1*2/*2 and GSTM1*2/*2 was 12.8 (4.09-49.7). The calculated OR for the combined genotype was 13.4 (2.92-97.7) for CYP1B1*2/*2 and GSTT1*2/*2, and 24.1 (9.36-70.5) for the combination of CYP1B1*2/*2 and GSTT1-expressors. The impact of the polymorphic variants of the CYP1B1 gene on HNSCC risk is reflected by the strong association with the frequency of somatic mutations of the p53 gene. Smokers with susceptible genotype CYP1B1*2/*2 were 20 times more likely to show evidence of p53 mutations than were those with CYP1B1 wild type. Combined genotype analysis of CYP1B1 and GSTM1 or GSTT1 revealed interactive effects on the occurrence of p53 gene mutations. The results of the present study indicate that polymorphic variants of CYP1B1 relate significantly to the individual susceptibility of smokers to HNSCC.

Influence of polymorphisms of the human glutathione transferases and cytochrome P450 2E1 enzyme on the metabolism and toxicity of ethylene oxide and acrylonitrile

A cohort of 59 persons with industrial handling of low levels of acrylonitrile is being studied as part of a medical surveillance programme. Previously, an extended haemoglobin adduct monitoring (N-(cyanoethyl)valine and N-(hydroxyethyl)-valine) was performed regarding the glutathione transferases hGSTM1 and hGSTT1 polymorphisms but no influence of hGSTM1 or hGSTT1 polymorphisms on specific adduct levels was found. A compilation of case reports of human accidental poisonings had pointed to significant individual differences in human acrylonitrile metabolism and toxicity. Therefore, a re-evaluation of the industrial cohort included known polymorphisms of the glutathione transferases hGSTM3 and hGSTP1 as well as of the cytochrome P450 CYP2E1. A detailed statistical analysis revealed that exposed carriers of the allelic variants of hGSTP1, hGSTP1*B/hGSTP1*C, characterized by a single nucleotide polymorphism at nucleotide 313 which results in a change from Ile to Val at codon 104, had higher levels of the acrylonitrile-specific haemoglobin adduct N-(cyanoethyl)valine compared to the carriers of the codon 113 alleles hGSTP1*A and hGSTP1*D. The single nucleotide polymorphism at codon 113 of hGSTP1 (hGSTP1*A/hGSTP1*B versus hGSTP1*C/hGSTP1*D) did not show an effect, and also no influence was seen on specific haemoglobin adduct levels of the polymorphisms of hGSTM3 or CYP2E1. The data, therefore, point to a possible influence of a human enzyme polymorphism of the GSTP1 gene at codon 104 on the detoxication of acrylonitrile which calls for experimental toxicological investigation. The study also confirmed the impact of GSTT1 polymorphism on background N-(hydroxyethyl)-valine adduct levels in haemoglobin which are caused by endogenous ethylene oxide.

Cytochrome p450 1b1, a new keystone in gene-environment interactions related to human head and neck cancer?

Alcohol consumption and tobacco smoking are major causes of head and neck cancers, and regional differences point to the importance of research into gene-environment interactions. Much interest has been focused on polymorphisms of CYP1A1 and of GSTM1 and GSTT1, but a number of studies have not demonstrated significant effects. This has mostly been ascribed to small sample sizes. In general, the impact of polymorphisms of metabolic enzymes appears inconsistent, with some reports of weak-to-moderate associations, and with others of no elevation of risks. The classical cytochrome P450 isoenzyme considered for metabolic activation of polycyclic aromatic hydrocarbons (PAH) is CYP1A1. A new member of the CYP1 family, CYP1B1, was cloned in 1994, currently representing the only member of the CYP1B subfamily. A number of single nucleotide polymorphisms of the CYP1B1 gene have been reported. The amino acid substitutions Val432Leu (CYP1B1*3) and Asn453Ser (CYP1B1*4), located in the heme binding domain of CYP1B1, appear as likely candidates to be linked with biological effects. CYP1B1 activates a wide range of PAH, aromatic and heterocyclic amines. Very recently, the CYP1B1 codon 432 polymorphism (CYP1B1*3) has been identified as a susceptibility factor in smoking-related head-and-neck squamous cell cancer. The impact of this polymorphic variant of CYP1B1 on cancer risk was also reflected by an association with the frequency of somatic mutations of the p53 gene. Combined genotype analysis of CYP1B1 and the glutathione transferases GSTM1 or GSTT1 has pointed to interactive effects. This provides new molecular evidence that tobacco smoke-specific compounds relevant to head and neck carcinogenesis are metabolically activated through CYP1B1 and is consistent with a major pathogenetic relevance of PAH as ingredients of tobacco smoke.

Evidence of associations between cytokine gene polymorphisms and quality of life in patients with cancer and their family caregivers

PURPOSE/OBJECTIVES: To identify latent classes of individuals with distinct quality-of-life (QOL) trajectories, to evaluate for differences in demographic characteristics between the latent classes, and to evaluate for variations in pro- and anti-inflammatory cytokine genes between the latent classes. DESIGN: Descriptive, longitudinal study. SETTING: Two radiation therapy departments located in a comprehensive cancer center and a community-based oncology program in northern California. SAMPLE: 168 outpatients with prostate, breast, brain, or lung cancer and 85 of their family caregivers (FCs). METHODS: Growth mixture modeling (GMM) was employed to identify latent classes of individuals based on QOL scores measured prior to, during, and for four months following completion of radiation therapy. Single nucleotide polymorphisms (SNPs) and haplotypes in 16 candidate cytokine genes were tested between the latent classes. Logistic regression was used to evaluate the relationships among genotypic and phenotypic characteristics and QOL GMM group membership. MAIN RESEARCH VARIABLES: QOL latent class membership and variations in cytokine genes. FINDINGS: Two latent QOL classes were found: higher and lower. Patients and FCs who were younger, identified with an ethnic minority group, had poorer functional status, or had children living at home were more likely to belong to the lower QOL class. After controlling for significant covariates, between-group differences were found in SNPs in interleukin 1 receptor 2 (IL1R2) and nuclear factor kappa beta 2 (NFKB2). For IL1R2, carrying one or two doses of the rare C allele was associated with decreased odds of belonging to the lower QOL class. For NFKB2, carriers with two doses of the rare G allele were more likely to belong to the lower QOL class. CONCLUSIONS: Unique genetic markers in cytokine genes may partially explain interindividual variability in QOL. IMPLICATIONS FOR NURSING: Determination of high-risk characteristics and unique genetic markers would allow for earlier identification of patients with cancer and FCs at higher risk for poorer QOL. Knowledge of these risk factors could assist in the development of more targeted clinical or supportive care interventions for those identified.

Subclinical psychotic experiences in healthy young adults : associations with stress and genetic predisposition

Stress has been identified as a common trigger for psychosis. Dopamine pathways are suggested to be affected by chronic and severe stress and to play an important role in psychosis. This pilot study investigates the potential relationship of stress and psychosis in subclinical psychotic experiences. It was hypothesized that single-nucleotide polymorphisms (SNPs) previously found to be associated with psychiatric disorders would be associated with both stress and subclinical psychotic experiences. University students (N=182) were genotyped for 17 SNPs across 11 genes. Higher stress reporting was associated with rs4680 COMT, rs13211507 HLA region, and rs13107325 SLC39A8. Reports of higher subclinical psychotic experiences were associated with DRD2 SNPs rs17601612 and rs658986 and an AKT1 SNP rs2494732. Replication studies are recommended to further pursue this line of research for identification of markers of psychosis for early diagnosis and intervention.

Association study of MTHFD1 coding polymorphisms R134K and R653Q with migraine susceptibility

Objective There is evidence that folate metabolism has a role in migraine pathophysiology, particularly in the migraine with aura subtype. In this study we investigate whether two non-synonymous single nucleotide polymorphisms (SNPs), rs1950902 (C401T; R134K) and rs2236225 (G1958A; R653Q), in MTHDF1 are associated with migraine in an Australian case-control population. Background Increased plasma levels of homocysteine (HCy), one of the metabolites produced in the folate pathway, has been found to be a risk factor for migraine. There is also a genetic link, as a common polymorphism (C667T) that reduces the catalytic activity of MTHFR, the enzyme that catalyses the formation of HCy, is associated with an increase in risk of the migraine with aura (MA) subtype. MTHFD1 is a crucial multifunctional enzyme that catalyses three separate reactions of the folate pathway and therefore variants in MTHFD1 may also influence migraine susceptibility. Methods The R134K and R653Q variants in MTHFD1 were genotyped in an Australian cohort of 520 unrelated migraineurs (162 were diagnosed with migraine without aura [MO] and 358 with MA) and 520 matched controls. Data were analysed for association with migraine and for interaction with the MTHFR C667T polymorphism. Results We find no significant differences in genotype or allele frequencies for either SNP between migraineurs and controls, or when either MO or MA cases were compared to controls. In addition these MTHFD1 polymorphisms did not appear to influence the risk of MA conferred by the MTHFR 667T allele. Conclusions We find no evidence for association of the MTHFD1 R134K and R653Q polymorphisms with migraine in our Australian case-control population. However, as folate metabolism appears to be important in migraine, particularly with respect to the aura component, future studies using high throughput methods to expand the number of SNPs in folate-related genes genotyped and investigation of interactions between SNPs may be justified.

Environmental, personal, and genetic determinants of response to vitamin D supplementation in older adults

CONTEXT AND OBJECTIVE: Suboptimal vitamin D status can be corrected by vitamin D supplementation, but individual responses to supplementation vary. We aimed to examine genetic and nongenetic determinants of change in serum 25-hydroxyvitamin D (25(OH)D) after supplementation. DESIGN AND PARTICIPANTS: We used data from a pilot randomized controlled trial in which 644 adults aged 60 to 84 years were randomly assigned to monthly doses of placebo, 30 000 IU, or 60 000 IU vitamin D3 for 12 months. Baseline characteristics were obtained from a self-administered questionnaire. Eighty-eight single-nucleotide polymorphisms (SNPs) in 41 candidate genes were genotyped using Sequenom MassArray technology. Serum 25(OH)D levels before and after the intervention were measured using the Diasorin Liaison platform immunoassay. We used linear regression models to examine associations between genetic and nongenetic factors and change in serum 25(OH)D levels. RESULTS: Supplement dose and baseline 25(OH)D level explained 24% of the variability in response to supplementation. Body mass index, self-reported health status, and ambient UV radiation made a small additional contribution. SNPs in CYP2R1, IRF4, MC1R, CYP27B1, VDR, TYRP1, MCM6, and HERC2 were associated with change in 25(OH)D level, although only CYP2R1 was significant after adjustment for multiple testing. Models including SNPs explained a similar proportion of variability in response to supplementation as models that included personal and environmental factors. CONCLUSION: Stepwise regression analyses suggest that genetic variability may be associated with response to supplementation, perhaps suggesting that some people might need higher doses to reach optimal 25(OH)D levels or that there is variability in the physiologically normal level of 25(OH)D.

The plasticity of NBS resistance genes in sorghum is driven by multiple evolutionary processes

Background Increased disease resistance is a key target of cereal breeding programs, with disease outbreaks continuing to threaten global food production, particularly in Africa. Of the disease resistance gene families, the nucleotide-binding site plus leucine-rich repeat (NBS-LRR) family is the most prevalent and ancient and is also one of the largest gene families known in plants. The sequence diversity in NBS-encoding genes was explored in sorghum, a critical food staple in Africa, with comparisons to rice and maize and with comparisons to fungal pathogen resistance QTL. Results In sorghum, NBS-encoding genes had significantly higher diversity in comparison to non NBS-encoding genes and were significantly enriched in regions of the genome under purifying and balancing selection, both through domestication and improvement. Ancestral genes, pre-dating species divergence, were more abundant in regions with signatures of selection than in regions not under selection. Sorghum NBS-encoding genes were also significantly enriched in the regions of the genome containing fungal pathogen disease resistance QTL; with the diversity of the NBS-encoding genes influenced by the type of co-locating biotic stress resistance QTL. Conclusions NBS-encoding genes are under strong selection pressure in sorghum, through the contrasting evolutionary processes of purifying and balancing selection. Such contrasting evolutionary processes have impacted ancestral genes more than species-specific genes. Fungal disease resistance hot-spots in the genome, with resistance against multiple pathogens, provides further insight into the mechanisms that cereals use in the “arms race” with rapidly evolving pathogens in addition to providing plant breeders with selection targets for fast-tracking the development of high performing varieties with more durable pathogen resistance.

A transcriptome resource for the koala (Phascolarctos cinereus) : insights into koala retrovirus transcription and sequence diversity

Background The koala, Phascolarctos cinereus, is a biologically unique and evolutionarily distinct Australian arboreal marsupial. The goal of this study was to sequence the transcriptome from several tissues of two geographically separate koalas, and to create the first comprehensive catalog of annotated transcripts for this species, enabling detailed analysis of the unique attributes of this threatened native marsupial, including infection by the koala retrovirus. Results RNA-Seq data was generated from a range of tissues from one male and one female koala and assembled de novo into transcripts using Velvet-Oases. Transcript abundance in each tissue was estimated. Transcripts were searched for likely protein-coding regions and a non-redundant set of 117,563 putative protein sequences was produced. In similarity searches there were 84,907 (72%) sequences that aligned to at least one sequence in the NCBI nr protein database. The best alignments were to sequences from other marsupials. After applying a reciprocal best hit requirement of koala sequences to those from tammar wallaby, Tasmanian devil and the gray short-tailed opossum, we estimate that our transcriptome dataset represents approximately 15,000 koala genes. The marsupial alignment information was used to look for potential gene duplications and we report evidence for copy number expansion of the alpha amylase gene, and of an aldehyde reductase gene. Koala retrovirus (KoRV) transcripts were detected in the transcriptomes. These were analysed in detail and the structure of the spliced envelope gene transcript was determined. There was appreciable sequence diversity within KoRV, with 233 sites in the KoRV genome showing small insertions/deletions or single nucleotide polymorphisms. Both koalas had sequences from the KoRV-A subtype, but the male koala transcriptome has, in addition, sequences more closely related to the KoRV-B subtype. This is the first report of a KoRV-B-like sequence in a wild population. Conclusions This transcriptomic dataset is a useful resource for molecular genetic studies of the koala, for evolutionary genetic studies of marsupials, for validation and annotation of the koala genome sequence, and for investigation of koala retrovirus. Annotated transcripts can be browsed and queried at http://koalagenome.org

Centralspindlin and α-catenin regulate Rho signalling at the epithelial zonula adherens

The biological impact of Rho depends critically on the precise subcellular localization of its active, GTP-loaded form. This can potentially be determined by the balance between molecules that promote nucleotide exchange or GTP hydrolysis. However, how these activities may be coordinated is poorly understood. We now report a molecular pathway that achieves exactly this coordination at the epithelial zonula adherens. We identify an extramitotic activity of the centralspindlin complex, better understood as a cytokinetic regulator, which localizes to the interphase zonula adherens by interacting with the cadherin-associated protein, α-catenin. Centralspindlin recruits the RhoGEF, ECT2, to activate Rho and support junctional integrity through myosin IIA. Centralspindlin also inhibits the junctional localization of p190 B RhoGAP, which can inactivate Rho. Thus, a conserved molecular ensemble that governs Rho activation during cytokinesis is used in interphase cells to control the Rho GTPase cycle at the zonula adherens

Comparative genomics of koala, cattle and sheep strains of Chlamydia pecorum

Background Chlamydia pecorum is an important pathogen of domesticated livestock including sheep, cattle and pigs. This pathogen is also a key factor in the decline of the koala in Australia. We sequenced the genomes of three koala C. pecorum strains, isolated from the urogenital tracts and conjunctiva of diseased koalas. The genome of the C. pecorum VR629 (IPA) strain, isolated from a sheep with polyarthritis, was also sequenced. Results Comparisons of the draft C. pecorum genomes against the complete genomes of livestock C. pecorum isolates revealed that these strains have a conserved gene content and order, sharing a nucleotide sequence similarity > 98%. Single nucleotide polymorphisms (SNPs) appear to be key factors in understanding the adaptive process. Two regions of the chromosome were found to be accumulating a large number of SNPs within the koala strains. These regions include the Chlamydia plasticity zone, which contains two cytotoxin genes (toxA and toxB), and a 77 kbp region that codes for putative type III effector proteins. In one koala strain (MC/MarsBar), the toxB gene was truncated by a premature stop codon but is full-length in IPTaLE and DBDeUG. Another five pseudogenes were also identified, two unique to the urogenital strains C. pecorum MC/MarsBar and C. pecorum DBDeUG, respectively, while three were unique to the koala C. pecorum conjunctival isolate IPTaLE. An examination of the distribution of these pseudogenes in C. pecorum strains from a variety of koala populations, alongside a number of sheep and cattle C. pecorum positive samples from Australian livestock, confirmed the presence of four predicted pseudogenes in koala C. pecorum clinical samples. Consistent with our genomics analyses, none of these pseudogenes were observed in the livestock C. pecorum samples examined. Interestingly, three SNPs resulting in pseudogenes identified in the IPTaLE isolate were not found in any other C. pecorum strain analysed, raising questions over the origin of these point mutations. Conclusions The genomic data revealed that variation between C. pecorum strains were mainly due to the accumulation of SNPs, some of which cause gene inactivation. The identification of these genetic differences will provide the basis for further studies to understand the biology and evolution of this important animal pathogen. Keywords: Chlamydia pecorum; Single nucleotide polymorphism; Pseudogene; Cytotoxin