Dynamic growth QTL action in diverse light environments: characterization of light regime-specific and stable QTL in Arabidopsis.

Plant growth is a complex process affected by a multitude of genetic and environmental factors and their interactions. To identify genetic factors influencing plant performance under different environmental conditions, vegetative growth was assessed in Arabidopsis thaliana cultivated under constant or fluctuating light intensities, using high-throughput phenotyping and genome-wide association studies. Daily automated non-invasive phenotyping of a collection of 382 Arabidopsis accessions provided growth data during developmental progression under different light regimes at high temporal resolution. Quantitative trait loci (QTL) for projected leaf area, relative growth rate, and PSII operating efficiency detected under the two light regimes were predominantly condition-specific and displayed distinct temporal activity patterns, with active phases ranging from 2 d to 9 d. Eighteen protein-coding genes and one miRNA gene were identified as potential candidate genes at 10 QTL regions consistently found under both light regimes. Expression patterns of three candidate genes affecting projected leaf area were analysed in time-series experiments in accessions with contrasting vegetative leaf growth. These observations highlight the importance of considering both environmental and temporal patterns of QTL/allele actions and emphasize the need for detailed time-resolved analyses under diverse well-defined environmental conditions to effectively unravel the complex and stage-specific contributions of genes affecting plant growth processes.

A Similar Genetic Architecture Underlies the Convergent Evolution of the Selfing Syndrome in Capsella.

Whether, and to what extent, phenotypic evolution follows predictable genetic paths remains an important question in evolutionary biology. Convergent evolution of similar characters provides a unique opportunity to address this question. The transition to selfing and the associated changes in flower morphology are among the most prominent examples of repeated evolution in plants. In this study, we take advantage of the independent transitions to self-fertilization in the genus Capsella to compare the similarities between parallel modifications of floral traits and test for genetic and developmental constraints imposed on flower evolution in the context of the selfing syndrome. Capsella rubella and Capsella orientalis emerged independently but evolved almost identical flower characters. Not only is the evolutionary outcome identical but the same developmental strategies underlie the convergent reduction of flower size. This has been associated with convergent evolution of gene expression changes. The transcriptomic changes common to both selfing lineages are enriched in genes with low network connectivity and with organ-specific expression patterns. Comparative genetic mapping also suggests that, at least in the case of petal size evolution, these similarities have a similar genetic basis. Based on these results, we hypothesize that the limited availability of low-pleiotropy paths predetermines closely related species to similar evolutionary outcomes.

Combining next-generation sequencing and progeny testing for rapid identification of induced recessive and dominant mutations in maize M2 individuals.

Molecular identification of mutant alleles responsible for certain phenotypic alterations is a central goal of genetic analyses. In this study we describe a rapid procedure suitable for the identification of induced recessive and dominant mutations applied to two Zea mays mutants expressing a dwarf and a pale green phenotype, respectively, which were obtained through pollen ethyl methanesulfonate (EMS) mutagenesis. First, without prior backcrossing, induced mutations (single nucleotide polymorphisms, SNPs) segregating in a (M2 ) family derived from a heterozygous (M1 ) parent were identified using whole-genome shotgun (WGS) sequencing of a small number of (M2 ) individuals with mutant and wild-type phenotypes. Second, the state of zygosity of the mutation causing the phenotype was determined for each sequenced individual by phenotypic segregation analysis of the self-pollinated (M3 ) offspring. Finally, we filtered for segregating EMS-induced SNPs whose state of zygosity matched the determined state of zygosity of the mutant locus in each sequenced (M2 ) individuals. Through this procedure, combining sequencing of individuals and Mendelian inheritance, three and four SNPs in linkage passed our zygosity filter for the homozygous dwarf and heterozygous pale green mutation, respectively. The dwarf mutation was found to be allelic to the an1 locus and caused by an insertion in the largest exon of the AN1 gene. The pale green mutation affected the nuclear W2 gene and was caused by a non-synonymous amino acid exchange in encoded chloroplast DNA polymerase with a predicted deleterious effect. This coincided with lower cpDNA levels in pale green plants.

The poly(A) polymerase PAPS1 interacts with the RNA-directed DNA-methylation pathway in sporophyte and pollen development.

RNA-based processes play key roles in the regulation of eukaryotic gene expression. This includes both the processing of pre-mRNAs into mature mRNAs ready for translation and RNA-based silencing processes, such as RNA-directed DNA methylation (RdDM). Polyadenylation of pre-mRNAs is one important step in their processing and is carried out by three functionally specialized canonical nuclear poly(A) polymerases in Arabidopsis thaliana. Null mutations in one of these, termed PAPS1, result in a male gametophytic defect. Using a fluorescence-labelling strategy, we have characterized this defect in more detail using RNA and small-RNA sequencing. In addition to global defects in the expression of pollen-differentiation genes, paps1 null-mutant pollen shows a strong overaccumulation of transposable element (TE) transcripts, yet a depletion of 21- and particularly 24-nucleotide-long short interfering RNAs (siRNAs) and microRNAs (miRNAs) targeting the corresponding TEs. Double-mutant analyses support a specific functional interaction between PAPS1 and components of the RdDM pathway, as evident from strong synergistic phenotypes in mutant combinations involving paps1, but not paps2 paps4, mutations. In particular, the double-mutant of paps1 and rna-dependent rna polymerase 6 (rdr6) shows a synergistic developmental phenotype disrupting the formation of the transmitting tract in the female gynoecium. Thus, our findings in A. thaliana uncover a potentially general link between canonical poly(A) polymerases as components of mRNA processing and RdDM, reflecting an analogous interaction in fission yeast.

Discovery of key regulators of dark gland development and hypericin biosynthesis in St. John's Wort (Hypericum perforatum).

Hypericin is a molecule of high pharmaceutical importance that is synthesized and stored in dark glands (DGs) of St. John's Wort (Hypericum perforatum). Understanding which genes are involved in dark gland development and hypericin biosynthesis is important for the development of new Hypericum extracts that are highly demanded for medical applications. We identified two transcription factors whose expression is strictly synchronized with the differentiation of DGs. We correlated the content of hypericin, pseudohypericin, endocrocin, skyrin glycosides and several flavonoids with gene expression and DG development to obtain a revised model for hypericin biosynthesis. Here, we report for the first time genotypes which are polymorphic for the presence/total absence (G+/G-) of DGs in their placental tissues (PTs). DG development was characterized in PTs using several microscopy techniques. Fourier transform infrared microscopy was established as a novel method to precisely locate polyaromatic compounds, such as hypericin, in plant tissues. In addition, we obtained transcriptome and metabolome profiles of unprecedented resolution in Hypericum. This study addresses for the first time the development of dark glands and identifies genes that constitute strong building blocks for the further elucidation of hypericin synthesis, its manipulation in plants, its engineering in microbial systems and its applications in medical research.

EFFECTOR OF TRANSCRIPTION factors are novel plant-specific regulators associated with genomic DNA methylation in Arabidopsis.

Plant-specific EFFECTORS OF TRANSCRIPTION (ET) are characterised by a variable number of highly conserved ET repeats, which are involved in zinc and DNA binding. In addition, ETs share a GIY-YIG domain, involved in DNA nicking activity. It was hypothesised that ETs might act as epigenetic regulators. Here, methylome, transcriptome and phenotypic analyses were performed to investigate the role of ET factors and their involvement in DNA methylation in Arabidopsis thaliana. Comparative DNA methylation and transcriptome analyses in flowers and seedlings of et mutants revealed ET-specific differentially expressed genes and mostly independently characteristic, ET-specific differentially methylated regions. Loss of ET function results in pleiotropic developmental defects. The accumulation of cyclobutane pyrimidine dimers after ultraviolet stress in et mutants suggests an ET function in DNA repair.

A Poly(A) Ribonuclease Controls the Cellotriose-Based Interaction between Piriformospora indica and Its Host Arabidopsis.

Piriformospora indica, an endophytic root-colonizing fungus, efficiently promotes plant growth and induces resistance to abiotic stress and biotic diseases. P. indica fungal cell wall extract induces cytoplasmic calcium elevation in host plant roots. Here, we show that cellotriose (CT) is an elicitor-active cell wall moiety released by P. indica into the medium. CT induces a mild defense-like response, including the production of reactive oxygen species, changes in membrane potential, and the expression of genes involved in growth regulation and root development. CT-based cytoplasmic calcium elevation in Arabidopsis (Arabidopsis thaliana) roots does not require the BAK1 coreceptor or the putative Ca2+ channels TPC1, GLR3.3, GLR2.4, and GLR2.5 and operates synergistically with the elicitor chitin. We identified an ethyl methanesulfonate-induced mutant (cytoplasmiccalcium elevation mutant) impaired in the response to CT and various other cellooligomers (n = 2-7), but not to chitooligomers (n = 4-8), in roots. The mutant contains a single nucleotide exchange in the gene encoding a poly(A) ribonuclease (AtPARN; At1g55870) that degrades the poly(A) tails of specific mRNAs. The wild-type PARN cDNA, expressed under the control of a 35S promoter, complements the mutant phenotype. Our identification of cellotriose as a novel chemical mediator casts light on the complex P. indica-plant mutualistic relationship.

Transcriptomic basis for reinforcement of elm antiherbivore defence mediated by insect egg deposition.

Plant responses to insect egg depositions are known to shape subsequent defensive responses to larvae hatching from the eggs. Elm (Ulmus minor) leaves, on which elm leaf beetles laid their eggs, mount a more efficient defence against larvae hatching from the eggs. However, the molecular mechanisms of this egg-mediated, improved defence are insufficiently understood and have so far only been studied in annual plants. We analysed the dynamics of transcriptomic changes in larval feeding-damaged elm leaves with and without prior egg deposition using de novo assembled RNA-seq data. Compared to egg-free leaves, egg deposition-treated leaves showed earlier and/or faster transcriptional regulations, as well as slightly enhanced differential transcriptional regulation after the onset of larval feeding. These early responding transcripts were overrepresented in gene ontology terms associated with post-translational protein modification, signalling and stress (defence) responses. We found evidence of transcriptional memory in initially egg deposition-induced transcripts whose differential expression was reset prior to larval hatching, but was more rapidly induced again by subsequent larval feeding. This potential memory effect of prior egg deposition, as well as the earlier/faster and enhanced feeding-induced differential regulation of transcripts in egg deposition-treated leaves, may contribute to the egg-mediated reinforcing effect on the elm's defence against larvae. Hence, our study shows that a plant's experience of a stress-indicating environmental cue (here: insect eggs) can push the dynamics of the plant's transcriptomic response to subsequent stress (here: larval feeding). Such experience-mediated acceleration of a stress-induced plant response may result in improved stress resistance.

RWP-RK domain-containing transcription factors control cell differentiation during female gametophyte development in Arabidopsis.

The formation of gametes is a prerequisite for any sexually reproducing organism in order to complete its life cycle. In plants, female gametes are formed in a multicellular tissue, the female gametophyte or embryo sac. Although the events leading to the formation of the female gametophyte have been morphologically characterized, the molecular control of embryo sac development remains elusive. We used single and double mutants as well as cell-specific marker lines to characterize a novel class of gene regulators in Arabidopsis thaliana, the RWP-RK domain-containing (RKD) transcription factors. Morphological and histological analyses were conducted using confocal laser scanning and differential interference contrast microscopy. Gene expression and transcriptome analyses were performed using quantitative reverse transcription-PCR and RNA sequencing, respectively. Our results showed that RKD genes are expressed during distinct stages of embryo sac development. Morphological analysis of the mutants revealed severe distortions in gametophyte polarity and cell differentiation. Transcriptome analysis revealed changes in the expression of several gametophyte-specific gene families (RKD2 and RKD3) and ovule development-specific genes (RKD3), and identified pleiotropic effects on phytohormone pathways (RKD5). Our data provide novel insight into the regulatory control of female gametophyte development. RKDs are involved in the control of cell differentiation and are required for normal gametophytic development.

Rye B chromosomes encode a functional Argonaute-like protein with in vitro slicer activities similar to its A chromosome paralog.

B chromosomes (Bs) are supernumerary, dispensable parts of the nuclear genome, which appear in many different species of eukaryote. So far, Bs have been considered to be genetically inert elements without any functional genes. Our comparative transcriptome analysis and the detection of active RNA polymerase II (RNAPII) in the proximity of B chromatin demonstrate that the Bs of rye (Secale cereale) contribute to the transcriptome. In total, 1954 and 1218 B-derived transcripts with an open reading frame were expressed in generative and vegetative tissues, respectively. In addition to B-derived transposable element transcripts, a high percentage of short transcripts without detectable similarity to known proteins and gene fragments from A chromosomes (As) were found, suggesting an ongoing gene erosion process. In vitro analysis of the A- and B-encoded AGO4B protein variants demonstrated that both possess RNA slicer activity. These data demonstrate unambiguously the presence of a functional AGO4B gene on Bs and that these Bs carry both functional protein coding genes and pseudogene copies. Thus, B-encoded genes may provide an additional level of gene control and complexity in combination with their related A-located genes. Hence, physiological effects, associated with the presence of Bs, may partly be explained by the activity of B-located (pseudo)genes.

Genes on B chromosomes: old questions revisited with new tools.

BACKGROUND: B chromosomes are supernumerary dispensable parts of the karyotype which appear in some individuals of some populations in some species. Often, they have been considered as 'junk DNA' or genomic parasites without functional genes. SCOPE OF REVIEW: Due to recent advances in sequencing technologies, it became possible to investigate their DNA composition, transcriptional activity and effects on the host transcriptome profile in detail. Here, we review the most recent findings regarding the gene content of B chromosomes and their transcriptional activities and discuss these findings in the context of comparable biological phenomena, like sex chromosomes, aneuploidy and pseudogenes. MAJOR CONCLUSIONS: Recent data suggest that B chromosomes carry transcriptionally active genic sequences which could affect the transcriptome profile of their host genome. GENERAL SIGNIFICANCE: These findings are gradually changing our view that B chromosomes are solely genetically inert selfish elements without any functional genes. This at one side could partly explain the deleterious effects which are associated with their presence. On the other hand it makes B chromosome a nice model for studying regulatory mechanisms of duplicated genes and their evolutionary consequences.

Sequencing of 15 622 gene-bearing BACs clarifies the gene-dense regions of the barley genome.

Barley (Hordeum vulgare L.) possesses a large and highly repetitive genome of 5.1 Gb that has hindered the development of a complete sequence. In 2012, the International Barley Sequencing Consortium released a resource integrating whole-genome shotgun sequences with a physical and genetic framework. However, because only 6278 bacterial artificial chromosome (BACs) in the physical map were sequenced, fine structure was limited. To gain access to the gene-containing portion of the barley genome at high resolution, we identified and sequenced 15 622 BACs representing the minimal tiling path of 72 052 physical-mapped gene-bearing BACs. This generated ~1.7 Gb of genomic sequence containing an estimated 2/3 of all Morex barley genes. Exploration of these sequenced BACs revealed that although distal ends of chromosomes contain most of the gene-enriched BACs and are characterized by high recombination rates, there are also gene-dense regions with suppressed recombination. We made use of published map-anchored sequence data from Aegilops tauschii to develop a synteny viewer between barley and the ancestor of the wheat D-genome. Except for some notable inversions, there is a high level of collinearity between the two species. The software HarvEST:Barley provides facile access to BAC sequences and their annotations, along with the barley-Ae. tauschii synteny viewer. These BAC sequences constitute a resource to improve the efficiency of marker development, map-based cloning, and comparative genomics in barley and related crops. Additional knowledge about regions of the barley genome that are gene-dense but low recombination is particularly relevant.

Plantago lagopus B Chromosome Is Enriched in 5S rDNA-Derived Satellite DNA.

B chromosomes are supernumerary dispensable parts of the karyotype which appear in some individuals of some populations in some species. Using advanced sequencing technology, we in silico characterized the high-copy DNA composition of Plantago lagopus with and without B chromosomes. The nuclear genome (2.46 pg/2C) was found to be relatively rich in repetitive sequences, with highly and moderately repeated elements making up 68% of the genome. Besides a centromere-specific marker, we identified a B-specific satellite and a repeat enriched in polymorphic A chromosome segments. The B-specific tandem repeat PLsatB originated from sequence amplification including 5S rDNA fragments.

Combination of transcriptomic and proteomic approaches helps to unravel the protein composition of Chelidonium majus L. milky sap.

MAIN CONCLUSION: A novel annotated Chelidonium majus L. transcriptome database composed of 23,004 unique coding sequences allowed to significantly improve the sensitivity of proteomic C. majus assessments, which showed novel defense-related proteins characteristic to its latex. To date, the composition of Chelidonium majus L. milky sap and biosynthesis of its components are poorly characterized. We, therefore, performed de novo sequencing and assembly of C. majus transcriptome using Illumina technology. Approximately, 119 Mb of raw sequence data was obtained. Assembly resulted in 107,088 contigs, with N50 of 1913 bp and N90 of 450 bp. Among 34,965 unique coding sequences (CDS), 23,004 obtained CDS database served as a basis for further proteomic analyses. The database was then used for the identification of proteins from C. majus milky sap, and whole plant extracts analyzed using liquid chromatography-electrospray ionization-tandem mass spectrometry (LC-ESI-MS/MS) approach. Of about 334 different putative proteins were identified in C. majus milky sap and 1155 in C. majus whole plant extract. The quantitative comparative analysis confirmed that C. majus latex contains proteins connected with response to stress conditions and generation of precursor metabolites and energy. Notable proteins characteristic to latex include major latex protein (MLP, presumably belonging to Bet v1-like superfamily), polyphenol oxidase (PPO, which could be responsible for browning of the sap after exposure to air), and enzymes responsible for anthocyanidin, phenylpropanoid, and alkaloid biosynthesis.

Toward the identification and regulation of the Arabidopsis thaliana ABI3 regulon.

The plant-specific, B3 domain-containing transcription factor ABSCISIC ACID INSENSITIVE3 (ABI3) is an essential component of the regulatory network controlling the development and maturation of the Arabidopsis thaliana seed. Genome-wide chromatin immunoprecipitation (ChIP-chip), transcriptome analysis, quantitative reverse transcriptase-polymerase chain reaction and a transient promoter activation assay have been combined to identify a set of 98 ABI3 target genes. Most of these presumptive ABI3 targets require the presence of abscisic acid for their activation and are specifically expressed during seed maturation. ABI3 target promoters are enriched for G-box-like and RY-like elements. The general occurrence of these cis motifs in non-ABI3 target promoters suggests the existence of as yet unidentified regulatory signals, some of which may be associated with epigenetic control. Several members of the ABI3 regulon are also regulated by other transcription factors, including the seed-specific, B3 domain-containing FUS3 and LEC2. The data strengthen and extend the notion that ABI3 is essential for the protection of embryonic structures from desiccation and raise pertinent questions regarding the specificity of promoter recognition.

Elongation-related functions of LEAFY COTYLEDON1 during the development of Arabidopsis thaliana.

The transcription factor LEAFY COTYLEDON1 (LEC1) controls aspects of early embryogenesis and seed maturation in Arabidopsis thaliana. To identify components of the LEC1 regulon, transgenic plants were derived in which LEC1 expression was inducible by dexamethasone treatment. The cotyledon-like leaves and swollen root tips developed by these plants contained seed-storage compounds and resemble the phenotypes produced by increased auxin levels. In agreement with this, LEC1 was found to mediate up-regulation of the auxin synthesis gene YUCCA10. Auxin accumulated primarily in the elongation zone at the root-hypocotyl junction (collet). This accumulation correlates with hypocotyl growth, which is either inhibited in LEC1-induced embryonic seedlings or stimulated in the LEC1-induced long-hypocotyl phenotype, therefore resembling etiolated seedlings. Chromatin immunoprecipitation analysis revealed a number of phytohormone- and elongation-related genes among the putative LEC1 target genes. LEC1 appears to be an integrator of various regulatory events, involving the transcription factor itself as well as light and hormone signalling, especially during somatic and early zygotic embryogenesis. Furthermore, the data suggest non-embryonic functions for LEC1 during post-germinative etiolation.

Promoters of the barley germin-like GER4 gene cluster enable strong transgene expression in response to pathogen attack.

Immunity of plants triggered by pathogen-associated molecular patterns (PAMPs) is based on the execution of an evolutionarily conserved defense response that includes the accumulation of pathogenesis-related (PR) proteins as well as multiple other defenses. The most abundant PR transcript of barley (Hordeum vulgare) leaf epidermis attacked by the powdery mildew fungus Blumeria graminis f. sp hordei encodes the germin-like protein GER4, which has superoxide dismutase activity and functions in PAMP-triggered immunity. Here, we show that barley GER4 is encoded by a dense cluster of tandemly duplicated genes (GER4a-h) that underwent several cycles of duplication. The genomic organization of the GER4 locus also provides evidence for repeated gene birth and death cycles. The GER4 promoters contain multiple WRKY factor binding sites (W-boxes) preferentially located in promoter fragments that were exchanged between subfamily members by gene conversion. Mutational analysis of TATA-box proximal W-boxes used GER4c promoter-beta-glucuronidase fusions to reveal their enhancing effects and functional redundancy on pathogen-induced promoter activity. The data suggest enhanced transcript dosage as an evolutionary driving force for the local expansion and functional redundancy of the GER4 locus. In addition, the GER4c promoter provides a tool to study signal transduction of PAMP-triggered immunity and to engineer strictly localized and pathogen-regulated disease resistance in transgenic cereal crops.

Members of the RKD transcription factor family induce an egg cell-like gene expression program.

In contrast to animals, the life cycle of higher plants alternates between a gamete-producing (gametophyte) and a spore-producing (sporophyte) generation. The female gametophyte of angiosperms consists of four distinct cell types, including two gametes, the egg and the central cell, which give rise to embryo and endosperm, respectively. Based on a combined subtractive hybridization and virtual subtraction approach in wheat (Triticum aestivum L.), we have isolated a class of transcription factors not found in animal genomes, the RKD (RWP-RK domain-containing) factors, which share a highly conserved RWP-RK domain. Single-cell RT-PCR revealed that the genes TaRKD1 and TaRKD2 are preferentially expressed in the egg cell of wheat. The Arabidopsis genome contains five RKD genes, at least two of them, AtRKD1 and AtRKD2, are preferentially expressed in the egg cell of Arabidopsis. Ectopic expression of the AtRKD1 and AtRKD2 genes induces cell proliferation and the expression of an egg cell marker. Analyses of RKD-induced proliferating cells exhibit a shift of gene expression towards an egg cell-like transcriptome. Promoters of selected RKD-induced genes were shown to be predominantly active in the egg cell and can be activated by RKD in a transient protoplast expression assay. The data show that egg cell-specific RKD factors control a transcriptional program, which is characteristic for plant egg cells.

Identification and genetic analysis of the APOSPORY locus in Hypericum perforatum L.

The introduction of apomixis - seed formation without fertilization - into crop plants is a long-held goal of breeding research, since it would allow for the ready fixation of heterozygosity. The genetic basis of apomixis, whether of the aposporous or the diplosporous type, is still only poorly understood. Hypericum perforatum (St John's wort), a plant with a small genome and a short generation time, can be aposporous and/or parthenogenetic, and so represents an interesting model dicot for apomixis research. Here we describe a genetic analysis which first defined and then isolated a locus (designated HAPPY for Hypericum APOSPORY) associated with apospory. Amplified fragment length polymorphism (AFLP) profiling was used to generate a cleaved amplified polymorphic sequence (CAPS) marker for HAPPY which co-segregated with apospory but not with parthenogenesis, showing that these two components of apomixis are independently controlled. Apospory was inherited as a dominant simplex gene at the tetraploid level. Part of the HAPPY sequence is homologous to the Arabidopsis thaliana gene ARI7 encoding the ring finger protein ARIADNE7. This protein is predicted to be involved in various regulatory processes, including ubiquitin-mediated protein degradation. While the aposporous and sexual alleles of the HAPPY component HpARI were co-expressed in many parts of the plant, the gene product of the apomict's allele is truncated. Cloning HpARI represents the first step towards the full characterization of HAPPY and the elucidation of the molecular mechanisms underlying apomixis in H. perforatum.

No Evidence of Unexpected Transgenic Insertions in T1190 - A Transgenic Apple Used in Rapid Cycle Breeding - Following Whole Genome Sequencing.

Rapid cycle breeding uses transgenic early flowering plants as crossbreed parents to facilitate the shortening of breeding programs for perennial crops with long-lasting juvenility. Rapid cycle breeding in apple was established using the transgenic genotype T1190 expressing the BpMADS4 gene of silver birch. In this study, the genomes of T1190 and its non-transgenic wild-type PinS (F1-offspring of 'Pinova' and 'Idared') were sequenced by Illumina short-read sequencing in two separate experiments resulting in a mean sequencing depth of 182× for T1190 and 167× for PinS. The sequencing revealed 8,450 reads, which contain sequences of ≥20 bp identical to the plant transformation vector. These reads were assembled into 125 contigs, which were examined to see whether they contained transgenic insertions or if they are not using a five-step procedure. The sequence of one contig represents the known T-DNA insertion on chromosome 4 of T1190. The sequences of the remaining contigs were either equally present in T1190 and PinS, their part with sequence identity to the vector was equally present in apple reference genomes, or they seem to result from endophytic contaminations rather than from additional transgenic insertions. Therefore, we conclude that the transgenic apple plant T1190 contains only one transgenic insertion, located on chromosome 4, and shows no further partial insertions of the transformation vector. Accession Numbers: JQ974028.1.