RNA-Seq analysis of mung bean (Vigna radiata L.) roots shows differential gene expression and predicts regulatory pathways responding to taxonomically different rhizobia.

Symbiotic interaction among legume and rhizobia is a complex phenomenon which results in the formation of nitrogen-fixing nodules. Mung bean is promiscuous host however expression profile of this important legume plant in response to rhizobial infection was particularly lacking and urgently needed. We have demonstrated the pattern of gene expression of mung bean roots inoculated with two symbionts Bradyrhizobium yuanmingense Vr50 and Sinorhizobium (Ensifer) aridi Vr33 and non-inoculated control (CK). The RNA-Seq data analyzed at two growth stages i.e., 1-3 h and 10-16 days post inoculation revealed significantly higher number of differentially expressed genes (DEGs) at nodulation stage. The DEGs encoding receptor kinases identified at early stage might be involved in perception of Nod factors produced by different rhizobia. At nodulation stage important genes involved in plant hormone signal transduction, nitrogen and sulfur metabolism were identified. KEGG pathway enrichment analysis showed that metabolic pathways were most prominent in both groups (Group 1: Vr33 vs CK; Group 2: Vr50 vs CK), followed by biosynthesis of secondary metabolites, plant hormone signal transduction and biosynthesis of amino acids. Furthermore, DEGs involved in cell communication and plant hormone signal transduction were found to be different among two symbiotic systems while DEGs involved in carbon, nitrogen and sulfur metabolism were similar but their expression varied in response to two rhizobial strains. This study provides the first insight into the mechanisms underlying interactions of mung bean host with two taxonomically different symbionts (Bradyrhizobium and Sinorhizobium) and the candidate genes for better understanding the mechanisms of symbiotic host-specificity.

Exploring biocontrol and growth-promoting potential of multifaceted PGPR isolated from natural suppressive soil against the causal agent of chickpea wilt.

Chickpea is an important nutritive food crop both for humans and animals. Chickpea wilt caused by Fusarium oxysporum f.sp. ciceris (Foc) results in huge yield losses every year. Chickpea being a food crop requires the development of an eco-friendly bio-pesticide to effectively control the chickpea wilt disease. In this study, more than 50 bacterial stains isolated from the rhizosphere of healthy plants growing in wilt sick soil were examined for their Foc antagonist activities. Out of these, 17 strains showing > 90% growth inhibition of Foc were then characterized for their plant growth-promoting (PGP) and biocontrol traits. The biocontrol and PGP traits identified include amylase, hydrogen cyanide, protease, cellulase, chitinase activities, p-solubilization, nitrogen-fixing, and indole-3-acetic acid production. Two bacterial strains, IR-27 and IR-57, exhibiting the highest Foc proliferation inhibition and the PGP potential along with a consortium of four different strains (Serratia sp. IN-1, Serratia sp. IS-1, Enterobacter sp. IN-2, Enterobacter sp. IN-6) were used for controlling the chickpea wilt disease and growth promotion of the chickpea plants. Confocal laser scanning microscopy revealed their root colonization ability with partial or complete elimination of broken Foc mycelia and hyphae from roots. The bacterial inoculations particularly the consortium significantly suppressed the disease and improved the overall root morphology traits (root length, root surface area, root volume, forks, tips, and crossings), resulting in enhanced growth of the chickpea plants. Significant changes in growth (107% increase in root length, 23% increase in shoot length, and 54% increase in branches) in Foc-challenged plants were observed when inoculated with the consortium. Further investigations revealed that the chickpea plants inoculated with bacterial strains induced the expression of a number of key defence enzymes, including the phenylalanine ammonia lyase, peroxidase, polyphenol peroxidase, ß-1,3 glucanase, which might have helped the plants to thwart the pathogen attack. These findings indicate the potential of our identified bacterial strains to be used as a natural biopesticide for controlling the chickpea wilt disease.

Phylogenetic diversity analysis reveals Bradyrhizobium yuanmingense and Ensifer aridi as major symbionts of mung bean (Vigna radiata L.) in Pakistan.

The present study was carried out to evaluate the diversity of rhizobia associated with nodules of mung bean in Pakistan, because this information is necessary for inoculum development. Based on sequence analysis of 16S rRNA gene of thirty-one bacteria, 11 were assigned to genus Bradyrhizobium, 17 to Ensifer, and 3 to Rhizobium. Phylogenetic analyses on the basis of 16S-23S ITS region, atpD, recA, nifH, and nodA of representative strains revealed that B. yuanmingense is the predominant species distributed throughout different mung bean-growing areas. Among the fast-growing rhizobia, Ensifer aridi was predominant in Faisalabad, Layyah, and Rawalpindi, while E. meliloti in Thal desert. Sequence variations and phylogeny of nifH and nodA genes suggested that these genes might have been co-evolved with the housekeeping genes and maintained by vertical gene transfer in rhizobia detected in the present study. Host infectivity assay revealed the successful nodulation of host by rhizobia related to genera Bradyrhizobium, Ensifer and Rhizobium. Among all, Bradyrhizobium and Ensifer spp. inoculation exhibited a significantly higher number of nodules (11-34 nodules plant-1) and nitrogenase activity (nodule ARA 60-110 µmol g-1 h-1). Contrary to the previous studies, our data reveal that B. yuanmingense and E. aridi are predominant species forming effective nodules in mung bean in Pakistan. Furthermore, to the best of our knowledge, this is the first report showing the effective symbiosis of E. aridi, E. meliloti, and Rhizobium pusense with mung bean. The diversity of rhizobia in different habitats revealed in the present study will contribute towards designing site-specific inocula for mung bean.

Illumina sequencing of 16S rRNA tag shows disparity in rhizobial and non-rhizobial diversity associated with root nodules of mung bean (Vigna radiata L.) growing in different habitats in Pakistan.

In Rhizobium-legume symbiosis, the nodule is the most frequently studied compartment, where the endophytic/symbiotic microbiota demands critical investigation for development of specific inocula. We identified the bacterial diversity within root nodules of mung bean from different growing areas of Pakistan using Illumina sequencing of 16S rRNA gene. We observed specific OTUs related to specific site where Bradyrhizobium was found to be the dominant genus comprising of 82-94% of total rhizobia in nodules with very minor fraction of sequences from other rhizobia at three sites. In contrast, Ensifer (Sinorhizobium) was single dominant genus comprising 99.9% of total rhizobial sequences at site four. Among non-rhizobial sequences, the genus Acinetobacter was abundant (7-18% of total sequences), particularly in Bradyrhizobium-dominated nodule samples. Rhizobia and non-rhizobial PGPR isolated from nodule samples include Ensifer, Bradyrhizobium, Acinetobacter, Microbacterium and Pseudomonas strains. Co-inoculation of multi-trait PGPR Acinetobacter sp. VrB1 with either of the two rhizobia in field exhibited more positive effect on nodulation and plant growth than single-strain inoculation which favors the use of Acinetobacter as an essential component for development of mung bean inoculum. Furthermore, site-specific dominance of rhizobia and non-rhizobia revealed in this study may contribute towards decision making for development and application of specific inocula in different habitats.

Glucose dehydrogenase gene containing phosphobacteria for biofortification of Phosphorus with growth promotion of rice.

Phosphorus (P) is an essential plant nutrient, but often limited in soils for plant uptake. A major economic constraint in the rice production is excessive use of chemical fertilizers to meet the P requirement. Bioaugmentation of phosphate solubilizing rhizobacteria (PSB) can be used as promising alternative. In the present study 11 mineral PSB were isolated from Basmati rice growing areas of Pakistan. In broth medium, PSB solubilized tricalcium phosphate (27-354 µg mL-1) with concomitant decrease in pH up to 3.6 due to the production of different organic acids, predominantly gluconic acid. Of these, 4 strains also have ability to mineralize phytate (245-412 µg mL-1). Principle component analysis showed that the gluconic acid producing PSB strains (Acinetobacter sp. MR5 and Pseudomonas sp. MR7) have pronounced effect on grain yield (up to 55%), plant P (up to 67%) and soil available P (up to 67%), with 20% reduced fertilization. For simultaneous validation of gluconic acid production by MR5 and MR7 through PCR, new specific primers were designed to amplify gcd, pqqE, pqqC genes responsible for glucose dehydrogenase (gcd) mediated phosphate solubilization. These findings for the first time demonstrated Acinetobacter soli as potent P solubilizer for rice and expands our knowledge about genus specific pqq and gcd primers. These two gcd containing PSB Acinetobacter sp. MR5 (DSM 106631) and Pseudomonas sp. MR7 (DSM 106634) submitted to German culture collection (DSMZ), serve as global valuable pool to significantly increase the P uptake, growth and yield of Basmati rice with decreased dependence on chemical fertilizer in P deficit agricultural soils.

Association of plant growth-promoting Serratia spp. with the root nodules of chickpea.

Serratia species-affiliated DNA sequences have recently been discovered in the root nodules of two chickpea cultivars; however, little is known about their potential influence on chickpea plant growth. All Serratia-affiliated sequences (1136) could be grouped into two clusters at 98% DNA similarity. The major cluster, represented by 96% of sequences, was closely associated with Serratia marcescens sequences from GenBank. In the current study, we isolated two Serratia strains, 5D and RTL100, from root nodules of a field-grown Desi cultivar from Faisalabad and Thal areas, respectively. In vitro, strain 5D showed significantly higher phosphate (P) solubilization and lactic acid production than RTL100, whereas a comparable concentration of phytohormone was produced by both isolates. The application of Serratia strain 5D as an inoculum resulted in 25.55% and 30.85% increases in the grain yield of crops grown on fertile soil in irrigated areas and nutrient-deficient soil in rainfed areas, respectively, compared to the non-inoculated control. Results of plant inoculations indicated that Serratia sp. 5D and RTL100 can serve as effective microbial inoculants, particularly in nutrient-deficient soils in rainfed areas, where chickpea is the only major crop grown during the entire year.

Cultivation-Based and Molecular Assessment of Bacterial Diversity in the Rhizosheath of Wheat under Different Crop Rotations.

A field study was conducted to compare the formationand bacterial communities of rhizosheaths of wheat grown under wheat-cotton and wheat-rice rotation and to study the effects of bacterial inoculation on plant growth. Inoculation of Azospirillum sp. WS-1 and Bacillus sp. T-34 to wheat plants increased root length, root and shoot dry weight and dry weight of rhizosheathsoil when compared to non-inoculated control plants, and under both crop rotations. Comparing both crop rotations, root length, root and shoot dry weight and dry weight of soil attached with roots were higher under wheat-cotton rotation. Organic acids (citric acid, malic acid, acetic acid and oxalic acid) were detected in rhizosheaths from both rotations, with malic acid being most abundant with 24.8±2 and 21.3±1.5 µg g(-1) dry soil in wheat-cotton and wheat-rice rotation, respectively. Two sugars (sucrose, glucose) were detected in wheat rhizosheath under both rotations, with highest concentrations of sucrose (4.08±0.5 µg g(-1) and 7.36±1.0 µg g(-1)) and glucose (3.12±0.5 µg g(-1) and 3.01± µg g(-1)) being detected in rhizosheaths of non-inoculated control plants under both rotations. Diversity of rhizosheath-associated bacteria was evaluated by cultivation, as well as by 454-pyrosequencing of PCR-tagged 16S rRNA gene amplicons. A total of 14 and 12 bacterial isolates predominantly belonging to the genera Arthrobacter, Azospirillum, Bacillus, Enterobacter and Pseudomonaswere obtained from the rhizosheath of wheat grown under wheat-cotton and wheat-rice rotation, respectively. Analysis of pyrosequencing data revealed Proteobacteria, Bacteriodetes and Verrucomicrobia as the most abundant phyla in wheat-rice rotation, whereas Actinobacteria, Firmicutes, Chloroflexi, Acidobacteria, Planctomycetes and Cyanobacteria were predominant in wheat-cotton rotation. From a total of 46,971 sequences, 10.9% showed ≥97% similarity with 16S rRNA genes of 32 genera previously shown to include isolates with plant growth promoting activity (nitrogen fixation, phosphate-solubilization, IAA production). Among these, the most predominant genera were Arthrobacter, Azoarcus, Azospirillum, Bacillus, Cyanobacterium, Paenibacillus, Pseudomonas and Rhizobium.

DNA barcoding of Bemisia tabaci complex (Hemiptera: Aleyrodidae) reveals southerly expansion of the dominant whitefly species on cotton in Pakistan.

BACKGROUND: Although whiteflies (Bemisia tabaci complex) are an important pest of cotton in Pakistan, its taxonomic diversity is poorly understood. As DNA barcoding is an effective tool for resolving species complexes and analyzing species distributions, we used this approach to analyze genetic diversity in the B. tabaci complex and map the distribution of B. tabaci lineages in cotton growing areas of Pakistan. METHODS/PRINCIPAL FINDINGS: Sequence diversity in the DNA barcode region (mtCOI-5') was examined in 593 whiteflies from Pakistan to determine the number of whitefly species and their distributions in the cotton-growing areas of Punjab and Sindh provinces. These new records were integrated with another 173 barcode sequences for B. tabaci, most from India, to better understand regional whitefly diversity. The Barcode Index Number (BIN) System assigned the 766 sequences to 15 BINs, including nine from Pakistan. Representative specimens of each Pakistan BIN were analyzed for mtCOI-3' to allow their assignment to one of the putative species in the B. tabaci complex recognized on the basis of sequence variation in this gene region. This analysis revealed the presence of Asia II 1, Middle East-Asia Minor 1, Asia 1, Asia II 5, Asia II 7, and a new lineage "Pakistan". The first two taxa were found in both Punjab and Sindh, but Asia 1 was only detected in Sindh, while Asia II 5, Asia II 7 and "Pakistan" were only present in Punjab. The haplotype networks showed that most haplotypes of Asia II 1, a species implicated in transmission of the cotton leaf curl virus, occurred in both India and Pakistan. CONCLUSIONS: DNA barcodes successfully discriminated cryptic species in B. tabaci complex. The dominant haplotypes in the B. tabaci complex were shared by India and Pakistan. Asia II 1 was previously restricted to Punjab, but is now the dominant lineage in southern Sindh; its southward spread may have serious implications for cotton plantations in this region.

Analyzing mosquito (Diptera: culicidae) diversity in Pakistan by DNA barcoding.

BACKGROUND: Although they are important disease vectors mosquito biodiversity in Pakistan is poorly known. Recent epidemics of dengue fever have revealed the need for more detailed understanding of the diversity and distributions of mosquito species in this region. DNA barcoding improves the accuracy of mosquito inventories because morphological differences between many species are subtle, leading to misidentifications. METHODOLOGY/PRINCIPAL FINDINGS: Sequence variation in the barcode region of the mitochondrial COI gene was used to identify mosquito species, reveal genetic diversity, and map the distribution of the dengue-vector species in Pakistan. Analysis of 1684 mosquitoes from 491 sites in Punjab and Khyber Pakhtunkhwa during 2010-2013 revealed 32 species with the assemblage dominated by Culex quinquefasciatus (61% of the collection). The genus Aedes (Stegomyia) comprised 15% of the specimens, and was represented by six taxa with the two dengue vector species, Ae. albopictus and Ae. aegypti, dominant and broadly distributed. Anopheles made up another 6% of the catch with An. subpictus dominating. Barcode sequence divergence in conspecific specimens ranged from 0-2.4%, while congeneric species showed from 2.3-17.8% divergence. A global haplotype analysis of disease-vectors showed the presence of multiple haplotypes, although a single haplotype of each dengue-vector species was dominant in most countries. Geographic distribution of Ae. aegypti and Ae. albopictus showed the later species was dominant and found in both rural and urban environments. CONCLUSIONS: As the first DNA-based analysis of mosquitoes in Pakistan, this study has begun the construction of a barcode reference library for the mosquitoes of this region. Levels of genetic diversity varied among species. Because of its capacity to differentiate species, even those with subtle morphological differences, DNA barcoding aids accurate tracking of vector populations.

Effects of Azospirillum brasilense with genetically modified auxin biosynthesis gene ipdC upon the diversity of the indigenous microbiota of the wheat rhizosphere.

The phytostimulatory properties of Azospirillum inoculants, which entail production of the phytohormone indole-3-acetic acid (IAA), can be enhanced by genetic means. However, it is not known whether this could affect their interactions with indigenous soil microbes. Here, wheat seeds were inoculated with the wild-type strain Azospirillum brasilense Sp245 or one of three genetically modified (GM) derivatives and grown for one month. The GM derivatives contained a plasmid vector harboring the indole-3-pyruvate/phenylpyruvate decarboxylase gene ipdC (IAA production) controlled either by the constitutive promoter PnptII or the root exudate-responsive promoter PsbpA, or by an empty vector (GM control). All inoculants displayed equal rhizosphere population densities. Only inoculation with either ipdC construct increased shoot biomass compared with the non-inoculated control. At one month after inoculation, automated ribosomal intergenic spacer analysis (ARISA) revealed that the effect of the PsbpA construct on bacterial community structure differed from that of the GM control, which was confirmed by 16S rDNA-based denaturing gradient gel electrophoresis (DGGE). The fungal community was sensitive to inoculation with the PsbpA construct and especially the GM control, based on ARISA data. Overall, fungal and bacterial communities displayed distinct responses to inoculation of GM A. brasilense phytostimulators, whose effects could differ from those of the wild-type.