Intronic T-DNA insertion in Arabidopsis NBR1 conditionally affects wild-type transcript level

Abstract

The SALK_135513 line of Arabidopsis thaliana is annotated by GenBank to have the T-DNA insertion in the fourth exon of NBR1 (At4g24690). Careful molecular analyses of the homozygous plants of SALK_135513 line indicated the place of T-DNA insertion in the fourth intron. Unexpectedly, 2 kinds of NBR1 transcripts, the wild-type and the mutated, resulting from alternative splicing events, were detected in those plants. Our findings explain the problems encountered by us with phenotypic evaluation of this line and emphasize the necessity for independent verification of the exact insertion site followed by careful expression studies when working with Arabidopsis T-DNA insertional mutants.

Keywords:

• arabidopsis
• At4g24690
• intronic insertion
• NBR1
• SALK_135513
• T-DNA mutant
• transcription

Abbreviations

NBR1	=	Neighbor of BRCA1
WT	=	wild-type
T-DNA	=	transferred DNA
KO	=	knock-out
MV	=	methyl viologen
sqRT-PCR	=	semiquantitative RT-PCR
qRT-PCR	=	quantitative RT-PCR

Introduction

Transferred DNA (T-DNA) from Agrobacterium tumefaciens Ti plasmids is nowadays popular tool for plant genetic manipulation.^1,2 T-DNA insertional mutagenesis helps in identification of the gene responsible for an observed phenotype and allows to establish its function. Currently, there are several collections of the Arabidopsis T-DNA mutants freely distributed to all researchers.³

The protein product of the NBR1 (At4g24690) gene was recently identified as a hybrid ortholog of the mammalian NBR1 (Neighbor of BRCA1) and p62 proteins which are selective autophagy cargo receptors participating in recycling of cellular components.⁴ The selective autophagy process involving NBR1-like proteins in plants is now very intensively studied.^4–6 The phenotypic and genetic studies of the effective knockout (KO) mutants could be very helpful in understanding the in vivo functions of NBR1. Although there are several insertion mutant lines available for NBR1, there are reports that they might be ineffective in knocking out the transcription.⁴ Only recently, a mutant line SALK_135513 was characterized as showing clear phenotypic differences in several stress-generating conditions.⁵ The aim of our study was to further characterize the NBR1 function based on the genetic and phenotypic observations of the SALK_135513 mutant. However, the precise mapping of the insertion site led us to conclude that T-DNA locates in the fourth intron instead of previously annotated fourth exon of NBR1 (GenBank Acc. No BZ384415.1). Additionally, we proved that the wild-type (WT) mRNA of NBR1 is present in SALK_135513 making the mutant ineffective in aborting the transcription. This explains why in our laboratory the mutant plants exhibited no obvious phenotypic and developmental defects under normal and unfavorable growth conditions.

Materials and Methods

The Arabidopsis thaliana SALK_135513 T-DNA insertion line was obtained from NASC. The plants (accession Columbia (WT) and SALK_135513 (homozygous)) were grown in soil. Total DNA and RNA were isolated from rosette leaves of 4- and 8-week-old plants, respectively. Genomic DNA was extracted as described by Edwards et al.⁷ Total RNA (5 μg) isolated independently from 2 leaves per plant with TRI Reagent (Molecular Research Center, Inc.) was used for cDNA synthesis as described earlier.^8,9 The list of oligonucleotides used as primers is shown in Table S1. The PCR conditions were adjusted depending on the primers and are shown in Table S2.

For methyl viologen (MV) treatment, heat, drought and salt stresses the plants were 5 weeks old grown in soil in short-day conditions (8 h light-22°C/16 h dark-19°C; in greenhouse or growth chamber with a light intensity of 157.67 ± 19.78 μE m⁻² s⁻¹). For the oxidative stress the plants were sprayed with 20 μM MV and kept well watered for 2 d under constant light in a growth chamber (181.33 ± 13.32 μE m⁻² s⁻¹ light intensity). The heat stress was conducted by placing the well watered plants for 24 h under 4 Dual 2700 K/6500 K Red/Blue spectrum lamps generating a temperature of 45°C; the pictures were taken after 3 d of recovery. The dark stress was performed by keeping the plants in a complete darkness for 6 d followed by 4 d of recovery. For drought stress plants received no water for 4 d and for salt stress they were kept on 0.6 M NaCl for 4 d and next the pictures were taken.

Results and Discussion

Four PCR primers applied in 3 PCR combinations were used for initial genotyping of the SALK_135513 mutant. The details of PCR1, 2 and 3, each performed with genomic DNA of the WT and 3 individual SALK_135513 plants, are listed in Table S2, while the overall scheme and results of the analysis is shown in Figure 1. The products of the expected size with the primers flanking the T-DNA insertion (PCR1 and PCR3) were obtained only in the case of the WT genomic DNA indicating homozygosity of the SALK_135513 plants. The small amounts of the shorter than expected PCR product in PCR1 in the case of SALK_135513 result from an unspecific reaction as proved by sequencing (data not shown). Furthermore, presence of the PCR2 product in SALK_135513 confirms that the analyzed plants indeed contain T-DNA insert within NBR1.

Figure 1. Mapping of the T-DNA insert in SALK_135513. (A) Scheme of the expected PCR products and their localization in the At4g24690 chromosomal region. (B) PCR results with the pairs of S13LP/S13RP (lanes 1; PCR1), LBb1.3/S13RP (lanes 2; PCR2) and S13LP/S05LP (lanes 3; PCR3) oligonucleotides as primers and genomic DNA isolated from the WT or 3 independent SALK_135513 plants (S13.2, S13.3 and S13.7) as templates. Bands marked with a black arrowhead (S13.7 plant, PCR1) and a white arrowhead (3 SALK_135513 plants, PCR2), were sequenced. Notice, that a non-specific band of 0.6 kb (PCR1) was amplified also in WT; L denotes 1kb DNA Ladder GeneRuler (Fermentas, cat No SMO311), with sizes of the selected bands indicated.

All PCR fragments were sequenced and aligned to the corresponding region of the WT genomic DNA allowing the precise insertion site of the T-DNA (GenBank Acc. No KC820135) (Fig. S1 and Fig. S2). Surprisingly, the insertion mapped to the fourth intron but not to the fourth exon of NBR1 as annotated (GenBank Acc. No BZ384415.1). This incorrect annotation might be a result of the rapid sequencing method used for indexing the SALK mutant collection.

The next step was to investigate if the T-DNA insert affects NBR1 transcription level. It has been reported that only 0.7% of the T-DNA intronic insertions do not have any effect on transcript level.¹⁰ The semi-quantitative RT-PCR (sqRT-PCR) was carried out using 3 combinations of the NBR1-specific oligonucleotides with cDNA prepared from WT and SALK_135513 plants. The expected NBR1-specific PCR products for 2 kinds of the transcripts (WT and truncated) are schematically marked in Figure 2A. Expression of the PetC (photosynthetic electron transfer c; At4g03280) served as a quantitative reference.¹¹ Reactions were performed separately and pooled before loading on to a gel to facilitate analysis. The PCR products of the size expected for the cDNA corresponding to the WT transcript (1522, 1155 and 145 bp) were obtained for all plants indicating the presence of the WT transcript in the homozygous SALK_135513 line (Fig. 2B). Sequencing of those PCR fragments (GenBank Acc No KC848420) confirmed their 100% identity with relevant fragments of WT NBR1 mRNA. The subsequent quantitative RT-PCR (qRT-PCR) using slightly different primers combinations (Fig. 3A) than the combination used for sqRT-PCR indicated the similar amounts of transcripts corresponding to the region upstream to the mapped insertion and much reduced amount of transcripts corresponding to the region downstream the T-DNA insert in the SALK_135513 line (Fig. 3B).

Figure 2. Two variants of NBR1 transcripts in SALK_135513. (A) Scheme of the expected PCR fragments and their localization on the expected mRNAs. The exons of NBR1 (represented as boxes) and the T-DNA position is marked. The size of the expected PCR products is indicated. (B) Results of semiquantitative RT-PCR indicating the presence of the WT transcript in the homozygous SALK_135513 plants. The PCR was performed in duplicate using 2 pools of independently prepared cDNA for the indicated plants (WT or the mutant). The expected sizes of the PCR products in WT are (1522, 1155 and 145 bp). (C) Results of semiquantitative RT-PCR showing the presence of mutated mRNA in SALK_135513, spliced within the T-DNA insert. The oligonucleotide pairs LBb/S05LP (lanes labeled as (1) and TUBF/TUBR (lanes labeled as (2) were used as primers. The PCR products marked with black arrowheads were cut off from the gel and sequenced. L, 1kb DNA Ladder GeneRuler (Fermentas, cat No SMO311).

Figure 3. Quantitative RT-PCR for detection of 3 different parts of NBR1 transcript (Nt, N-terminal region of the cDNA, 79 bp; IR, intermediate region of the cDNA, 116 bp, Ct, C-terminal region of the cDNA, 113 bp) in non-stressed WT (Col-0) and SALK_135513 plants. (A) Schematic location of the expected PCR fragments on the WT and mutated mRNA of NBR1 (boxes represent individual exons); (B) Results of qRT-PCR in WT (Col-0) and SALK_135513 (S13) plants indicating the similar level of Nt but reduced amounts of IR and Ct in SALK_135513 as compared to WT. Data bars represent the mean (±SD) of 2 biological repeats performed in technical triplicate for WT and 4 biological repeats performed in technical triplicate for SALK_135513.

Further work aimed to check if the chimeric transcripts containing truncated NBR1 and fragments of T-DNA insert can be detected in SALK_135513 plants. The results of this analysis are shown in Figure 2C and the expected product of sqRT-PCR (1176 bp) is schematically marked in Figure 2A. The expression of tubulin α3 served as a quantitative reference. The results of the sqRT-PCR confirmed the presence of the chimeric transcripts in SALK_135513. Sequencing of the 1176 bp product indicated that splicing occurred within the T-DNA (removing its 97-nt fragment). The remaining T-DNA stayed inserted into the fourth intron of the NBR1 transcript (GenBank Acc. No KC848419), however the exact length of such truncated mRNA is not known because its sequence downstream of the T-DNA insertion point has not been verified.

We have demonstrated that mixture of the WT and the truncated transcripts exists in homozygous SALK_135513. Such cases have been reported before in T-DNA insertion mutants.¹⁰ The mechanisms of splicing in plants are not completely understood but from preliminary studies¹² and comparison to other organisms it is possible to distinguish some common features which may serve as recognition sites for the plant splicing machinery. There are several features required for splicing within the intron: a donor site (5′ end of the intron), a branch site (near the 3′ end of the intron) and an acceptor site (3′ end of the intron) composed of specific nucleotides. The alignment of the acceptor sites of the introns present in At4g24690 indicated a limited similarity and allowed the distinguishing of 2 groups. The putative site of the splicing in the T-DNA insert resembles the acceptor sites in introns 2, 3 and 6 of NBR1 (Table 1). These sequences coincide with the less restrictive definition criteria of the acceptor site which identified AG as acceptor sites for plants.¹³

Table 1. Alignment of acceptor sites of the introns in NBR1 gene and a sequence of T-DNA insertion functioning as an acceptor site of an intron

3ss	−18	−17	−16	−15	−14	−13	−12	−11	−10	−9	−8	−7	−6	−5	−4	−3	−2	−1	+1	+2	+3
Intron 1	T	T	G	T	T	T	G	T	G	A	T	T	A	T	G	T	A	G	G	T	G
Intron 2	C	T	T	G	T	G	A	A	T	A	T	T	G	T	G	C	A	G	T	A	A
Intron 3	C	T	T	T	T	T	T	T	T	T	T	G	G	T	A	T	A	G	T	T	T
Intron 4	T	C	T	C	T	T	T	A	C	A	C	A	T	G	G	C	A	G	A	T	T
Intron 5	T	T	C	T	T	T	T	T	G	C	A	T	G	G	T	T	A	G	G	T	T
Intron 6	A	A	T	A	T	T	G	T	G	A	T	T	T	T	G	C	A	G	G	G	C
T-DNA	G	G	T	T	T	T	T	C	T	T	T	T	C	A	C	C	A	G	T	G	A

Zhou et al.⁵ have reported that SALK_135513 mutant has increased sensitivity to various abiotic stresses. However, in our laboratory, no visible differences between the WT and SALK_135513 plants were noticed in any of the conditions of applied oxidative, heat, drought, salt or dark stresses (Fig. S3). We followed the expression of both types of NBR1 transcripts (WT and mutated) in SALK_135513 plants challenged by salt and drought stresses in comparison to the control conditions. Interestingly, the ratio between both types of NBR1 mRNA appeared to be differentially affected by stresses (Fig. 4). In plants challenged with the salt stress there was more of the mutated NBR1 produced while in drought stress the splicing machinery favored the WT NBR1 assembly.

Figure 4. Stress-dependent changes of the relative abundance of 2 variants of NBR1 transcript in SALK_135513 plants. Localization of the PCR product specific for SALK_135513 WT (Nt, 145 bp; IR, 116 and T-DNA, 145 bp) transcripts is as in Fig. 3A. The amounts of the indicated PCR products in SALK_135513 plants treated by salt stress (SB) and drought stress (SD) is shown in relation to the amount of the same products in SALK_135513 plants grown in the control conditions (CNTR). Data bars represent the mean (±SD) of 2 biological repeats performed in technical duplicate.

The observations of the phenotypes of SALK_135513 described in this work disagree with phenotypes reported by Zhou et al.⁵ These discrepancies might be attributed to differences in growing conditions that might affect the splicing machinery. Similarly to long intronic non-coding RNA causing epigenetic changes in the level of gene expression, T-DNA insertion may affect the proper splicing.¹⁴

Recently, it was demonstrated for the intronic T-DNA insertion mutant opr3 that it shows clear KO phenotype after wounding while upon fungal infection there is no effect on the normal transcription of OPR3 probably due to the splicing of the T-DNA-containing intron.¹⁵ Ulker et al. 2008¹⁶ reported that introns containing T-DNA inserts are less efficiently spliced than WT introns in the various mutant lines. These findings corroborates our data showing the influence of different conditions on the WT copy of the transcript in the SALK_135512 (Fig. 4).

Conclusions

Our results emphasize the need for detailed molecular analyses of Arabidopsis T-DNA insertion lines. This characterization should include the precise mapping of the T-DNA insert to confirm annotated position within the gene of interest as well as expression analysis with gene-specific primers to ensure the absence of WT mRNA copies in T-DNA mutant line. An extra care has to be taken when working with intronic T-DNA insertions because in some of them the expression may undergo environmental changes complicating interpretations of the mutant phenotype.

Disclosure of Potential Conflicts of Interest

No potential conflicts of interest were disclosed.

Supplemental Material

Supplemental data for this article can be accessed on thepublisher's website.

Additional information

Funding

This work was supported by the Marie Curie Initial Training Network, project No 264296, BIONUT.