miRDeepFinder is a software package developed to identify and functionally analyze plant microRNAs (miRNAs) and their targets from small RNA datasets obtained from deep sequencing. The functions available in miRDeepFinder include pre-processing of raw data, identifying conserved miRNAs, mining and classifying novel miRNAs, miRNA expression profiling, predicting miRNA targets, and gene pathway and gene network analysis involving miRNAs. The fundamental design of miRDeepFinder is based on miRNA biogenesis, miRNA-mediated gene regulation and target recognition, such as perfect or near perfect hairpin structures, different read abundances of miRNA and miRNA*, and targeting patterns of plant miRNAs. To test the accuracy and robustness of miRDeepFinder, we analyzed a small RNA deep sequencing dataset of Arabidopsis thaliana published in the GEO database of NCBI. Our test retrieved 128 of 131 (97.7%) known miRNAs that have a more than 3 read count in Arabidopsis. Because many known miRNAs are not associated with miRNA*s in small RNA datasets, miRDeepFinder was also designed to recover miRNA candidates without the presence of miRNA*. To mine as many miRNAs as possible, miRDeepFinder allows users to compare mature miRNAs and their miRNA*s with other small RNA datasets from the same species. Cleaveland software package was also incorporated into miRDeepFinder for miRNA target identification using degradome sequencing analysis. Using this new computational tool, we identified 13 novel miRNA candidates with miRNA*s from Arabidopsis and validated 12 of them experimentally. Interestingly, of the 12 verified novel miRNAs, a miRNA named AC1 spans the exons of two genes (UTG71C4 and UGT71C3). Both the mature AC1 miRNA and its miRNA* were also found in four other small RNA datasets. We also developed a tool, "miRNA primer designer" to design primers for any type of miRNAs. miRDeepFinder provides a powerful tool for analyzing small RNA datasets from all species, with or without the availability of genome information. miRDeepFinder and miRNA primer designer are freely available at http://www.leonxie.com/DeepFinder.php and at http://www.leonxie.com/miRNAprimerDesigner.php , respectively. A program (called RefFinder: http://www.leonxie.com/referencegene.php ) was also developed for assessing the reliable reference genes for gene expression analysis, including miRNAs.

No abstract

No abstract

Background Heavy metal toxicity has become a major threat to sustainable crop production worldwide. Thus, considerable interest has been placed on deciphering the mechanisms that allow plants to combat heavy metal stress. Strategies to deal with heavy metals are largely focused on detoxification, transport and/or sequestration. The P 1B subfamily of the H eavy M etal-transporting P-type A TPases (HMAs) was shown to play a crucial role in the uptake and translocation of heavy metals in plants. Here, we report the locus-specific expression changes in the rice HMA genes together with several low-copy cellular genes and transposable elements upon the heavy metal treatment and monitored the transgenerational inheritance of the altered expression states. We reveal that plants cope with heavy metal stress by making heritable changes in gene expression and further determined gene-specific responses to heavy metal stress. Results We found most HMA genes were upregulated in response to heavy metal stress, and furthermore found evidence of transgenerational memory via changes in gene regulation even after the removal of heavy metals. To explore whether DNA methylation was also altered in response to the heavy metal stress, we selected a Tos17 retrotransposon for bisulfite sequencing and studied its methylation state across three generations. We found the DNA methylation state of Tos17 was altered in response to the heavy metal stress and showed transgenerational inheritance. Conclusions Collectively, the present study elucidates heritable changes in gene expression and DNA methylation in rice upon exposure to heavy metal stress and discusses implications of this knowledge in breeding for heavy metal tolerant crops. Electronic supplementary material The online version of this article (10.1186/s12870-019-1887-7) contains supplementary material, which is available to authorized users.

No abstract

Daphnia on the Tibetan Plateau has been little studied, and information on species diversity and biogeography is lacking. Here, we conducted a 4‐year survey using the barcoding fragment of the mitochondrial COI gene to determine the distribution and diversity of Daphnia species found across the Plateau. Our results show that species richness is higher than previously thought, with total described and provisional species number doubling from 5 to 10. Six of the taxonomic units recovered by DNA taxonomy agreed well with morphology, but DNA barcoding distinguished three clades each for the D. longispina (D. galeata, D. dentifera, and D. longispina) and D. pulex (D. pulex, D. cf. tenebrosa, and D. pulicaria) complexes. The sequence divergence between congeneric species varied within a large range, from 9.25% to 30.71%. The endemic D. tibetana was the most common and widespread species, occurring in 12 hyposaline to mesosaline lakes. The lineage of D. longispina is the first confirmed occurrence in west Tibet.

Seashore paspalum (Paspalum vaginatum) is among the most salt- and cadmium-tolerant warm-season perennial grass species widely used as turf or forage. The objective of this study was to select stable reference genes for quantitative real-time polymerase chain reaction (qRT-PCR) analysis of seashore paspalum in response to four abiotic stresses. The stability of 12 potential reference genes was evaluated by four programs (geNorm, NormFinder, BestKeeper, and RefFinder). U2AF combined with glyceraldehyde-3-phosphate dehydrogenase (GAPDH) showed stable expression in Cd-treated leaves and cold-treated roots. U2AF and FBOX were the most stable reference genes in Cd-treated roots and cold-treated leaves. In Polyethylene Glycol (PEG)- or salt-treated roots, the reference gene U2AF paired with either ACT or CYP were stable. SAND and CACS exhibited the most stability in salt-treated leaves, and combining UPL, PP2A, and EF1a was most suitable for PEG-treated leaves. The stability of U2AF and instability of UPL and TUB was validated by analyzing the expression levels of four target genes (MT2a, VP1, PIP1, and Cor413), and were shown to be capable of detecting subtle changes in expression levels of the target genes in seashore paspalum. This study demonstrated that FBOX, U2AF, and PP2A could be used in future molecular studies that aim to understand the mechanisms of abiotic stress tolerance in seashore paspalum.

BackgroundFlorfenicol is a derivative of chloramphenicol that is used only for the treatment of animal diseases. A key resistance gene for florfenicol, floR, can spread among bacteria of the same and different species or genera through horizontal gene transfer. To analyze the potential transmission of resistance genes between animal and human pathogens, we investigated floR in Klebsiella pneumoniae isolates from patient samples. floR in human pathogens may originate from animal pathogens and would reflect the risk to human health of using antimicrobial agents in animals.MethodsPCR was used to identify floR-positive strains. The floR genes were cloned, and the minimum inhibitory concentrations (MICs) were determined to assess the relative resistance levels of the genes and strains. Sequencing and comparative genomics methods were used to analyze floR gene-related sequence structure as well as the molecular mechanism of resistance dissemination.ResultsOf the strains evaluated, 20.42% (67/328) were resistant to florfenicol, and 86.96% (20/23) of the floR-positive strains demonstrated high resistance to florfenicol with MICs ≥512 μg/mL. Conjugation experiments showed that transferrable plasmids carried the floR gene in three isolates. Sequencing analysis of a plasmid approximately 125 kb in size (pKP18–125) indicated that the floR gene was flanked by multiple copies of mobile genetic elements. Comparative genomics analysis of a 9-kb transposon-like fragment of pKP18–125 showed that an approximately 2-kb sequence encoding lysR-floR-virD2 was conserved in the majority (79.01%, 83/105) of floR sequences collected from NCBI nucleotide database. Interestingly, the most similar sequence was a 7-kb fragment of plasmid pEC012 from an Escherichia coli strain isolated from a chicken.ConclusionsIdentified on a transferable plasmid in the human pathogen K. pneumoniae, the floR gene may be disseminated through horizontal gene transfer from animal pathogens. Studies on the molecular mechanism of resistance gene dissemination in different bacterial species of animal origin could provide useful information for preventing or controlling the spread of resistance between animal and human pathogens.Electronic supplementary materialThe online version of this article (10.1186/s13756-018-0415-0) contains supplementary material, which is available to authorized users.