An Arabidopsis hAT-like transposase is essential for plant development

Divergent molecular pathways govern temperature-dependent wheat stem rust resistance genes

The wheat stem rust pathogen Puccinia graminis f. sp. tritici (Pgt) causes severe crop losses worldwide. Several stem rust resistance (Sr) genes exhibit temperature-dependent immune responses. Sr6-mediated resistance is enhanced at lower temperatures, whereas Sr13 and Sr21 resistances are enhanced at higher temperatures. Here, we clone Sr6 using mutagenesis and resistance gene enrichment and sequencing (MutRenSeq), identifying it to encode a nucleotide-binding leucine-rich repeat (NLR) protein with an integrated BED domain. Sr6 temperature sensitivity is also transferred to wheat plants transformed with the Sr6 gene. Differential gene expression analysis of near-isogenic lines inoculated with Pgt at varying temperatures reveals that genes upregulated in the low-temperature-effective Sr6 response differ from those upregulated in the high-temperature-effective responses associated with Sr13 and Sr21. These findings highlight divergent molecular pathways involved in temperature-sensitive immunity and inform future strategies for deployment and engineering of genetic resistance in response to a changing climate.

The Silene latifolia genome and its giant Y chromosome

In many species with sex chromosomes, the Y is a tiny chromosome. However, the dioecious plant Silene latifolia has a giant ~550-megabase Y chromosome, which has remained unsequenced so far. We used a long- and short-read hybrid approach to obtain a high-quality male genome. Comparative analysis of the sex chromosomes with their homologs in outgroups showed that the Y is highly rearranged and degenerated. Recombination suppression between X and Y extended in several steps and triggered a massive accumulation of repeats on the Y as well as in the nonrecombining pericentromeric region of the X, leading to giant sex chromosomes. Using sex phenotype mutants, we identified candidate sex-determining genes on the Y in locations consistent with their favoring recombination suppression events 11 and 5 million years ago.

From the cauldron of conflict: Endogenous gene regulation by piRNA and other modes of adaptation enabled by selfish transposable elements

Transposable elements (TEs) provide a prime example of genetic conflict because they can proliferate in genomes and populations even if they harm the host. However, numerous studies have shown that TEs, though typically harmful, can also provide fuel for adaptation. This is because they code functional sequences that can be useful for the host in which they reside. In this review, I summarize the "how" and "why" of adaptation enabled by the genetic conflict between TEs and hosts. In addition, focusing on mechanisms of TE control by small piwi-interacting RNAs (piRNAs), I highlight an indirect form of adaptation enabled by conflict. In this case, mechanisms of host defense that regulate TEs have been redeployed for endogenous gene regulation. I propose that the genetic conflict released by meiosis in early eukaryotes may have been important because, among other reasons, it spurred evolutionary innovation on multiple interwoven trajectories - on the part of hosts and also embedded genetic parasites. This form of evolution may function as a complexity generating engine that was a critical player in eukaryotic evolution.

Genome-wide identification and characterization of SLEEPER, a transposon-derived gene family and their expression pattern in Brassica napus L

BACKGROUND

The transposons of the hAT superfamily are the most widespread transposons ever known. SLEEPER genes encode domesticated transposases from the hAT superfamily, which may have lost their transposable functions during long-term evolution and transformed into host proteins that regulate plant growth and development.

RESULTS

This study identified 162 members of the SLEEPER gene family from Brassica napus. These members are widely distributed on 19 chromosomes, mainly in the Cn subgenome, and have promoters with various cis-acting elements related to hormone regulation, abiotic stress, and growth and development regulation. Most of the genes in this family contain similar conserved domains and motifs, and the closer the genes are distributed on evolutionary branches, the more similar their structures are. Transcriptome sequencing performed on tissues at different growth stages from B. napus line 3529 indicated that these genes had different expression patterns, and nearly half of the genes were not detectably expressed in all samples.

CONCLUSIONS

This study investigated the gene structure, expression patterns, evolutionary features, and gene localization of the SLEEPER family members to confirm the significance of these genes in the growth of B. napus, providing a reference for the study of transposon domestication and outstanding genetic resources for the genetic improvement of B. napus.

Kinetochore and ionomic adaptation to whole-genome duplication in Cochlearia shows evolutionary convergence in three autopolyploids

Whole-genome duplication (WGD) occurs in all kingdoms and impacts speciation, domestication, and cancer outcome. However, doubled DNA management can be challenging for nascent polyploids. The study of within-species polyploidy (autopolyploidy) permits focus on this DNA management aspect, decoupling it from the confounding effects of hybridization (in allopolyploid hybrids). How is autopolyploidy tolerated, and how do young polyploids stabilize? Here, we introduce a powerful model to address this: the genus Cochlearia, which has experienced many polyploidization events. We assess meiosis and other polyploid-relevant phenotypes, generate a chromosome-scale genome, and sequence 113 individuals from 33 ploidy-contrasting populations. We detect an obvious autopolyploidy-associated selection signal at kinetochore components and ion transporters. Modeling the selected alleles, we detail evidence of the kinetochore complex mediating adaptation to polyploidy. We compare candidates in independent autopolyploids across three genera separated by 40 million years, highlighting a common function at the process and gene levels, indicating evolutionary flexibility in response to polyploidy.

Transcriptomic Analyses Reveal Insights into the Shared Regulatory Network of Phenolic Compounds and Steviol Glycosides in Stevia rebaudiana

Stevia rebaudiana (Bertoni) is a highly valuable crop for the steviol glycoside content in its leaves, which are no-calorie sweeteners hundreds of times more potent than sucrose. The presence of health-promoting phenolic compounds, particularly flavonoids, in the leaf of S. rebaudiana adds further nutritional value to this crop. Although all these secondary metabolites are highly desirable in S. rebaudiana leaves, the genes regulating the biosynthesis of phenolic compounds and the shared gene network between the regulation of biosynthesis of steviol glycosides and phenolic compounds still need to be investigated in this species. To identify putative candidate genes involved in the synergistic regulation of steviol glycosides and phenolic compounds, four genotypes with different contents of these compounds were selected for a pairwise comparison RNA-seq analysis, yielding 1136 differentially expressed genes. Genes that highly correlate with both steviol glycosides and phenolic compound accumulation in the four genotypes of S. rebaudiana were identified using the weighted gene co-expression network analysis. The presence of UDP-glycosyltransferases 76G1, 76H1, 85C1, and 91A1, and several genes associated with the phenylpropanoid pathway, including peroxidase, caffeoyl-CoA O-methyltransferase, and malonyl-coenzyme A:anthocyanin 3-O-glucoside-6″-O-malonyltransferase, along with 21 transcription factors like SCL3, WRK11, and MYB111, implied an extensive and synergistic regulatory network involved in enhancing the production of such compounds in S. rebaudiana leaves. In conclusion, this work identified a variety of putative candidate genes involved in the biosynthesis and regulation of particular steviol glycosides and phenolic compounds that will be useful in gene editing strategies for increasing and steering the production of such compounds in S. rebaudiana as well as in other species.

Genetic mapping and transcriptome profiling of a chickpea (Cicer arietinum L.) mutant identifies a novel locus (CaEl) regulating organ size and early vigor

Chickpea is among the top three legumes produced and consumed worldwide. Early plant vigor, characterized by good germination and rapid seedling growth, is an important agronomic trait in many crops including chickpea, and shows a positive correlation with seed size. In this study, we report a gamma-ray-induced chickpea mutant with a larger organ and seed size. The mutant (elm) exhibits increased early vigor and contains higher proline that contributes to a better tolerance under salt stress at germination, seedling, and early vegetative phase. The trait is governed as monogenic recessive, with wild-type allele being incompletely dominant over the mutant. Genetic mapping of this locus (CaEl) identified it as a previously uncharacterized gene (101503252) in chromosome 1 of the chickpea genome. There is a deletion of this gene in the mutant with a complete loss of expression. In silico analysis suggests that the gene is present as a single copy in chickpea and related legumes of the galegoid clade. In the mutant, cell division and expansion are affected. Transcriptome profiling identified differentially regulated transcripts related to cell division, expansion, cell wall organization, and metabolism in the mutant. The mutant can be exploited in chickpea breeding programs for increasing plant vigor and seed size.

Identification and characterization of RuvBL DNA helicase genes for tolerance against abiotic stresses in bread wheat (Triticum aestivum L.) and related species

Recombination UVB (sensitivity) like (RuvBL) helicase genes represent a conserved family of genes, which are known to be involved in providing tolerance against abiotic stresses like heat and drought. We identified nine wheat RuvBL genes, one each on nine different chromosomes, belonging to homoeologous groups 2, 3, and 4. The lengths of genes ranged from 1647 to 2197 bp and exhibited synteny with corresponding genes in related species including Ae. tauschii, Z. mays, O. sativa, H. vulgare, and B. distachyon. The gene sequences were associated with regulatory cis-elements and transposable elements. Two genes, namely TaRuvBL1a-4A and TaRuvBL1a-4B, also carried targets for a widely known miRNA, tae-miR164. Gene ontology revealed that these genes were closely associated with ATP-dependent formation of histone acetyltransferase complex. Analysis of the structure and function of RuvBL proteins revealed that the proteins were localized mainly in the cytoplasm. A representative gene, namely TaRuvBL1a-4A, was also shown to be involved in protein-protein interactions with ten other proteins. On the basis of phylogeny, RuvBL proteins were placed in two sub-divisions, namely RuvBL1 and RuvBL2, which were further classified into clusters and sub-clusters. In silico studies suggested that these genes were differentially expressed under heat/drought. The qRT-PCR analysis confirmed that expression of TaRuvBL genes differed among wheat cultivars, which differed in the level of thermotolerance. The present study advances our understanding of the biological role of wheat RuvBL genes and should help in planning future studies on RuvBL genes in wheat including use of RuvBL genes in breeding thermotolerant wheat cultivars.

Comparative Transcriptome Analysis Reveals Novel Candidate Resistance Genes Involved in Defence against Phytophthora cactorum in Strawberry

Crown rot, caused by Phytophthora cactorum, is a devastating disease of strawberry. While most commercial octoploid strawberry cultivars (Fragaria × ananassa Duch) are generally susceptible, the diploid species Fragaria vesca is a potential source of resistance genes to P. cactorum. We previously reported several F. vesca genotypes with varying degrees of resistance to P. cactorum. To gain insights into the strawberry defence mechanisms, comparative transcriptome profiles of two resistant genotypes (NCGR1603 and Bukammen) and a susceptible genotype (NCGR1218) of F. vesca were analysed by RNA-Seq after wounding and subsequent inoculation with P. cactorum. Differential gene expression analysis identified several defence-related genes that are highly expressed in the resistant genotypes relative to the susceptible genotype in response to P. cactorum after wounding. These included putative disease resistance (R) genes encoding receptor-like proteins, receptor-like kinases, nucleotide-binding sites, leucine-rich repeat proteins, RPW8-type disease resistance proteins, and 'pathogenesis-related protein 1'. Seven of these R-genes were expressed only in the resistant genotypes and not in the susceptible genotype, and these appeared to be present only in the genomes of the resistant genotypes, as confirmed by PCR analysis. We previously reported a single major gene locus RPc-1 (Resistance to Phytophthora cactorum 1) in F. vesca that contributed resistance to P. cactorum. Here, we report that 4-5% of the genes (35-38 of ca 800 genes) in the RPc-1 locus are differentially expressed in the resistant genotypes compared to the susceptible genotype after inoculation with P. cactorum. In particular, we identified three defence-related genes encoding wall-associated receptor-like kinase 3, receptor-like protein 12, and non-specific lipid-transfer protein 1-like that were highly expressed in the resistant genotypes compared to the susceptible one. The present study reports several novel candidate disease resistance genes that warrant further investigation for their role in plant defence against P. cactorum.

PANTOTHENATE KINASE4, LOSS OF GDU2, and TRANSPOSON PROTEIN1 affect the canalization of tomato fruit metabolism

Most studies investigating quantitative traits focus on mean levels per genotype rather than the variation between different individuals of one genotype or the variation elicited by different environments. Consequently, the genes that govern this effect are not well understood. The concept, named canalization, which describes a lack of variation, is well known in the context of developmental processes but is poorly studied for quantitative traits such as metabolism. In this study, we selected 8 putative candidate genes from previously identified canalized metabolic quantitative trait loci and created genome-edited tomato (Solanum lycopersicum) mutants of these genes for experimental validation. Most lines showed wild-type morphology, except for an ARF-like protein mutant showing aberrant phenotypes in the form of scarred fruit cuticles. In greenhouse trials with different irrigation conditions, whole-plant traits showed a general increase of their level toward the more optimal irrigation conditions, whereas most metabolic traits showed an increase toward the other end of the gradient. Mutants of a PANTOTHENATE KINASE 4, the AIRP ubiquitin gene LOSS OF GDU2, and TRANSPOSON PROTEIN 1 grown under these conditions showed an overall improved plant performance. Additional effects, on both target and other metabolites in tomato fruits, regarding the mean level at specific conditions and, ergo, the cross-environment coefficient of variation, were observed. However, variation between individuals remained unaffected. In conclusion, this study supports the idea of distinct sets of genes regulating different types of variation.

Genome-wide identification and characterization of exapted transposable elements in the large genome of sunflower (Helianthus annuus L.)

Transposable elements (TEs) are an important source of genome variability, playing many roles in the evolution of eukaryotic species. Besides well-known phenomena, TEs may undergo the exaptation process and generate the so-called exapted transposable element genes (ETEs). Here we present a genome-wide survey of ETEs in the large genome of sunflower (Helianthus annuus L.), in which the massive amount of TEs, provides a significant source for exaptation. A library of sunflower TEs was used to build TE-specific Hidden Markov Model profiles, to search for all available sunflower gene products. In doing so, 20 016 putative ETEs were identified and further investigated for the characteristics that distinguish TEs from genes, leading to the validation of 3530 ETEs. The analysis of ETEs transcription patterns under different stress conditions showed a differential regulation triggered by treatments mimicking biotic and abiotic stress; furthermore, the distribution of functional domains of differentially regulated ETEs revealed a relevant presence of domains involved in many aspects of cellular functions. A comparative genomic investigation was performed including species representative of Asterids and appropriate outgroups: the bulk of ETEs that resulted were specific to the sunflower, while few ETEs presented orthologues in the genome of all analyzed species, making the hypothesis of a conserved function. This study highlights the crucial role played by exaptation, actively contributing to species evolution.

Genomic analyses provide insights into the evolution and salinity adaptation of halophyte Tamarix chinensis

BACKGROUND

The woody halophyte Tamarix chinensis is a pioneer tree species in the coastal wetland ecosystem of northern China, exhibiting high resistance to salt stress. However, the genetic information underlying salt tolerance in T. chinensis remains to be seen. Here we present a genomic investigation of T. chinensis to elucidate the underlying mechanism of its high resistance to salinity.

RESULTS

Using a combination of PacBio and high-throughput chromosome conformation capture data, a chromosome-level T. chinensis genome was assembled with a size of 1.32 Gb and scaffold N50 of 110.03 Mb. Genome evolution analyses revealed that T. chinensis significantly expanded families of HAT and LIMYB genes. Whole-genome and tandem duplications contributed to the expansion of genes associated with the salinity adaptation of T. chinensis. Transcriptome analyses were performed on root and shoot tissues during salt stress and recovery, and several hub genes responding to salt stress were identified. WRKY33/40, MPK3/4, and XBAT31 were critical in responding to salt stress during early exposure, while WRKY40, ZAT10, AHK4, IRX9, and CESA4/8 were involved in responding to salt stress during late stress and recovery. In addition, PER7/27/57/73 encoding class III peroxidase and MCM3/4/5/7 encoding DNA replication licensing factor maintained up/downregulation during salt stress and recovery stages.

CONCLUSIONS

The results presented here reveal the genetic mechanisms underlying salt adaptation in T. chinensis, thus providing important genomic resources for evolutionary studies on tamarisk and plant salt tolerance genetic improvement.

Multiomic Data Integration in the Analysis of Drought-Responsive Mechanisms in Quercus ilex Seedlings

The integrated analysis of different omic layers can provide new knowledge not provided by their individual analysis. This approach is also necessary to validate data and reveal post-transcriptional and post-translational mechanisms of gene expression regulation. In this work, we validated the possibility of applying this approach to non-model species such as Quercus ilex. Transcriptomics, proteomics, and metabolomics from Q. ilex seedlings subjected to drought-like conditions under the typical summer conditions in southern Spain were integrated using a non-targeted approach. Two integrative approaches, PCA and DIABLO, were used and compared. Both approaches seek to reduce dimensionality, preserving the maximum information. DIABLO also allows one to infer interconnections between the different omic layers. For easy visualization and analysis, these interconnections were analyzed using functional and statistical networks. We were able to validate results obtained by analyzing the omic layers separately. We identified the importance of protein homeostasis with numerous protease and chaperones in the networks. We also discovered new key processes, such as transcriptional control, and identified the key function of transcription factors, such as DREB2A, WRKY65, and CONSTANS, in the early response to drought.

Integrated Metabolites and Transcriptomics at Different Growth Stages Reveal Polysaccharide and Flavonoid Biosynthesis in Cynomorium songaricum

Cynomorium songaricum is a perennial parasitic herb, and its stem is widely used as a traditional Chinese medicine, which largely relies on bioactive compounds (e.g., polysaccharides, flavonoids, and triterpenes). To date, although the optimum harvest time of stems has been demonstrated at the unearthed stage (namely the early flowering stage, EFS), the accumulation mechanism of polysaccharides and flavonoids during growth stages is still limited. In this study, the physiological characteristics (stem fresh weight, contents of soluble sugar and flavonoids, and antioxidant capacity) at four different growth stages (germination stage (GS), vegetative growth stage (VGS), EFS, and flowering stage (FS)) were determined, transcriptomics were analyzed by illumina sequencing, and expression levels of key genes were validated by qRT-PCR at the GS, VGS, and EFS. The results show that the stem biomass, soluble sugar and total flavonoids contents, and antioxidant capacity peaked at EFS compared with GS, VGS, and FS. A total of 6098 and 13,023 differentially expressed genes (DEGs) were observed at VGS and EFS vs. GS, respectively, with 367 genes co-expressed. Based on their biological functions, 109 genes were directly involved in polysaccharide and flavonoid biosynthesis as well as growth and development. The expression levels of key genes involved in polysaccharides (e.g., GLCs, XTHs and PMEs), flavonoids (e.g., 4CLLs, CYPs and UGTs), growth and development (e.g., AC58, TCPs and AP1), hormones biosynthesis and signaling (e.g., YUC8, AIPT and ACO1), and transcription factors (e.g., MYBs, bHLHs and WRKYs) were in accordance with changes of physiological characteristics. The combinational analysis of metabolites with transcriptomics provides insight into the mechanism of polysaccharide and flavonoid biosynthesis in C. songaricum during growth stages.

Genome-wide association analysis of 101 accessions dissects the genetic basis of shell thickness for genetic improvement in Persian walnut (Juglans regia L.)

BACKGROUND

Understanding the underlying genetic mechanisms that drive phenotypic variations is essential for enhancing the efficacy of crop improvement. Persian walnut (Juglans regia L.), which is grown extensively worldwide, is an important economic tree fruit due to its horticultural, medicinal, and material value. The quality of the walnut fruit is related to the selection of traits such as thinner shells, larger filling rates, and better taste, which is very important for breeding in China. The complex quantitative fruit-related traits are influenced by a variety of physiological and environmental factors, which can vary widely between walnut genotypes.

RESULTS

For this study, a set of 101 Persian walnut accessions were re-sequenced, which generated a total of 906.2 Gb of Illumina sequence data with an average read depth of 13.8× for each accession. We performed the genome-wide association study (GWAS) using 10.9 Mb of high-quality single-nucleotide polymorphisms (SNPs) and 10 agronomic traits to explore the underlying genetic basis of the walnut fruit. Several candidate genes are proposed to be involved in walnut characteristics, including JrPXC1, JrWAKL8, JrGAMYB, and JrFRK1. Specifically, the JrPXC1 gene was confirmed to participate in the regulation of secondary wall cellulose thickening in the walnut shell.

CONCLUSION

In addition to providing considerable available genetic resources for walnut trees, this study revealed the underlying genetic basis involved in important walnut agronomic traits, particularly shell thickness, as well as providing clues for the improvement of genetic breeding and domestication in other perennial economic crops.

The genome of Aechmea fasciata provides insights into the evolution of tank epiphytic habits and ethylene-induced flowering

Aechmea fasciata is one of the most popular bromeliads and bears a water-impounding tank with a vase-like rosette. The tank habit is a key innovation that has promoted diversity among bromeliads. To reveal the genomic basis of tank habit formation and ethylene-induced flowering, we sequenced the genome of A. fasciata and assembled 352 Mb of sequences into 24 chromosomes. Comparative genomic analysis showed that the chromosomes experienced at least two fissions and two fusions from the ancestral genome of A. fasciata and Ananas comosus. The gibberellin receptor gene GID1C-like was duplicated by a segmental duplication event. This duplication may affect GA signalling and promote rosette expansion, which may permit water-impounding tank formation. During ethylene-induced flowering, AfFTL2 expression is induced and targets the EIN3 binding site 'ATGTAC' by AfEIL1-like. The data provided here will serve as an important resource for studying the evolution and mechanisms underlying flowering time regulation in bromeliads.

A genome-wide comparative evolutionary analysis of zinc finger-BED transcription factor genes in land plants

Zinc finger (Zf)-BED proteins are a novel superfamily of transcription factors that controls numerous activities in plants including growth, development, and cellular responses to biotic and abiotic stresses. Despite their important roles in gene regulation, little is known about the specific functions of Zf-BEDs in land plants. The current study identified a total of 750 Zf-BED-encoding genes in 35 land plant species including mosses, bryophytes, lycophytes, gymnosperms, and angiosperms. The gene family size was somewhat proportional to genome size. All identified genes were categorized into 22 classes based on their specific domain architectures. Of these, class I (Zf-BED_DUF-domain_Dimer_Tnp_hAT) was the most common in the majority of the land plants. However, some classes were family-specific, while the others were species-specific, demonstrating diversity at different classification levels. In addition, several novel functional domains were also predicated including WRKY and nucleotide-binding site (NBS). Comparative genomics, transcriptomics, and proteomics provided insights into the evolutionary history, duplication, divergence, gene gain and loss, species relationship, expression profiling, and structural diversity of Zf-BEDs in land plants. The comprehensive study of Zf-BEDs in Gossypium sp., (cotton) also demonstrated a clear footprint of polyploidization. Overall, this comprehensive evolutionary study of Zf-BEDs in land plants highlighted significant diversity among plant species.

High-Density Mapping of Quantitative Trait Loci Controlling Agronomically Important Traits in Quinoa (Chenopodium quinoa Willd.)

Quinoa is a pseudocereal originating from the Andean regions. Despite quinoa's long cultivation history, genetic analysis of this crop is still in its infancy. We aimed to localize quantitative trait loci (QTL) contributing to the phenotypic variation of agronomically important traits. We crossed the Chilean accession PI-614889 and the Peruvian accession CHEN-109, which depicted significant differences in days to flowering, days to maturity, plant height, panicle length, and thousand kernel weight (TKW), saponin content, and mildew susceptibility. We observed sizeable phenotypic variation across F2 plants and F3 families grown in the greenhouse and the field, respectively. We used Skim-seq to genotype the F2 population and constructed a high-density genetic map with 133,923 single nucleotide polymorphism (SNPs). Fifteen QTL were found for ten traits. Two significant QTL, common in F2 and F3 generations, depicted pleiotropy for days to flowering, plant height, and TKW. The pleiotropic QTL harbored several putative candidate genes involved in photoperiod response and flowering time regulation. This study presents the first high-density genetic map of quinoa that incorporates QTL for several important agronomical traits. The pleiotropic loci can facilitate marker-assisted selection in quinoa breeding programs.

Recurrent mutations promote widespread structural and functional divergence of MULE-derived genes in plants

Transposable element (TE)-derived genes are increasingly recognized as major sources conferring essential traits in agriculturally important crops but underlying evolutionary mechanisms remain obscure. We updated previous annotations and constructed 18,744 FAR-RED IMPAIRED RESPONSE1 (FAR1) genes, a transcription factor family derived from Mutator-like elements (MULEs), from 80 plant species, including 15,546 genes omitted in previous annotations. In-depth sequence comparison of the updated gene repertoire revealed that FAR1 genes underwent continuous structural divergence via frameshift and nonsense mutations that caused premature translation termination or specific domain truncations. CRISPR/Cas9-based genome editing and transcriptome analysis determined a novel gene involved in fertility-regulating transcription of rice pollen, denoting the functional capacity of our re-annotated gene models especially in monocots which had the highest copy numbers. Genomic evidence showed that the functional gene adapted by obtaining a shortened form through a frameshift mutation caused by a tandem duplication of a 79-bp sequence resulting in premature translation termination. Our findings provide improved resources for comprehensive studies of FAR1 genes with beneficial agricultural traits and unveil novel evolutionary mechanisms generating structural divergence and subsequent adaptation of TE-derived genes in plants.

Impact of transposable elements on the evolution of complex living systems and their epigenetic control

Transposable elements (TEs) contribute to genomic innovations, as well as genome instability, across a wide variety of species. Popular designations such as 'selfish DNA' and 'junk DNA,' common in the 1980s, may be either inaccurate or misleading, while a more enlightened view of the TE-host relationship covers a range from parasitism to mutualism. Both plant and animal hosts have evolved epigenetic mechanisms to reduce the impact of TEs, both by directly silencing them and by reducing their ability to transpose in the genome. However, TEs have also been co-opted by both plant and animal genomes to perform a variety of physiological functions, ranging from TE-derived proteins acting directly in normal biological functions to innovations in transcription factor activity and also influencing gene expression. Their presence, in fact, can affect a range of features at genome, phenotype, and population levels. The impact TEs have had on evolution is multifaceted, and many aspects still remain unexplored. In this review, the epigenetic control of TEs is contextualized according to the evolution of complex living systems.

The Evolutionary Volte-Face of Transposable Elements: From Harmful Jumping Genes to Major Drivers of Genetic Innovation

Transposable elements (TEs) are self-replicating DNA elements that constitute major fractions of eukaryote genomes. Their ability to transpose can modify the genome structure with potentially deleterious effects. To repress TE activity, host cells have developed numerous strategies, including epigenetic pathways, such as DNA methylation or histone modifications. Although TE neo-insertions are mostly deleterious or neutral, they can become advantageous for the host under specific circumstances. The phenomenon leading to the appropriation of TE-derived sequences by the host is known as TE exaptation or co-option. TE exaptation can be of different natures, through the production of coding or non-coding DNA sequences with ultimately an adaptive benefit for the host. In this review, we first give new insights into the silencing pathways controlling TE activity. We then discuss a model to explain how, under specific environmental conditions, TEs are unleashed, leading to a TE burst and neo-insertions, with potential benefits for the host. Finally, we review our current knowledge of coding and non-coding TE exaptation by providing several examples in various organisms and describing a method to identify TE co-option events.

A domesticated Harbinger transposase forms a complex with HDA6 and promotes histone H3 deacetylation at genes but not TEs in Arabidopsis

In eukaryotes, histone acetylation is a major modification on histone N-terminal tails that is tightly connected to transcriptional activation. HDA6 is a histone deacetylase involved in the transcriptional regulation of genes and transposable elements (TEs) in Arabidopsis thaliana. HDA6 has been shown to participate in several complexes in plants, including a conserved SIN3 complex. Here, we uncover a novel protein complex containing HDA6, several Harbinger transposon-derived proteins (HHP1, SANT1, SANT2, SANT3, and SANT4), and MBD domain-containing proteins (MBD1, MBD2, and MBD4). We show that mutations of all four SANT genes in the sant-null mutant cause increased expression of the flowering repressors FLC, MAF4, and MAF5, resulting in a late flowering phenotype. Transcriptome deep sequencing reveals that while the SANT proteins and HDA6 regulate the expression of largely overlapping sets of genes, TE silencing is unaffected in sant-null mutants. Our global histone H3 acetylation profiling shows that SANT proteins and HDA6 modulate gene expression through deacetylation. Collectively, our findings suggest that Harbinger transposon-derived SANT domain-containing proteins are required for histone deacetylation and flowering time control in plants.

Unearthing LTR Retrotransposon gag Genes Co-opted in the Deep Evolution of Eukaryotes

LTR retrotransposons comprise a major component of the genomes of eukaryotes. On occasion, retrotransposon genes can be recruited by their hosts for diverse functions, a process formally referred to as co-option. However, a comprehensive picture of LTR retrotransposon gag gene co-option in eukaryotes is still lacking, with several documented cases exclusively involving Ty3/Gypsy retrotransposons in animals. Here, we use a phylogenomic approach to systemically unearth co-option of retrotransposon gag genes above the family level of taxonomy in 2,011 eukaryotes, namely co-option occurring during the deep evolution of eukaryotes. We identify a total of 14 independent gag gene co-option events across more than 740 eukaryote families, eight of which have not been reported previously. Among these retrotransposon gag gene co-option events, nine, four, and one involve gag genes of Ty3/Gypsy, Ty1/Copia, and Bel-Pao retrotransposons, respectively. Seven, four, and three co-option events occurred in animals, plants, and fungi, respectively. Interestingly, two co-option events took place in the early evolution of angiosperms. Both selective pressure and gene expression analyses further support that these co-opted gag genes might perform diverse cellular functions in their hosts, and several co-opted gag genes might be subject to positive selection. Taken together, our results provide a comprehensive picture of LTR retrotransposon gag gene co-option events that occurred during the deep evolution of eukaryotes and suggest paucity of LTR retrotransposon gag gene co-option during the deep evolution of eukaryotes.

Characteristics of the AT-Hook Motif Containing Nuclear Localized (AHL) Genes in Carrot Provides Insight into Their Role in Plant Growth and Storage Root Development

The AT-hook motif containing nuclear localized (AHL) gene family, controlling various developmental processes, is conserved in land plants. They comprise Plant and Prokaryote Conserved (PPC) domain and one or two AT-hook motifs. DcAHLc1 has been proposed as a candidate gene governing the formation of the carrot storage root. We identified and in-silico characterized carrot AHL proteins, performed phylogenetic analyses, investigated their expression profiles and constructed gene coexpression networks. We found 47 AHL genes in carrot and grouped them into two clades, A and B, comprising 29 and 18 genes, respectively. Within Clade-A, we distinguished three subclades, one of them grouping noncanonical AHLs differing in their structure (two PPC domains) and/or cellular localization (not nucleus). Coexpression network analysis attributed AHLs expressed in carrot roots into four of the 72 clusters, some of them showing a large number of interactions. Determination of expression profiles of AHL genes in various tissues and samples provided basis to hypothesize on their possible roles in the development of the carrot storage root. We identified a group of rapidly evolving noncanonical AHLs, possibly differing functionally from typical AHLs, as suggested by their expression profiles and their predicted cellular localization. We pointed at several AHLs likely involved in the development of the carrot storage root.

When junk DNA turns functional: transposon-derived non-coding RNAs in plants

Transposable elements (TEs) are major contributors to genome complexity in eukaryotes. TE mobilization may cause genome instability, although it can also drive genome diversity throughout evolution. TE transposition may influence the transcriptional activity of neighboring genes by modulating the epigenomic profile of the region or by altering the relative position of regulatory elements. Notably, TEs have emerged in the last few years as an important source of functional long and small non-coding RNAs. A plethora of small RNAs derived from TEs have been linked to the trans regulation of gene activity at the transcriptional and post-transcriptional levels. Furthermore, TE-derived long non-coding RNAs have been shown to modulate gene expression by interacting with protein partners, sequestering active small RNAs, and forming duplexes with DNA or other RNA molecules. In this review, we summarize our current knowledge of the functional and mechanistic paradigms of TE-derived long and small non-coding RNAs and discuss their role in plant development and evolution.

Genomic insights into the origin, domestication and genetic basis of agronomic traits of castor bean

BACKGROUND

Castor bean (Ricinus communis L.) is an important oil crop, which belongs to the Euphorbiaceae family. The seed oil of castor bean is currently the only commercial source of ricinoleic acid that can be used for producing about 2000 industrial products. However, it remains largely unknown regarding the origin, domestication, and the genetic basis of key traits of castor bean.

RESULTS

Here we perform a de novo chromosome-level genome assembly of the wild progenitor of castor bean. By resequencing and analyzing 505 worldwide accessions, we reveal that the accessions from East Africa are the extant wild progenitors of castor bean, and the domestication occurs ~ 3200 years ago. We demonstrate that significant genetic differentiation between wild populations in Kenya and Ethiopia is associated with past climate fluctuation in the Turkana depression ~ 7000 years ago. This dramatic change in climate may have caused the genetic bottleneck in wild castor bean populations. By a genome-wide association study, combined with quantitative trait locus analysis, we identify important candidate genes associated with plant architecture and seed size.

CONCLUSIONS

This study provides novel insights of domestication and genome evolution of castor bean, which facilitates genomics-based breeding of this important oilseed crop and potentially other tree-like crops in future.

Connecting the dots: from nanodomains to physiological functions of REMORINs

REMORINs (REMs) are a plant-specific protein family, proposed regulators of membrane-associated molecular assemblies and well-established markers of plasma membrane nanodomains. REMs play a diverse set of functions in plant interactions with pathogens and symbionts, responses to abiotic stresses, hormone signaling and cell-to-cell communication. In this review, we highlight the established and more putative roles of REMs throughout the literature. We discuss the physiological functions of REMs, the mechanisms underlying their nanodomain-organization and their putative role as regulators of nanodomain-associated molecular assemblies. Furthermore, we discuss how REM phosphorylation may regulate their functional versatility. Overall, through data-mining and comparative analysis of the literature, we suggest how to further study the molecular mechanisms underpinning the functions of REMs.

Novelty and convergence in adaptation to whole genome duplication

Whole genome duplication (WGD) can promote adaptation but is disruptive to conserved processes, especially meiosis. Studies in Arabidopsis arenosa revealed a coordinated evolutionary response to WGD involving interacting proteins controlling meiotic crossovers, which are minimised in an autotetraploid (within-species polyploid) to avoid mis-segregation. Here we test whether this surprising flexibility of a conserved essential process, meiosis, is recapitulated in an independent WGD system, Cardamine amara, 17 million years diverged from A. arenosa. We assess meiotic stability and perform population-based scans for positive selection, contrasting the genomic response to WGD in C. amara with that of A. arenosa. We found in C. amara the strongest selection signals at genes with predicted functions thought important to adaptation to WGD: meiosis, chromosome remodelling, cell cycle, and ion transport. However, genomic responses to WGD in the two species differ: minimal ortholog-level convergence emerged, with none of the meiosis genes found in A. arenosa exhibiting strong signal in C. amara. This is consistent with our observations of lower meiotic stability and occasional clonal spreading in diploid C. amara, suggesting that nascent C. amara autotetraploid lineages were preadapted by their diploid lifestyle to survive while enduring reduced meiotic fidelity. However, in contrast to a lack of ortholog convergence, we see process-level and network convergence in DNA management, chromosome organisation, stress signalling, and ion homeostasis processes. This gives the first insight into the salient adaptations required to meet the challenges of a WGD state and shows that autopolyploids can utilize multiple evolutionary trajectories to adapt to WGD.

Phenotype study of multifoliolate leaf formation in Trifolium alexandrinum L

Background

The genus Trifolium is characterized by typical trifoliolate leaves. Alterations in leaf formats from trifoliolate to multifoliolate, i.e., individual plants bearing trifoliolate, quadrifoliolate, pentafoliolate or more leaflets, were previously reported among many species of the genus. The study is an attempt to develop pure pentafoliolate plants of T. alexandrinum and to understand its genetic control.

Methods

The experimental material consisted of two populations of T. alexandrinum with multifoliolate leaf expression, i.e.,interspecific hybrid progenies of T. alexandrinum with T. apertum, and T. alexandrinum genotype Penta-1. Penetrance of the multifoliolate trait was observed among multifoliolate and trifoliolate plant progenies. In vitro culture and regeneration of plantlets from the axillary buds from different plant sources was also attempted.

Results

The inheritance among a large number of plant progenies together with in vitro micro-propagation results did not establish a definite pattern. The multifoliolate leaf formation was of chimeric nature, i.e., more than one leaf format appearing on individual branches. Reversal to normal trifoliolate from multifoliolate was also quite common. Penetrance and expression of multifoliolate leaf formation was higher among the plants raised from multifoliolate plants. Multifoliolate and pure pentafoliolate plants were observed in the progenies of pure trifoliolate plants and vice-versa. There was an apparent increase in the pentafoliolate leaf formation frequency over the years due to targeted selection. A few progenies of the complete pentafoliolate plants in the first year were true breeding in the second year. Frequency of plantlets with multifoliolate leaf formation was also higher in in vitro axillary bud multiplication when the explant bud was excised from the multifoliolate leaf node.

Conclusion

Number of leaflets being a discrete variable, occurrence of multifoliolate leaves on individual branches, reversal of leaf formats on branches and developing true breeding pentafoliolates were the factors leading to a hypothesis beyond normal Mendelian inheritance. Transposable elements (TEs) involved in leaf development in combination with epigenetics were probably responsible for alterations in the expression of leaflet number. Putative TE’s movement owing to chromosomal rearrangements possibly resulted in homozygous pentafoliolate trait with evolutionary significance. The hypothesis provides a new insight into understanding the genetic control of this trait in T. alexandrinum and may also be useful in other Trifolium species where such observations are reported.

Genome-Wide Association Study Dissects Resistance Loci against Bacterial Blight in a Diverse Rice Panel from the 3000 Rice Genomes Project

BACKGROUND

Bacterial blight (BB), caused by Xanthomonas oryzae pv. oryzae (Xoo) is one of the most devastating bacterial diseases of rice in temperate and tropical regions. Breeding and deployment of resistant cultivars carrying major resistance (R) genes has been the most effective approach for BB management. However, because of specific interaction of each R gene with the product of the corresponding pathogen avirulence or effector gene, new pathogen strains that can overcome the deployed resistance often emerge rapidly. To deal with ever-evolving Xoo, it is necessary to identify novel R genes and resistance quantitative trait loci (QTL).

RESULTS

BB resistance of a diverse panel of 340 accessions from the 3000 Rice Genomes Project (3 K RGP) was evaluated by artificial inoculation with four representative Xoo strains, namely Z173 (C4), GD1358 (C5), V from China and PXO339 (P9a) from Philippines. Using the 3 K RG 4.8mio filtered SNP Dataset, a total of 11 QTL associated with BB resistance on chromosomes 4, 5, 11 and 12 were identified through a genome-wide association study (GWAS). Among them, eight resistance loci, which were narrowed down to relatively small genomic intervals, coincided with previously reported QTL or R genes, e.g. xa5, xa25, xa44(t). The other three QTL were putative novel loci associated with BB resistance. Linear regression analysis showed a dependence of BB lesion length on the number of favorable alleles, suggesting that pyramiding QTL using marker-assisted selection would be an effective approach for improving resistance. In addition, the Hap2 allele of LOC_Os11g46250 underlying qC5-11.1 was validated as positively regulating resistance against strain C5.

CONCLUSIONS

Our findings provide valuable information for the genetic improvement of BB resistance and application of germplasm resources in rice breeding programs.

In Silico Identification of Stress-Associated Transposable Elements in Arabidopsis thaliana Using Public Transcriptome Data

Transposable elements (TEs) have been associated with stress response in many plants, making them a key target of study. However, the high variability, genomic repeat-heavy nature, and widely noncoding character of TEs have made them difficult to study using non-specialized methods, whether experimental or computational. In this chapter, we introduce two computational workflows to analyze transposable elements using publicly available transcriptome data. In the first of these methods, we identify TEs, which show differential expression under salt stress using sample transcriptome libraries that includes noncoding transcripts. In the second, we identify protein-coding genes with differential expression under the same conditions, and determine which TEs are enriched in the promoter regions of these stress-related genes.

Comparative Genomics and Functional Studies of Wheat BED-NLR Loci

Nucleotide-binding leucine-rich-repeat (LRR) receptors (NLRs) with non-canonical integrated domains (NLR-IDs) are widespread in plant genomes. Zinc-finger BED (named after the Drosophila proteins Boundary Element-Associated Factor and DNA Replication-related Element binding Factor, named BED hereafter) are among the most frequently found IDs. Five BED-NLRs conferring resistance against bacterial and fungal pathogens have been characterized. However, it is unknown whether BED-NLRs function in a manner similar to other NLR-IDs. Here, we used chromosome-level assemblies of wheat to explore the Yr7 and Yr5a genomic regions and show that, unlike known NLR-ID loci, there is no evidence for a NLR-partner in their vicinity. Using neighbor-network analyses, we observed that BED domains from BED-NLRs share more similarities with BED domains from single-BED proteins and from BED-containing proteins harboring domains that are conserved in transposases. We identified a nuclear localization signal (NLS) in Yr7, Yr5, and the other characterized BED-NLRs. We thus propose that this is a feature of BED-NLRs that confer resistance to plant pathogens. We show that the NLS was functional in truncated versions of the Yr7 protein when expressed in N. benthamiana. We did not observe cell-death upon the overexpression of Yr7 full-length, truncated, and 'MHD' variants in N. benthamiana. This suggests that either this system is not suitable to study BED-NLR signaling or that BED-NLRs require additional components to trigger cell death. These results define novel future directions to further understand the role of BED domains in BED-NLR mediated resistance.

The Rice DNA-Binding Protein ZBED Controls Stress Regulators and Maintains Disease Resistance After a Mild Drought

BACKGROUND

Identifying new sources of disease resistance and the corresponding underlying resistance mechanisms remains very challenging, particularly in Monocots. Moreover, the modification of most disease resistance pathways made so far is detrimental to tolerance to abiotic stresses such as drought. This is largely due to negative cross-talks between disease resistance and abiotic stress tolerance signaling pathways. We have previously described the role of the rice ZBED protein containing three Zn-finger BED domains in disease resistance against the fungal pathogen Magnaporthe oryzae. The molecular and biological functions of such BED domains in plant proteins remain elusive.

RESULTS

Using Nicotiana benthamiana as a heterologous system, we show that ZBED localizes in the nucleus, binds DNA, and triggers basal immunity. These activities require conserved cysteine residues of the Zn-finger BED domains that are involved in DNA binding. Interestingly, ZBED overexpressor rice lines show increased drought tolerance. More importantly, the disease resistance response conferred by ZBED is not compromised by drought-induced stress.

CONCLUSIONS

Together our data indicate that ZBED might represent a new type of transcriptional regulator playing simultaneously a positive role in both disease resistance and drought tolerance. We demonstrate that it is possible to provide disease resistance and drought resistance simultaneously.

Twenty years of transposable element analysis in the Arabidopsis thaliana genome

Transposable elements (TEs) are mobile repetitive DNA sequences shown to be major drivers of genome evolution. As the first plant to have its genome sequenced and analyzed at the genomic scale, Arabidopsis thaliana has largely contributed to our TE knowledge.

The present report describes 20 years of accumulated TE knowledge gained through the study of the Arabidopsis genome and covers the known TE families, their relative abundance, and their genomic distribution. It presents our knowledge of the different TE family activities, mobility, population and long-term evolutionary dynamics. Finally, the role of TE as substrates for new genes and their impact on gene expression is illustrated through a few selected demonstrative cases. Promising future directions for TE studies in this species conclude the review.

Characterization of B-Genome Specific High Copy hAT MITE Families in Brassica nigra Genome

Miniature inverted-repeat transposable elements (MITEs) are non-autonomous class II transposons which have been shown to influence genome evolution. Brassica nigra L. (B-genome) is one of three Brassica diploids cultivated primarily as an oil crop, which harbors novel alleles important for breeding. Two new high copy hAT MITE families (BniHAT-1 and BniHAT-2) from the B-genome were characterized and their prevalence assessed in the genomes of the related diploids, rapa L. (A) and Brassica oleracea L. (C). Both novel MITE families were present at high copy numbers in the B-genome with 434 and 331 copies of BniHAT-1 and BniHAT-2, respectively. Yet less than 20 elements were identified in the genome assemblies of the A, and C -genomes, supporting B-genome specific proliferation of these MITE families. Although apparently randomly distributed across the genome, 68 and 70% of the B-genome MITEs were present within 2 kb flanking regions of annotated genes suggesting they might influence gene expression and/or function. In addition, MITE derived microRNAs and transcription factor binding sites suggested a putative role in gene regulation. Age of insertion analysis revealed that the major proliferation of these elements occurred during 2-3 million years ago. Additionally, site-specific polymorphism analyses showed that 44% MITEs were undergoing active amplification into the B-genome. Overall, this study provides a comprehensive analysis of two high copy MITE families, which were specifically amplified in the B-genome, suggesting a potential role in shaping the Brassica B-genome.

Genome assembly and characterization of a complex zfBED-NLR gene-containing disease resistance locus in Carolina Gold Select rice with Nanopore sequencing

Long-read sequencing facilitates assembly of complex genomic regions. In plants, loci containing nucleotide-binding, leucine-rich repeat (NLR) disease resistance genes are an important example of such regions. NLR genes constitute one of the largest gene families in plants and are often clustered, evolving via duplication, contraction, and transposition. We recently mapped the Xo1 locus for resistance to bacterial blight and bacterial leaf streak, found in the American heirloom rice variety Carolina Gold Select, to a region that in the Nipponbare reference genome is NLR gene-rich. Here, toward identification of the Xo1 gene, we combined Nanopore and Illumina reads and generated a high-quality Carolina Gold Select genome assembly. We identified 529 complete or partial NLR genes and discovered, relative to Nipponbare, an expansion of NLR genes at the Xo1 locus. One of these has high sequence similarity to the cloned, functionally similar Xa1 gene. Both harbor an integrated zfBED domain, and the repeats within each protein are nearly perfect. Across diverse Oryzeae, we identified two sub-clades of NLR genes with these features, varying in the presence of the zfBED domain and the number of repeats. The Carolina Gold Select genome assembly also uncovered at the Xo1 locus a rice blast resistance gene and a gene encoding a polyphenol oxidase (PPO). PPO activity has been used as a marker for blast resistance at the locus in some varieties; however, the Carolina Gold Select sequence revealed a loss-of-function mutation in the PPO gene that breaks this association. Our results demonstrate that whole genome sequencing combining Nanopore and Illumina reads effectively resolves NLR gene loci. Our identification of an Xo1 candidate is an important step toward mechanistic characterization, including the role(s) of the zfBED domain. Finally, the Carolina Gold Select genome assembly will facilitate identification of other useful traits in this historically important variety.

Plants lack adaptive immunity, and instead contain repeat-rich, disease resistance genes that evolve rapidly through duplication, recombination, and transposition. The number, variation, and often clustered arrangement of these genes make them challenging to sequence and catalog. The US heirloom rice variety Carolina Gold Select has resistance to two important bacterial diseases. Toward identifying the responsible gene(s), we combined long- and short-read sequencing technologies to assemble the whole genome and identify the resistance gene repertoire. We previously narrowed the location of the gene(s) to a region on chromosome four. The region in Carolina Gold Select is larger than in the rice reference genome (Nipponbare) and contains twice as many resistance genes. One shares unusual features with a known bacterial disease resistance gene, suggesting that it confers the resistance. Across diverse varieties and related species, we identified two widely-distributed groups of such genes. The results are an important step toward mechanistic characterization and deployment of the bacterial disease resistance. The genome assembly also identified a resistance gene for a fungal disease and predicted a marker phenotype used in breeding for resistance. Thus, the Carolina Gold Select genome assembly can be expected to aid in the identification and deployment of other valuable traits.

Regulation of Cell Wall Thickening by a Medley of Mechanisms

To provide physical support for developing structures and to withstand the pressures associated with water and nutrient transport, some cells deposit a secondary cell wall, a rigid matrix of polysaccharide and phenolic biopolymers. The biosynthesis and deposition of these materials and the patterning of secondary wall-forming cells is controlled by a network of transcription factors. However, recent work suggests that this network forms the core of a more complex, multilevel regulatory system. This expanded system includes epigenetic, post-transcriptional, and post-translational regulation, and is coordinated with other pathways controlling primary growth and responses to environmental stimuli. New findings expand the set of transcription factors identified as secondary cell wall regulators and reveal novel regulatory processes that further govern secondary wall biogenesis.

The unique epigenetic features of Pack-MULEs and their impact on chromosomal base composition and expression spectrum

Acquisition and rearrangement of host genes by transposable elements (TEs) is an important mechanism to increase gene diversity as exemplified by the ∼3000 Pack-Mutator-like TEs in the rice genome which have acquired gene sequences (Pack-MULEs), yet remain enigmatic. To identify signatures of functioning Pack-MULEs and Pack-MULE evolution, we generated transcriptome, translatome, and epigenome datasets and compared Pack-MULEs to genes and other TE families. Approximately 40% of Pack-MULEs were transcribed with 9% having translation evidence, clearly distinguishing them from other TEs. Pack-MULEs exhibited a unique expression profile associated with specificity in reproductive tissues that may be associated with seed traits. Expressed Pack-MULEs resemble regular protein-coding genes as exhibited by a low level of DNA methylation, association with active histone marks and DNase I hypersensitive sites, and an absence of repressive histone marks, suggesting that a substantial fraction of Pack-MULEs are potentially functional in vivo. Interestingly, the expression capacity of Pack-MULEs is independent of the local genomic environment, and the insertion and expression of Pack-MULEs may have altered the local chromosomal expression pattern as well as counteracted the impact of recombination on chromosomal base composition, which has profound consequences on the evolution of chromosome structure.

Diversification of Transposable Elements in Arthropods and Its Impact on Genome Evolution

Transposable elements (TEs) are ubiquitous in arthropods. However, analyses of large-scale and long-term coevolution between TEs and host genomes remain scarce in arthropods. Here, we choose 14 representative Arthropoda species from eight orders spanning more than 500 million years of evolution. By developing an unbiased TE annotation pipeline, we obtained 87 to 2266 TE reference sequences in a species, which is a considerable improvement compared to the reference TEs previously annotated in Repbase. We find that TE loads are diversified among species and were previously underestimated. The highly species- and time-specific expansions and contractions, and intraspecific sequence diversification are the leading driver of long terminal repeat (LTR) dynamics in Lepidoptera. Terminal inverted repeats (TIRs) proliferated substantially in five species with large genomes. A phylogenetic comparison reveals that the loads of multiple TE subfamilies are positively correlated with genome sizes. We also identified a few horizontally transferred TE candidates across nine species. In addition, we set up the Arthropod Transposable Elements database (ArTEdb) to provide TE references and annotations. Collectively, our results provide high-quality TE references and uncover that TE loads and expansion histories vary greatly among arthropods, which implies that TEs are an important driving force shaping the evolution of genomes through gain and loss.

A novel variant of Gh_D02G0276 is required for root-knot nematode resistance on chromosome 14 (D02) in Upland cotton

MAGIC population sequencing and virus-induced gene silencing identify Gh_D02G0276 as a novel root-knot nematode resistance gene on chromosome 14 in Upland cotton. The southern root-knot nematode [RKN; Meloidogyne incognita (Kofoid & White)] remains the primary yield-limiting biotic stress to Upland cotton (Gossypium hirsutum L.) throughout the southeastern USA. While useful genetic markers have been developed for two major RKN resistance loci on chromosomes 11 (A11) and 14 (D02), these markers are not completely effective because the causative genes have not been identified. Here, we sequenced 550 recombinant inbred lines (RILs) from a multi-parent advanced generation intercross (MAGIC) population to identify five RILs that had informative recombinations near the D02-RKN resistance locus. The RKN resistance phenotypes of these five RILs narrowed the D02-RKN locus to a 30-kb region with four candidate genes. We conducted virus-induced gene silencing (VIGS) on each of these genes and found that Gh_D02G0276 was required for suppression of RKN egg production conferred by the Chr. D02 resistance gene. The resistant lines all possessed an allele of Gh_D02G0276 that showed non-synonymous mutations and was prematurely truncated. Furthermore, a Gh_D02G0276-specific marker for the resistance allele variant was able to identify RKN-resistant germplasm from a collection of 367 cotton accessions. The Gh_D02G0276 peptide shares similarity with domesticated hAT-like transposases with additional novel N- and C-terminal domains that resemble the target of known RKN effector molecules and a prokaryotic motif, respectively. The truncation in the resistance allele results in a loss of a plant nuclear gene-specific C-terminal motif, potentially rendering this domain antigenic due to its high homology with bacterial proteins. The conclusive identification of this RKN resistance gene opens new avenues for understanding plant resistance mechanisms to RKN as well as opportunities to develop more efficient marker-assisted selection in cotton breeding programs.

BED-domain-containing immune receptors confer diverse resistance spectra to yellow rust

Crop diseases reduce wheat yields by ~25% globally and thus pose a major threat to global food security¹. Genetic resistance can reduce crop losses in the field and can be selected through the use of molecular markers. However, genetic resistance often breaks down following changes in pathogen virulence, as experienced with the wheat yellow (stripe) rust fungus Puccinia striiformis f. sp. tritici (Pst)². This highlights the need to (1) identify genes that, alone or in combination, provide broad-spectrum resistance, and (2) increase our understanding of the underlying molecular modes of action. Here we report the isolation and characterization of three major yellow rust resistance genes (Yr7, Yr5 and YrSP) from hexaploid wheat (Triticum aestivum), each having a distinct recognition specificity. We show that Yr5, which remains effective to a broad range of Pst isolates worldwide, is closely related yet distinct from Yr7, whereas YrSP is a truncated version of Yr5 with 99.8% sequence identity. All three Yr genes belong to a complex resistance gene cluster on chromosome 2B encoding nucleotide-binding and leucine-rich repeat proteins (NLRs) with a non-canonical N-terminal zinc-finger BED domain³ that is distinct from those found in non-NLR wheat proteins. We developed diagnostic markers to accelerate haplotype analysis and for marker-assisted selection to expedite the stacking of the non-allelic Yr genes. Our results provide evidence that the BED-NLR gene architecture can provide effective field-based resistance to important fungal diseases such as wheat yellow rust.

Abiotic Stress Phenotypes Are Associated with Conserved Genes Derived from Transposable Elements

Plant phenomics offers unique opportunities to accelerate our understanding of gene function and plant response to different environments, and may be particularly useful for studying previously uncharacterized genes. One important type of poorly characterized genes is those derived from transposable elements (TEs), which have departed from a mobility-driven lifestyle to attain new adaptive roles for the host (exapted TEs). We used phenomics approaches, coupled with reverse genetics, to analyze T-DNA insertion mutants of both previously reported and novel protein-coding exapted TEs in the model plant Arabidopsis thaliana. We show that mutations in most of these exapted TEs result in phenotypes, particularly when challenged by abiotic stress. We built statistical multi-dimensional phenotypic profiles and compared them to wild-type and known stress responsive mutant lines for each particular stress condition. We found that these exapted TEs may play roles in responses to phosphate limitation, tolerance to high salt concentration, freezing temperatures, and arsenic toxicity. These results not only experimentally validate a large set of putative functional exapted TEs recently discovered through computational analysis, but also uncover additional novel phenotypes for previously well-characterized exapted TEs in A. thaliana.

Regulation of rice root development by a retrotransposon acting as a microRNA sponge

It is well documented that transposable elements (TEs) can regulate the expression of neighbouring genes. However, their ability to act in trans and influence ectopic loci has been reported rarely. We searched in rice transcriptomes for tissue-specific expression of TEs and found them to be regulated developmentally. They often shared sequence homology with co-expressed genes and contained potential microRNA-binding sites, which suggested possible contributions to gene regulation. In fact, we have identified a retrotransposon that is highly transcribed in roots and whose spliced transcript constitutes a target mimic for miR171. miR171 destabilizes mRNAs encoding the root-specific family of SCARECROW-Like transcription factors. We demonstrate that retrotransposon-derived transcripts act as decoys for miR171, triggering its degradation and thus results in the root-specific accumulation of SCARECROW-Like mRNAs. Such transposon-mediated post-transcriptional control of miR171 levels is conserved in diverse rice species.

Regulation of rice root development by a retrotransposon acting as a microRNA sponge

It is well documented that transposable elements (TEs) can regulate the expression of neighbouring genes. However, their ability to act in trans and influence ectopic loci has been reported rarely. We searched in rice transcriptomes for tissue-specific expression of TEs and found them to be regulated developmentally. They often shared sequence homology with co-expressed genes and contained potential microRNA-binding sites, which suggested possible contributions to gene regulation. In fact, we have identified a retrotransposon that is highly transcribed in roots and whose spliced transcript constitutes a target mimic for miR171. miR171 destabilizes mRNAs encoding the root-specific family of SCARECROW-Like transcription factors. We demonstrate that retrotransposon-derived transcripts act as decoys for miR171, triggering its degradation and thus results in the root-specific accumulation of SCARECROW-Like mRNAs. Such transposon-mediated post-transcriptional control of miR171 levels is conserved in diverse rice species.

LTR-retrotransposons in plants: Engines of evolution

LTR retrotransposons are the most abundant group of transposable elements (TEs) in plants. These elements can fall inside or close to genes, and therefore influence their expression and evolution. This review aims to examine how LTR retrotransposons, especially Ty1-copia elements, mediate gene regulation and evolution. Various stimuli, including polyploidization and biotic and abiotic elicitors, result in the transcription and movement of these retrotransposons, and can facilitate adaptation. The presence of cis-regulatory motifs in the LTRs are central to their stress-mediated responses and are shared with host stress-responsive genes, showing a complex evolutionary history in which TEs provide new regulatory units to genes. The presence of retrotransposon remnants in genes that are necessary for normal gene function, demonstrates the importance of exaptation and co-option, and is also a consequence of the abundance of these elements in plant genomes. Furthermore, insertions of LTR retrotransposons in and around genes provide potential for alternative splicing, epigenetic control, transduction, duplication and recombination. These characteristics can become an active part of the evolution of gene families as in the case of resistance genes (R-genes). The character of TEs as exclusively selfish is now being re-evaluated. Since genome-wide reprogramming via TEs is a long evolutionary process, the changes we can examine are case-specific and their fitness advantage may not be evident until TE-derived motifs and domains have been completely co-opted and fixed. Nevertheless, the presence of LTR retrotransposons inside genes and as part of gene promoter regions is consistent with their roles as engines of plant genome evolution.