A family of human microRNA genes from miniature inverted-repeat transposable elements

While hundreds of novel microRNA (miRNA) genes have been discovered in the last few years alone, the origin and evolution of these non-coding regulatory sequences remain largely obscure. In this report, we demonstrate that members of a recently discovered family of human miRNA genes, hsa-mir-548, are derived from Made1 transposable elements. Made1 elements are short miniature inverted-repeat transposable elements (MITEs), which consist of two 37 base pair (bp) terminal inverted repeats that flank 6 bp of internal sequence. Thus, Made1 elements are nearly perfect palindromes, and when expressed as RNA they form highly stable hairpin loops. Apparently, these Made1-related structures are recognized by the RNA interference enzymatic machinery and processed to form 22 bp mature miRNA sequences. Consistent with their origin from MITEs, hsa-mir-548 genes are primate-specific and have many potential paralogs in the human genome. There are more than 3,500 putative hsa-mir-548 target genes; analysis of their expression profiles and functional affinities suggests cancer-related regulatory roles for hsa-mir-548. Taken together, the characteristics of Made1 elements, and MITEs in general, point to a specific mechanism for the generation of numerous small regulatory RNAs and target sites throughout the genome. The evolutionary lineage-specific nature of MITEs could also provide for the generation of novel regulatory phenotypes related to species diversification. Finally, we propose that MITEs may represent an evolutionary link between siRNAs and miRNAs.     

Evolution of microRNA genes in Oryza sativa and Arabidopsis thaliana: an update of the inverted duplication model

The origin and evolution of microRNA (miRNA) genes, which are of significance in tuning and buffering gene expressions in a number of critical cellular processes, have long attracted evolutionary biologists. However, genome-wide perspectives on their origins, potential mechanisms of their de novo generation and subsequent evolution remain largely unsolved in flowering plants. Here, genome-wide analyses of Oryza sativa and Arabidopsis thaliana revealed apparently divergent patterns of miRNA gene origins. A large proportion of miRNA genes in O. sativa were TE-related and MITE-related miRNAs in particular, whereas the fraction of these miRNA genes much decreased in A. thaliana. Our results show that the majority of TE-related and pseudogene-related miRNA genes have originated through inverted duplication instead of segmental or tandem duplication events. Based on the presented findings, we hypothesize and illustrate the four likely molecular mechanisms to de novo generate novel miRNA genes from TEs and pseudogenes. Our rice genome analysis demonstrates that non-MITEs and MITEs mediated inverted duplications have played different roles in de novo generating miRNA genes. It is confirmed that the previously proposed inverted duplication model may give explanations for non-MITEs mediated duplication events. However, many other miRNA genes, known from the earlier proposed model, were rather arisen from MITE transpositions into target genes to yield binding sites. We further investigated evolutionary processes spawned from de novo generated to maturely-formed miRNA genes and their regulatory systems. We found that miRNAs increase the tunability of some gene regulatory systems with low gene copy numbers. The results also suggest that gene balance effects may have largely contributed to the evolution of miRNA regulatory systems.     

Phylogeny and origin of pearl millet (Pennisetum glaucum [L.] R. Br) as revealed by microsatellite loci

During the last 12,000 years, different cultures around the world have domesticated cereal crops. Several studies investigated the evolutionary history and domestication of cereals such as wheat in the Middle East, rice in Asia or maize in America. The domestication process in Africa has led to the emergence of important cereal crops like pearl millet in Sahelian Africa. In this study, we used 27 microsatellite loci to analyze 84 wild accessions and 355 cultivated accessions originating from the whole pearl millet distribution area in Africa and Asia. We found significantly higher diversity in the wild pearl millet group. The cultivated pearl millet sample possessed 81% of the alleles and 83% of the genetic diversity of the wild pearl millet sample. Using Bayesian approaches, we identified intermediate genotypes between the cultivated and wild groups. We then analyzed the phylogenetic relationship among accessions not showing introgression and found that a monophyletic origin of cultivated pearl millet in West Africa is the most likely scenario supported by our data set.     

Origins of weedy rice revealed by polymorphisms of chloroplast DNA sequences and nuclear microsatellites

Conspecific weeds that permanently infest worldwide agroecosystems are evolved from their crop species. These weeds cause substantial problems for crop production by competing for resources in agricultural fields. Weedy rice represents such a conspecific weed infesting rice ecosystems, and causing tremendous rice yield losses owing to its strong competitiveness and abundant genetic diversity, likely resulted from its complex origins. Here, we report the use of chloroplast DNA (cpDNA) fingerprints to determine whether weedy rice is evolved from its wild (exo‐feral) or cultivated (endo‐feral) rice progenitor as the maternal donor in recent hybridization events. In addition, we also applied nuclear simple sequence repeat (SSR) markers to confirm the exo‐feral or endo‐feral origins of weedy rice accessions determined by the cpDNA fingerprints. We found that the studied weedy rice accessions evolved either from their wild or cultivated rice progenitor, as the maternal donor, based on the cpDNA network and structure analyses. Combined analyses of cpDNA and nuclear SSR markers indicated that a much greater proportion of weedy rice accessions had the endo‐feral origin. In addition, results from the genetic structure of nuclear SSR markers indicated that weedy rice accessions from the endo‐feral pathway are distinctly associated with either indica or japonica rice cultivars, suggesting their complex origins through crop–weed introgression. The complex pathways of origin and evolution could greatly promote genetic diversity of weedy rice. Therefore, innovative methods should be developed for effective weedy rice control.     

Complex regulation of two target genes encoding SPX-MFS proteins by rice miR827 in response to phosphate starvation

Here we report on the characterization of rice osa-miR827 and its two target genes, OsSPX-MFS1 and OsSPX-MFS2, which encode SPX-MFS proteins predicted to be implicated in phosphate (Pi) sensing or transport. We first show by Northern blot analysis that osa-miR827 is strongly induced by Pi starvation in both shoots and roots. Hybridization of osa-miR827 in situ confirms its strong induction by Pi starvation, with signals concentrated in mesophyll, epidermis and ground tissues of roots. In parallel, we analyzed the responses of the two OsSPX-MFS1 and OsSPX-MFS2 gene targets to Pi starvation. OsSPX-MFS1 mRNA is mainly expressed in shoots under sufficient Pi supply while its expression is reduced on Pi starvation, revealing a direct relationship between induction of osa-miR827 and down-regulation of OsSPX-MFS1. In contrast, OsSPX-MFS2 responds in a diametrically opposed manner to Pi starvation. The accumulation of OsSPX-MFS2 mRNA is dramatically enhanced under Pi starvation, suggesting the involvement of complex regulation of osa-miR827 and its two target genes. We further produced transgenic rice lines overexpressing osa-miR827 and T-DNA knockout mutant lines in which the expression of osa-miR827 is abolished. Compared with wild-type controls, both target mRNAs exhibit similar changes, their expression being reduced and increased in overexpressing and knockout lines, respectively. This suggests that OsSPX-MFS1 and OsSPX-MFS2 are both negatively regulated by osa-miR827 abundance although they respond differently to external Pi conditions. We propose that this is a complex mechanism comprising fine tuning of spatial or temporal regulation of both targets by osa-miR827.     

PIF- and Pong-like transposable elements: distribution, evolution and relationship with Tourist-like miniature inverted-repeat transposable elements

Miniature inverted-repeat transposable elements (MITEs) are short, nonautonomous DNA elements that are widespread and abundant in plant genomes. Most of the hundreds of thousands of MITEs identified to date have been divided into two major groups on the basis of shared structural and sequence characteristics: Tourist-like and Stowaway-like. Since MITEs have no coding capacity, they must rely on transposases encoded by other elements. Two active transposons, the maize P Instability Factor (PIF) and the rice Pong element, have recently been implicated as sources of transposase for Tourist-like MITEs. Here we report that PIF- and Pong-like elements are widespread, diverse, and abundant in eukaryotes with hundreds of element-associated transposases found in a variety of plant, animal, and fungal genomes. The availability of virtually the entire rice genome sequence facilitated the identification of all the PIF/Pong-like elements in this organism and permitted a comprehensive analysis of their relationship with Tourist-like MITEs. Taken together, our results indicate that PIF and Pong are founding members of a large eukaryotic transposon superfamily and that members of this superfamily are responsible for the origin and amplification of Tourist-like MITEs.     

MITE-Hunter: a program for discovering miniature inverted-repeat transposable elements from genomic sequences

Miniature inverted-repeat transposable elements (MITEs) are a special type of Class 2 non-autonomous transposable element (TE) that are abundant in the non-coding regions of the genes of many plant and animal species. The accurate identification of MITEs has been a challenge for existing programs because they lack coding sequences and, as such, evolve very rapidly. Because of their importance to gene and genome evolution, we developed MITE-Hunter, a program pipeline that can identify MITEs as well as other small Class 2 non-autonomous TEs from genomic DNA data sets. The output of MITE-Hunter is composed of consensus TE sequences grouped into families that can be used as a library file for homology-based TE detection programs such as RepeatMasker. MITE-Hunter was evaluated by searching the rice genomic database and comparing the output with known rice TEs. It discovered most of the previously reported rice MITEs (97.6%), and found sixteen new elements. MITE-Hunter was also compared with two other MITE discovery programs, FINDMITE and MUST. Unlike MITE-Hunter, neither of these programs can search large genomic data sets including whole genome sequences. More importantly, MITE-Hunter is significantly more accurate than either FINDMITE or MUST as the vast majority of their outputs are false-positives.     

Amplification dynamics of miniature inverted-repeat transposable elements and their impact on rice trait variability

Transposable elements (TEs) are a rich source of genetic variability. Among TEs, miniature inverted-repeat TEs (MITEs) are of particular interest as they are present in high copy numbers in plant genomes and are closely associated with genes. MITEs are deletion derivatives of class II transposons, and can be mobilized by the transposases encoded by the latter through a typical cut-and-paste mechanism. However, MITEs are typically present at much higher copy numbers than class II transposons. We present here an analysis of 103 109 transposon insertion polymorphisms (TIPs) in 738 Oryza sativa genomes representing the main rice population groups. We show that an important fraction of MITE insertions has been fixed in rice concomitantly with its domestication. However, another fraction of MITE insertions is present at low frequencies. We performed MITE TIP-genome-wide association studies (TIP-GWAS) to study the impact of these elements on agronomically important traits and found that these elements uncover more trait associations than single nucleotide polymorphisms (SNPs) on important phenotypes such as grain width. Finally, using SNP-GWAS and TIP-GWAS we provide evidence of the replicative amplification of MITEs.     

Insertion Preference of Maize and Rice Miniature Inverted Repeat Transposable Elements as Revealed by the Analysis of Nested Elements

A 128-bp insertion into the maize waxy-B2 allele led to the discovery of Tourist, a family of miniature inverted repeat transposable elements (MITEs). As a special category of nonautonomous elements, MITEs are distinguished by their high copy number, small size, and close association with plant genes. In maize, some Tourist elements (named Tourist-Zm) are present as adjacent or nested insertions. To determine whether the formation of multimers is a common feature of MITEs, we performed a more thorough survey, including an estimation of the proportion of multimers, with 30.2 Mb of publicly available rice genome sequence. Among the 6600 MITEs identified, >10% were present as multimers. The proportion of multimers differs for different MITE families. For some MITE families, a high frequency of self-insertions was found. The fact that all 340 multimers are unique indicates that the multimers are not capable of further amplification.     

Characterization of Stowaway MITEs in pea (Pisum sativum L.) and identification of their potential master elements.

We have investigated miniature inverted-repeat transposable elements (MITEs) of the Stowaway family and corresponding Mariner-like master elements that could potentially facilitate their mobilization in the genome of the garden pea (Pisum sativum L.). The population of pea Stowaway MITEs consists of 103-104 copies dispersed in the genome. Judging from a sequence analysis of 17 isolated Stowaway elements and their flanking genomic regions, the elements are relatively uniform in size and sequence and occur in the vicinity of genes as well as within repetitive sequences. Insertional polymorphism of several elements was detected among various Pisum accessions, suggesting they were still transpositionally active during diversification of these taxa. The identification of several Mariner-like elements (MLEs) harboring intact open reading frames, capable of encoding a transposase, further supports a recent mobilization of the Stowaway elements. Using transposase-coding sequences as a hybridization probe, we estimated that there are about 50 MLE sequences in the pea genome. Among the 5 elements sequenced, 3 distinct subfamilies showing mutual similarities within their transposase-coding regions, but otherwise diverged in sequence, were distinguished and designated as Psmar-1 to Psmar-3. The terminal inverted repeats (TIRs) of these MLE subfamilies differed in their homology to the TIRs of Stowaway MITEs. The homlogy ranged from 9 bp in Psmar-3 to 30 bp in Psmar-1, which corresponds to the complete Stowaway TIR sequence. Based on this feature, the Psmar-1 elements are believed to be the most likely candidates for the master elements of the Stowaway MITEs in pea.     

Evolutionary dynamics of wheat mitochondrial gene structure with special remarks on the origin and effects of RNA editing in cereals

We investigated the evolutionary dynamics of wheat mitochondrial genes with respect to their structural differentiation during organellar evolution, and to mutations that occurred during cereal evolution. First, we compared the nucleotide sequences of three wheat mitochondrial genes to those of wheat chloroplast, alpha-proteobacterium and cyanobacterium orthologs. As a result, we were able to (1) differentiate the conserved and variable segments of the orthologs, (2) reveal the functional importance of the conserved segments, and (3) provide a corroborative support for the alpha-proteobacterial and cyanobacterial origins of those mitochondrial and chloroplast genes, respectively. Second, we compared the nucleotide sequences of wheat mitochondrial genes to those of rice and maize to determine the types and frequencies of base changes and indels occurred in cereal evolution. Our analyses showed that both the evolutionary speed, in terms of number of base substitutions per site, and the transition/transversion ratio of the cereal mitochondrial genes were less than two-fifths of those of the chloroplast genes. Eight mitochondrial gene groups differed in their evolutionary variability, RNA and Complex I (nad) genes being most stable whereas Complex V (atp) and ribosomal protein genes most variable. C-to-T transition was the most frequent type of base change; C-to-G and G-to-C transversions occurred at lower rates than all other changes. The excess of C-to-T transitions was attributed to C-to-U RNA editing that developed in early stage of vascular plant evolution. On the contrary, the editing of C residues at cereal T-to-C transition sites developed mostly during cereal divergence. Most indels were associated with short direct repeats, suggesting intra- and intermolecular recombination as an important mechanism for their origin. Most of the repeats associated with indels were di- or trinucleotides, although no preference was noticed for their sequences. The maize mt genome was characterized by a high incidence of indels, comparing to the wheat and rice mt genomes.     

The CRISPR/Cas9 system produces specific and homozygous targeted gene editing in rice in one generation

The CRISPR/Cas9 system has been demonstrated to efficiently induce targeted gene editing in a variety of organisms including plants. Recent work showed that CRISPR/Cas9‐induced gene mutations in Arabidopsis were mostly somatic mutations in the early generation, although some mutations could be stably inherited in later generations. However, it remains unclear whether this system will work similarly in crops such as rice. In this study, we tested in two rice subspecies 11 target genes for their amenability to CRISPR/Cas9‐induced editing and determined the patterns, specificity and heritability of the gene modifications. Analysis of the genotypes and frequency of edited genes in the first generation of transformed plants (T0) showed that the CRISPR/Cas9 system was highly efficient in rice, with target genes edited in nearly half of the transformed embryogenic cells before their first cell division. Homozygotes of edited target genes were readily found in T0 plants. The gene mutations were passed to the next generation (T1) following classic Mendelian law, without any detectable new mutation or reversion. Even with extensive searches including whole genome resequencing, we could not find any evidence of large‐scale off‐targeting in rice for any of the many targets tested in this study. By specifically sequencing the putative off‐target sites of a large number of T0 plants, low‐frequency mutations were found in only one off‐target site where the sequence had 1‐bp difference from the intended target. Overall, the data in this study point to the CRISPR/Cas9 system being a powerful tool in crop genome engineering.     

Matrix attachment regions and structural colinearity in the genomes of two grass species.

In order to gain insights into the relationship between spatial organization of the genome and genome function we have initiated studies of the co-linear Sh2/A1- homologous regions of rice (30 kb) and sorghum (50 kb). We have identified the locations of matrix attachment regions (MARs) in these homologous chromosome segments, which could serve as anchors for individual structural units or loops. Despite the fact that the nucleotide sequences serving as MARs were not detectably conserved, the general organizational patterns of MARs relative to the neighboring genes were preserved. All identified genes were placed in individual loops that were of comparable size for homologous genes. Hence, gene composition, gene orientation, gene order and the placement of genes into structural units has been evolutionarily conserved in this region. Our analysis demonstrated that the occurrence of various 'MAR motifs' is not indicative of MAR location. However, most of the MARs discovered in the two genomic regions were found to co-localize with miniature inverted repeat transposable elements (MITEs), suggesting that MITEs preferentially insert near MARs and/or that they can serve as MARs.     

Hybridisation‐based target enrichment of phenology genes to dissect the genetic basis of yield and adaptation in barley

Barley (Hordeum vulgare L.) is a major cereal grain widely used for livestock feed, brewing malts and human food. Grain yield is the most important breeding target for genetic improvement and largely depends on optimal timing of flowering. Little is known about the allelic diversity of genes that underlie flowering time in domesticated barley, the genetic changes that have occurred during breeding, and their impact on yield and adaptation. Here, we report a comprehensive genomic assessment of a worldwide collection of 895 barley accessions based on the targeted resequencing of phenology genes. A versatile target‐capture method was used to detect genome‐wide polymorphisms in a panel of 174 flowering time‐related genes, chosen based on prior knowledge from barley, rice and Arabidopsis thaliana. Association studies identified novel polymorphisms that accounted for observed phenotypic variation in phenology and grain yield, and explained improvements in adaptation as a result of historical breeding of Australian barley cultivars. We found that 50% of genetic variants associated with grain yield, and 67% of the plant height variation was also associated with phenology. The precise identification of favourable alleles provides a genomic basis to improve barley yield traits and to enhance adaptation for specific production areas.     

Evolutionary relationship of plant catalase genes inferred from exon-intron structures: isozyme divergence after the separation of monocots and dicots

In order to understand the molecular evolution of catalase genes in higher plants, we compared the exon-intron structures of 12 genomic sequences from six plant species. It was assumed that the putative single primordial catalase gene had seven introns, because only those catalase genes having this structure are found in the monocotyledonae and dicotyledonae classes. After the evolutionary divergence of monocots from dicots, consecutive duplication of the primordial gene followed by the differential loss of introns occurred in each class to form three (or possibly four in dicots) diverse isozyme genes. In monocots, three ancestral isozyme genes were formed before the divergence of ancestral rice and maize. One of the rice genes, CatA, has an entirely new short intron which was not found in any other plant catalase gene examined. We have investigated the existence of the intron in the CatA homolog in other rice species by polymerase chain reaction (PCR) analysis. One major PCR product was found with the genomic DNAs from O. sativa (indica and japonica types), O. rufipogon and O. glaberrima. DNAs from several accessions of O. longistaminata showed variation in both the number and size of the DNA fragments amplified. PCR analyses and sequencing of the PCR products revealed that there are several CatA homologs having different sequences in some accessions of O. longistaminata. We have extended our study to other species in the Poaceae. The results suggest that the gain of the intron, most likely by insertion of a retroposon, took place in the ancestral genome of rice after its evolutionary divergence from other ancestral cereals such as barley, wheat and oat. Received: 20 November 1997 / Accepted: 5 January 1998