Domestication of rice has reduced the occurrence of transposable elements within gene coding regions

A chromosome-scale genome sequence of sudangrass (Sorghum sudanense) highlights the genome evolution and regulation of dhurrin biosynthesis

Sudangrass is more similar to US commercial sorghums than to cultivated sorghums from Africa sequence-wise and contain significantly lower dhurrin than sorghums. CYP79A1 is linked to dhurrin content in sorghum. Sudangrass [Sorghum sudanense (Piper) Stapf] is a hybrid between grain sorghum and its wild relative S. bicolor ssp. verticilliflorum and is grown as a forage crop due to its high biomass production and low dhurrin content compared to sorghum. In this study, we sequenced the sudangrass genome and showed that the assembled genome was 715.95 Mb with 35,243 protein-coding genes. Phylogenetic analysis with whole genome proteomes demonstrated that the sudangrass genome was more similar to US commercial sorghums than to its wild relatives and cultivated sorghums from Africa. We confirmed that at seedling stage, sudangrass accessions contained significantly lower dhurrin as measured by hydrocyanic acid potential (HCN-p) than cultivated sorghum accessions. Genome-wide association study identified a QTL most tightly associated with HCN-p and the linked SNPs were located in the 3' UTR of Sobic.001G012300 which encodes CYP79A1, the enzyme that catalyzes the first step of dhurrin biosynthesis. As in other grasses such as maize and rice, we also found that copia/gypsy long terminal repeat (LTR) retrotransposons were more abundant in cultivated than in wild sorghums, implying that crop domestication in the grasses was accompanied by increased copia/gypsy LTR retrotransposon insertions in the genomes.

Miniature Inverted-Repeat Transposable Elements: Small DNA Transposons That Have Contributed to Plant MICRORNA Gene Evolution

Angiosperms form the largest phylum within the Plantae kingdom and show remarkable genetic variation due to the considerable difference in the nuclear genome size of each species. Transposable elements (TEs), mobile DNA sequences that can amplify and change their chromosome position, account for much of the difference in nuclear genome size between individual angiosperm species. Considering the dramatic consequences of TE movement, including the complete loss of gene function, it is unsurprising that the angiosperms have developed elegant molecular strategies to control TE amplification and movement. Specifically, the RNA-directed DNA methylation (RdDM) pathway, directed by the repeat-associated small-interfering RNA (rasiRNA) class of small regulatory RNA, forms the primary line of defense to control TE activity in the angiosperms. However, the miniature inverted-repeat transposable element (MITE) species of TE has at times avoided the repressive effects imposed by the rasiRNA-directed RdDM pathway. MITE proliferation in angiosperm nuclear genomes is due to their preference to transpose within gene-rich regions, a pattern of transposition that has enabled MITEs to gain further transcriptional activity. The sequence-based properties of a MITE results in the synthesis of a noncoding RNA (ncRNA), which, after transcription, folds to form a structure that closely resembles those of the precursor transcripts of the microRNA (miRNA) class of small regulatory RNA. This shared folding structure results in a MITE-derived miRNA being processed from the MITE-transcribed ncRNA, and post-maturation, the MITE-derived miRNA can be used by the core protein machinery of the miRNA pathway to regulate the expression of protein-coding genes that harbor homologous MITE insertions. Here, we outline the considerable contribution that the MITE species of TE have made to expanding the miRNA repertoire of the angiosperms.

Polymorphic inverted repeats near coding genes impact chromatin topology and phenotypic traits in Arabidopsis thaliana

Transposons are mobile elements that are commonly silenced to protect eukaryotic genome integrity. In plants, transposable element (TE)-derived inverted repeats (IRs) are commonly found near genes, where they affect host gene expression. However, the molecular mechanisms of such regulation are unclear in most cases. Expression of these IRs is associated with production of 24-nt small RNAs, methylation of the IRs, and drastic changes in local 3D chromatin organization. Notably, many of these IRs differ between Arabidopsis thaliana accessions, causing variation in short-range chromatin interactions and gene expression. CRISPR-Cas9-mediated disruption of two IRs leads to a switch in genome topology and gene expression with phenotypic consequences. Our data show that insertion of an IR near a gene provides an anchor point for chromatin interactions that profoundly impact the activity of neighboring loci. This turns IRs into powerful evolutionary agents that can contribute to rapid adaptation.

Genomic-environmental associations in wild cranberry (Vaccinium macrocarpon Ait.)

Understanding the genetic basis of local adaptation in natural plant populations, particularly crop wild relatives, may be highly useful for plant breeding. By characterizing genetic variation for adaptation to potentially stressful environmental conditions, breeders can make targeted use of crop wild relatives to develop cultivars for novel or changing environments. This is especially appealing for improving long-lived woody perennial crops such as the American cranberry (Vaccinium macrocarpon Ait.), the cultivation of which is challenged by biotic and abiotic stresses. In this study, we used environmental association analyses in a collection of 111 wild cranberry accessions to identify potentially adaptive genomic regions for a range of bioclimatic and soil conditions. We detected 126 significant associations between SNP marker loci and environmental variables describing temperature, precipitation, and soil attributes. Many of these markers tagged genes with functional annotations strongly suggesting a role in adaptation to biotic or abiotic conditions. Despite relatively low genetic variation in cranberry, our results suggest that local adaptation to divergent environments is indeed present, and the identification of potentially adaptive genetic variation may enable a selective use of this germplasm for breeding more stress-tolerant cultivars.

Genetic Conservation of CBS Domain Containing Protein Family in Oryza Species and Their Association with Abiotic Stress Responses

Crop Wild Relatives (CWRs) form a comprehensive gene pool that can answer the queries related to plant domestication, speciation, and ecological adaptation. The genus ‘Oryza’ comprises about 27 species, of which two are cultivated, while the remaining are wild. Here, we have attempted to understand the conservation and diversification of the genes encoding Cystathionine β-synthase (CBS) domain-containing proteins (CDCPs) in domesticated and CWRs of rice. Few members of CDCPs were previously identified to be stress-responsive and associated with multiple stress tolerance in rice. Through genome-wide analysis of eleven rice genomes, we identified a total of 36 genes encoding CDCPs in O. longistaminata, 38 in O. glaberrima, 39 each in O. rufipogon, O. glumaepatula, O. brachyantha, O. punctata, and O. sativa subsp. japonica, 40 each in O. barthii and O. meridionalis, 41 in O. nivara, and 42 in O. sativa subsp. indica. Gene duplication analysis as well as non-synonymous and synonymous substitutions in the duplicated gene pairs indicated that this family is shaped majorly by the negative or purifying selection pressure through the long-term evolution process. We identified the presence of two additional hetero-domains, namely TerCH and CoatomerE (specifically in O. sativa subsp. indica), which were not reported previously in plant CDCPs. The in silico expression analysis revealed some of the members to be responsive to various abiotic stresses. Furthermore, the qRT-PCR based analysis identified some members to be highly inducive specifically in salt-tolerant genotype in response to salinity. The cis-regulatory element analysis predicted the presence of numerous stress as well as a few phytohormone-responsive elements in their promoter region. The data presented in this study would be helpful in the characterization of these CDCPs from rice, particularly in relation to abiotic stress tolerance.

Development of molecular markers based on LTR retrotransposon in the Cleistogenes songorica genome

Long terminal repeat retrotransposons (LTR-RTs) contribute a large fraction of many sequenced plant genomes and play important roles in genomic diversity and phenotypic variations. LTR-RTs are abundantly distributed in plant genomes, facilitating the development of markers based on LTR-RTs for a variety of genotyping purposes. Whole-genome analysis of LTR-RTs was performed in Cleistogenes songorica. A total of 299,079 LTR-RTs were identified and classified as Gypsy type, Copia type, or other type. LTR-RTs were widely distributed in the genome, enriched in the heterochromatic region of the chromosome, and negatively correlated with gene distribution. However, approximately one-fifth of genes were still interrupted by LTR-RTs, and these genes are annotated. Furthermore, four types of primer pairs (PPs) were designed, namely, retrotransposon-based insertion polymorphisms, inter-retrotransposon amplified polymorphisms, insertion site-based polymorphisms, and retrotransposon-microsatellite amplified polymorphisms. A total of 350 PPs were screened in 23 accessions of the genus Cleistogenes, of which 80 PPs showed polymorphism, and 72 PPs showed transferability among Gramineae and non-Gramineae species. In addition, a comparative analysis of homologous LTR-RTs was performed with other related grasses. Taken together, the study will serve as a valuable resource for genotyping applications for C. songorica and related grasses.

Gene co-expression analysis reveals transcriptome divergence between wild and cultivated chickpea under drought stress

Ancestral adaptations in crop wild relatives can provide a genetic reservoir for crop improvement. Here we document physiological changes to mild and severe drought stress, and the associated transcriptome dynamics in both wild and cultivated chickpea. Over 60% of transcriptional changes were related to metabolism, indicating that metabolic plasticity is a core and conserved drought response. In addition, changes in RNA processing and protein turnover were predominant in the data, suggestive of broad restructuring of the chickpea proteome in response to drought. While 12% of the drought-responsive transcripts have similar dynamics in cultivated and wild accessions, numerous transcripts had expression patterns unique to particular genotypes, or that distinguished wild from cultivated genotypes and whose divergence may be a consequence of domestication. These and other comparisons provide a transcriptional correlate of previously described species' genetic diversity, with wild accessions well differentiated from each other and from cultivars, and cultivars essentially indistinguishable at the broad transcriptome level. We identified metabolic pathways such as phenylpropanoid metabolism, and biological processes such as stomatal development, which are differentially regulated across genotypes with potential consequences on drought tolerance. These data indicate that wild Cicer reticulatum may provide both conserved and divergent mechanisms as a resource in breeding for drought tolerance in cultivated chickpea.

Comparative Study of Pine Reference Genomes Reveals Transposable Element Interconnected Gene Networks

Sequencing the giga-genomes of several pine species has enabled comparative genomic analyses of these outcrossing tree species. Previous studies have revealed the wide distribution and extraordinary diversity of transposable elements (TEs) that occupy the large intergenic spaces in conifer genomes. In this study, we analyzed the distribution of TEs in gene regions of the assembled genomes of Pinus taeda and Pinus lambertiana using high-performance computing resources. The quality of draft genomes and the genome annotation have significant consequences for the investigation of TEs and these aspects are discussed. Several TE families frequently inserted into genes or their flanks were identified in both species’ genomes. Potentially important sequence motifs were identified in TEs that could bind additional regulatory factors, promoting gene network formation with faster or enhanced transcription initiation. Node genes that contain many TEs were observed in multiple potential transposable element-associated networks. This study demonstrated the increased accumulation of TEs in the introns of stress-responsive genes of pines and suggests the possibility of rewiring them into responsive networks and sub-networks interconnected with node genes containing multiple TEs. Many such regulatory influences could lead to the adaptive environmental response clines that are characteristic of naturally spread pine populations.

DNA methylation is reconfigured at the onset of reproduction in rice shoot apical meristem

DNA methylation is an epigenetic modification that specifies the basic state of pluripotent stem cells and regulates the developmental transition from stem cells to various cell types. In flowering plants, the shoot apical meristem (SAM) contains a pluripotent stem cell population which generates the aerial part of plants including the germ cells. Under appropriate conditions, the SAM undergoes a developmental transition from a leaf-forming vegetative SAM to an inflorescence- and flower-forming reproductive SAM. While SAM characteristics are largely altered in this transition, the complete picture of DNA methylation remains elusive. Here, by analyzing whole-genome DNA methylation of isolated rice SAMs in the vegetative and reproductive stages, we show that methylation at CHH sites is kept high, particularly at transposable elements (TEs), in the vegetative SAM relative to the differentiated leaf, and increases in the reproductive SAM via the RNA-dependent DNA methylation pathway. We also show that half of the TEs that were highly methylated in gametes had already undergone CHH hypermethylation in the SAM. Our results indicate that changes in DNA methylation begin in the SAM long before germ cell differentiation to protect the genome from harmful TEs.

The shoot apical meristem of flowering plants transitions from forming leaves to floral organs. Here Higo et al. show that DNA methylation of many transposons that are hypermethylated in gametes is established in the SAM before flowering, suggesting it protects against harmful transposition long before germ cell differentiation.

Transcriptional regulation of genes bearing intronic heterochromatin in the rice genome

Intronic regions of eukaryotic genomes accumulate many Transposable Elements (TEs). Intronic TEs often trigger the formation of transcriptionally repressive heterochromatin, even within transcription-permissive chromatin environments. Although TE-bearing introns are widely observed in eukaryotic genomes, their epigenetic states, impacts on gene regulation and function, and their contributions to genetic diversity and evolution, remain poorly understood. In this study, we investigated the genome-wide distribution of intronic TEs and their epigenetic states in the Oryza sativa genome, where TEs comprise 35% of the genome. We found that over 10% of rice genes contain intronic heterochromatin, most of which are associated with TEs and repetitive sequences. These heterochromatic introns are longer and highly enriched in promoter-proximal positions. On the other hand, introns also accumulate hypomethylated short TEs. Genes with heterochromatic introns are implicated in various biological functions. Transcription of genes bearing intronic heterochromatin is regulated by an epigenetic mechanism involving the conserved factor OsIBM2, mutation of which results in severe developmental and reproductive defects. Furthermore, we found that heterochromatic introns evolve rapidly compared to non-heterochromatic introns. Our study demonstrates that heterochromatin is a common epigenetic feature associated with actively transcribed genes in the rice genome.

Chromosome-Scale Assembly of Winter Oilseed Rape Brassica napus

Rapeseed (Brassica napus), the second most important oilseed crop globally, originated from an interspecific hybridization between B. rapa and B. oleracea. After this genome collision, B. napus underwent extensive genome restructuring, via homoeologous chromosome exchanges, resulting in widespread segmental deletions and duplications. Illicit pairing among genetically similar homoeologous chromosomes during meiosis is common in recent allopolyploids like B. napus, and post-polyploidization restructuring compounds the difficulties of assembling a complex polyploid plant genome. Specifically, genomic rearrangements between highly similar chromosomes are challenging to detect due to the limitation of sequencing read length and ambiguous alignment of reads. Recent advances in long read sequencing technologies provide promising new opportunities to unravel the genome complexities of B. napus by encompassing breakpoints of genomic rearrangements with high specificity. Moreover, recent evidence revealed ongoing genomic exchanges in natural B. napus, highlighting the need for multiple reference genomes to capture structural variants between accessions. Here we report the first long-read genome assembly of a winter B. napus cultivar. We sequenced the German winter oilseed rape accession 'Express 617' using 54.5x of long reads. Short reads, linked reads, optical map data and high-density genetic maps were used to further correct and scaffold the assembly to form pseudochromosomes. The assembled Express 617 genome provides another valuable resource for Brassica genomics in understanding the genetic consequences of polyploidization, crop domestication, and breeding of recently-formed crop species.

Chromosome-Scale Assembly of Winter Oilseed Rape Brassica napus

Rapeseed (Brassica napus), the second most important oilseed crop globally, originated from an interspecific hybridization between B. rapa and B. oleracea. After this genome collision, B. napus underwent extensive genome restructuring, via homoeologous chromosome exchanges, resulting in widespread segmental deletions and duplications. Illicit pairing among genetically similar homoeologous chromosomes during meiosis is common in recent allopolyploids like B. napus, and post-polyploidization restructuring compounds the difficulties of assembling a complex polyploid plant genome. Specifically, genomic rearrangements between highly similar chromosomes are challenging to detect due to the limitation of sequencing read length and ambiguous alignment of reads. Recent advances in long read sequencing technologies provide promising new opportunities to unravel the genome complexities of B. napus by encompassing breakpoints of genomic rearrangements with high specificity. Moreover, recent evidence revealed ongoing genomic exchanges in natural B. napus, highlighting the need for multiple reference genomes to capture structural variants between accessions. Here we report the first long-read genome assembly of a winter B. napus cultivar. We sequenced the German winter oilseed rape accession 'Express 617' using 54.5x of long reads. Short reads, linked reads, optical map data and high-density genetic maps were used to further correct and scaffold the assembly to form pseudochromosomes. The assembled Express 617 genome provides another valuable resource for Brassica genomics in understanding the genetic consequences of polyploidization, crop domestication, and breeding of recently-formed crop species.

New clues into the mechanisms of rice domestication

Domestication of rice involved incorporation of specific yield-related changes in wild species of rice. This agricultural process has been of significant interest for plant biologists. The recent advance in genomics has provided new tools to investigate the genetic basis and consequences of domestication. Several genes involved in domestication and diversification process have been characterized, and as expected, this list is over-represented by transcription factors and their cofactors. Most often the modification orchestrated expression levels of genes such as those coding for transcription factors. It has been proposed that transcriptional regulators and their regulation is likely a major theme controlling morphological differences between crops and their progenitors. However, recent data indicate that single amino acid changes in genes coding for key proteins as well as epigenetic and small RNA-mediated pathways also contributed towards domesticationassociated phenotypes.

Identification of an active miniature inverted-repeat transposable element mJing in rice

Miniature inverted-repeat transposable elements (MITEs) are structurally homogeneous non-autonomous DNA transposons with high copy numbers that play important roles in genome evolution and diversification. Here, we analyzed the rice high-tillering dwarf (htd) mutant in an advanced backcross population between cultivated and wild rice, and identified an active MITE named miniature Jing (mJing). The mJing element belongs to the PIF/Harbinger superfamily. japonica rice var. Nipponbare and indica var. 93-11 harbor 72 and 79 mJing family members, respectively, have undergone multiple rounds of amplification bursts during the evolution of Asian cultivated rice (Oryza sativa L.). A heterologous transposition experiment in Arabidopsis thaliana indicated that the autonomous element Jing is likely to have provides the transposase needed for mJing mobilization. We identified 297 mJing insertion sites and their presence/absence polymorphism among 71 rice samples through targeted high-throughput sequencing. The results showed that the copy number of mJing varies dramatically among Asian cultivated rice (O. sativa), its wild ancestor (O. rufipogon), and African cultivated rice (O. glaberrima) and that some mJing insertions are subject to directional selection. These findings suggest that the amplification and removal of mJing elements have played an important role in rice genome evolution and species diversification.

Major Domestication-Related Phenotypes in Indica Rice Are Due to Loss of miRNA-Mediated Laccase Silencing

Domestication of rice (Oryza sativa) included conversion of perennial wild species with few seeds to short plants that produced abundant seeds. Most domestication-associated changes were due to variations in transcription factors and other key proteins such as enzymes. Here, we show that multiple yield-related traits associated with indica rice domestication are linked to micro (mi) RNA-mediated regulation. Analysis of small (s) RNA data sets from cultivated indica rice lines, a few landraces, and two wild relatives of rice revealed the presence of abundant 22-nucleotide (nt) reads in wild relatives that mapped to miR397 precursors. miR397 was expressed at very high levels in wild relatives and at negligible levels in high-yielding cultivated lines. In its genera-specific form of 22-nt, miR397 targeted mRNAs encoding laccases that decayed and induced robust secondary cascade silencing in wild species that required RNA-dependent RNA polymerase 6. In wild species of rice, reduced expression of laccases resulted in low lignification. As expected, overexpression of miR397 induced de-domestication phenotypes. At least 26 uncharacterized QTLs previously implicated in rice yield overlapped with laccases and miR397 genes. These results suggest that miRNAs contribute to rice domestication-associated phenotypes.

Transposons: a blessing curse

The genomes of most plant species are dominated by transposable elements (TEs). Once considered as 'junk DNA', TEs are now known to have a major role in driving genome evolution. Over the last decade, it has become apparent that some stress conditions and other environmental stimuli can drive bursts of activity of certain TE families and consequently new TE insertions. These can give rise to altered gene expression patterns and phenotypes, with new TE insertions sometimes causing flanking genes to become transcriptionally responsive to the same stress conditions that activated the TE in the first place. Such connections between TE-mediated increases in diversity and an accelerated rate of genome evolution provide powerful mechanisms for plants to adapt more rapidly to new environmental conditions. This review will focus on environmentally induced transposition, the mechanisms by which it alters gene expression, and the consequences for plant genome evolution and breeding.