Large-scale insertional mutagenesis using the Tnt1 retrotransposon in the model legume Medicago truncatula

Medicago truncatula is a fast-emerging model for the study of legume functional biology. We used the tobacco retrotransposon Tnt1 to tag the Medicago genome and generated over 7600 independent lines representing an estimated 190 000 insertion events. Tnt1 inserted on average at 25 different locations per genome during tissue culture, and insertions were stable during subsequent generations in soil. Analysis of 2461 Tnt1 flanking sequence tags (FSTs) revealed that Tnt1 appears to prefer gene-rich regions. The proportion of Tnt1 insertion in coding sequences was 34.1%, compared to the expected 15.9% if random insertions were to occur. However, Tnt1 showed neither unique target site specificity nor strong insertion hot spots, although some genes were more frequently tagged than others. Forward-genetic screening of 3237 R₁ lines resulted in identification of visible mutant phenotypes in approximately 30% of the regenerated lines. Tagging efficiency appears to be high, as all of the 20 mutants examined so far were found to be tagged. Taking the properties of Tnt1 into account and assuming 1.7 kb for the average M. truncatula gene size, we estimate that approximately 14 000–16 000 lines would be sufficient for 90% gene tagging coverage in M. truncatula . This is in contrast to more than 500 000 lines required to achieve the same saturation level using T-DNA tagging. Our data demonstrate that Tnt1 is an efficient insertional mutagen in M. truncatula , and could be a primary choice for other plant species with large genomes.     

Osmotic shock improves <Emphasis Type="Italic">Tnt1</Emphasis> transposition frequency in <Emphasis Type="Italic">Medicago truncatula</Emphasis> cv Jemalong during in vitro regeneration

Insertion mutant collections are powerful tools for genetic studies in plants. Although large-scale insertional mutagenesis using T-DNA is not feasible in legumes, the Tnt1 tobacco retrotransposon can be used as a very efficient mutagen in the Medicago truncatula R108 genotype. In this article, we show that Tnt1 can also be exploited to create insertional mutants via transformation and/or regeneration in the reference cultivar Jemalong. Tnt1 insertional mutagenesis in Jemalong following Agrobacterium tumefaciens-mediated transformation was found to be very efficient, with an average of greater than 15 insertions/line. In contrast, regeneration using low-copy transgenic starter lines resulted in a highly variable rate of new Tnt1 insertions. With the goal of increasing the number of additional Tnt1 insertions during regeneration of starter lines, we have compared the insertion frequencies for a number of different regeneration protocols. In addition, we have been able to show that sucrose-mediated osmotic shock preceding regeneration significantly increases the transposition frequency. Under optimal conditions, 95% of the regenerated Jemalong plants possess new insertions.     

FLAGdb/FST: a database of mapped flanking insertion sites (FSTs) of Arabidopsis thaliana T-DNA transformants

A large collection of T-DNA insertion transformants of Arabidopsis thaliana has been generated at the Institute of Agronomic Research, Versailles, France. The molecular characterisation of the insertion sites is currently performed by sequencing genomic regions flanking the inserted T-DNA (FST). The almost complete sequence of the nuclear genome of A.thaliana provides the framework for organising FSTs in a genome oriented database, FLAGdb/FST (http://genoplante-info.infobiogen.fr). The main scope of FLAGdb/FST is to help biologists to find the FSTs that interrupt the genes in which they are interested. FSTs are anchored to the genome sequences of A.thaliana and positions of both predicted genes and FSTs are shown graphically on sequences. Requests to locate the genomic position of a query sequence are made using BLAST programs. The response delivered by FLAGdb/FST is a graphical representation of the putative FSTs and of predicted genes in a 20 kb region.     

Molecular mechanism underlying pseudopeloria in Habenaria radiata (Orchidaceae)

Habenaria radiata (Orchidaceae) has two whorls of perianth, comprising three greenish sepals, two white petals and one lip (labellum). By contrast, the pseudopeloric (with a decreased degree of zygomorphy) mutant cultivar of H. radiata , ‘Hishou’, has changes in the identities of the dorsal sepal to a petaloid organ and the two ventral sepals to lip‐like organs. Here, we isolated four DEFICIENS ‐ like and two AGL 6 ‐like genes from H. radiata , and characterized their expression. Most of these genes revealed similar expression patterns in the wild type and in the ‘Hishou’ cultivar, except Hr DEF ‐C3. The Hr DEF ‐C3 gene was expressed in petals and lip in the wild type but was ectopically expressed in sepal, petals, lip, leaf, root and bulb in ‘Hishou’. Sequence analysis of the Hr DEF ‐C3 loci revealed that the ‘Hishou’ genome harbored two types of Hr DEF ‐C3 genes: one identical to wild‐type Hr DEF ‐C3 and the other carrying a retrotransposon insertion in its promoter. Genetic linkage analysis of the progeny derived from an intraspecific cross between ‘Hishou’ and the wild type demonstrated that the mutant pseudopeloric trait was dominantly inherited and was linked to the Hr DEF ‐C3 gene carrying the retrotransposon. These results indicate that the pseudopeloric phenotype is caused by retrotransposon insertion in the Hr DEF ‐C3 promoter, resulting in the ectopic expression of Hr DEF ‐C3 . As the expression of Hr AGL 6‐C2 was limited to lateral sepals and lip, the overlapping expression of Hr DEF ‐C3 and Hr AGL 6‐C2 is likely to be responsible for the sepal to lip‐like identity in the lateral sepals of the ‘Hishou’ cultivar.     

The high‐quality genome of Brassica napus cultivar ‘ZS11’ reveals the introgression history in semi‐winter morphotype

Allotetraploid oilseed rape (Brassica napus L.) is an agriculturally important crop. Cultivation and breeding of B. napus by humans has resulted in numerous genetically diverse morphotypes with optimized agronomic traits and ecophysiological adaptation. To further understand the genetic basis of diversification and adaptation, we report a draft genome of an Asian semi‐winter oilseed rape cultivar ‘ZS11’ and its comprehensive genomic comparison with the genomes of the winter‐type cultivar ‘Darmor‐bzh’ as well as two progenitors. The integrated BAC‐to‐BAC and whole‐genome shotgun sequencing strategies were effective in the assembly of repetitive regions (especially young long terminal repeats) and resulted in a high‐quality genome assembly of B. napus ‘ZS11’. Within a short evolutionary period (~6700 years ago), semi‐winter‐type ‘ZS11’ and the winter‐type ‘Darmor‐bzh’ maintained highly genomic collinearity. Even so, certain genetic differences were also detected in two morphotypes. Relative to ‘Darmor‐bzh’, both two subgenomes of ‘ZS11’ are closely related to its progenitors, and the ‘ZS11’ genome harbored several specific segmental homoeologous exchanges (HEs). Furthermore, the semi‐winter‐type ‘ZS11’ underwent potential genomic introgressions with B. rapa (A_r). Some of these genetic differences were associated with key agronomic traits. A key gene of A03.FLC3 regulating vernalization‐responsive flowering time in ‘ZS11’ was first experienced HE, and then underwent genomic introgression event with A_r, which potentially has led to genetic differences in controlling vernalization in the semi‐winter types. Our observations improved our understanding of the genetic diversity of different B. napus morphotypes and the cultivation history of semi‐winter oilseed rape in Asia.     

Rapid evolutionary dynamics in a 2.8‐Mb chromosomal region containing multiple prolamin and resistance gene families in Aegilops tauschii

Prolamin and resistance gene families are important in wheat food use and in defense against pathogen attacks, respectively. To better understand the evolution of these multi‐gene families, the DNA sequence of a 2.8‐Mb genomic region, representing an 8.8 cM genetic interval and harboring multiple prolamin and resistance‐like gene families, was analyzed in the diploid grass Aegilops tauschii, the D‐genome donor of bread wheat. Comparison with orthologous regions from rice, Brachypodium, and sorghum showed that the Ae. tauschii region has undergone dramatic changes; it has acquired more than 80 non‐syntenic genes and only 13 ancestral genes are shared among these grass species. These non‐syntenic genes, including prolamin and resistance‐like genes, originated from various genomic regions and likely moved to their present locations via sequence evolution processes involving gene duplication and translocation. Local duplication of non‐syntenic genes contributed significantly to the expansion of gene families. Our analysis indicates that the insertion of prolamin‐related genes occurred prior to the separation of the Brachypodieae and Triticeae lineages. Unlike in Brachypodium, inserted prolamin genes have rapidly evolved and expanded to encode different classes of major seed storage proteins in Triticeae species. Phylogenetic analyses also showed that the multiple insertions of resistance‐like genes and subsequent differential expansion of each R gene family. The high frequency of non‐syntenic genes and rapid local gene evolution correlate with the high recombination rate in the 2.8‐Mb region with nine‐fold higher than the genome‐wide average. Our results demonstrate complex evolutionary dynamics in this agronomically important region of Triticeae species.     

Improving Nelumbo nucifera genome assemblies using high‐resolution genetic maps and BioNano genome mapping reveals ancient chromosome rearrangements

Genetic and physical maps are powerful tools to anchor fragmented draft genome assemblies generated from next‐generation sequencing. Currently, two draft assemblies of Nelumbo nucifera, the genomes of ‘China Antique’ and ‘Chinese Tai‐zi’, have been released. However, there is presently no information on how the sequences are assembled into chromosomes in N. nucifera. The lack of physical maps and inadequate resolution of available genetic maps hindered the assembly of N. nucifera chromosomes. Here, a linkage map of N. nucifera containing 2371 bin markers [217 577 single nucleotide polymorphisms (SNPs)] was constructed using restriction‐site associated DNA sequencing data of 181 F₂ individuals and validated by adding 197 simple sequence repeat (SSR) markers. Additionally, a BioNano optical map covering 86.20% of the ‘Chinese Tai‐zi’ genome was constructed. The draft assembly of ‘Chinese Tai‐zi’ was improved based on the BioNano optical map, showing an increase of the scaffold N50 from 0.989 to 1.48 Mb. Using a combination of multiple maps, 97.9% of the scaffolds in the ‘Chinese Tai‐zi’ draft assembly and 97.6% of the scaffolds in the ‘China Antique’ draft assembly were anchored into pseudo‐chromosomes, and the centromere regions along the pseudo‐chromosomes were identified. An evolutionary scenario was proposed to reach the modern N. nucifera karyotype from the seven ancestral eudicot chromosomes. The present study provides the highest‐resolution linkage map, the optical map and chromosome level genome assemblies for N. nucifera, which are valuable for the breeding and cultivation of N. nucifera and future studies of comparative and evolutionary genomics in angiosperms.     

Successful gene tagging in lettuce using the Tnt1 retrotransposon from tobacco

The tobacco (Nicotiana tabacum) element Tnt1 is one of the few identified active retrotransposons in plants. These elements possess unique properties that make them ideal genetic tools for gene tagging. Here, we demonstrate the feasibility of gene tagging using the retrotransposon Tnt1 in lettuce (Lactuca sativa), which is the largest genome tested for retrotransposon mutagenesis so far. Of 10 different transgenic bushes carrying a complete Tnt1 containing T-DNA, eight contained multiple transposed copies of Tnt1. The number of transposed copies of the element per plant was particularly high, the smallest number being 28. Tnt1 transposition in lettuce can be induced by a very simple in vitro culture protocol. Tnt1 insertions were stable in the progeny of the primary transformants and could be segregated genetically. Characterization of the sequences flanking some insertion sites revealed that Tnt1 often inserted into genes. The progeny of some primary transformants showed phenotypic alterations due to recessive mutations. One of these mutations was due to Tnt1 insertion in the gibberellin 3beta-hydroxylase gene. Taken together, these results indicate that Tnt1 is a powerful tool for insertion mutagenesis especially in plants with a large genome.     

Development of iFOX‐hunting as a functional genomic tool and demonstration of its use to identify early senescence‐related genes in the polyploid Brassica napus

Functional genomic studies of many polyploid crops, including rapeseed (Brassica napus), are constrained by limited tool sets. Here we report development of a gain‐of‐function platform, termed ‘iFOX (inducible Full‐length cDNA OvereXpressor gene)‐Hunting’, for inducible expression of B. napus seed cDNAs in Arabidopsis. A Gateway‐compatible plant gene expression vector containing a methoxyfenozide‐inducible constitutive promoter for transgene expression was developed. This vector was used for cloning of random cDNAs from developing B. napus seeds and subsequent Agrobacterium‐mediated transformation of Arabidopsis. The inducible promoter of this vector enabled identification of genes upon induction that are otherwise lethal when constitutively overexpressed and to control developmental timing of transgene expression. Evaluation of a subset of the resulting ~6000 Arabidopsis transformants revealed a high percentage of lines with full‐length B. napus transgene insertions. Upon induction, numerous iFOX lines with visible phenotypes were identified, including one that displayed early leaf senescence. Phenotypic analysis of this line (rsl‐1327) after methoxyfenozide induction indicated high degree of leaf chlorosis. The integrated B. napuscDNA was identified as a homolog of an Arabidopsis acyl‐CoA binding protein (ACBP) gene designated BnACBP1‐like. The early senescence phenotype conferred by BnACBP1‐like was confirmed by constitutive expression of this gene in Arabidopsis and B. napus. Use of the inducible promoter in the iFOX line coupled with RNA‐Seq analyses allowed mechanistic clues and a working model for the phenotype associated with BnACBP1‐like expression. Our results demonstrate the utility of iFOX‐Hunting as a tool for gene discovery and functional characterization of Brassica napus genome.     

Evolution and selection of Rhg1, a copy‐number variant nematode‐resistance locus

The soybean cyst nematode (SCN) resistance locus Rhg1 is a tandem repeat of a 31.2 kb unit of the soybean genome. Each 31.2‐kb unit contains four genes. One allele of Rhg1, Rhg1‐b, is responsible for protecting most US soybean production from SCN. Whole‐genome sequencing was performed, and PCR assays were developed to investigate allelic variation in sequence and copy number of the Rhg1 locus across a population of soybean germplasm accessions. Four distinct sequences of the 31.2‐kb repeat unit were identified, and some Rhg1 alleles carry up to three different types of repeat unit. The total number of copies of the repeat varies from 1 to 10 per haploid genome. Both copy number and sequence of the repeat correlate with the resistance phenotype, and the Rhg1 locus shows strong signatures of selection. Significant linkage disequilibrium in the genome outside the boundaries of the repeat allowed the Rhg1 genotype to be inferred using high‐density single nucleotide polymorphism genotyping of 15 996 accessions. Over 860 germplasm accessions were found likely to possess Rhg1 alleles. The regions surrounding the repeat show indications of non‐neutral evolution and high genetic variability in populations from different geographic locations, but without evidence of fixation of the resistant genotype. A compelling explanation of these results is that balancing selection is in operation at Rhg1.     

Resistance gene enrichment sequencing (RenSeq) enables reannotation of the NB‐LRR gene family from sequenced plant genomes and rapid mapping of resistance loci in segregating populations

RenSeq is a NB‐LRR (nucleotide binding‐site leucine‐rich repeat) gene‐targeted, Resistance gene enrichment and sequencing method that enables discovery and annotation of pathogen resistance gene family members in plant genome sequences. We successfully applied RenSeq to the sequenced potato Solanum tuberosum clone DM, and increased the number of identified NB‐LRRs from 438 to 755. The majority of these identified R gene loci reside in poorly or previously unannotated regions of the genome. Sequence and positional details on the 12 chromosomes have been established for 704 NB‐LRRs and can be accessed through a genome browser that we provide. We compared these NB‐LRR genes and the corresponding oligonucleotide baits with the highest sequence similarity and demonstrated that ~80% sequence identity is sufficient for enrichment. Analysis of the sequenced tomato S. lycopersicum ‘Heinz 1706’ extended the NB‐LRR complement to 394 loci. We further describe a methodology that applies RenSeq to rapidly identify molecular markers that co‐segregate with a pathogen resistance trait of interest. In two independent segregating populations involving the wild Solanum species S. berthaultii (Rpi‐ber2) and S. ruiz‐ceballosii (Rpi‐rzc1), we were able to apply RenSeq successfully to identify markers that co‐segregate with resistance towards the late blight pathogen Phytophthora infestans. These SNP identification workflows were designed as easy‐to‐adapt Galaxy pipelines.     

Sequence‐tagged high‐density genetic maps of Zoysia japonica provide insights into genome evolution in Chloridoideae

Zoysiagrass (Zoysia spp.), belonging to the genus Zoysia in the subfamily Chloridoideae, is widely used in domestic lawns, sports fields and as forage. We constructed high‐density genetic maps of Zoysia japonica using a restriction site‐associated DNA sequencing (RAD‐Seq) approach and an F₁ mapping population derived from a cross between ‘Carrizo’ and ‘El Toro’. Two linkage maps were constructed, one for each of the parents. A map consisting of 2408 RAD markers distributed on 21 linkage groups was constructed for ‘Carrizo’. Another map with 1230 RAD markers mapped on 20 linkage groups was constructed for ‘El Toro’. The average distance between adjacent markers of the two maps was at 0.56 and 1.4 cM, respectively. Comparative genomics analysis was carried out among zoysiagrass, rice and sorghum genomes and a highly conserved collinearity in the gene order was observed among the three genomes. Chromosome collinearity was disrupted at centromeric regions for each chromosome pair between zoysiagrass and sorghum genomes. However, no obvious synteny gaps were observed across the centromeric regions between zoysiagrass and rice genomes. Two homologous chromosomes for each of the 10 sorghum chromosomes were found in the zoysiagrass genome, indicating an allotetraploid origin for zoysiagrass. The reduction of the basic chromosome number from 12 to 10 in chloridoids and panicoids took place via independent single‐step nested chromosome fusion events after the two subfamilies diverged from a common ancestor. The genetic maps will assist in genome sequence assembly, targeted gene isolation and comparative genomic analyses among grasses.     

Induction and recovery of copy number variation in banana through gamma irradiation and low‐coverage whole‐genome sequencing

Traditional breeding methods are hindered in bananas due to the fact that major cultivars are sterile, parthenocarpic, triploid and thus clonally propagated. This has resulted in a narrow genetic base and limited resilience to biotic and abiotic stresses. Mutagenesis of in vitro propagated bananas is one method to introduce novel alleles and broaden genetic diversity. We previously established a method for the induction and recovery of single nucleotide mutations generated with the chemical mutagen EMS. However, officially released mutant banana varieties have been created using gamma rays, a mutagen that can produce large genomic insertions and deletions (indels). Such dosage mutations may be important for generating observable phenotypes in polyploids. In this study, we establish a low‐coverage whole‐genome sequencing approach in triploid bananas to recover large genomic indels caused by treatment with gamma irradiation. We first evaluated the commercially released mutant cultivar ‘Novaria’ and found that it harbours multiple predicted deletions, ranging from 0.3 to 3.8 million base pairs (Mbp). In total, predicted deletions span 189 coding regions. To evaluate the feasibility of generating and maintaining new mutations, we developed a pipeline for mutagenesis and screening for copy number variation in Cavendish bananas using the cultivar ‘Williams’. Putative mutations were recovered in 70% of lines treated with 20 Gy and 60% of the lines treated with 40 Gy. While deletion events predominate, insertions were identified in 20 Gy‐treated material. Based on these results, we believe this approach can be scaled up to support large breeding projects.     

Population genomic analyses from low‐coverage RAD‐Seq data: a case study on the non‐model cucurbit bottle gourd

Restriction site‐associated DNA sequencing (RAD‐Seq), a next‐generation sequencing‐based genome ‘complexity reduction’ protocol, has been useful in population genomics in species with a reference genome. However, the application of this protocol to natural populations of genomically underinvestigated species, particularly under low‐to‐medium sequencing depth, has not been well justified. In this study, a Bayesian method was developed for calling genotypes from an F₂ population of bottle gourd [Lagenaria siceraria (Mol.) Standl.] to construct a high‐density genetic map. Low‐depth genome shotgun sequencing allowed the assembly of scaffolds/contigs comprising approximately 50% of the estimated genome, of which 922 were anchored for identifying syntenic regions between species. RAD‐Seq genotyping of a natural population comprising 80 accessions identified 3226 single nuclear polymorphisms (SNPs), based on which two sub‐gene pools were suggested for association with fruit shape. The two sub‐gene pools were moderately differentiated, as reflected by the Hudson's F_ST value of 0.14, and they represent regions on LG7 with strikingly elevated F_ST values. Seven‐fold reduction in heterozygosity and two times increase in LD (r²) were observed in the same region for the round‐fruited sub‐gene pool. Outlier test suggested the locus LX3405 on LG7 to be a candidate site under selection. Comparative genomic analysis revealed that the cucumber genome region syntenic to the high F_ST island on LG7 harbors an ortholog of the tomato fruit shape gene OVATE. Our results point to a bright future of applying RAD‐Seq to population genomic studies for non‐model species even under low‐to‐medium sequencing efforts. The genomic resources provide valuable information for cucurbit genome research.     

Functional specialization of duplicated AP3‐like genes in Medicago truncatula

The B–class of MADS box genes has been studied in a wide range of plant species, but has remained largely uncharacterized in legumes. Here we investigate the evolutionary fate of the duplicated AP3‐like genes of a legume species. To obtain insight into the extent to which B‐class MADS box gene functions are conserved or have diversified in legumes, we isolated and characterized the two members of the AP3 lineage in Medicago truncatula: MtNMH7 and MtTM6 (euAP3 and paleoAP3 genes, respectively). A non‐overlapping and complementary expression pattern of both genes was observed in petals and stamens. MtTM6 was expressed predominantly in the outer cell layers of both floral organs, and MtNMH7 in the inner cell layers of petals and stamens. Functional analyses by reverse genetics approaches (RNAi and Tnt1 mutagenesis) showed that the contribution of MtNMH7 to petal identity is more important than that of MtTM6, whereas MtTM6 plays a more important role in stamen identity than its paralog MtNMH7. Our results suggest that the M. truncatula AP3‐like genes have undergone a functional specialization process associated with complete partitioning of gene expression patterns of the ancestral gene lineage. We provide information regarding the similarities and differences in petal and stamen development among core eudicots.     

Molecular characterization of cisgenic lines of apple ‘Gala’ carrying the Rvi6 scab resistance gene

Using resistance genes from a crossable donor to obtain cultivars resistant to diseases and the use of such cultivars in production appears an economically and environmentally advantageous approach. In apple, introgression of resistance genes by classical breeding results in new cultivars, while introducing cisgenes by biotechnological methods maintains the original cultivar characteristics. Recently, plants of the popular apple ‘Gala’ were genetically modified by inserting the apple scab resistance gene Rvi6 (formerly HcrVf2) under control of its own regulatory sequences. This gene is derived from the scab‐resistant apple ‘Florina’ (originally from the wild apple accession Malus floribunda 821). The vector used for genetic modification allowed a postselection marker gene elimination to achieve cisgenesis. In this work, three cisgenic lines were analysed to assess copy number, integration site, expression level and resistance to apple scab. For two of these lines, a single insertion was observed and, despite a very low expression of 0.07‐ and 0.002‐fold compared with the natural expression of ‘Florina’, this was sufficient to induce plant reaction and reduce fungal growth by 80% compared with the scab‐susceptible ‘Gala’. Similar results for resistance and expression analysis were obtained also for the third line, although it was impossible to determine the copy number and TDNA integration site–such molecular characterization is requested by the (EC) Regulation No. 1829/2003, but may become unnecessary if cisgenic crops become exempt from GMO regulation.