Background: The MYB gene family comprises one of the richest groups of transcription factors in plants. Plant MYB proteins are characterized by a highly conserved MYB DNA-binding domain. MYB proteins are classified into four major groups namely, 1R-MYB, 2R-MYB, 3R-MYB and 4R-MYB based on the number and position of MYB repeats. MYB transcription factors are involved in plant development, secondary metabolism, hormone signal transduction, disease resistance and abiotic stress tolerance. A comparative analysis of MYB family genes in rice and Arabidopsis will help reveal the evolution and function of MYB genes in plants. Results: A genome-wide analysis identified at least 155 and 197 MYB genes in rice and Arabidopsis, respectively. Gene structure analysis revealed that MYB family genes possess relatively more number of introns in the middle as compared with C- and N-terminal regions of the predicted genes. Intronless MYB-genes are highly conserved both in rice and Arabidopsis. MYB genes encoding R2R3 repeat MYB proteins retained conserved gene structure with three exons and two introns, whereas genes encoding R1R2R3 repeat containing proteins consist of six exons and five introns. The splicing pattern is similar among R1R2R3 MYB genes in Arabidopsis. In contrast, variation in splicing pattern was observed among R1R2R3 MYB members of rice. Consensus motif analysis of 1kb upstream region (5′ to translation initiation codon) of MYB gene ORFs led to the identification of conserved and over-represented cis-motifs in both rice and Arabidopsis. Real-time quantitative RT-PCR analysis showed that several members of MYBs are up-regulated by various abiotic stresses both in rice and Arabidopsis. Conclusion: A comprehensive genome-wide analysis of chromosomal distribution, tandem repeats and phylogenetic relationship of MYB family genes in rice and Arabidopsis suggested their evolution via duplication. Genome-wide comparative analysis of MYB genes and their expression analysis identified several MYBs with potential role in development and stress response of plants.

Small non-coding RNAs (sRNAs) namely microRNAs (miRNAs) and trans-acting small interfering RNAs (tasi-RNAs) play a crucial role in post-transcriptional regulation of gene expression and thus the control plant development and stress responses. In order to identify drought-responsive miRNAs and tasi-RNAs in sorghum, we constructed small RNA libraries from a drought tolerant (M35-1) and susceptible (C43) sorghum genotypes grown under control and drought stress conditions, and sequenced by Illumina Genome Analyzer IIx. Ninety seven conserved and 526 novel miRNAs representing 472 unique miRNA families were identified from sorghum. Ninety six unique miRNAs were found to be regulated by drought stress, of which 32 were up- and 49 were down-regulated (fold change ≥ 2 or ≤ -2) at least in one genotype, while the remaining 15 miRNAs showed contrasting drought-regulated expression pattern between genotypes. A maximum of 17 and 18 miRNAs was differentially regulated under drought stress condition in the sensitive and tolerant genotypes, respectively. These results suggest that genotype dependent stress responsive regulation of miRNAs may contribute, at least in part, to the differential drought tolerance of sorghum genotypes. We also identified two miR390-directed TAS3 gene homologs and the auxin response factors as tasi-RNA targets. We predicted more than 1,300 unique target genes for the novel and conserved miRNAs. These target genes were predicted to be involved in different cellular, metabolic, response to stimulus, biological regulation, and developmental processes. Genome-wide identification of stress-responsive miRNAs, tasi-RNAs and their targets identified in this study will be useful in unravelling the molecular mechanisms underlying drought stress responses and genetic improvement of biomass production and stress tolerance in sorghum.

MicroRNAs (miRNAs) regulate gene expression mainly by post-transcriptional gene silencing (PTGS) and in some cases by transcriptional genes silencing (TGS). miRNAs play critical roles in developmental processes, nutrient homeostasis, abiotic stress and pathogen responses of plants. In contrast to the large number of miRNAs predicted in cereal model plant rice, only 148 miRNAs were predicted in sorghum till date (miRBase release 17). This suggested that miRNAs identified in sorghum is far from saturation. Hence, we developed a bioinformatics pipeline using an in-house PERL script and publicly available structure prediction tools to identify miRNAs and their target genes from publically available Expressed Sequence Tags (EST) and Genomic Survey Sequence (GSS). About 1,379 known and unique plant miRNAs from 33 different crops were used to predict new miRNAs in sorghum. We identified 31 new miRNAs belonging to 10 different miRNA families. We predicted 72 potential target genes for 31 miRNAs, and most of these target genes are predicted to be involved in plant growth and development. These newly identified miRNAs add to the growing database of miRNA and lay the foundation for further understanding of miRNA function in sorghum plant development.

The GNE-associated V727M mutation is one of the most prevalent ethnic founder mutations in the Asian HIBM cohort; however, its role in inducing disease phenotype remains largely elusive. In this study, the function of this hotspot mutation was profoundly investigated. For this, V727M mutation-specific altered expression profile and potential networks were explored. The relevant muscular disorder-specific in vivo studies and patient data were further analyzed, and the key altered molecular pathways were identified. Our study found that the GNEV727M mutation resulted in a deregulated lincRNA profile, the majority of which (91%) were associated with a down-regulation trend. Further, in silico analysis of associated targets showed their active role in regulating Wnt, TGF-β, and apoptotic signaling. Interestingly, COL6a3 was found as a key target of these lincRNAs. Further, GSEA analysis showed HIBM patients with variable COL6A3 transcript levels have significant alteration in many critical pathways, including epithelial-mesenchymal-transition, myogenesis, and apoptotic signaling. Interestingly, 12 of the COL6A3 coexpressed genes also showed a similar altered expression profile in HIBM. A similar altered trend in COL6A3 and coexpressed genes were found in in vivo HIBM disease models as well as in multiple other skeletal disorders. Thus, the COL6A3-specific 13 gene signature seems to be altered in multiple muscular disorders. Such deregulation could play a pivotal role in regulating many critical processes such as extracellular matrix organization, cell adhesion, and skeletal muscle development. Thus, investigating this novel COL6A3-specific 13 gene signature provides valuable information for understanding the molecular cause of HIBM and may also pave the way for better diagnosis and effective therapeutic strategies for many muscular disorders.