Distribution and characterization of over 1000 T-DNA tags in rice genome

We generated T-DNA insertions throughout the rice genome for saturation mutagenesis. More than 1,000 flanking sequences were mapped on 12 rice chromosomes. Our results showed that T-DNA tags were not randomly spread on rice chromosomes and were preferentially inserted in gene-rich regions. Few insertions (2.4%) were found in repetitive regions. T-DNA insertions in genic (58.1%) and intergenic regions (41.9%) showed a good correlation with the predicted size distribution of these sequences in the rice genome. Whereas, obvious biases were found for the insertions in the 5'- and 3'-regulatory regions outside the coding regions both at 500-bp size and in introns rather than in exons. Such distribution patterns and biases for T-DNA integration in rice are similar to that of the previous report in Arabidopsis, which may result from T-DNA integration mechanism itself. Rice will require approximately the same number of T-DNA insertions for saturation mutagenesis as will Arabidopsis. A database of the T-DNA insertion sites in rice is publicly available at our web site (http://www.genomics.zju.edu.cn/ricetdna).     

Characterization of T-DNA insertion sites in Arabidopsis thaliana and the implications for saturation mutagenesis

A key component of a sound functional genomics infrastructure is the availability of a knockout mutant for every gene in the genome. A fruitful approach to systematically knockingout genes in the plant Arabidopsis thaliana has been the use of transferred-DNA (T-DNA) from Agrobacterium tumefaciens as an insertional mutagen. One of the assumptions underlying the use of T-DNA as a mutagen is that the insertion of these DNA elements into the Arabidopsis genome occurs at randomly selected locations. We have directly investigated the distribution of T-DNA insertions sites in populations of transformed Arabidopsis using two different approaches. To begin with, we utilized a polymerase chain reaction (PCR) procedure to systematically catalog the precise locations of all the T-DNA elements inserted within a 65 kb segment of chromosome IV. Of the 47 T-DNA insertions identified, 30% were found within the coding regions of genes. We also documented the insertion of T-DNA elements within the centromeric region of chromosome IV. In addition to these targeted T-DNA screens, we also mapped the genomic locations of 583 randomly chosen T-DNA elements by sequencing the genomic DNA flanking the insertion sites from individual T-DNA-transformed lines. 35% of these randomly chosen T-DNA insertions were located within the coding regions of genes. For comparison, coding sequences account for 44% of the Arabidopsis genome. Our results demonstrate that there is a small bias towards recovering T-DNA insertions within intergenic regions. However, this bias does not limit the utility of T-DNA as an effective insertional mutagen for use in reverse-genetic strategies.     

A high-throughput Arabidopsis reverse genetics system

A collection of Arabidopsis lines with T-DNA insertions in known sites was generated to increase the efficiency of functional genomics. A high-throughput modified thermal asymmetric interlaced (TAIL)-PCR protocol was developed and used to amplify DNA fragments flanking the T-DNA left borders from approximately 100000 transformed lines. A total of 85108 TAIL-PCR products from 52964 T-DNA lines were sequenced and compared with the Arabidopsis genome to determine the positions of T-DNAs in each line. Predicted T-DNA insertion sites, when mapped, showed a bias against predicted coding sequences. Predicted insertion mutations in genes of interest can be identified using Arabidopsis Gene Index name searches or by BLAST (Basic Local Alignment Search Tool) search. Insertions can be confirmed by simple PCR assays on individual lines. Predicted insertions were confirmed in 257 of 340 lines tested (76%). This resource has been named SAIL (Syngenta Arabidopsis Insertion Library) and is available to the scientific community at www.tmri.org.     

Generation and analysis of end sequence database for T-DNA tagging lines in rice

We analyzed 6749 lines tagged by the gene trap vector pGA2707. This resulted in the isolation of 3793 genomic sequences flanking the T-DNA. Among the insertions, 1846 T-DNAs were integrated into genic regions, and 1864 were located in intergenic regions. Frequencies were also higher at the beginning and end of the coding regions and upstream near the ATG start codon. The overall GC content at the insertion sites was close to that measured from the entire rice (Oryza sativa) genome. Functional classification of these 1846 tagged genes showed a distribution similar to that observed for all the genes in the rice chromosomes. This indicates that T-DNA insertion is not biased toward a particular class of genes. There were 764, 327, and 346 T-DNA insertions in chromosomes 1, 4 and 10, respectively. Insertions were not evenly distributed; frequencies were higher at the ends of the chromosomes and lower near the centromere. At certain sites, the frequency was higher than in the surrounding regions. This sequence database will be valuable in identifying knockout mutants for elucidating gene function in rice. This resource is available to the scientific community at http://www.postech.ac.kr/life/pfg/risd.     

Non-random distribution of T-DNA insertions at various levels of the genome hierarchy as revealed by analyzing 13 804 T-DNA flanking sequences from an enhancer-trap mutant library

We isolated 13 804 T-DNA flanking sequence tags (FSTs) from a T-DNA insertion library of rice. A comprehensive analysis of the 13 804 FSTs revealed a number of features demonstrating a highly non-random distribution of the T-DNA insertions in the rice genome: T-DNA insertions were biased towards large chromosomes, not only in the absolute number of insertions but also in the relative density; within chromosomes the insertions occurred more densely in the distal ends, and less densely in the centromeric regions; the distribution of the T-DNA insertions was highly correlated with that of full-length cDNAs, but the correlations were highly heterogeneous among the chromosomes; T-DNA insertions strongly disfavored transposable element (TE)-related sequences, but favored genic sequences with a strong bias toward the 5' upstream and 3' downstream regions of the genes; T-DNA insertions preferentially occurred among the various classes of functional genes, such that the numbers of insertions were in excess in certain functional categories but were deficient in other categories. The analysis of DNA sequence compositions around the T-DNA insertion sites also revealed several prominent features, including an elevated bendability from -200 to 200 bp relative to the insertion sites, an inverse relationship between the GC and TA skews, and reversed GC and TA skews in sequences upstream and downstream of the insertion sites, with both GC and TA skews equal to zero at the insertion sites. It was estimated that 365 380 insertions are needed to saturate the genome with P = 0.95, and that the 45 441 FSTs that have been isolated so far by various groups tagged 14 287 of the 42 653 non-TE related genes.     

An Arabidopsis thaliana T-DNA mutagenized population (GABI-Kat) for flanking sequence tag-based reverse genetics

The GABI-Kat population of T-DNA mutagenized Arabidopsis thaliana lines with sequence-characterized insertion sites is used extensively for efficient progress in plant functional genomics. Here we provide details about the establishment of the material, demonstrate the population's functionality and discuss results from quality control studies. T-DNA insertion mutants of the accession Columbia (Col-0) were created by Agrobacterium tumefaciens-mediated transformation. To allow selection of transformed plants under greenhouse conditions, a sulfadiazine resistance marker was employed. DNA from leaves of T1 plants was extracted and used as a template for PCR-based amplification of DNA fragments spanning insertion site borders. After sequencing, the data were placed in a flanking sequence tag (FST) database describing which mutant allele was present in which line. Analysis of the distribution of T-DNA insertions revealed a clear bias towards intergenic regions. Insertion sites appeared more frequent in regions in front of the ATG and after STOP codons of predicted genes. Segregation analysis for sulfadiazine resistance showed that 62% of the transformants contain an insertion at only one genetic locus. In quality control studies with gene-specific primers in combination with T-DNA primers, 76% of insertions could be confirmed. Finally, the functionality of the GABI-Kat population was demonstrated by exemplary confirmation of several new transparent testa alleles, as well as a number of other mutants, which were identified on the basis of the FST data.     

Natural selection maintains the transcribed LTR retrotransposons in Nosema bombycis

Eight intact LTR retrotransposons(Nbr1?Nbr8)have been previously characterized from the genome of Nosema bombycis,a eu-karyotic parasite with a compact and reduced genome.Here we describe six novel transcribed Nbr elements(Nbr9?Nbr14)identified through either cDNA library or RT-PCR.Like previously determined ones,all of them belong to the Ty3/Gypsy superfamily.Retrotransposon diversity and incomplete domains with insertions(Nbr12),deletions(Nbr11)and in-frame stop codons in coding regions(Nbr9)were detected...     

Roles of the sequence encoding tobacco etch virus capsid protein in genome amplification: requirements for the translation process and a cis-active element.

The roles of the capsid protein (CP) and the CP coding sequence of tobacco etch potyvirus (TEV) in genome amplification were analyzed. A series of frameshift-stop codon mutations that interrupted translation of the CP coding sequence at various positions were introduced into the TEV genome. A series of 3' deletion mutants that lacked the CP coding sequence beyond each of the frameshift-stop codon mutations were also produced. In addition, a series of 5' CP deletion mutants were generated. Amplification of genomes containing either frameshift-stop codon insertions after codons 1, 59, 103, and 138 or genomes containing the corresponding 3' deletions of the CP coding sequence was reduced by 100- to 1,000-fold relative to that of the parental genome in inoculated protoplasts. In contrast, a mutant containing a frameshift-stop codon after CP position 189 was amplified to 27% of the level of the parental virus, but the corresponding 3' deletion mutant lacking codons 190 to 261 was nonviable. Deletion mutants lacking CP codons 2 to 100, 2 to 150, 2 to 189, and 2 to 210 were amplified relatively efficiently in protoplasts, but a deletion mutant lacking codons 2 to 230 was nonviable. None of the amplification-defective frameshift-stop codon or deletion mutants was rescued in transgenic cells expressing TEV CP, although the transgenic CP was able to rescue intercellular movement defects of replication-competent CP mutants. Coupled with previous results, these data led to the conclusions that (i) TEV genome amplification requires translation to a position between CP codons 138 and 189 but does not require the CP product and (ii) the TEV CP coding sequence contains a cis-active RNA element between codons 211 and 246. The implications of these findings on mechanisms of RNA replication and genome evolution are discussed.     

A systematic analysis of T-DNA insertion events in Magnaporthe oryzae

We describe here the analysis of random T-DNA insertions that were generated as part of a large-scale insertional mutagenesis project for Magnaporthe oryzae. Chromosomal regions flanking T-DNA insertions were rescued by inverse PCR, sequenced and used to search the M. oryzae genome assembly. Among the 175 insertions for which at least one flank was rescued, 137 had integrated in single-copy regions of the genome, 17 were in repeated sequences, one had no match to the genome, and the remainder were unassigned due to illegitimate T-DNA integration events. These included in order of abundance: head-to-tail tandem insertions, right border excision failures, left border excision failures and insertion of one T-DNA into another. The left borders of the T-DNA were frequently truncated and inserted in sequences with micro-homology to the left terminus. By contrast the right borders were less prone to degradation and appeared to have been integrated in a homology-independent manner. Gross genome rearrangements rarely occurred when the T-DNAs integrated in single-copy regions, although most insertions did cause small deletions at the target site. Significant insertion bias was detected, with promoters receiving two times more T-DNA hits than expected, and open reading frames receiving three times fewer. In addition, we found that the distribution of T-DNA inserts among the M. oryzae chromosomes was not random. The implications of these findings with regard to saturation mutagenesis of the M. oryzae genome are discussed.     

Analysis of Genomic G + C Content, Codon Usage, Initiator Codon Context and Translation Termination Sites In Tetrahymena Thermophila

In recent years, the amount of molecular sequencing data from Tetrahymena thermophila has dramatically increased. We analyzed G + C content, codon usage, initiator codon context and stop codon sites in the extremely A + T rich genome of this ciliate. Average G + C content was 38% for protein coding regions, 21% for 5' non-coding sequences, 19% for 3' non-coding sequences, 15% for introns, 19% for micronuclear limited sequences and 17% for macronuclear retained sequences flanking micronuclear specific regions. The 75 available T. thermophila protein coding sequences favored codons ending in T and, where possible, avoided those with G in the third position. Highly expressed genes were relatively G + C-rich and exhibited an extremely biased pattern of codon usage while developmentally regulated genes were more A + T-rich and showed less codon usage bias. Regions immediately preceding Tetrahymena translation initiator codons were generally A-rich. For the 60 stop codons examined, the frequency of G in the end + 1 site was much higher than expected whereas C never occupied this position.     

A comprehensive characterization of single-copy T-DNA insertions in the Arabidopsis thaliana genome

T-DNA flanking sequences were isolated from 112 Arabidopsis thaliana single-copy T-DNA lines and sequence mapped to the chromosomes. Even though two T-DNA insertions mapped to a heterochromatic domain located in the pericentromeric region of chromosome I, expression of reporter genes was detected in these transgenic lines. T-DNA insertion did not seem to be biased toward any of Arabidopsis' five chromosomes. The observed distribution of T-DNA copies in intergenic sequence versus gene sequence (i.e. 5-upstream regions, coding sequences and 3-downstream regions) appeared randomly. An evaluation of T-DNA insertion frequencies within gene sequence revealed that integration into 5-upstream regions occurred more frequently than expected, whereas insertions in coding sequences (exons and introns) were found less frequently than expected based on random distribution predictions. In the majority of cases, single-copy T-DNA insertions were associated with small or large rearrangements such as deletions and/or duplications of target site sequences, deletions and/or duplications of T-DNA sequences, and gross chromosomal rearrangements such as translocations. The accuracy of integration was similarly high for both left- and right-border sequences. These results may be called upon when making detailed molecular analyses of transgenic plants or T-DNA induced mutants.     

Partial shotgun sequencing of the Boechera stricta genome reveals extensive microsynteny and promoter conservation with Arabidopsis

Comparative genomics provides insight into the evolutionary dynamics that shape discrete sequences as well as whole genomes. To advance comparative genomics within the Brassicaceae, we have end sequenced 23,136 medium-sized insert clones from Boechera stricta, a wild relative of Arabidopsis (Arabidopsis thaliana). A significant proportion of these sequences, 18,797, are nonredundant and display highly significant similarity (BLASTn e-value < or = 10(-30)) to low copy number Arabidopsis genomic regions, including more than 9,000 annotated coding sequences. We have used this dataset to identify orthologous gene pairs in the two species and to perform a global comparison of DNA regions 5' to annotated coding regions. On average, the 500 nucleotides upstream to coding sequences display 71.4% identity between the two species. In a similar analysis, 61.4% identity was observed between 5' noncoding sequences of Brassica oleracea and Arabidopsis, indicating that regulatory regions are not as diverged among these lineages as previously anticipated. By mapping the B. stricta end sequences onto the Arabidopsis genome, we have identified nearly 2,000 conserved blocks of microsynteny (bracketing 26% of the Arabidopsis genome). A comparison of fully sequenced B. stricta inserts to their homologous Arabidopsis genomic regions indicates that indel polymorphisms >5 kb contribute substantially to the genome size difference observed between the two species. Further, we demonstrate that microsynteny inferred from end-sequence data can be applied to the rapid identification and cloning of genomic regions of interest from nonmodel species. These results suggest that among diploid relatives of Arabidopsis, small- to medium-scale shotgun sequencing approaches can provide rapid and cost-effective benefits to evolutionary and/or functional comparative genomic frameworks.     

Efficient Identification of Arabidopsis Knock-Out Mutants Using DNA-Arrays of Transposon Flanking Sequences

Abstract: Mutants obtained by insertional mutagenesis are widely used for determining gene-phenotype relationships. In Arabidopsis thaliana, several populations mutagenized either by T-DNA or transposon insertion are available for screening for knockout mutants in genes of interest. We have so far screened our Arabidopsis population mutagenized with the Zea mays transposon En-1/Spm for insertion mutations in 718 genes, using PCR on DNA pools. Although successful, this common approach is too time consuming for use in systematic screening of all 25 498 predicted genes of the Arabidopsis genome. We therefore investigated the use of DNA arrays for the direct identification of mutants in our population. All transposon-flanking regions from individual plants are amplified by PCR and subsequently spotted at high density onto nylon membranes. A single hybridization experiment with a gene-specific probe then allows one to identify candidate mutant plants. The efficiency of each separate step was determined and optimized. Screening of filters representing 2880 plants for insertions in 144 genes and subsequent investigation of some of the potential insertion mutants suggest that an overall screening efficiency of 50 % is attained.     

Patterns of selection against transposons inferred from the distribution of Tc1, Tc3 and Tc5 insertions in the mut-7 line of the nematode Caenorhabditis elegans

To identify the factors (selective or mutational) that affect the distribution of transposable elements (TEs) within a genome, it is necessary to compare the pattern of newly arising element insertions to the pattern of element insertions that have been fixed in a population. To do this, we analyzed the distribution of recent mutant insertions of the Tc1, Tc3, and Tc5 elements in a mut-7 background of the nematode Caenorhabditis elegans and compared it to the distribution of element insertions (presumably fixed) within the sequenced genome. Tc1 elements preferentially insert in regions with high recombination rates, whereas Tc3 and Tc5 do not. Although Tc1 and Tc3 both insert in TA dinucleotides, there is no clear relationship between the frequency of insertions and the TA dinucleotide density. There is a strong selection against TE insertions within coding regions: the probability that a TE will be fixed is at least 31 times lower in coding regions than in noncoding regions. Contrary to the prediction of theoretical models, we found that the selective pressure against TE insertions does not increase with the recombination rate. These findings indicate that the distribution of these three transposon families in the genome of C. elegans is determined essentially by just two factors: the pattern of insertions, which is a characteristic of each family, and the selection against insertions within coding regions.     

Establishing an efficient Ac/Ds tagging system in rice: large-scale analysis of Ds flanking sequences

A two-element Activator/Dissociation (Ac/Ds) gene trap system was successfully established in rice (Oryza sativa ssp. japonica cv. Nipponbare) to generate a collection of stable, unlinked and single-copy Ds transposants. The germinal transposition frequency of Ds was estimated as an average of 51% by analyzing 4413 families. Study of Ds transposition pattern in siblings revealed that 79% had at least two different insertions, suggesting late transposition during rice development. Analysis of 2057 Ds flanking sequences showed that 88% of them were unique, whereas the rest within T-DNA. The insertions were distributed randomly throughout the genome; however, there was a bias toward chromosomes 4 and 7, which had two times as many insertions as that expected. A hot spot for Ds insertions was identified on chromosome 7 within a 40-kbp region. One-third of Ds flanking sequences was homologous to either proteins or rice expressed sequence tags (ESTs), confirming a preference for Ds transposition into coding regions. Analysis of 200 Ds lines on chromosome 1 revealed that 72% insertions were found in genic region. Anchoring of more than 800 insertions to yeast artificial chromosome (YAC)-based EST map showed that Ds transposes preferentially into regions rich in expressed sequences. High germinal transposition frequency and independent transpositions among siblings show that the efficiency of this system is suitable for large-scale transposon mutagenesis in rice.