Genome-wide identification and characterization of novel microRNAs in seed development of soybean

MicroRNAs (miRNAs) are important and ubiquitous regulators of gene expression in eukaryotes. However, the information about miRNAs population and their regulatory functions involving in soybean seed development remains incomplete. Base on the Dicer-like1-mediated cleavage signals during miRNA processing could be employed for novel miRNA discovery, a genome-wide search for miRNA candidates involved in seed development was carried out. As a result, 17 novel miRNAs, 14 isoforms of miRNA (isomiRs) and 31 previously validated miRNAs were discovered. These novel miRNAs and isomiRs represented tissue-specific expression and the isomiRs showed significantly higher abundance than that of their miRNA counterparts in different tissues. After target prediction and degradome sequencing data-based validation, 13 novel miRNA–target pairs were further identified. Besides, five targets of 22-nt iso-gma-miR393h were found to be triggered to produce secondary trans-acting siRNA (ta-siRNAs). Summarily, our results could expand the repertoire of miRNAs with potentially important functions in soybean.     

The evolution of plant microRNAs: insights from a basal eudicot sacred lotus

microRNAs (miRNAs) are important noncoding small RNAs that regulate mRNAs in eukaryotes. However, under which circumstances different miRNAs/miRNA families exhibit different evolutionary trajectories in plants remains unclear. In this study, we sequenced the small RNAs and degradome from a basal eudicot, sacred lotus (Nelumbo nucifera or lotus), to identify miRNAs and their targets. Combining with public miRNAs, we predicted 57 pre‐eudicot miRNA families from different evolutionary stages. We found that miRNA families featuring older age, higher copy and target number tend to show lower propensity for miRNA family loss (PGL) and stronger signature of purifying selection during divergence of temperate and tropical lotus. Further analyses of lotus genome revealed that there is an association between loss of miRNA families in descendent plants and in duplicated genomes. Gene dosage balance is crucial in maintaining those preferentially retained MIRNA duplicates by imposing stronger purifying selection. However, these factors and selection influencing miRNA family evolution are not applicable to the putative MIRNA‐likes. Additionally, the MIRNAs participating in lotus pollen–pistil interaction, a conserved process in angiosperms, also have a strong signature of purifying selection. Functionally, sequence divergence in MIRNAs escalates expression divergence of their target genes between temperate and tropical lotus during rhizome and leaf growth. Overall, our study unravels several important factors and selection that determine the miRNA family distribution in plants and duplicated genomes, and provides evidence for functional impact of MIRNA sequence evolution.     

Analyses of a Glycine max Degradome Library Identify microRNA Targets and MicroRNAs that Trigger Secondary SiRNA Biogenesis

Plant microRNAs (miRNAs) regulate gene expression mainly by guiding cleavage of target mRNAs. In this study, a degradome library constructed from different soybean (Glycine max (L.) Merr.) tissues was deep-sequenced. 428 potential targets of small interfering RNAs and 25 novel miRNA families were identified. A total of 211 potential miRNA targets, including 174 conserved miRNA targets and 37 soybean-specific miRNA targets, were identified. Among them, 121 targets were first discovered in soybean. The signature distribution of soybean primary miRNAs (pri-miRNAs) showed that most pri-miRNAs had the characteristic pattern of Dicer processing. The biogenesis of TAS3 small interfering RNAs (siRNAs) was conserved in soybean, and nine Auxin Response Factors were identified as TAS3 siRNA targets. Twenty-three miRNA targets produced secondary small interfering RNAs (siRNAs) in soybean. These targets were guided by five miRNAs: gma-miR393, gma-miR1508, gma-miR1510, gma-miR1514, and novel-11. Multiple targets of these secondary siRNAs were detected. These 23 miRNA targets may be the putative novel TAS genes in soybean. Global identification of miRNA targets and potential novel TAS genes will contribute to research on the functions of miRNAs in soybean.     

Construction and comparison of three reference‐quality genome assemblies for soybean

We report reference‐quality genome assemblies and annotations for two accessions of soybean (Glycine max) and for one accession of Glycine soja, the closest wild relative of G. max. The G. max assemblies provided are for widely used US cultivars: the northern line Williams 82 (Wm82) and the southern line Lee. The Wm82 assembly improves the prior published assembly, and the Lee and G. soja assemblies are new for these accessions. Comparisons among the three accessions show generally high structural conservation, but nucleotide difference of 1.7 single‐nucleotide polymorphisms (snps) per kb between Wm82 and Lee, and 4.7 snps per kb between these lines and G. soja. snp distributions and comparisons with genotypes of the Lee and Wm82 parents highlight patterns of introgression and haplotype structure. Comparisons against the US germplasm collection show placement of the sequenced accessions relative to global soybean diversity. Analysis of a pan‐gene collection shows generally high conservation, with variation occurring primarily in genomically clustered gene families. We found approximately 40–42 inversions per chromosome between either Lee or Wm82v4 and G. soja, and approximately 32 inversions per chromosome between Wm82 and Lee. We also investigated five domestication loci. For each locus, we found two different alleles with functional differences between G. soja and the two domesticated accessions. The genome assemblies for multiple cultivated accessions and for the closest wild ancestor of soybean provides a valuable set of resources for identifying causal variants that underlie traits for the domestication and improvement of soybean, serving as a basis for future research and crop improvement efforts for this important crop species.     

The role of Bh4 in parallel evolution of hull colour in domesticated and weedy rice

The two independent domestication events in the genus Oryza that led to African and Asian rice offer an extremely useful system for studying the genetic basis of parallel evolution. This system is also characterized by parallel de‐domestication events, with two genetically distinct weedy rice biotypes in the US derived from the Asian domesticate. One important trait that has been altered by rice domestication and de‐domestication is hull colour. The wild progenitors of the two cultivated rice species have predominantly black‐coloured hulls, as does one of the two U.S. weed biotypes; both cultivated species and one of the US weedy biotypes are characterized by straw‐coloured hulls. Using Black hull 4 (Bh4) as a hull colour candidate gene, we examined DNA sequence variation at this locus to study the parallel evolution of hull colour variation in the domesticated and weedy rice system. We find that independent Bh4‐coding mutations have arisen in African and Asian rice that are correlated with the straw hull phenotype, suggesting that the same gene is responsible for parallel trait evolution. For the U.S. weeds, Bh4 haplotype sequences support current hypotheses on the phylogenetic relationship between the two biotypes and domesticated Asian rice; straw hull weeds are most similar to indica crops, and black hull weeds are most similar to aus crops. Tests for selection indicate that Asian crops and straw hull weeds deviate from neutrality at this gene, suggesting possible selection on Bh4 during both rice domestication and de‐domestication.     

A draft genome of sweet cherry (Prunus avium L.) reveals genome‐wide and local effects of domestication

Sweet cherry (Prunus avium L.) trees are both economically important fruit crops but also important components of natural forest ecosystems in Europe, Asia and Africa. Wild and domesticated trees currently coexist in the same geographic areas with important questions arising on their historical relationships. Little is known about the effects of the domestication process on the evolution of the sweet cherry genome. We assembled and annotated the genome of the cultivated variety “Big Star*” and assessed the genetic diversity among 97 sweet cherry accessions representing three different stages in the domestication and breeding process (wild trees, landraces and modern varieties). The genetic diversity analysis revealed significant genome‐wide losses of variation among the three stages and supports a clear distinction between wild and domesticated trees, with only limited gene flow being detected between wild trees and domesticated landraces. We identified 11 domestication sweeps and five breeding sweeps covering, respectively, 11.0 and 2.4 Mb of the P. avium genome. A considerable fraction of the domestication sweeps overlaps with those detected in the related species, Prunus persica (peach), indicating that artificial selection during domestication may have acted independently on the same regions and genes in the two species. We detected 104 candidate genes in sweep regions involved in different processes, such as the determination of fruit texture, the regulation of flowering and fruit ripening and the resistance to pathogens. The signatures of selection identified will enable future evolutionary studies and provide a valuable resource for genetic improvement and conservation programs in sweet cherry.     

Comprehensive analyses of microRNA gene evolution in paleopolyploid soybean genome

miRNA genes are thought to undergo quick birth and death processes in genomes and the emergence of a MIRNA‐like hairpin provides the base for functional miRNA gene formation. However, the factors affecting the formation of an active miRNA gene from a MIRNA‐like hairpin within a genome remain unclear. We performed a genome‐wide investigation of MIRNA‐like hairpin accumulation, expression, structural changes and relationships with annotated genomic features in the paleopolyploid soybean genome. Our results showed that adjacent gene and transposable element content, rates of genetic recombination at location of emergence, along with its own gene structure divergence greatly affected miRNA gene evolution. Further investigation suggested that miRNA genes from different duplication sources followed distinct evolutionary trajectories and that the accumulation of MIRNA‐like hairpins might be an important factor in causing long terminal repeat retrotransposons to lose activity during genome evolution.     

Genome-wide signatures of the geographic expansion and breeding of soybean

Soybean is a leguminous crop that provides oil and protein. Exploring the genomic signatures of soybean evolution is crucial for breeding varieties with improved adaptability to environmental extremes. We analyzed the genome sequences of 2,214 soybeans and proposed a soybean evolutionary route, i.e., the expansion of annual wild soybean (Glycine soja Sieb. & Zucc.) from southern China and its domestication in central China, followed by the expansion and local breeding selection of its landraces (G. max (L.) Merr.). We observed that the genetic introgression in soybean landraces was mostly derived from sympatric rather than allopatric wild populations during the geographic expansion. Soybean expansion and breeding were accompanied by the positive selection of flowering time genes, including GmSPA3c. Our study sheds light on the evolutionary history of soybean and provides valuable genetic resources for its future breeding.

     

The Ghost of Selection Past: Rates of Evolution and Functional Divergence of Anciently Duplicated Genes

The duplication of genes and even complete genomes may be a prerequisite for major evolutionary transitions and the origin of evolutionary novelties. However, the evolutionary mechanisms of gene evolution and the origin of novel gene functions after gene duplication have been a subject of many debates. Recently, we compiled 26 groups of orthologous genes, which included one gene from human, mouse, and chicken, one or two genes from the tetraploid Xenopus and two genes from zebrafish. Comparative analysis and mapping data showed that these pairs of zebrafish genes were probably produced during a fish-specific genome duplication that occurred between 300 and 450 Mya, before the teleost radiation (Taylor et al. 2001). As discussed here, many of these retained duplicated genes code for DNA binding proteins. Different models have been developed to explain the retention of duplicated genes and in particular the subfunctionalization model of Force et al. (1999) could explain why so many developmental control genes have been retained. Other models are harder to reconcile with this particular set of duplicated genes. Most genes seem to have been subjected to strong purifying selection, keeping properties such as charge and polarity the same in both duplicates, although some evidence was found for positive Darwinian selection, in particular for Hox genes. However, since only the cumulative pattern of nucleotide substitutions can be studied, clear indications of positive Darwinian selection or neutrality may be hard to find for such anciently duplicated genes. Nevertheless, an increase in evolutionary rate in about half of the duplicated genes seems to suggest that either positive Darwinian selection has occurred or that functional constraints have been relaxed at one point in time during functional divergence. Received: 4 January 2001 / Accepted: 29 March 2001     

Significant sequence similarities in promoters and precursors of Arabidopsis thaliana non-conserved microRNAs

Some plant microRNAs have been shown to be de novo generated by inverted duplication from their target genes. Subsequent duplication events potentially generate multigene microRNA families. Within this article we provide supportive evidence for the inverted duplication model of plant microRNA evolution. First, we report that the precursors of four Arabidopsis thaliana microRNA families, miR157, miR158, miR405 and miR447 share nearly identical nucleotide sequences throughout the whole miRNA precursor between the family members. The extent and degree of sequence conservation is suggestive of recent evolutionary duplication events. Furthermore we found that sequence similarities are not restricted to the transcribed part but extend into the promoter regions. Thus the duplication event most probably included the promoter regions as well. Conserved elements in upstream regions of miR163 and its targets were also detected. This implies that the inverted duplication of target genes, at least in certain cases, had included the promoters of the target genes. Sequence conservation within promoters of miRNA families as well as between miRNA and its potential progenitor gene can be exploited for understanding the regulation of microRNA genes.     

Evolution of a secondary metabolic pathway from primary metabolism: shikimate and quinate biosynthesis in plants

The shikimate pathway synthesizes aromatic amino acids essential for protein biosynthesis. Shikimate dehydrogenase (SDH) is a central enzyme of this primary metabolic pathway, producing shikimate. The structurally similar quinate is a secondary metabolite synthesized by quinate dehydrogenase (QDH). SDH and QDH belong to the same gene family, which diverged into two phylogenetic clades after a defining gene duplication just prior to the angiosperm/gymnosperm split. Non‐seed plants that diverged before this duplication harbour only a single gene of this family. Extant representatives from the chlorophytes (Chlamydomonas reinhardtii), bryophytes (Physcomitrella patens) and lycophytes (Selaginella moellendorfii) encoded almost exclusively SDH activity in vitro. A reconstructed ancestral sequence representing the node just prior to the gene duplication also encoded SDH activity. Quinate dehydrogenase activity was gained only in seed plants following gene duplication. Quinate dehydrogenases of gymnosperms, represented here by Pinus taeda, may be reminiscent of an evolutionary intermediate since they encode equal SDH and QDH activities. The second copy in P. taeda maintained specificity for shikimate similar to the activity found in the angiosperm SDH sister clade. The codon for a tyrosine residue within the active site displayed a signature of positive selection at the node defining the QDH clade, where it changed to a glycine. Replacing the tyrosine with a glycine in a highly shikimate‐specific angiosperm SDH was sufficient to gain some QDH function. Thus, very few mutations were necessary to facilitate the evolution of QDH genes.     

Roles of small RNAs in soybean defense against Phytophthora sojae infection

The genus Phytophthora consists of many notorious pathogens of crops and forestry trees. At present, battling Phytophthora diseases is challenging due to a lack of understanding of their pathogenesis. We investigated the role of small RNAs in regulating soybean defense in response to infection by Phytophthora sojae, the second most destructive pathogen of soybean. Small RNAs, including microRNAs (miRNAs) and small interfering RNAs (siRNAs), are universal regulators that repress target gene expression in eukaryotes. We identified known and novel small RNAs that differentially accumulated during P. sojae infection in soybean roots. Among them, miR393 and miR166 were induced by heat‐inactivated P. sojae hyphae, indicating that they may be involved in soybean basal defense. Indeed, knocking down the level of mature miR393 led to enhanced susceptibility of soybean to P. sojae; furthermore, the expression of isoflavonoid biosynthetic genes was drastically reduced in miR393 knockdown roots. These data suggest that miR393 promotes soybean defense against P. sojae. In addition to miRNAs, P. sojae infection also resulted in increased accumulation of phased siRNAs (phasiRNAs) that are predominantly generated from canonical resistance genes encoding nucleotide binding‐leucine rich repeat proteins and genes encoding pentatricopeptide repeat‐containing proteins. This work identifies specific miRNAs and phasiRNAs that regulate defense‐associated genes in soybean during Phytophthora infection.     

Plasticity and innovation of regulatory mechanisms underlying seed oil content mediated by duplicated genes in the palaeopolyploid soybean

Many plants have undergone whole genome duplication (WGD). However, how regulatory networks underlying a particular trait are reshaped in polyploids has not been experimentally investigated. Here we show that the regulatory pathways modulating seed oil content, which involve WRINKLED1 (WRI1), LEAFY COTYLEDON1 (LEC1), and LEC2 in Arabidopsis, have been modified in the palaeopolyploid soybean. Such modifications include functional reduction of GmWRI1b of the GmWRI1a/GmWRI1b homoeologous pair relevant to WRI1, complementary non‐allelic dosage effects of the GmLEC1a/GmLEC1b homoeologous pair relevant to LEC1, pseudogenization of the singleton GmLEC2 relevant to LEC2, and the rise of the LEC2‐like function of GmABI3b, contrasting to its homoeolog GmABI3a, which maintains the ABSCISIC ACID INSENSITIVE 3 (ABI3)‐like function in modulating seed maturation and dormancy. The function of GmABI3b in modulating seed oil biosynthesis was fulfilled by direct binding to a RY (CATGCA) cis‐regulatory element in the GmWRI1a promoter, which was absent in the GmWRI1b promoter, resulting in reduction of the GmWRI1b expression. Nevertheless, the three regulators each exhibited similar intensities of purifying selection to their respective duplicates since these pairs were formed by a WGD event that is proposed to have occurred approximately 13 million years ago (mya), suggesting that the differentiation in spatiotemporal expression between the duplicated genes is more likely to be the outcome of neutral variation in regulatory sequences. This study thus exemplifies the plasticity, dynamics, and novelty of regulatory networks mediated by WGD.     

miRTar Hunter: A prediction system for identifying human microRNA target sites

MicroRNAs (miRNAs) are important regulators of gene expression and play crucial roles in many biological processes including apoptosis, differentiation, development, and tumorigenesis. Recent estimates suggest that more than 50% of human protein coding genes may be regulated by miRNAs and that each miRNA may bind to 300–400 target genes. Approximately 1,000 human miRNAs have been identified so far with each having up to hundreds of unique target mRNAs. However, the targets for a majority of these miRNAs have not been identified due to the lack of large-scale experimental detection techniques. Experimental detection of miRNA target sites is a costly and time-consuming process, even though identification of miRNA targets is critical to unraveling their functions in various biological processes. To identify miRNA targets, we developed miRTar Hunter, a novel computational approach for predicting target sites regardless of the presence or absence of a seed match or evolutionary sequence conservation. Our approach is based on a dynamic programming algorithm that incorporates more sequence-specific features and reflects the properties of various types of target sites that determine diverse aspects of complementarities between miRNAs and their targets. We evaluated the performance of our algorithm on 532 known human miRNA:target pairs and 59 experimentally-verified negative miRNA:target pairs, and also compared our method with three popular programs for 481 miRNA:target pairs. miRTar Hunter outperformed three popular existing algorithms in terms of recall and precision, indicating that our unique scheme to quantify the determinants of complementary sites is effective at detecting miRNA targets. miRTar Hunter is now available at http://203.230.194.162/~kbkim.     

Largely unlinked gene sets targeted by selection for domestication syndrome phenotypes in maize and sorghum

The domestication of diverse grain crops from wild grasses was a result of artificial selection for a suite of overlapping traits producing changes referred to in aggregate as ‘domestication syndrome’. Parallel phenotypic change can be accomplished by either selection on orthologous genes or selection on non‐orthologous genes with parallel phenotypic effects. To determine how often artificial selection for domestication traits in the grasses targeted orthologous genes, we employed resequencing data from wild and domesticated accessions of Zea (maize) and Sorghum (sorghum). Many ‘classic’ domestication genes identified through quantitative trait locus mapping in populations resulting from wild/domesticated crosses indeed show signatures of parallel selection in both maize and sorghum. However, the overall number of genes showing signatures of parallel selection in both species is not significantly different from that expected by chance. This suggests that while a small number of genes will extremely large phenotypic effects have been targeted repeatedly by artificial selection during domestication, the optimization part of domestication targeted small and largely non‐overlapping subsets of all possible genes which could produce equivalent phenotypic alterations.