Genetic Characterization of Cancer of Unknown Primary Using Liquid Biopsy Approaches.

Cancers of unknown primary (CUPs) comprise a heterogeneous group of rare metastatic tumors whose primary site cannot be identified after extensive clinical-pathological investigations. CUP patients are generally treated with empirical chemotherapy and have dismal prognosis. As recently reported, CUP genome presents potentially druggable alterations for which targeted therapies could be proposed. The paucity of tumor tissue, as well as the difficult DNA testing and the lack of dedicated panels for target gene sequencing are further relevant limitations. Here, we propose that circulating tumor cells (CTCs) and circulating tumor DNA (ctDNA) could be used to identify actionable mutations in CUP patients. Blood was longitudinally collected from two CUP patients. CTCs were isolated with CELLSEARCH® and DEPArrayTM NxT and Parsortix systems, immunophenotypically characterized and used for single-cell genomic characterization with Ampli1TM kits. Circulating cell-free DNA (ccfDNA), purified from plasma at different time points, was tested for tumor mutations with a CUP-dedicated, 92-gene custom panel using SureSelect Target Enrichment technology. In parallel, FFPE tumor tissue was analyzed with three different assays: FoundationOne CDx assay, DEPArray LibPrep and OncoSeek Panel, and the SureSelect custom panel. These approaches identified the same mutations, when the gene was covered by the panel, with the exception of an insertion in APC gene. which was detected by OncoSeek and SureSelect panels but not FoundationOne. FGFR2 and CCNE1 gene amplifications were detected in single CTCs, tumor tissue, and ccfDNAs in one patient. A somatic variant in ARID1A gene (p.R1276∗) was detected in the tumor tissue and ccfDNAs. The alterations were validated by Droplet Digital PCR in all ccfDNA samples collected during tumor evolution. CTCs from a second patient presented a pattern of recurrent amplifications in ASPM and SEPT9 genes and loss of FANCC. The 92-gene custom panel identified 16 non-synonymous somatic alterations in ccfDNA, including a deletion (I1485Rfs∗19) and a somatic mutation (p. A1487V) in ARID1A gene and a point mutation in FGFR2 gene (p.G384R). Our results support the feasibility of non-invasive liquid biopsy testing in CUP cases, either using ctDNA or CTCs, to identify CUP genetic alterations with broad NGS panels covering the most frequently mutated genes.

Whole-Transcriptome Analysis Unveils the Synchronized Activities of Genes for Fructans in Developing Tubers of the Jerusalem Artichoke.

Helianthus tuberosus L., known as the Jerusalem artichoke, is a hexaploid plant species, adapted to low-nutrient soils, that accumulates high levels of inulin in its tubers. Inulin is a fructose-based polysaccharide used either as dietary fiber or for the production of bioethanol. Key enzymes involved in inulin biosynthesis are well known. However, the gene networks underpinning tuber development and inulin accumulation in H. tuberous remain elusive. To fill this gap, we selected 6,365 expressed sequence tags (ESTs) from an H. tuberosus library to set up a microarray platform and record their expression across three tuber developmental stages, when rhizomes start enlarging (T0), at maximum tuber elongation rate (T3), and at tuber physiological maturity (Tm), in "VR" and "K8-HS142"clones. The former was selected as an early tuberizing and the latter as a late-tuberizing clone. We quantified inulin and starch levels, and qRT-PCR confirmed the expression of critical genes accounting for inulin biosynthesis. The microarray analysis revealed that the differences in morphological and physiological traits between tubers of the two clones are genetically determined since T0 and that is relatively low the number of differentially expressed ESTs across the stages shared between the clones (93). The expression of ESTs for sucrose:sucrose 1-fructosyltransferase (1-SST) and fructan:fructan 1-fructosyltransferase (1-FFT), the two critical genes for fructans polymerization, resulted to be temporarily synchronized and mirror the progress of inulin accumulation and stretching. The expression of ESTs for starch biosynthesis was insignificant throughout the developmental stages of the clones in line with the negligible level of starch into their mature tubers, where inulin was the dominant polysaccharide. Overall, our study disclosed candidate genes underpinning the development and storage of carbohydrates in the tubers of two H. tuberosus clones. A model according to which the steady-state levels of 1-SST and 1-FFT transcripts are developmentally controlled and might represent a limiting factor for inulin accumulation has been provided. Our finding may have significant repercussions for breeding clones with improved levels of inulin for food and chemical industry.

A chromosome-anchored eggplant genome sequence reveals key events in Solanaceae evolution.

With approximately 450 species, spiny Solanum species constitute the largest monophyletic group in the Solanaceae family, but a high-quality genome assembly from this group is presently missing. We obtained a chromosome-anchored genome assembly of eggplant (Solanum melongena), containing 34,916 genes, confirming that the diploid gene number in the Solanaceae is around 35,000. Comparative genomic studies with tomato (S. lycopersicum), potato (S. tuberosum) and pepper (Capsicum annuum) highlighted the rapid evolution of miRNA:mRNA regulatory pairs and R-type defense genes in the Solanaceae, and provided a genomic basis for the lack of steroidal glycoalkaloid compounds in the Capsicum genus. Using parsimony methods, we reconstructed the putative chromosomal complements of the key founders of the main Solanaceae clades and the rearrangements that led to the karyotypes of extant species and their ancestors. From 10% to 15% of the genes present in the four genomes were syntenic paralogs (ohnologs) generated by the pre-γ, γ and T paleopolyploidy events, and were enriched in transcription factors. Our data suggest that the basic gene network controlling fruit ripening is conserved in different Solanaceae clades, and that climacteric fruit ripening involves a differential regulation of relatively few components of this network, including CNR and ethylene biosynthetic genes.

A streamlined workflow for single-cells genome-wide copy-number profiling by low-pass sequencing of LM-PCR whole-genome amplification products.

Chromosomal instability and associated chromosomal aberrations are hallmarks of cancer and play a critical role in disease progression and development of resistance to drugs. Single-cell genome analysis has gained interest in latest years as a source of biomarkers for targeted-therapy selection and drug resistance, and several methods have been developed to amplify the genomic DNA and to produce libraries suitable for Whole Genome Sequencing (WGS). However, most protocols require several enzymatic and cleanup steps, thus increasing the complexity and length of protocols, while robustness and speed are key factors for clinical applications. To tackle this issue, we developed a single-tube, single-step, streamlined protocol, exploiting ligation mediated PCR (LM-PCR) Whole Genome Amplification (WGA) method, for low-pass genome sequencing with the Ion Torrent™ platform and copy number alterations (CNAs) calling from single cells. The method was evaluated on single cells isolated from 6 aberrant cell lines of the NCI-H series. In addition, to demonstrate the feasibility of the workflow on clinical samples, we analyzed single circulating tumor cells (CTCs) and white blood cells (WBCs) isolated from the blood of patients affected by prostate cancer or lung adenocarcinoma. The results obtained show that the developed workflow generates data accurately representing whole genome absolute copy number profiles of single cell and allows alterations calling at resolutions down to 100 Kbp with as few as 200,000 reads. The presented data demonstrate the feasibility of the Ampli1™ WGA-based low-pass workflow for detection of CNAs in single tumor cells which would be of particular interest for genome-driven targeted therapy selection and for monitoring of disease progression.

Whole-genome sequencing and SNV genotyping of 'Nebbiolo' (Vitis vinifera L.) clones.

'Nebbiolo' (Vitis vinifera) is among the most ancient and prestigious wine grape varieties characterised by a wide genetic variability exhibited by a high number of clones (vegetatively propagated lines of selected mother plants). However, limited information is available for this cultivar at the molecular and genomic levels. The whole-genomes of three 'Nebbiolo' clones (CVT 71, CVT 185 and CVT 423) were re-sequenced and a de novo transcriptome assembly was produced. Important remarks about the genetic peculiarities of 'Nebbiolo' and its intra-varietal variability useful for clonal identification were reported. In particular, several varietal transcripts identified for the first time in 'Nebbiolo' were disease resistance genes and single-nucleotide variants (SNVs) identified in 'Nebbiolo', but not in other cultivars, were associated with genes involved in the stress response. Ten newly discovered SNVs were successfully employed to identify some periclinal chimeras and to classify 98 'Nebbiolo' clones in seven main genotypes, which resulted to be linked to the geographical origin of accessions. In addition, for the first time it was possible to discriminate some 'Nebbiolo' clones from the others.

Interference with ethylene perception at receptor level sheds light on auxin and transcriptional circuits associated with the climacteric ripening of apple fruit (Malus x domestica Borkh.).

Apple (Malus x domestica Borkh.) is a model species for studying the metabolic changes that occur at the onset of ripening in fruit crops, and the physiological mechanisms that are governed by the hormone ethylene. In this study, to dissect the climacteric interplay in apple, a multidisciplinary approach was employed. To this end, a comprehensive analysis of gene expression together with the investigation of several physiological entities (texture, volatilome and content of polyphenolic compounds) was performed throughout fruit development and ripening. The transcriptomic profiling was conducted with two microarray platforms: a dedicated custom array (iRIPE) and a whole genome array specifically enriched with ripening-related genes for apple (WGAA). The transcriptomic and phenotypic changes following the application of 1-methylcyclopropene (1-MCP), an ethylene inhibitor leading to important modifications in overall fruit physiology, were also highlighted. The integrative comparative network analysis showed both negative and positive correlations between ripening-related transcripts and the accumulation of specific metabolites or texture components. The ripening distortion caused by the inhibition of ethylene perception, in addition to affecting the ethylene pathway, stimulated the de-repression of auxin-related genes, transcription factors and photosynthetic genes. Overall, the comprehensive repertoire of results obtained here advances the elucidation of the multi-layered climacteric mechanism of fruit ripening, thus suggesting a possible transcriptional circuit governed by hormones and transcription factors.

DNA methylation and gene expression profiles show novel regulatory pathways in hepatocellular carcinoma.

BACKGROUND: Alcohol is a well-known risk factor for hepatocellular carcinoma (HCC), but the mechanisms underlying the alcohol-related hepatocarcinogenesis are still poorly understood. Alcohol alters the provision of methyl groups within the hepatic one-carbon metabolism, possibly inducing aberrant DNA methylation. Whether specific pathways are epigenetically regulated in alcohol-associated HCC is, however, unknown. The aim of the present study was to investigate the genome-wide promoter DNA methylation and gene expression profiles in non-viral, alcohol-associated HCC. From eight HCC patients undergoing curative surgery, array-based DNA methylation and gene expression data of all annotated genes were analyzed by comparing HCC tissue and homologous cancer-free liver tissue. RESULTS: After merging the DNA methylation with gene expression data, we identified 159 hypermethylated-repressed, 30 hypomethylated-induced, 49 hypermethylated-induced, and 56 hypomethylated-repressed genes. Notably, promoter DNA methylation emerged as a novel regulatory mechanism for the transcriptional repression of genes controlling the retinol metabolism (ADH1A, ADH1B, ADH6, CYP3A43, CYP4A22, RDH16), iron homeostasis (HAMP), one-carbon metabolism (SHMT1), and genes with a putative, newly identified function as tumor suppressors (FAM107A, IGFALS, MT1G, MT1H, RNF180). CONCLUSIONS: A genome-wide DNA methylation approach merged with array-based gene expression profiles allowed identifying a number of novel, epigenetically regulated candidate tumor-suppressor genes in alcohol-associated hepatocarcinogenesis. Retinol metabolism genes and SHMT1 are also epigenetically regulated through promoter DNA methylation in alcohol-associated HCC. Due to the reversibility of epigenetic mechanisms by environmental/nutritional factors, these findings may open up to novel interventional strategies for hepatocarcinogenesis prevention in HCC related to alcohol, a modifiable dietary component.

Physical Mapping of Bread Wheat Chromosome 5A: An Integrated Approach.

The huge size, redundancy, and highly repetitive nature of the bread wheat [Triticum aestivum (L.)] genome, makes it among the most difficult species to be sequenced. To overcome these limitations, a strategy based on the separation of individual chromosomes or chromosome arms and the subsequent production of physical maps was established within the frame of the International Wheat Genome Sequence Consortium (IWGSC). A total of 95,812 bacterial artificial chromosome (BAC) clones of short-arm chromosome 5A (5AS) and long-arm chromosome 5A (5AL) arm-specific BAC libraries were fingerprinted and assembled into contigs by complementary analytical approaches based on the FingerPrinted Contig (FPC) and Linear Topological Contig (LTC) tools. Combined anchoring approaches based on polymerase chain reaction (PCR) marker screening, microarray, and sequence homology searches applied to several genomic tools (i.e., genetic maps, deletion bin map, neighbor maps, BAC end sequences (BESs), genome zipper, and chromosome survey sequences) allowed the development of a high-quality physical map with an anchored physical coverage of 75% for 5AS and 53% for 5AL with high portions (64 and 48%, respectively) of contigs ordered along the chromosome. In the genome of grasses, Brachypodium [Brachypodium distachyon (L.) Beauv.], rice (Oryza sativa L.), and sorghum [Sorghum bicolor (L.) Moench] homologs of genes on wheat chromosome 5A were separated into syntenic blocks on different chromosomes as a result of translocations and inversions during evolution. The physical map presented represents an essential resource for fine genetic mapping and map-based cloning of agronomically relevant traits and a reference for the 5A sequencing projects.

Transcriptional Reprogramming of the Mycoparasitic Fungus Ampelomyces quisqualis During the Powdery Mildew Host-Induced Germination.

Ampelomyces quisqualis is a mycoparasite of a diverse range of phytopathogenic fungi associated with the powdery mildew disease. Among them are several Erysiphaceae species with great economic impact on high-value crops such as grape. Due to its ability to parasitize and prevent the spread of powdery mildews, A. quisqualis has received considerable attention for its biocontrol potential. However, and in sharp contrast to the extensively studied biocontrol species belonging to the genus Trichoderma, little is known about the biology of A. quisqualis at the molecular and genetic levels. We present the first genome-wide transcription profiling in A. quisqualis during host-induced germination. A total of 1,536 putative genes showed significant changes in transcription during the germination of A. quisqualis. This finding denotes an extensive transcriptional reprogramming of A. quisqualis induced by the presence of the host. Several upregulated genes were predicted to encode for putative mycoparasitism-related proteins such as secreted proteases, virulence factors, and proteins related to toxin biosynthesis. Our data provide the most comprehensive sequence resource currently available for A. quisqualis in addition to offering valuable insights into the biology of A. quisqualis and its mycoparasitic lifestyle. Eventually, this may improve the biocontrol capacity of this mycoparasite.

Nitrate induction triggers different transcriptional changes in a high and a low nitrogen use efficiency maize inbred line.

In higher plants, NO3(-) can induce its own uptake and the magnitude of this induction is positively related to the external anion concentration. This phenomenon has been characterized in both herbaceous and woody plants. Here, different adaptation strategies of roots from two maize (Zea mays L., ZmAGOs) inbred lines differing in nitrogen use efficiency (NUE) and exhibiting different timing of induction were discussed by investigating NO3(-) -induced changes in their transcriptome. Lo5 line (high NUE) showing the maximum rate of NO3(-) uptake 4 h after the provision of 200 µmol/L NO3(-) treatment modulated a higher number of transcripts relative to T250 (low NUE) that peaked after 12 h. The two inbred lines share only 368 transcripts that are modulated by the treatment with NO3(-) and behaved differently when transcripts involved in anion uptake and assimilation were analyzed. T250 line responded to the NO3(-) induction modulating this group of genes as reported for several plant species. On the contrary, the Lo5 line did not exhibit during the induction changes in this set of genes. Obtained data suggest the importance of exploring the physiological and molecular variations among different maize genotypes in response to environmental clues like NO3(-) provision, in order to understand mechanisms underlying NUE.

De novo genome assembly of the soil-borne fungus and tomato pathogen Pyrenochaeta lycopersici.

BACKGROUND: Pyrenochaeta lycopersici is a soil-dwelling ascomycete pathogen that causes corky root rot disease in tomato (Solanum lycopersicum) and other Solanaceous crops, reducing fruit yields by up to 75%. Fungal pathogens that infect roots receive less attention than those infecting the aerial parts of crops despite their significant impact on plant growth and fruit production. RESULTS: We assembled a 54.9Mb P. lycopersici draft genome sequence based on Illumina short reads, and annotated approximately 17,000 genes. The P. lycopersici genome is closely related to hemibiotrophs and necrotrophs, in agreement with the phenotypic characteristics of the fungus and its lifestyle. Several gene families related to host-pathogen interactions are strongly represented, including those responsible for nutrient absorption, the detoxification of fungicides and plant cell wall degradation, the latter confirming that much of the genome is devoted to the pathogenic activity of the fungus. We did not find a MAT gene, which is consistent with the classification of P. lycopersici as an imperfect fungus, but we observed a significant expansion of the gene families associated with heterokaryon incompatibility (HI). CONCLUSIONS: The P. lycopersici draft genome sequence provided insight into the molecular and genetic basis of the fungal lifestyle, characterizing previously unknown pathogenic behaviors and defining strategies that allow this asexual fungus to increase genetic diversity and to acquire new pathogenic traits.

Extreme sensitivity of gene expression in human SH-SY5Y neurocytes to ultra-low doses of Gelsemium sempervirens.

BACKGROUND: Gelsemium sempervirens L. (Gelsemium s.) is a traditional medicinal plant, employed as an anxiolytic at ultra-low doses and animal models recently confirmed this activity. However the mechanisms by which it might operate on the nervous system are largely unknown. This work investigates the gene expression of a human neurocyte cell line treated with increasing dilutions of Gelsemium s. extract. METHODS: Starting from the crude extract, six 100 × (centesimal, c) dilutions of Gelsemium s. (2c, 3c, 4c, 5c, 9c and 30c) were prepared according to the French homeopathic pharmacopoeia. Human SH-SY5Y neuroblastoma cells were exposed for 24 h to test dilutions, and their transcriptome compared by microarray to that of cells treated with control vehicle solutions. RESULTS: Exposure to the Gelsemium s. 2c dilution (the highest dose employed, corresponding to a gelsemine concentration of 6.5 × 10(-9) M) significantly changed the expression of 56 genes, of which 49 were down-regulated and 7 were overexpressed. Several of the down-regulated genes belonged to G-protein coupled receptor signaling pathways, calcium homeostasis, inflammatory response and neuropeptide receptors. Fisher exact test, applied to the group of 49 genes down-regulated by Gelsemium s. 2c, showed that the direction of effects was significantly maintained across the treatment with high homeopathic dilutions, even though the size of the differences was distributed in a small range. CONCLUSIONS: The study shows that Gelsemium s., a medicinal plant used in traditional remedies and homeopathy, modulates a series of genes involved in neuronal function. A small, but statistically significant, response was detected even to very low doses/high dilutions (up to 30c), indicating that the human neurocyte genome is extremely sensitive to this regulation.

Patchwork sequencing of tomato San Marzano and Vesuviano varieties highlights genome-wide variations.

BACKGROUND: Investigation of tomato genetic resources is a crucial issue for better straight evolution and genetic studies as well as tomato breeding strategies. Traditional Vesuviano and San Marzano varieties grown in Campania region (Southern Italy) are famous for their remarkable fruit quality. Owing to their economic and social importance is crucial to understand the genetic basis of their unique traits. RESULTS: Here, we present the draft genome sequences of tomato Vesuviano and San Marzano genome. A 40x genome coverage was obtained from a hybrid Illumina paired-end reads assembling that combines de novo assembly with iterative mapping to the reference S. lycopersicum genome (SL2.40). Insertions, deletions and SNP variants were carefully measured. When assessed on the basis of the reference annotation, 30% of protein-coding genes are predicted to have variants in both varieties. Copy genes number and gene location were assessed by mRNA transcripts mapping, showing a closer relationship of San Marzano with reference genome. Distinctive variations in key genes and transcription/regulation factors related to fruit quality have been revealed for both cultivars. CONCLUSIONS: The effort performed highlighted varieties relationships and important variants in fruit key processes useful to dissect the path from sequence variant to phenotype.

The high polyphenol content of grapevine cultivar tannat berries is conferred primarily by genes that are not shared with the reference genome.

The grapevine (Vitis vinifera) cultivar Tannat is cultivated mainly in Uruguay for the production of high-quality red wines. Tannat berries have unusually high levels of polyphenolic compounds, producing wines with an intense purple color and remarkable antioxidant properties. We investigated the genetic basis of these important characteristics by sequencing the genome of the Uruguayan Tannat clone UY11 using Illumina technology, followed by a mixture of de novo assembly and iterative mapping onto the PN40024 reference genome. RNA sequencing data for genome reannotation were processed using a combination of reference-guided annotation and de novo transcript assembly, allowing 5901 previously unannotated or unassembled genes to be defined and resulting in the discovery of 1873 genes that were not shared with PN40024. Expression analysis showed that these cultivar-specific genes contributed substantially (up to 81.24%) to the overall expression of enzymes involved in the synthesis of phenolic and polyphenolic compounds that contribute to the unique characteristics of the Tannat berries. The characterization of the Tannat genome therefore indicated that the grapevine reference genome lacks many genes that appear to be relevant for the varietal phenotype.

De novo transcriptome characterization of Vitis vinifera cv. Corvina unveils varietal diversity.

BACKGROUND: Plants such as grapevine (Vitis spp.) display significant inter-cultivar genetic and phenotypic variation. The genetic components underlying phenotypic diversity in grapevine must be understood in order to disentangle genetic and environmental factors. RESULTS: We have shown that cDNA sequencing by RNA-seq is a robust approach for the characterization of varietal diversity between a local grapevine cultivar (Corvina) and the PN40024 reference genome. We detected 15,161 known genes including 9463 with novel splice isoforms, and identified 2321 potentially novel protein-coding genes in non-annotated or unassembled regions of the reference genome. We also discovered 180 apparent private genes in the Corvina genome which were missing from the reference genome. CONCLUSIONS: The de novo assembly approach allowed a substantial amount of the Corvina transcriptome to be reconstructed, improving known gene annotations by robustly defining gene structures, annotating splice isoforms and detecting genes without annotations. The private genes we discovered are likely to be nonessential but could influence certain cultivar-specific characteristics. Therefore, the application of de novo transcriptome assembly should not be restricted to species lacking a reference genome because it can also improve existing reference genome annotations and identify novel, cultivar-specific genes.

Identification of putative stage-specific grapevine berry biomarkers and omics data integration into networks.

The analysis of grapevine (Vitis vinifera) berries at the transcriptomic, proteomic, and metabolomic levels can provide great insight into the molecular events underlying berry development and postharvest drying (withering). However, the large and very different data sets produced by such investigations are difficult to integrate. Here, we report the identification of putative stage-specific biomarkers for berry development and withering and, to our knowledge, the first integrated systems-level study of these processes. Transcriptomic, proteomic, and metabolomic data were integrated using two different strategies, one hypothesis free and the other hypothesis driven. A multistep hypothesis-free approach was applied to data from four developmental stages and three withering intervals, with integration achieved using a hierarchical clustering strategy based on the multivariate bidirectional orthogonal projections to latent structures technique. This identified stage-specific functional networks of linked transcripts, proteins, and metabolites, providing important insights into the key molecular processes that determine the quality characteristics of wine. The hypothesis-driven approach was used to integrate data from three withering intervals, starting with subdata sets of transcripts, proteins, and metabolites. We identified transcripts and proteins that were modulated during withering as well as specific classes of metabolites that accumulated at the same time and used these to select subdata sets of variables. The multivariate bidirectional orthogonal projections to latent structures technique was then used to integrate the subdata sets, identifying variables representing selected molecular processes that take place specifically during berry withering. The impact of this holistic approach on our knowledge of grapevine berry development and withering is discussed.

GRAMINIFOLIA homolog expression in Streptocarpus rexii is associated with the basal meristems in phyllomorphs, a morphological novelty in Gesneriaceae.

Streptocarpus is a genus showing great variation in vegetative plant architecture and hence provides an attractive system to study the evolution of morphological diversity. Besides species showing an orthodox caulescent plant organization, producing leaves from a conventional shoot apical meristem (SAM), there are species whose body plan is composed of units (phyllomorphs) consisting of a petiole-like structure and a lamina that has the ability of continued growth. The first of these units is the macrocotyledon, derived from the continued growth of one of the two cotyledons by the activity of a basal meristem (BM), whereas further phyllomorphs develop from a SAM-like meristem. We carried out anatomical and morphological studies on the macrocotyledon of Streptocarpus rexii showing that the lamina has a bifacial structure, whereas the petiolode is partially unifacial. YABBY transcription factors are known to be involved in organ polarity and also promote lamina growth. We characterized the expression of SrGRAM, an ortholog of the YABBY genes GRAMINIFOLIA (GRAM) and FILAMENTOUS FLOWER (FIL), in S. rexii by in situ hybridization and RT-PCR. Gene expression pattern during embryogenesis was found to be conserved between SrGRAM and FIL from Arabidopsis. During subsequent seedling development SrGRAM expression in S. rexii was closely associated with the activity of the BM of the macrocotyledon and consecutively produced phyllomorphs, whereas it was excluded from the SAM-like meristem. Our results suggest that SrGRAM acts in intercalary growth and that an altered regulation of SrGRAM may underlay the evolution of the BM in S. rexii.

WUS and STM homologs are linked to the expression of lateral dominance in the acaulescent Streptocarpus rexii (Gesneriaceae).

Acaulescent species of Streptocarpus Lindl. show unusual patterns of growth, characterized by anisocotyly (i.e. the unequal growth of cotyledons after germination) and lack of a conventional embryonic shoot apical meristem (SAM). A SAM-like structure appears during post-embryonic development on the axis of the continuously growing cotyledon. Since we have shown previously that KNOX genes are involved in this unusual morphology of Streptocarpus rexii, here we investigated the expression pattern of WUSCHEL (WUS), which is also required for the indeterminacy of the SAM, but is expressed independently from KNOX in Arabidopsis thaliana. In A. thaliana WUSCHEL is involved in the maintenance of the stem cell fate in the organizing centre. The expression pattern of the WUS ortholog in S. rexii (SrWUS) strongly deviates from that of the model plant, suggesting a fundamentally different spatial and temporal regulation of signalling involved in meristem initiation and maintenance. In S. rexii, exogenous application of growth regulators, i.e. gibberellin (GA(3)), cytokinin (CK) and a gibberellin biosynthesis inhibitor (PAC), prevents anisocotyly and relocates meristematic cells to a position of conventional SAMs; this coincides with a re-localization of the two main pathways controlling meristem formation, the SrWUS and the KNOX pathways. Our results suggest that the establishment of a hormone imbalance in the seedlings is the basis of anisocotyly, causing a lateral dominance of the macrocotyledon over the microcotyledon. The peculiar morphogenetic program in S. rexii is linked to this delicate hormone balance and is the result of crosstalk between endogenous hormones and regulatory genes.