A membrane associated tandem kinase from wild emmer wheat confers broad-spectrum resistance to powdery mildew.

Crop wild relatives offer natural variations of disease resistance for crop improvement. Here, we report the isolation of broad-spectrum powdery mildew resistance gene Pm36, originated from wild emmer wheat, that encodes a tandem kinase with a transmembrane domain (WTK7-TM) through the combination of map-based cloning, PacBio SMRT long-read genome sequencing, mutagenesis, and transformation. Mutagenesis assay reveals that the two kinase domains and the transmembrane domain of WTK7-TM are critical for the powdery mildew resistance function. Consistently, in vitro phosphorylation assay shows that two kinase domains are indispensable for the kinase activity of WTK7-TM. Haplotype analysis uncovers that Pm36 is an orphan gene only present in a few wild emmer wheat, indicating its single ancient origin and potential contribution to the current wheat gene pool. Overall, our findings not only provide a powdery mildew resistance gene with great potential in wheat breeding but also sheds light into the mechanism underlying broad-spectrum resistance.

Development of a multiplex PCR assay for the detection of tomato wilt caused by co-infection of Fusarium brachygibbosum, Fusarium oxysporum, and Ralstonia solanacearum based on comparative genomics.

Tomato is widely consumed worldwide as fresh or processed food products. However, soil-borne diseases of tomato plants caused by co-infection of various pathogens result in great economic losses to the tomato industry. It is difficult to accurately identify and diagnose soil-borne diseases of tomato plants caused by pathogen complexes. In this study, we investigated field diseases of tomato plants by pathogen isolation and molecular identification and found that tomato wilt was caused by co-infection of Fusarium brachygibbosum, Fusarium oxysporum, and Ralstonia solanacearum. Therefore, the development of a method for simultaneous detection of DNA from Fusarium brachygibbosum, Fusarium oxysporum, and Ralstonia solanacearum can efficiently and accurately monitor disease development at different growth stages of tomato plants. In this study, we performed a comparative genomic analysis of Fusarium brachygibbosum, Fusarium oxysporum, and Ralstonia solanacearum, and determined the primer sets for simultaneous detection of DNA from these target pathogens. Then, we tested the reagent and condition parameters of multiplex PCR, including primers, dNTP and Mg2+ concentrations, and the annealing temperatures, to determine the optimal parameters of a multiplex PCR system. We evaluated the specificity, sensitivity and stability of the multiplex PCR system based on the optimized reaction conditions. The multiplex PCR system can specifically identify 13 target pathogens from 57 different fungal and bacterial pathogens, at the lower detection limit of the three target pathogens at concentrations of 100pg/ul. In addition, we can accurately identify the three pathogens in tomato plants using the optimized multiplex PCR method. These results demonstrated that the multiplex PCR method developed in this study can simultaneously detect DNA from Fusarium brachygibbosum, Fusarium oxysporum, and Ralstonia solanacearum in a single PCR to accurately identify and diagnose the pathogen causing tomato wilt.

Indole-3-Carboxylic Acid From the Endophytic Fungus Lasiodiplodia pseudotheobromae LPS-1 as a Synergist Enhancing the Antagonism of Jasmonic Acid Against Blumeria graminis on Wheat.

The discovery of natural bioactive compounds from endophytes or medicinal plants against plant diseases is an attractive option for reducing the use of chemical fungicides. In this study, three compounds, indole-3-carbaldehyde, indole-3-carboxylic acid (3-ICA), and jasmonic acid (JA), were isolated from the EtOAc extract of the culture filtrate of the endophytic fungus Lasiodiplodia pseudotheobromae LPS-1, which was previously isolated from the medicinal plant, Ilex cornuta. Some experiments were conducted to further determine the antifungal activity of these compounds on wheat powdery mildew. The results showed that JA was much more bioactive than indole-3-carbaldehyde and 3-ICA against Blumeria graminis, and the disease severity caused by B. graminis decreased significantly with the concentration increase of JA treatment. The assay of the interaction of 3-ICA and JA indicated that there was a significant synergistic effect between the two compounds on B. graminis in each of the ratios of 3-ICA to JA (3-ICA:JA) ranging from 1:9 to 9:1. When the compound ratio of 3-ICA to JA was 2:8, the synergistic coefficient was the highest as 22.95. Meanwhile, a histological investigation indicated that, under the treatment of JA at 500 µg/ml or 3-ICA:JA (2:8) at 40 µg/ml, the appressorium development and haustorium formation of B. graminis were significantly inhibited. Taken together, we concluded that JA plays an important role in the infection process of B. graminis and that 3-ICA as a synergist of JA enhances the antagonism against wheat powdery mildew.

Diversity and similarity of wheat powdery mildew resistance among three allelic functional genes at the Pm60 locus.

Cultivated wheat is continually exposed to various pathogens. Blumeria graminis f. sp. tritici (Bgt) causes powdery mildew disease and significant yield loss. Pm60 was cloned from Triticum urartu and confers race-specific powdery mildew resistance in wheat. Pm60a and Pm60b are allelic variants of Pm60 and have two leucine-rich repeat motifs deletions and insertions, respectively, which were detected in other T. urartu accessions. Through map-based cloning, virus-induced gene silencing, and stable transformation assays, we demonstrated that Pm60a and Pm60b conferred Bgt E09 resistance resembling that provided by Pm60. However, the homozygous Pm60a (but not Pm60 or Pm60b) transformants driven by the native promoters lacked race-specific resistance when they were inoculated with Bgt E18. As all three T. urartu accessions contained the three foregoing alleles, they had high resistance to Bgt E18. Pyramiding Pm60a with either of the allelic genes in F1 plants did not cause mutual allele suppression or interference with Bgt E18 resistance. Deletion (but not insertion) of the two leucine-rich repeat motifs in Pm60a substantially narrowed the resistance spectrum. In T. urartu accession PI428210, we identified another locus adjacent to Pm60a and resistant to Bgt E18. Characterization of the alleles at the Pm60 locus revealed their diversity and similarity and may facilitate wheat breeding for resistance to powdery mildew disease caused by B. graminis f. sp. tritici.

Complete Genome Sequence of Bacillus amyloliquefaciens EA19, an Endophytic Bacterium with Biocontrol Potential Isolated from Erigeron annuus.

Bacillus amyloliquefaciens strain EA19 is an endophyte isolated from Erigeron annuus with antifungal activity against Blumeria graminis f. sp. tritici, Magnaporthe oryzae, and Fusarium graminearum. The genome sequence of this strain is 3.96 Mb and contains 3,421 coding sequences, which will facilitate an understanding of the mechanisms of biocontrol.

Pepper root rot resistance and pepper yield are enhanced through biological agent G15 soil amelioration.

Pepper root rot is a serious soil-borne disease that hinders pepper production, and efforts are being made to identify biological agents that can prevent and control pepper root rot. Our group recently discovered and produced a biological agent, named G15, which reduces the diversity and richness of fungi and bacteria when applied to pepper fields. In the soil of the G15-treatment condition, the pathogenic fungus Fusarium was inhibited, while the richness of beneficial bacteria Rhodanobacter was increased. Also, the ammonia nitrogen level was decreased in the G15-treatment soil, and the pH, total carbon, and total potassium levels were increased. Compared to the control condition, pepper yield was increased in the treatment group (by 16,680 kg acre-1). We found that G15 could alter the microbial community structure of the pepper rhizosphere. These changes alter the physical and chemical properties of the soil and, ultimately, improve resistance to pepper root rot and increase pepper yield.

The influence of rhizosphere soil fungal diversity and complex community structure on wheat root rot disease.

Wheat root rot disease due to soil-borne fungal pathogens leads to tremendous yield losses worth billions of dollars worldwide every year. It is very important to study the relationship between rhizosphere soil fungal diversity and wheat roots to understand the occurrence and development of wheat root rot disease. A significant difference in fungal diversity was observed in the rhizosphere soil of healthy and diseased wheat roots in the heading stage, but the trend was the opposite in the filling stage. The abundance of most genera with high richness decreased significantly from the heading to the filling stage in the diseased groups; the richness of approximately one-third of all genera remained unchanged, and only a few low-richness genera, such as Fusarium and Ceratobasidium, had a very significant increase from the heading to the filling stage. In the healthy groups, the abundance of most genera increased significantly from the heading to filling stage; the abundance of some genera did not change markedly, or the abundance of very few genera increased significantly. Physical and chemical soil indicators showed that low soil pH and density, increases in ammonium nitrogen, nitrate nitrogen and total nitrogen contributed to the occurrence of wheat root rot disease. Our results revealed that in the early stages of disease, highly diverse rhizosphere soil fungi and a complex community structure can easily cause wheat root rot disease. The existence of pathogenic fungi is a necessary condition for wheat root rot disease, but the richness of pathogenic fungi is not necessarily important. The increases in ammonium nitrogen, nitrate nitrogen and total nitrogen contributed to the occurrence of wheat root rot disease. Low soil pH and soil density are beneficial to the occurrence of wheat root rot disease.

A CNL protein in wild emmer wheat confers powdery mildew resistance.

Powdery mildew, a fungal disease caused by Blumeria graminis f. sp. tritici (Bgt), has a serious impact on wheat production. Loss of resistance in cultivars prompts a continuing search for new sources of resistance. Wild emmer wheat (Triticum turgidum ssp. dicoccoides, WEW), the progenitor of both modern tetraploid and hexaploid wheats, harbors many powdery mildew resistance genes. We report here the positional cloning and functional characterization of Pm41, a powdery mildew resistance gene derived from WEW, which encodes a coiled-coil, nucleotide-binding site and leucine-rich repeat protein (CNL). Mutagenesis and stable genetic transformation confirmed the function of Pm41 against Bgt infection in wheat. We demonstrated that Pm41 was present at a very low frequency (1.81%) only in southern WEW populations. It was absent in other WEW populations, domesticated emmer, durum, and common wheat, suggesting that the ancestral Pm41 was restricted to its place of origin and was not incorporated into domesticated wheat. Our findings emphasize the importance of conservation and exploitation of the primary WEW gene pool, as a valuable resource for discovery of resistance genes for improvement of modern wheat cultivars.

A rare gain of function mutation in a wheat tandem kinase confers resistance to powdery mildew.

Powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt), is one of the most destructive diseases that pose a great threat to wheat production. Wheat landraces represent a rich source of powdery mildew resistance. Here, we report the map-based cloning of powdery mildew resistance gene Pm24 from Chinese wheat landrace Hulutou. It encodes a tandem kinase protein (TKP) with putative kinase-pseudokinase domains, designated WHEAT TANDEM KINASE 3 (WTK3). The resistance function of Pm24 was validated by transgenic assay, independent mutants, and allelic association analyses. Haplotype analysis revealed that a rare 6-bp natural deletion of lysine-glycine codons, endemic to wheat landraces of Shaanxi Province, China, in the kinase I domain (Kin I) of WTK3 is critical for the resistance function. Transgenic assay of WTK3 chimeric variants revealed that only the specific two amino acid deletion, rather than any of the single or more amino acid deletions, in the Kin I of WTK3 is responsible for gaining the resistance function of WTK3 against the Bgt fungus.

Comprehensive transcriptome analysis reveals genes potentially involved in isoflavone biosynthesis in Pueraria thomsonii Benth.

Pueraria thomsonii Benth is an important medicinal plant. Transcriptome sequencing, unigene assembly, the annotation of transcripts and the study of gene expression profiles play vital roles in gene function research. However, the full-length transcriptome of P. thomsonii remains unknown. Here, we obtained 44,339 nonredundant transcripts of P. thomsonii by using the PacBio RS II Isoform and Illumina sequencing platforms, of which 43,195 were annotated genes. Compared with the expression levels in the plant roots, those of transcripts with a |fold change| ≥ 4 and FDR < 0.01 in the leaves or stems were assigned as differentially expressed transcripts (DETs). In total, we found 9,225 DETs, 32 of which came from structural genes that were potentially involved in isoflavone biosynthesis. The expression profiles of 8 structural genes from the RNA-Seq data were validated by qRT-PCR. We identified 437 transcription factors (TFs) that were positively or negatively correlated with at least 1 of the structural genes involved in isoflavone biosynthesis using Pearson correlation coefficients (r) (r > 0.8 or r < -0.8). We also identified a total of 32 microRNAs (miRNAs), which targeted 805 transcripts. These miRNAs caused enriched function in 'ATP binding', 'defense response', 'ADP binding', and 'signal transduction'. Interestingly, MIR156a potentially promoted isoflavone biosynthesis by repressing SBP, and MIR319 promoted isoflavone biosynthesis by repressing TCP and HB-HD-ZIP. Finally, we identified 2,690 alternative splicing events, including that of the structural genes of trans-cinnamate 4-monooxygenase and pullulanase, which are potentially involved in the biosynthesis of isoflavone and starch, respectively, and of three TFs potentially involved in isoflavone biosynthesis. Together, these results provide us with comprehensive insight into the gene expression and regulation of P. thomsonii.

The AvrPm3-Pm3 effector-NLR interactions control both race-specific resistance and host-specificity of cereal mildews on wheat.

The wheat Pm3 resistance gene against the powdery mildew pathogen occurs as an allelic series encoding functionally different immune receptors which induce resistance upon recognition of isolate-specific avirulence (AVR) effectors from the pathogen. Here, we describe the identification of five effector proteins from the mildew pathogens of wheat, rye, and the wild grass Dactylis glomerata, specifically recognized by the PM3B, PM3C and PM3D receptors. Together with the earlier identified AVRPM3A2/F2, the recognized AVRs of PM3B/C, (AVRPM3B2/C2), and PM3D (AVRPM3D3) belong to a large group of proteins with low sequence homology but predicted structural similarities. AvrPm3b2/c2 and AvrPm3d3 are conserved in all tested isolates of wheat and rye mildew, and non-host infection assays demonstrate that Pm3b, Pm3c, and Pm3d are also restricting the growth of rye mildew on wheat. Furthermore, divergent AVR homologues from non-adapted rye and Dactylis mildews are recognized by PM3B, PM3C, or PM3D, demonstrating their involvement in host specificity.

Transcriptome analysis reveals downregulation of virulence-associated genes expression in a low virulence Verticillium dahliae strain.

Verticillium dahliae causes wilt diseases and early senescence in numerous plants, including agricultural crops such as cotton. In this study, we studied two closely related V. dahliae strains, and found that V991w showed significantly reduced virulence on cotton than V991b. Comprehensive transcriptome analysis revealed various differentially expressed genes between the two strains, with more genes repressed in V991w. The downregulated genes in V991w were involved in production of hydrophobins, melanin, predicted aflatoxin, and membrane proteins, most of which are related to pathogenesis and multidrug resistance. Consistently, melanin production in V991w in vitro was compromised. We next obtained genomic variations between the two strains, demonstrating that transcription factor genes containing fungi specific transcription factor domain and fungal Zn2-Cys6 binuclear cluster domain were enriched in V991w, which might be related to pathogenicity-related genes downregulation. Thus, this study supports a model in which some virulence factors involved in V. dahliae pathogenicity were pre-expressed during in vitro growth before host interaction.

Distinct domains of the AVRPM3A2/F2 avirulence protein from wheat powdery mildew are involved in immune receptor recognition and putative effector function.

Recognition of the AVRPM3A2/F2 avirulence protein from powdery mildew by the wheat PM3A/F immune receptor induces a hypersensitive response after co-expression in Nicotiana benthamiana. The molecular determinants of this interaction and how they shape natural AvrPm3a2/f2 allelic diversity are unknown. We sequenced the AvrPm3a2/f2 gene in a worldwide collection of 272 mildew isolates. Using the natural polymorphisms of AvrPm3a2/f2 as well as sequence information from related gene family members, we tested 85 single-residue-altered AVRPM3A2/F2 variants with PM3A, PM3F and PM3FL456P/Y458H (modified for improved signaling) in Nicotiana benthamiana for effects on recognition. An intact AvrPm3a2/f2 gene was found in all analyzed isolates and the protein variant recognized by PM3A/F occurred globally at high frequencies. Single-residue alterations in AVRPM3A2/F2 mostly disrupted, but occasionally enhanced, the recognition response by PM3A, PM3F and PM3FL456P/Y458H . Residues enhancing hypersensitive responses constituted a protein domain separate from both naturally occurring polymorphisms and positively selected residues of the gene family. These results demonstrate the utility of using gene family sequence diversity to screen residues for their role in recognition. This approach identified a putative interaction surface in AVRPM3A2/F2 not polymorphic in natural alleles. We conclude that molecular mechanisms besides recognition drive AvrPm3a2/f2 diversification.

A Mucin-Like Protein of Planthopper Is Required for Feeding and Induces Immunity Response in Plants.

The brown planthopper, Nilaparvata lugens, is a pest that threatens rice (Oryza sativa) production worldwide. While feeding on rice plants, planthoppers secrete saliva, which plays crucial roles in nutrient ingestion and modulating plant defense responses, although the specific functions of salivary proteins remain largely unknown. We identified an N. lugens-secreted mucin-like protein (NlMLP) by transcriptome and proteome analyses and characterized its function, both in brown planthopper and in plants. NlMLP is highly expressed in salivary glands and is secreted into rice during feeding. Inhibition of NlMLP expression in planthoppers disturbs the formation of salivary sheaths, thereby reducing their performance. In plants, NlMLP induces cell death, the expression of defense-related genes, and callose deposition. These defense responses are related to Ca2+ mobilization and the MEK2 MAP kinase and jasmonic acid signaling pathways. The active region of NlMLP that elicits plant responses is located in its carboxyl terminus. Our work provides a detailed characterization of a salivary protein from a piercing-sucking insect other than aphids. Our finding that the protein functions in plant immune responses offers new insights into the mechanism underlying interactions between plants and herbivorous insects.

Transcriptome Analyses Shed New Insights into Primary Metabolism and Regulation of Blumeria graminis f. sp. tritici during Conidiation.

Conidia of the obligate biotrophic fungal pathogen Blumeria graminis f. sp. tritici (Bgt) play a vital role in its survival and rapid dispersal. However, little is known about the genetic basis for its asexual reproduction. To uncover the primary metabolic and regulatory events during conidiation, we sequenced the transcriptome of Bgt epiphytic structures at 3 (vegetative hyphae growth), 4 (foot cells initiation), and 5 (conidiophore erection) days post-inoculation (dpi). RNA-seq analyses identified 556 and 404 (combined 685) differentially expressed genes (DEGs) at 4 and 5 dpi compared with their expression levels at 3 dpi, respectively. We found that several genes involved in the conversion from a variety of sugars to glucose, glycolysis, the tricarboxylic acid cycle (TAC), the electron transport chain (ETC), and unsaturated fatty acid oxidation were activated during conidiation, suggesting that more energy supply is required during this process. Moreover, we found that glucose was converted into glycogen, which was accumulated in developing conidiophores, indicating that it could be the primary energy storage molecule in Bgt conidia. Clustering for the expression profiles of 91 regulatory genes showed that calcium (Ca2+), H2O2, and phosphoinositide (PIP) signaling were involved in Bgt conidiation. Furthermore, a strong accumulation of H2O2 in developing conidiophores was detected. Application of EGTA, a Ca2+ chelator, and trifluoperazine dihydrochloride (TFP), a calmodulin (CaM) antagonist, markedly suppressed the generation of H2O2, affected foot cell and conidiophore development and reduced conidia production significantly. These results suggest that Ca2+ and H2O2 signaling play important roles in conidiogenesis and a crosslink between them is present. In addition to some conidiation-related orthologs known in other fungi, such as the velvet complex components, we identified several other novel B. graminis-specific genes that have not been previously found to be implicated in fungal conidiation, reflecting a unique molecular mechanism underlying asexual development of cereal powdery mildews.

Identification and evaluation of resistance to powdery mildew and yellow rust in a wheat mapping population.

Deployment of cultivars with genetic resistance is an effective approach to control the diseases of powdery mildew (PM) and yellow rust (YR). Chinese wheat cultivar XK0106 exhibits high levels of resistance to both diseases, while cultivar E07901 has partial, adult plant resistance (APR). The aim of this study was to map resistance loci derived from the two cultivars and analyze their effects against PM and YR in a range of environments. A doubled haploid population (388 lines) was used to develop a framework map consisting of 117 SSR markers, while a much higher density map using the 90K Illumina iSelect SNP array was produced with a subset of 80 randomly selected lines. Seedling resistance was characterized against a range of PM and YR isolates, while field scores in multiple environments were used to characterize APR. Composite interval mapping (CIM) of seedling PM scores identified two QTLs (QPm.haas-6A and QPm.haas-2A), the former being located at the Pm21 locus. These QTLs were also significant in field scores, as were Qpm.haas-3A and QPm.haas-5A. QYr.haas-1B-1 and QYr.haas-2A were identified in field scores of YR and were located at the Yr24/26 and Yr17 chromosomal regions respectively. A second 1B QTL, QYr.haas-1B-2 was also identified. QPm.haas-2A and QYr.haas-1B-2 are likely to be new QTLs that have not been previously identified. Effects of the QTLs were further investigated in multiple environments through the testing of selected lines predicted to contain various QTL combinations. Significant additive interactions between the PM QTLs highlighted the ability to pyramid these loci to provide higher level of resistance. Interactions between the YR QTLs gave insights into the pathogen populations in the different locations as well as showing genetic interactions between these loci.

Transcriptome analysis reveals the complexity of alternative splicing regulation in the fungus Verticillium dahliae.

BACKGROUND: Alternative splicing (AS) regulation is extensive and shapes the functional complexity of higher organisms. However, the contribution of alternative splicing to fungal biology is not well studied. RESULTS: This study provides sequences of the transcriptomes of the plant wilt pathogen Verticillium dahliae, using two different strains and multiple methods for cDNA library preparations. We identified alternatively spliced mRNA isoforms in over a half of the multi-exonic fungal genes. Over one-thousand isoforms involve TopHat novel splice junction; multiple types of combinatory alternative splicing patterns were identified. We showed that one Verticillium gene could use four different 5' splice sites and two different 3' donor sites to produce up to five mature mRNAs, representing one of the most sophisticated alternative splicing model in eukaryotes other than animals. Hundreds of novel intron types involving a pair of new splice sites were identified in the V. dahliae genome. All the types of AS events were validated by using RT-PCR. Functional enrichment analysis showed that AS genes are involved in most known biological functions and enriched in ATP biosynthesis, sexual/asexual reproduction, morphogenesis, signal transduction etc., predicting that the AS regulation modulates mRNA isoform output and shapes the V. dahliae proteome plasticity of the pathogen in response to the environmental and developmental changes. CONCLUSIONS: These findings demonstrate the comprehensive alternative splicing mechanisms in a fungal plant pathogen, which argues the importance of this fungus in developing complicate genome regulation strategies in eukaryotes.

AvrPm2 encodes an RNase-like avirulence effector which is conserved in the two different specialized forms of wheat and rye powdery mildew fungus.

There is a large diversity of genetically defined resistance genes in bread wheat against the powdery mildew pathogen Blumeria graminis (B. g.) f. sp. tritici. Many confer race-specific resistance to this pathogen, but until now only the mildew avirulence gene AvrPm3a2/f2 that is recognized by Pm3a/f was known molecularly. We performed map-based cloning and genome-wide association studies to isolate a candidate for the mildew avirulence gene AvrPm2. We then used transient expression assays in Nicotiana benthamiana to demonstrate specific and strong recognition of AvrPm2 by Pm2. The virulent AvrPm2 allele arose from a conserved 12 kb deletion, while there is no protein sequence diversity in the gene pool of avirulent B. g. tritici isolates. We found one polymorphic AvrPm2 allele in B. g. triticale and one orthologue in B. g. secalis and both are recognized by Pm2. AvrPm2 belongs to a small gene family encoding structurally conserved RNase-like effectors, including Avra13 from B. g. hordei, the cognate Avr of the barley resistance gene Mla13. These results demonstrate the conservation of functional avirulence genes in two cereal powdery mildews specialized on different hosts, thus providing a possible explanation for successful introgression of resistance genes from rye or other grass relatives to wheat.

Candidate Resistant Genes of Sand Pear (Pyrus pyrifolia Nakai) to Alternaria alternata Revealed by Transcriptome Sequencing.

Pear black spot (PBS) disease, which is caused by Alternaria alternata (Aa), is one of the most serious diseases affecting sand pear (Pyrus pyrifolia Nakai) cultivation worldwide. To investigate the defense mechanisms of sand pear in response to Aa, the transcriptome of a sand pear germplasm with differential resistance to Aa was analyzed using Illumina paired-end sequencing. Four libraries derived from PBS-resistant and PBS-susceptible sand pear leaves were characterized through inoculation or mock-inoculation. In total, 20.5 Gbp of sequence data and 101,632,565 reads were generated, representing 44717 genes. Approximately 66% of the genes or sequenced reads could be aligned to the pear reference genome. A large number (5213) of differentially expressed genes related to PBS resistance were obtained; 34 microsatellites were detected in these genes, and 28 genes were found to be closely related to PBS resistance. Using a transcriptome analysis in response to PBS inoculation and comparison analysis to the PHI database, 4 genes (Pbr039001, Pbr001627, Pbr025080 and Pbr023112) were considered to be promising candidates for sand pear resistance to PBS. This study provides insight into changes in the transcriptome of sand pear in response to PBS infection, and the findings have improved our understanding of the resistance mechanism of sand pear to PBS and will facilitate future gene discovery and functional genome studies of sand pear.

Bioinformatics analysis of small RNAs in pima (Gossypium barbadense L.).

Small RNAs (sRNAs) are ~20 to 24 nucleotide single-stranded RNAs that play crucial roles in regulation of gene expression. In plants, sRNAs are classified into microRNAs (miRNAs), repeat-associated siRNAs (ra-siRNAs), phased siRNAs (pha-siRNAs), cis and trans natural antisense transcript siRNAs (cis- and trans-nat siRNAs). Pima (Gossypium barbadense L.) is one of the most economically important fiber crops, producing the best and longest spinnable fiber. Although some miRNAs are profiled in Pima, little is known about siRNAs, the largest subclass of plant sRNAs. In order to profile these gene regulators in Pima, a comprehensive analysis of sRNAs was conducted by mining publicly available sRNA data, leading to identification of 678 miRNAs, 3,559,126 ra-siRNAs, 627 pha-siRNAs, 136,600 cis-nat siRNAs and 79,994 trans-nat siRNAs. The 678 miRNAs, belonging to 98 conserved and 402 lineage-specific families, were produced from 2,138 precursors, of which 297 arose from introns, exons, or intron/UTR-exon junctions of protein-coding genes. Ra-siRNAs were produced from various repeat loci, while most (97%) were yielded from retrotransposons, especially LTRs (long terminal repeats). The genes encoding auxin-signaling-related proteins, NBS-LRRs and transcription factors were major sources of pha-siRNAs, while two conserved TAS3 homologs were found as well. Most cis-NATs in Pima overlapped in enclosed and convergent orientations, while a few hybridized in divergent and coincided orientations. Most cis- and trans-nat siRNAs were produced from overlapping regions. Additionally, characteristics of length and the 5'-first nucleotide of each sRNA class were analyzed as well. Results in this study created a valuable molecular resource that would facilitate studies on mechanism of controlling gene expression.