A Nest of LTR Retrotransposons Adjacent the Disease Resistance-Priming Gene NPR1 in Beta vulgaris L. U.S. Hybrid H20

Comparative transcriptome profiling of near isogenic lines PBW343 and FLW29 to unravel defense related genes and pathways contributing to stripe rust resistance in wheat

Stripe rust (Sr), caused by Puccinia striiformis f. sp. tritici (Pst), is the most devastating disease that poses serious threat to the wheat-growing nations across the globe. Developing resistant cultivars is the most challenging aspect in wheat breeding. The function of resistance genes (R genes) and the mechanisms by which they influence plant-host interactions are poorly understood. In the present investigation, comparative transcriptome analysis was carried out by involving two near-isogenic lines (NILs) PBW343 and FLW29. The seedlings of both the genotypes were inoculated with Pst pathotype 46S119. In total, 1106 differentially expressed genes (DEGs) were identified at early stage of infection (12 hpi), whereas expressions of 877 and 1737 DEGs were observed at later stages (48 and 72 hpi) in FLW29. The identified DEGs were comprised of defense-related genes including putative R genes, 7 WRKY transcriptional factors, calcium, and hormonal signaling associated genes. Moreover, pathways involved in signaling of receptor kinases, G protein, and light showed higher expression in resistant cultivar and were common across different time points. Quantitative real-time PCR was used to further confirm the transcriptional expression of eight critical genes involved in plant defense mechanism against stripe rust. The information about genes are likely to improve our knowledge of the genetic mechanism that controls the stripe rust resistance in wheat, and data on resistance response-linked genes and pathways will be a significant resource for future research.

New reference genome sequences of hot pepper reveal the massive evolution of plant disease-resistance genes by retroduplication

Background

Transposable elements are major evolutionary forces which can cause new genome structure and species diversification. The role of transposable elements in the expansion of nucleotide-binding and leucine-rich-repeat proteins (NLRs), the major disease-resistance gene families, has been unexplored in plants.

Results

We report two high-quality de novo genomes (Capsicum baccatum and C. chinense) and an improved reference genome (C. annuum) for peppers. Dynamic genome rearrangements involving translocations among chromosomes 3, 5, and 9 were detected in comparison between C. baccatum and the two other peppers. The amplification of athila LTR-retrotransposons, members of the gypsy superfamily, led to genome expansion in C. baccatum. In-depth genome-wide comparison of genes and repeats unveiled that the copy numbers of NLRs were greatly increased by LTR-retrotransposon-mediated retroduplication. Moreover, retroduplicated NLRs are abundant across the angiosperms and, in most cases, are lineage-specific.

Conclusions

Our study reveals that retroduplication has played key roles for the massive emergence of NLR genes including functional disease-resistance genes in pepper plants.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1341-9) contains supplementary material, which is available to authorized users.

New reference genome sequences of hot pepper reveal the massive evolution of plant disease-resistance genes by retroduplication

BACKGROUND

Transposable elements are major evolutionary forces which can cause new genome structure and species diversification. The role of transposable elements in the expansion of nucleotide-binding and leucine-rich-repeat proteins (NLRs), the major disease-resistance gene families, has been unexplored in plants.

RESULTS

We report two high-quality de novo genomes (Capsicum baccatum and C. chinense) and an improved reference genome (C. annuum) for peppers. Dynamic genome rearrangements involving translocations among chromosomes 3, 5, and 9 were detected in comparison between C. baccatum and the two other peppers. The amplification of athila LTR-retrotransposons, members of the gypsy superfamily, led to genome expansion in C. baccatum. In-depth genome-wide comparison of genes and repeats unveiled that the copy numbers of NLRs were greatly increased by LTR-retrotransposon-mediated retroduplication. Moreover, retroduplicated NLRs are abundant across the angiosperms and, in most cases, are lineage-specific.

CONCLUSIONS

Our study reveals that retroduplication has played key roles for the massive emergence of NLR genes including functional disease-resistance genes in pepper plants.

Genetic diversity of the conserved motifs of six bacterial leaf blight resistance genes in a set of rice landraces

BACKGROUND

Bacterial leaf blight (BLB) caused by the vascular pathogen Xanthomonas oryzae pv. oryzae (Xoo) is one of the most serious diseases leading to crop failure in rice growing countries. A total of 37 resistance genes against Xoo has been identified in rice. Of these, ten BLB resistance genes have been mapped on rice chromosomes, while 6 have been cloned, sequenced and characterized. Diversity analysis at the resistance gene level of this disease is scanty, and the landraces from West Bengal and North Eastern states of India have received little attention so far. The objective of this study was to assess the genetic diversity at conserved domains of 6 BLB resistance genes in a set of 22 rice accessions including landraces and check genotypes collected from the states of Assam, Nagaland, Mizoram and West Bengal.

RESULTS

In this study 34 pairs of primers were designed from conserved domains of 6 BLB resistance genes; Xa1, xa5, Xa21, Xa21(A1), Xa26 and Xa27. The designed primer pairs were used to generate PCR based polymorphic DNA profiles to detect and elucidate the genetic diversity of the six genes in the 22 diverse rice accessions of known disease phenotype. A total of 140 alleles were identified including 41 rare and 26 null alleles. The average polymorphism information content (PIC) value was 0.56/primer pair. The DNA profiles identified each of the rice landraces unequivocally. The amplified polymorphic DNA bands were used to calculate genetic similarity of the rice landraces in all possible pair combinations. The similarity among the rice accessions ranged from 18% to 89% and the dendrogram produced from the similarity values was divided into 2 major clusters. The conserved domains identified within the sequenced rare alleles include Leucine-Rich Repeat, BED-type zinc finger domain, sugar transferase domain and the domain of the carbohydrate esterase 4 superfamily.

CONCLUSIONS

This study revealed high genetic diversity at conserved domains of six BLB resistance genes in a set of 22 rice accessions. The inclusion of more genotypes from remote ecological niches and hotspots holds promise for identification of further genetic diversity at the BLB resistance genes.

Transposon variation by order during allopolyploidisation between Brassica oleracea and Brassica rapa

Although many studies have shown that transposable element (TE) activation is induced by hybridisation and polyploidisation in plants, much less is known on how different types of TE respond to hybridisation, and the impact of TE-associated sequences on gene function. We investigated the frequency and regularity of putative transposon activation for different types of TE, and determined the impact of TE-associated sequence variation on the genome during allopolyploidisation. We designed different types of TE primers and adopted the Inter-Retrotransposon Amplified Polymorphism (IRAP) method to detect variation in TE-associated sequences during the process of allopolyploidisation between Brassica rapa (AA) and Brassica oleracea (CC), and in successive generations of self-pollinated progeny. In addition, fragments with TE insertions were used to perform Blast2GO analysis to characterise the putative functions of the fragments with TE insertions. Ninety-two primers amplifying 548 loci were used to detect variation in sequences associated with four different orders of TE sequences. TEs could be classed in ascending frequency into LTR-REs, TIRs, LINEs, SINEs and unknown TEs. The frequency of novel variation (putative activation) detected for the four orders of TEs was highest from the F1 to F2 generations, and lowest from the F2 to F3 generations. Functional annotation of sequences with TE insertions showed that genes with TE insertions were mainly involved in metabolic processes and binding, and preferentially functioned in organelles. TE variation in our study severely disturbed the genetic compositions of the different generations, resulting in inconsistencies in genetic clustering. Different types of TE showed different patterns of variation during the process of allopolyploidisation.

Insights into the loblolly pine genome: characterization of BAC and fosmid sequences

Despite their prevalence and importance, the genome sequences of loblolly pine, Norway spruce, and white spruce, three ecologically and economically important conifer species, are just becoming available to the research community. Following the completion of these large assemblies, annotation efforts will be undertaken to characterize the reference sequences. Accurate annotation of these ancient genomes would be aided by a comprehensive repeat library; however, few studies have generated enough sequence to fully evaluate and catalog their non-genic content. In this paper, two sets of loblolly pine genomic sequence, 103 previously assembled BACs and 90,954 newly sequenced and assembled fosmid scaffolds, were analyzed. Together, this sequence represents 280 Mbp (roughly 1% of the loblolly pine genome) and one of the most comprehensive studies of repetitive elements and genes in a gymnosperm species. A combination of homology and de novo methodologies were applied to identify both conserved and novel repeats. Similarity analysis estimated a repetitive content of 27% that included both full and partial elements. When combined with the de novo investigation, the estimate increased to almost 86%. Over 60% of the repetitive sequence consists of full or partial LTR (long terminal repeat) retrotransposons. Through de novo approaches, 6,270 novel, full-length transposable element families and 9,415 sub-families were identified. Among those 6,270 families, 82% were annotated as single-copy. Several of the novel, high-copy families are described here, with the largest, PtPiedmont, comprising 133 full-length copies. In addition to repeats, analysis of the coding region reported 23 full-length eukaryotic orthologous proteins (KOGS) and another 29 novel or orthologous genes. These discoveries, along with other genomic resources, will be used to annotate conifer genomes and address long-standing questions about gymnosperm evolution.

New insights into nested long terminal repeat retrotransposons in Brassica species

Long terminal repeat (LTR) retrotransposons, one of the foremost types of transposons, continually change or modify gene function and reorganize the genome through bursts of dramatic proliferation. Many LTR-TEs preferentially insert within other LTR-TEs, but the cause and evolutionary significance of these nested LTR-TEs are not well understood. In this study, a total of 1.52Gb of Brassica sequence containing 2020 bacterial artificial chromosomes (BACs) was scanned, and six bacterial artificial chromosome (BAC) clones with extremely nested LTR-TEs (LTR-TEs density: 7.24/kb) were selected for further analysis. The majority of the LTR-TEs in four of the six BACs were found to be derived from the rapid proliferation of retrotransposons originating within the BAC regions, with only a few LTR-TEs originating from the proliferation and insertion of retrotransposons from outside the BAC regions approximately 5-23Mya. LTR-TEs also preferably inserted into TA-rich repeat regions. Gene prediction by Genescan identified 207 genes in the 0.84Mb of total BAC sequences. Only a few genes (3/207) could be matched to the Brassica expressed sequence tag (EST) database, indicating that most genes were inactive after retrotransposon insertion. Five of the six BACs were putatively centromeric. Hence, nested LTR-TEs in centromere regions are rapidly duplicated, repeatedly inserted, and act to suppress activity of genes and to reshuffle the structure of the centromeric sequences. Our results suggest that LTR-TEs burst and proliferate on a local scale to create nested LTR-TE regions, and that these nested LTR-TEs play a role in the formation of centromeres.

Highly diverse chromoviruses of Beta vulgaris are classified by chromodomains and chromosomal integration

Background

Chromoviruses are one of the three genera of Ty3-gypsy long terminal repeat (LTR) retrotransposons, and are present in high copy numbers in plant genomes. They are widely distributed within the plant kingdom, with representatives even in lower plants such as green and red algae. Their hallmark is the presence of a chromodomain at the C-terminus of the integrase. The chromodomain exhibits structural characteristics similar to proteins of the heterochromatin protein 1 (HP1) family, which mediate the binding of each chromovirus type to specific histone variants. A specific integration via the chromodomain has been shown for only a few chromoviruses. However, a detailed study of different chromoviral clades populating a single plant genome has not yet been carried out.

Results

We conducted a comprehensive survey of chromoviruses within the Beta vulgaris (sugar beet) genome, and found a highly diverse chromovirus population, with significant differences in element size, primarily caused by their flanking LTRs. In total, we identified and annotated full-length members of 16 families belonging to the four plant chromoviral clades: CRM, Tekay, Reina, and Galadriel. The families within each clade are structurally highly conserved; in particular, the position of the chromodomain coding region relative to the polypurine tract is clade-specific. Two distinct groups of chromodomains were identified. The group II chromodomain was present in three chromoviral clades, whereas families of the CRM clade contained a more divergent motif. Physical mapping using representatives of all four clades identified a clade-specific integration pattern. For some chromoviral families, we detected the presence of expressed sequence tags, indicating transcriptional activity.

Conclusions

We present a detailed study of chromoviruses, belonging to the four major clades, which populate a single plant genome. Our results illustrate the diversity and family structure of B. vulgaris chromoviruses, and emphasize the role of chromodomains in the targeted integration of these viruses. We suggest that the diverse sets of plant chromoviruses with their different localization patterns might help to facilitate plant-genome organization in a structural and functional manner.

Evolutionary history of LTR retrotransposon chromodomains in plants

Chromodomain-containing LTR retrotransposons are one of the most successful groups of mobile elements in plant genomes. Previously, we demonstrated that two types of chromodomains (CHDs) are carried by plant LTR retrotransposons. Chromodomains from group I (CHD_I) were detected only in Tcn1-like LTR retrotransposons from nonseed plants such as mosses (including the model moss species Physcomitrella) and lycophytes (the Selaginella species). LTR retrotransposon chromodomains from group II (CHD_II) have been described from a wide range of higher plants. In the present study, we performed computer-based mining of plant LTR retrotransposon CHDs from diverse plants with an emphasis on spike-moss Selaginella. Our extended comparative and phylogenetic analysis demonstrated that two types of CHDs are present only in the Selaginella genome, which puts this species in a unique position among plants. It appears that a transition from CHD_I to CHD_II and further diversification occurred in the evolutionary history of plant LTR retrotransposons at approximately 400 MYA and most probably was associated with the evolution of chromatin organization.

Genes Encoding Callose Synthase and Phytochrome A Are Adjacent to a MAP3Kα-Like Gene in Beta vulgaris US H20

MAP3Kα, a gene that encodes a key conserved protein kinase, is responsible for initiating a rapid cascade of cellular events leading to localized cell death. Hypersensitive response, as it is termed, enables genetically resistant plants to limit microbial invasion under the right environmental conditions. Since knowledge of close physically linked genes is important for genome analysis and possibly for improving disease resistance, systematic DNA sequence analysis, gene annotation, and protein BLASTs were performed to identify and characterize genes in close physical proximity to a MAP3Kα-like gene in Beta vulgaris L. US H20. On the same 125 Kb BAC, callose synthase (BvCS) and phytochrome A (PhyA) genes were within 50 Kb of MAP3Kα. The close physical linkage of these genes may result from selection for coordinated responses to disease pressure. Bert, a new chromodomain-carrying gypsy-like LTR retrotransposon, resides within an intron of the BvCS gene, where it is transcribed from the opposing strand.

A specific insertion of a solo-LTR characterizes the Y-chromosome of Bryonia dioica (Cucurbitaceae)

Background

Relatively few species of flowering plants are dioecious and even fewer are known to have sex chromosomes. Current theory posits that homomorphic sex chromosomes, such as found in Bryonia dioica (Cucurbitaceae), offer insight into the early stages in the evolution of sex chromosomes from autosomes. Little is known about these early steps, but an accumulation of transposable element sequences has been observed on the Y-chromosomes of some species with heteromorphic sex chromosomes. Recombination, by which transposable elements are removed, is suppressed on at least part of the emerging Y-chromosome, and this may explain the correlation between the emergence of sex chromosomes and transposable element enrichment.

Findings

We sequenced 2321 bp of the Y-chromosome in Bryonia dioica that flank a male-linked marker, BdY1, reported previously. Within this region, which should be suppressed for recombination, we observed a solo-LTR nested in a Copia-like transposable element. We also found other, presumably paralogous, solo-LTRs in a consensus sequence of the underlying Copia-like transposable element.

Conclusions

Given that solo-LTRs arise via recombination events, it is noteworthy that we find one in a genomic region where recombination should be suppressed. Although the solo-LTR could have arisen before recombination was suppressed, creating the male-linked marker BdY1, our previous study on B. dioica suggested that BdY1 may not lie in the recombination-suppressed region of the Y-chromosome in all populations. Presence of a solo-LTR near BdY1 therefore fits with the observed correlation between retrotransposon accumulation and the suppression of recombination early in the evolution of sex chromosomes. These findings further suggest that the homomorphic sex chromosomes of B. dioica, the first organism for which genetic XY sex-determination was inferred, are evolutionarily young and offer reference information for comparative studies of other plant sex chromosomes.