MicroRNA (miRNA) profiling of maize genotypes with differential response to Aspergillus flavus implies zma-miR156-squamosa promoter binding protein (SBP) and zma-miR398/zma-miR394-F -box combinations involved in resistance mechanisms.

Maize (Zea mays), a major food crop worldwide, is susceptible to infection by the saprophytic fungus Aspergillus flavus that can produce the carcinogenic metabolite aflatoxin (AF) especially under climate change induced abiotic stressors that favor mold growth. Several studies have used "-omics" approaches to identify genetic elements with potential roles in AF resistance, but there is a lack of research identifying the involvement of small RNAs such as microRNAs (miRNAs) in maize-A. flavus interaction. In this study, we compared the miRNA profiles of three maize lines (resistant TZAR102, moderately resistant MI82, and susceptible Va35) at 8 h, 3 d, and 7 d after A. flavus infection to investigate possible regulatory antifungal role of miRNAs. A total of 316 miRNAs (275 known and 41 putative novel) belonging to 115 miRNA families were identified in response to the fungal infection across all three maize lines. Eighty-two unique miRNAs were significantly differentially expressed with 39 miRNAs exhibiting temporal differential regulation irrespective of the maize genotype, which targeted 544 genes (mRNAs) involved in diverse molecular functions. The two most notable biological processes involved in plant immunity, namely cellular responses to oxidative stress (GO:00345990) and reactive oxygen species (GO:0034614) were significantly enriched in the resistant line TZAR102. Coexpression network analysis identified 34 hubs of miRNA-mRNA pairs where nine hubs had a node in the module connected to their target gene with potentially important roles in resistance/susceptible response of maize to A. flavus. The miRNA hubs in resistance modules (TZAR102 and MI82) were mostly connected to transcription factors and protein kinases. Specifically, the module of miRNA zma-miR156b-nb - squamosa promoter binding protein (SBP), zma-miR398a-3p - SKIP5, and zma-miR394a-5p - F-box protein 6 combinations in the resistance-associated modules were considered important candidates for future functional studies.

Exploring the sorghum race level diversity utilizing 272 sorghum accessions genomic resources.

Due to evolutionary divergence, sorghum race populations exhibit significant genetic and morphological variation. A k-mer-based sorghum race sequence comparison identified the conserved k-mers of all 272 accessions from sorghum and the race-specific genetic signatures identified the gene variability in 10,321 genes (PAVs). To understand sorghum race structure, diversity and domestication, a deep learning-based variant calling approach was employed in a set of genotypic data derived from a diverse panel of 272 sorghum accessions. The data resulted in 1.7 million high-quality genome-wide SNPs and identified selective signature (both positive and negative) regions through a genome-wide scan with different (iHS and XP-EHH) statistical methods. We discovered 2,370 genes associated with selection signatures including 179 selective sweep regions distributed over 10 chromosomes. Co-localization of these regions undergoing selective pressure with previously reported QTLs and genes revealed that the signatures of selection could be related to the domestication of important agronomic traits such as biomass and plant height. The developed k-mer signatures will be useful in the future to identify the sorghum race and for trait and SNP markers for assisting in plant breeding programs.

Construction of Practical Haplotype Graph (PHG) with the Whole-Genome Sequence Data.

With the emerging sequencing technologies and cost reduction, the sequence data generation has accelerated from a single individual to multiple (thousands of) individuals of a species. The terabytes of sequence data generated from thousands of individuals include the majority of the redundant sequence which depends on the level of sequence similarity within the population of individuals. Managing large datasets and creating the unique catalogue sequence from such a large population is challenging to analyze, store, and retrieve the information. In this chapter, we discuss the practical haplotype graph (PHG) which addresses the above said challenges and also able to retrieve required information such as variants and sequences more efficiently, which enable researchers to manage and assess large genomic data.

Sorghum Pan-Genome Explores the Functional Utility for Genomic-Assisted Breeding to Accelerate the Genetic Gain.

Sorghum (Sorghum bicolor L.) is a staple food crops in the arid and rainfed production ecologies. Sorghum plays a critical role in resilient farming and is projected as a smart crop to overcome the food and nutritional insecurity in the developing world. The development and characterisation of the sorghum pan-genome will provide insight into genome diversity and functionality, supporting sorghum improvement. We built a sorghum pan-genome using reference genomes as well as 354 genetically diverse sorghum accessions belonging to different races. We explored the structural and functional characteristics of the pan-genome and explain its utility in supporting genetic gain. The newly-developed pan-genome has a total of 35,719 genes, a core genome of 16,821 genes and an average of 32,795 genes in each cultivar. The variable genes are enriched with environment responsive genes and classify the sorghum accessions according to their race. We show that 53% of genes display presence-absence variation, and some of these variable genes are predicted to be functionally associated with drought adaptation traits. Using more than two million SNPs from the pan-genome, association analysis identified 398 SNPs significantly associated with important agronomic traits, of which, 92 were in genes. Drought gene expression analysis identified 1,788 genes that are functionally linked to different conditions, of which 79 were absent from the reference genome assembly. This study provides comprehensive genomic diversity resources in sorghum which can be used in genome assisted crop improvement.

Nitrogen Use Efficiency in Sorghum: Exploring Native Variability for Traits Under Variable N-Regimes.

Exploring the natural genetic variability and its exploitation for improved Nitrogen Use Efficiency (NUE) in sorghum is one of the primary goals in the modern crop improvement programs. The integrated strategies include high-throughput phenotyping, next generation sequencing (NGS)-based genotyping technologies, and a priori selected candidate gene studies that help understand the detailed physiological and molecular mechanisms underpinning this complex trait. A set of sixty diverse sorghum genotypes was evaluated for different vegetative, reproductive, and yield traits related to NUE in the field (under three N regimes) for two seasons. Significant variations for different yield and related traits under 0 and 50% N confirmed the availability of native genetic variability in sorghum under low N regimes. Sorghum genotypes with distinct genetic background had interestingly similar NUE associated traits. The Genotyping-By-Sequencing based SNPs (>89 K) were used to study the population structure, and phylogenetic groupings identified three distinct groups. The information of grain N and stalk N content of the individuals covered on the phylogenetic groups indicated randomness in the distribution for adaptation under variable N regimes. This study identified promising sorghum genotypes with consistent performance under varying environments, with buffer capacity for yield under low N conditions. We also report better performing genotypes for varied production use-grain, stover, and dual-purpose sorghum having differential adaptation response to NUE traits. Expression profiling of NUE associated genes in shoot and root tissues of contrasting lines (PVK801 and HDW703) grown in varying N conditions revealed interesting outcomes. Root tissues of contrasting lines exhibited differential expression profiles for transporter genes [ammonium transporter (SbAMT), nitrate transporters (SbNRT)]; primary assimilatory (glutamine synthetase (SbGS), glutamate synthase (SbGOGAT[NADH], SbGOGAT[Fd]), assimilatory genes [nitrite reductase (SbNiR[NADH]3)]; and amino acid biosynthesis associated gene [glutamate dehydrogenase (SbGDH)]. Identification and expression profiling of contrasting sorghum genotypes in varying N dosages will provide new information to understand the response of NUE genes toward adaptation to the differential N regimes in sorghum. High NUE genotypes identified from this study could be potential candidates for in-depth molecular analysis and contribute toward the development of N efficient sorghum cultivars.

Genome-wide comparative transcriptome analysis of the A4-CMS line ICPA 2043 and its maintainer ICPB 2043 during the floral bud development of pigeonpea.

Cytoplasmic male sterility (CMS) offers a unique system to understand cytoplasmic nuclear crosstalk, and is also employed for exploitation of hybrid vigor in various crops. Pigeonpea A4-CMS, a predominant source of male sterility, is being used for efficient hybrid seed production. The molecular mechanisms of CMS trait remain poorly studied in pigeonpea. We performed genome-wide transcriptome profiling of A4-CMS line ICPA 2043 and its isogenic maintainer ICPB 2043 at two different stages of floral bud development (stage S1 and stage S2). Consistent with the evidences from some other crops, we also observed significant difference in the expression levels of genes in the later stage, i.e., stage S2. Differential expression was observed for 143 and 55 genes within the two stages of ICPA 2043 and ICPB 2043, respectively. We obtained only 10 differentially expressed genes (DEGs) between the stage S1 of the two genotypes, whereas expression change was significant for 582 genes in the case of stage S2. The qRT-PCR assay of randomly selected six genes supported the differential expression of genes between ICPA 2043 and ICPB 2043. Further, GO and KEGG pathway mapping suggested a possible compromise in key bioprocesses during flower and pollen development. Besides providing novel insights into the functional genomics of CMS trait, our results were in strong agreement with the gene expression atlas of pigeonpea that implicated various candidate genes like sucrose-proton symporter 2 and an uncharacterized protein along with pectate lyase, pectinesterase inhibitors, L-ascorbate oxidase homolog, ATPase, ß-galactosidase, polygalacturonase, and aldose 1-epimerase for pollen development of pigeonpea. The dataset presented here provides a rich genomic resource to improve understanding of CMS trait and its deployment in heterosis breeding in pigeonpea.

Identification of microRNAs and their gene targets in cytoplasmic male sterile and fertile maintainer lines of pigeonpea.

MAIN CONCLUSION: Comparative analysis of genome-wide miRNAs and their gene targets between cytoplasmic male sterile (CMS) and fertile lines of pigeonpea suggests a possible role of miRNA-regulated pathways in reproductive development. Exploitation of hybrid vigor using CMS technology has delivered nearly 50% yield gain in pigeonpea. Among various sterility-inducing cytoplasms (A1-A9) reported so far in pigeonpea, A2 and A4 are the two major sources that facilitate hybrid seed production. Recent evidence suggests involvement of micro RNA in vast array of biological processes including plant reproductive development. In pigeonpea, information about the miRNAs is insufficient. In view of this, we sequenced six small RNA libraries of CMS line UPAS 120A and isogenic fertile line UPAS 120B using Illumina technology. Results revealed 316 miRNAs including 248 known and 68 novel types. A total of 637 gene targets were predicted for known miRNAs, while 324 genes were associated with novel miRNAs. Degradome analysis revealed 77 gene targets of predicted miRNAs, which included a variety of transcription factors playing key roles in plant reproduction such as F-box family proteins, apetala 2, auxin response factors, ethylene-responsive factors, homeodomain-leucine zipper proteins etc. Differential expression of both known and novel miRNAs implied roles for both conserved as well as species-specific players. We also obtained several miRNA families such as miR156, miR159, miR167 that are known to influence crucial aspects of plant fertility. Gene ontology and pathway level analyses of the target genes showed their possible implications for crucial events during male reproductive development such as tapetal degeneration, pollen wall formation, retrograde signaling etc. To the best of our knowledge, present study is first to combine deep sequencing of small RNA and degradome for elucidating the role of miRNAs in flower and male reproductive development in pigeonpea.

Genome-Wide Association Study for Major Biofuel Traits in Sorghum Using Minicore Collection.

BACKGROUND: Production of biofuels from lignocellulosic crop biomass is an alternative to reduce greenhouse gas emissions. The biofuel production involves collecting biomass, breaking down cell wall components followed by the conversion of sugars to ethanol. The lingo-cellulosic biomass comprises 40-50% cellulose, 20-30% hemicellulose, and 10-25% lignin. Sorghum is a widely adapted energy crop for biofuel production. Biomass with low lignin, high cellulose, and high hemicellulose contents are exploited to attain maximum biofuel production efficiency. Resistance to lodging, pest, disease, and abiotic stresses related to cell wall components is well documented, and quantitative trait loci were identified to understand these traits' genetic correlation. Selection for reduced lignin and increased cellulose content in stover can increase the ethanol yield. The Genome-Wide Association Studies (GWAS) is a complementary approach to evaluating the marker and phenotype associations among large diversity panels. Single nucleotide polymorphisms were scanned to identify loci associated with the traits of interest. In this study, the GWAS was performed on 245 sorghum minicore genotypes to analyze agronomic traits (days to 50%flowering, fresh biomass yield, dry biomass yield) and cell wall components (cellulose, hemicellulose, and lignin). Further, in-silico validation of the candidate genes was performed in a global gene expression data from large-scale RNA sequencing studies in sorghum available in the NCBI GEO database was used. OBJECTIVE: The objectives of this study are to evaluate native variations in biofuel related agronomic traits and stalk cell wall components and to identify significant SNPs or loci related to the cell wall components. METHODS: In this article, an association mapping panel, comprising of 245 sorghum minicore germplasm accessions, was evaluated during two post rainy seasons of 2013 and 2014, and observations were recorded on the whole plot- for days to 50% flowering, fresh biomass yield (tha-1, and dry biomass yield (tha-1). The biomass of sun-dried plants from both seasons was collected separately, chopped, dried, and ground to powder. The cellulose, hemicellulose, and lignin contents were determined in the powdered. The content of each of these three components in sorghum was expressed in percent of dry matter. The data on agronomic traits and composition analysis was subjected to Analysis of Variance. For the current study, we remapped the raw GBS data with the sorghum assembly version v3.1. A total of 27,589 SNPs were obtained with a minor allele frequency (MAF) >1% and missing data <50%. The GWAS was performed in a single minicore population using FarmCPU, in R software. The synteny positions of the identified significant SNPs between sorghum and other model crop species viz., maize, switchgrass, and Arabidopsis were represented using CIRCOS software for traits viz., dry biomass yield, cellulose, hemicellulose, and lignin. The transcriptome dataset from where sorghum gene atlas studies of grain, sweet, and bioenergy sorghums are available through NCBI's Gene Expression Omnibus (GEO) under accession number GSE49879, was used to cross-validate the identified SNPs for cellulose, hemicellulose, and lignin through GWAS. RESULTS: High broad-sense heritability was exhibited for all the traits in individual seasons along with significant genotype × environment interaction across seasons except lignin. Association mapping with a P < 1×10-4 revealed genomic regions associated with the- (i) agronomic traits (days to 50% flowering, fresh and dry biomass), and (ii) biochemical traits (cellulose, hemicellulose, and lignin) associated with biofuels production, in individual seasons. Twelve significant SNPs for flowering time, 30 fresh biomass yields, and 24 for dry biomass yield, 25 for cellulose, 7 for hemicellulose, and 21 for lignin were identified. CIRCOS plot was constructed to identify and analyze similarities and differences while comparing the sorghum genome with different crops. For cellulose high similarity of >80% was observed for all sorghum gene sequences with the maize homologs. The overall similarity of sorghum homologs with foxtail millet was >65%, for Arabidopsis from 30.6% to 48.6%, and rice from 28.2% to 92.8%. SNPs for hemicellulose displayed maximum similarity to foxtail millet followed by maize. The sequence similarity of lignin SNPs in sorghum was highest with the maize genome followed by Arabidopsis. Both rice and foxtail millet showed >55% similarity to the sorghum genome. CONCLUSION: This study reports large variability for agronomic and biofuel traits in the sorghum minicore collection with high heritability. The genetic architecture of cell wall components using the GWAS approach was studied and candidate genes for each component were annotated. These results give a better understanding of the genetic basis of the sorghum cell wall composition. The association analysis identified regions of the genome that could be targeted to enhance the quality of biomass and yield along with the desired composition promoting breeding efficiency for enhanced biofuel yield.

Apophysomyces variabilis: draft genome sequence and comparison of predictive virulence determinants with other medically important Mucorales.

BACKGROUND: Apophysomyces species are prevalent in tropical countries and A. variabilis is the second most frequent agent causing mucormycosis in India. Among Apophysomyces species, A. elegans, A. trapeziformis and A. variabilis are commonly incriminated in human infections. The genome sequences of A. elegans and A. trapeziformis are available in public database, but not A. variabilis. We, therefore, performed the whole genome sequence of A. variabilis to explore its genomic structure and possible genes determining the virulence of the organism. RESULTS: The whole genome of A. variabilis NCCPF 102052 was sequenced and the genomic structure of A. variabilis was compared with already available genome structures of A. elegans, A. trapeziformis and other medically important Mucorales. The total size of genome assembly of A. variabilis was 39.38 Mb with 12,764 protein-coding genes. The transposable elements (TEs) were low in Apophysomyces genome and the retrotransposon Ty3-gypsy was the common TE. Phylogenetically, Apophysomyces species were grouped closely with Phycomyces blakesleeanus. OrthoMCL analysis revealed 3025 orthologues proteins, which were common in those three pathogenic Apophysomyces species. Expansion of multiple gene families/duplication was observed in Apophysomyces genomes. Approximately 6% of Apophysomyces genes were predicted to be associated with virulence on PHIbase analysis. The virulence determinants included the protein families of CotH proteins (invasins), proteases, iron utilisation pathways, siderophores and signal transduction pathways. Serine proteases were the major group of proteases found in all Apophysomyces genomes. The carbohydrate active enzymes (CAZymes) constitute the majority of the secretory proteins. CONCLUSION: The present study is the maiden attempt to sequence and analyze the genomic structure of A. variabilis. Together with available genome sequence of A. elegans and A. trapeziformis, the study helped to indicate the possible virulence determinants of pathogenic Apophysomyces species. The presence of unique CAZymes in cell wall might be exploited in future for antifungal drug development.