Madhav P. Nepal

Identification, Nomenclature, and Evolutionary Relationships of Mitogen-Activated Protein Kinase (MAPK) Genes in Soybean

Mitogen-activated protein kinase (MAPK) genes in eukaryotes regulate various developmental and physiological processes including those associated with biotic and abiotic stresses. Although MAPKs in some plant species including Arabidopsis have been identified, they are yet to be identified in soybean. Major objectives of this study were to identify GmMAPKs, assess their evolutionary relationships, and analyze their functional divergence. We identified a total of 38 MAPKs, eleven MAPKKs, and 150 MAPKKKs in soybean. Within the GmMAPK family, we also identified a new clade of six genes: four genes with TEY and two genes with TQY motifs requiring further investigation into possible legume-specific functions. The results indicated the expansion of the GmMAPK families attributable to the ancestral polyploidy events followed by chromosomal rearrangements. The GmMAPK and GmMAPKKK families were substantially larger than those in other plant species. The duplicated GmMAPK members presented complex evolutionary relationships and functional divergence when compared to their counterparts in Arabidopsis. We also highlighted existing nomenclatural issues, stressing the need for nomenclatural consistency. GmMAPK identification is vital to soybean crop improvement, and novel insights into the evolutionary relationships will enhance our understanding about plant genome evolution.

Phylogenetics of Morus (Moraceae) Inferred from ITS and trnL-trnF Sequence Data

Abstract Morus (Tribe Moreae, Moraceae) consists of ca. 13 species of trees distributed in Asia, Africa, Europe, and North, Central, and South America. The broad geographical distribution of the genus, overlapping ranges of many taxa, and documented hybridization between some species present interesting questions of taxonomy, phylogeny, and biogeography. Phylogenetic data for Morus also contribute to higher level taxonomic work in the family. We used sequence data from ITS of the nrDNA and the chloroplast trnL-trnF intergenic spacer to study phylogenetic relationships of Morus. Phylogenies based on separate data sets were not statistically incongruent, and the combined tree reveals that Morus, as currently circumscribed, is non-monophyletic. Subgenus Morus (sensu Leroy) is resolved as a clade and consists of two well-supported clades: one of Asian taxa and one of North American taxa. Sampled members of the genus Trophis (two, including the type) form a clade sister to subgenus Morus. Morus mesozygia (Africa; subgenus Afromorus) and M. insignis (Neotropics; subgenus Gomphomorus), which have not been included to date in other phylogenetic studies of the family, are placed outside the subgenus Morus-Trophis clade. This work is an important step in elucidating relationships of Morus and along with other recent phylogenetic studies in Moraceae, underscores the need for further work within Tribe Moreae to clarify natural generic relationships.

CNL Disease Resistance Genes in Soybean and Their Evolutionary Divergence

Disease resistance genes (R-genes) encode proteins involved in detecting pathogen attack and activating downstream defense molecules. Recent availability of soybean genome sequences makes it possible to examine the diversity of gene families including disease-resistant genes. The objectives of this study were to identify coiled-coil NBS-LRR (= CNL) R-genes in soybean, infer their evolutionary relationships, and assess structural as well as functional divergence of the R-genes. Profile hidden Markov models were used for sequence identification and model-based maximum likelihood was used for phylogenetic analysis, and variation in chromosomal positioning, gene clustering, and functional divergence were assessed. We identified 188 soybean CNL genes nested into four clades consistent to their orthologs in Arabidopsis. Gene clustering analysis revealed the presence of 41 gene clusters located on 13 different chromosomes. Analyses of the Ks-values and chromosomal positioning suggest duplication events occurring at varying timescales, and an extrapericentromeric positioning may have facilitated their rapid evolution. Each of the four CNL clades exhibited distinct patterns of gene expression. Phylogenetic analysis further supported the extrapericentromeric positioning effect on the divergence and retention of the CNL genes. The results are important for understanding the diversity and divergence of CNL genes in soybean, which would have implication in soybean crop improvement in future.

Diversity and Evolution of Disease Resistance Genes in Barley (Hordeum vulgare L.)

Plant disease resistance genes (R-genes) play a critical role in the defense response to pathogens. Barley is one of the most important cereal crops, having a genome recently made available, for which the diversity and evolution of R-genes are not well understood. The main objectives of this research were to conduct a genome-wide identification of barley Coiled-coil, Nucleotide-binding site, Leucine-rich repeat (CNL) genes and elucidate their evolutionary history. We employed a Hidden Markov Model using 52 Arabidopsis thaliana CNL reference sequences and analyzed for phylogenetic relationships, structural variation, and gene clustering. We identified 175 barley CNL genes nested into three clades, showing (a) evidence of an expansion of the CNL-C clade, primarily due to tandem duplications; (b) very few members of clade CNL-A and CNL-B; and (c) a complete absence of clade CNL-D. Our results also showed that several of the previously identified mildew locus A (MLA) genes may be allelic variants of two barley CNL genes, MLOC_66581 and MLOC_10425, which respond to powdery mildew. Approximately 23% of the barley CNL genes formed 15 gene clusters located in the extra-pericentromeric regions on six of the seven chromosomes; more than half of the clustered genes were located on chromosomes 1H and 7H. Higher average numbers of exons and multiple splice variants in barley relative to those in Arabidopsis and rice may have contributed to a diversification of the CNL-C members. These results will help us understand the evolution of R-genes with potential implications for developing durable resistance in barley cultivars.

Distribution of Diverse Escherichia coli between Cattle and Pasture

Escherichia coli is widely considered to not survive for extended periods outside the intestines of warm-blooded animals; however, recent studies demonstrated that E. coli strains maintain populations in soil and water without any known fecal contamination. The objective of this study was to investigate whether the niche partitioning of E. coli occurs between cattle and their pasture. We attempted to clarify whether E. coli from bovine feces differs phenotypically and genotypically from isolates maintaining a population in pasture soil over winter. Soil, bovine fecal, and run-off samples were collected before and after the introduction of cattle to the pasture. Isolates (363) were genotyped by uidA and mutS sequences and phylogrouping, and evaluated for curli formation (Rough, Dry, And Red, or RDAR). Three types of clusters emerged, viz. bovine-associated, clusters devoid of cattle isolates and representing isolates endemic to the pasture environment, and clusters with both. All isolates clustered with strains of E. coli sensu stricto, distinct from the cryptic species Clades I, III, IV, and V. Pasture soil endemic and bovine fecal populations had very different phylogroup distributions, indicating niche partitioning. The soil endemic population was largely comprised of phylogroup B1 and had a higher average RDAR score than other isolates. These results indicate the existence of environmental E. coli strains that are phylogenetically distinct from bovine fecal isolates, and that have the ability to maintain populations in the soil environment.

Disease Resistance Mechanisms in Plants

Plants have developed a complex defense system against diverse pests and pathogens. Once pathogens overcome mechanical barriers to infection, plant receptors initiate signaling pathways driving the expression of defense response genes. Plant immune systems rely on their ability to recognize enemy molecules, carry out signal transduction, and respond defensively through pathways involving many genes and their products. Pathogens actively attempt to evade and interfere with response pathways, selecting for a decentralized, multicomponent immune system. Recent advances in molecular techniques have greatly expanded our understanding of plant immunity, largely driven by potential application to agricultural systems. Here, we review the major plant immune system components, state of the art knowledge, and future direction of research on plant–pathogen interactions. In our review, we will discuss how the decentralization of plant immune systems have provided both increased evolutionary opportunity for pathogen resistance, as well as additional mechanisms for pathogen inhibition of such defense responses. We conclude that the rapid advances in bioinformatics and molecular biology are driving an explosion of information that will advance agricultural production and illustrate how complex molecular interactions evolve.

Evolutionary Divergence of TNL Disease-Resistant Proteins in Soybean (Glycine max) and Common Bean (Phaseolus vulgaris)

Disease-resistant genes (R genes) encode proteins that are involved in protecting plants from their pathogens and pests. Availability of complete genome sequences from soybean and common bean allowed us to perform a genome-wide identification and analysis of the Toll interleukin-1 receptor-like nucleotide-binding site leucine-rich repeat (TNL) proteins. Hidden Markov model (HMM) profiling of all protein sequences resulted in the identification of 117 and 77 regular TNL genes in soybean and common bean, respectively. We also identified TNL gene homologs with unique domains, and signal peptides as well as nuclear localization signals. The TNL genes in soybean formed 28 clusters located on 10 of the 20 chromosomes, with the majority found on chromosome 3, 6 and 16. Similarly, the TNL genes in common bean formed 14 clusters located on five of the 11 chromosomes, with the majority found on chromosome 10. Phylogenetic analyses of the TNL genes from Arabidopsis, soybean and common bean revealed less divergence within legumes relative to the divergence between legumes and Arabidopsis. Syntenic blocks were found between chromosomes Pv10 and Gm03, Pv07 and Gm10, as well as Pv01 and Gm14. The gene expression data revealed basal level expression and tissue specificity, while analysis of available microRNA data showed 37 predicted microRNA families involved in targeting the identified TNL genes in soybean and common bean.

Comparative Genomics of Non-TNL Disease Resistance Genes from Six Plant Species

Disease resistance genes (R genes), as part of the plant defense system, have coevolved with corresponding pathogen molecules. The main objectives of this project were to identify non-Toll interleukin receptor, nucleotide-binding site, leucine-rich repeat (nTNL) genes and elucidate their evolutionary divergence across six plant genomes. Using reference sequences from Arabidopsis, we investigated nTNL orthologs in the genomes of common bean, Medicago, soybean, poplar, and rice. We used Hidden Markov Models for sequence identification, performed model-based phylogenetic analyses, visualized chromosomal positioning, inferred gene clustering, and assessed gene expression profiles. We analyzed 908 nTNL R genes in the genomes of the six plant species, and classified them into 12 subgroups based on the presence of coiled-coil (CC), nucleotide binding site (NBS), leucine rich repeat (LRR), resistance to Powdery mildew 8 (RPW8), and BED type zinc finger domains. Traditionally classified CC-NBS-LRR (CNL) genes were nested into four clades (CNL A-D) often with abundant, well-supported homogeneous subclades of Type-II R genes. CNL-D members were absent in rice, indicating a unique R gene retention pattern in the rice genome. Genomes from Arabidopsis, common bean, poplar and soybean had one chromosome without any CNL R genes. Medicago and Arabidopsis had the highest and lowest number of gene clusters, respectively. Gene expression analyses suggested unique patterns of expression for each of the CNL clades. Differential gene expression patterns of the nTNL genes were often found to correlate with number of introns and GC content, suggesting structural and functional divergence.

Genome-Wide Identification of NBS-Encoding Resistance Genes in Sunflower (Helianthus annuus L.)

Nucleotide Binding Site—Leucine-Rich Repeat (NBS-LRR) genes encode disease resistance proteins involved in plants’ defense against their pathogens. Although sunflower is affected by many diseases, only a few molecular details have been uncovered regarding pathogenesis and resistance mechanisms. Recent availability of sunflower whole genome sequences in publicly accessible databases allowed us to accomplish a genome-wide identification of Toll-interleukin-1 receptor-like Nucleotide-binding site Leucine-rich repeat (TNL), Coiled Coil (CC)-NBS-LRR (CNL), Resistance to powdery mildew8 (RPW8)-NBS-LRR (RNL) and NBS-LRR (NL) protein encoding genes. Hidden Markov Model (HMM) profiling of 52,243 putative protein sequences from sunflower resulted in 352 NBS-encoding genes, among which 100 genes belong to CNL group including 64 genes with RX_CC like domain, 77 to TNL, 13 to RNL, and 162 belong to NL group. We also identified signal peptides and nuclear localization signals present in the identified genes and their homologs. We found that NBS genes were located on all chromosomes and formed 75 gene clusters, one-third of which were located on chromosome 13. Phylogenetic analyses between sunflower and Arabidopsis NBS genes revealed a clade-specific nesting pattern in CNLs, with RNLs nested in the CNL-A clade, and species-specific nesting pattern for TNLs. Surprisingly, we found a moderate bootstrap support (BS = 50%) for CNL-A clade being nested within TNL clade making both the CNL and TNL clades paraphyletic. Arabidopsis and sunflower showed 87 syntenic blocks with 1049 high synteny hits between chromosome 5 of Arabidopsis and chromosome 6 of sunflower. Expression data revealed functional divergence of the NBS genes with basal level tissue-specific expression. This study represents the first genome-wide identification of NBS genes in sunflower paving avenues for functional characterization and potential crop improvement.

Data on the genome-wide identification of CNL R-genes in Setaria italica (L.) P. Beauv.

We report data associated with the identification of 242 disease resistance genes (R-genes) in the genome of Setaria italica as presented in “Genetic diversity of disease resistance genes in foxtail millet (Setaria italica L.)” (Andersen and Nepal, 2017) [1]. Our data describe the structure and evolution of the Coiled-coil, Nucleotide-binding site, Leucine-rich repeat (CNL) R-genes in foxtail millet. The CNL genes were identified through rigorous extraction and analysis of recently available plant genome sequences using cutting-edge analytical software. Data visualization includes gene structure diagrams, chromosomal syntenic maps, a chromosomal density plot, and a maximum-likelihood phylogenetic tree comparing Sorghum bicolor, Panicum virgatum, Setaria italica, and Arabidopsis thaliana. Compilation of InterProScan annotations, Gene Ontology (GO) annotations, and Basic Local Alignment Search Tool (BLAST) results for the 242 R-genes identified in the foxtail millet genome are also included in tabular format.