Magdalena Pawełkowicz

The Genome Sequence of the North-European Cucumber (Cucumis sativus L.) Unravels Evolutionary Adaptation Mechanisms in Plants

Cucumber (Cucumis sativus L.), a widely cultivated crop, has originated from Eastern Himalayas and secondary domestication regions includes highly divergent climate conditions e.g. temperate and subtropical. We wanted to uncover adaptive genome differences between the cucumber cultivars and what sort of evolutionary molecular mechanisms regulate genetic adaptation of plants to different ecosystems and organism biodiversity. Here we present the draft genome sequence of the Cucumis sativus genome of the North-European Borszczagowski cultivar (line B10) and comparative genomics studies with the known genomes of: C. sativus (Chinese cultivar – Chinese Long (line 9930)), Arabidopsis thaliana, Populus trichocarpa and Oryza sativa. Cucumber genomes show extensive chromosomal rearrangements, distinct differences in quantity of the particular genes (e.g. involved in photosynthesis, respiration, sugar metabolism, chlorophyll degradation, regulation of gene expression, photooxidative stress tolerance, higher non-optimal temperatures tolerance and ammonium ion assimilation) as well as in distributions of abscisic acid-, dehydration- and ethylene-responsive cis-regulatory elements (CREs) in promoters of orthologous group of genes, which lead to the specific adaptation features. Abscisic acid treatment of non-acclimated Arabidopsis and C. sativus seedlings induced moderate freezing tolerance in Arabidopsis but not in C. sativus. This experiment together with analysis of abscisic acid-specific CRE distributions give a clue why C. sativus is much more susceptible to moderate freezing stresses than A. thaliana. Comparative analysis of all the five genomes showed that, each species and/or cultivars has a specific profile of CRE content in promoters of orthologous genes. Our results constitute the substantial and original resource for the basic and applied research on environmental adaptations of plants, which could facilitate creation of new crops with improved growth and yield in divergent conditions.

Next generation sequencing and omics in cucumber (Cucumis sativus L.) breeding directed research.

In the post-genomic era the availability of genomic tools and resources is leading us to novel generation methods in plant breeding, as they facilitate the study of the genotype and its relationship with the phenotype, in particular for complex traits. In this study we have mainly concentrated on the Cucumis sativus and (but much less) Cucurbitaceae family several important vegetable crops. There are many reports on research conducted in Cucurbitaceae plant breeding programs on the ripening process, phloem transport, disease resistance, cold tolerance and fruit quality traits. This paper presents the role played by new omic technologies in the creation of knowledge on the mechanisms of the formation of the breeding features. The analysis of NGS (NGS-next generation sequencing) data allows the discovery of new genes and regulatory sequences, their positions, and makes available large collections of molecular markers. Genome-wide expression studies provide breeders with an understanding of the molecular basis of complex traits. Firstly a high density map should be created for the reference genome, then each re-sequencing data could be mapped and new markers brought out into breeding populations. The paper also presents methods that could be used in the future for the creation of variability and genomic modification of the species in question. It has been shown also the state and usefulness in breeding the chloroplastomic and mitochondriomic study.

The construction and characteristics of a BAC library for Cucumis sativus L. 'B10'.

Cloning using bacterial artificial chromosomes (BACs) can yield high quality genomic libraries, which are used for the physical mapping, identification and isolation of genes, and for gene sequencing. A BAC genomic library was constructed from high molecular weight DNA (HMW DNA) obtained from nuclei of the cucumber (Cucumis sativus L. cv. Borszczagowski; B10 line). The DNA was digested with the HindIII restriction enzyme and ligated into the pCC1BAC vector. The library consists of 34,560 BAC clones with an average insert size of 135 kb, and 12.7x genome coverage. Screening the library for chloroplast and mitochondrial DNA content indicated an exceptionally low 0.26% contamination with chloroplast DNA and 0.3% with mitochondrial DNA.

Molecular Cytogenetic Analysis of Cucumis Wild Species Distributed in Southern Africa: Physical Mapping of 5S and 45S rDNA with DAPI

Wild Cucumis species have been divided into Australian/Asian and African groups using morphological and phylogenetic characteristics, and new species have been described recently. No molecular cytogenetic information is available for most of these species. The crossability between 5 southern African Cucumis species (C. africanus, C. anguria, C. myriocarpus, C. zeyheri, and C. heptadactylus) has been reported; however, the evolutionary relationship among them is still unclear. Here, a molecular cytogenetic analysis using FISH with 5S and 45S ribosomal DNA (rDNA) was used to investigate these Cucumis species based on sets of rDNA-bearing chromosomes (rch) types I, II and III. The molecular cytogenetic and phylogenetic results suggested that at least 2 steps of chromosomal rearrangements may have occurred during the evolution of tetraploid C. heptadactylus. In step 1, an additional 45S rDNA site was observed in the chromosome (type III). In particular, C. myriocarpus had a variety of rch sets. Our results suggest that chromosomal rearrangements may have occurred in the 45S rDNA sites. We propose that polyploid evolution occurred in step 2. This study provides insights into the chromosomal characteristics of African Cucumis species and contributes to the understanding of chromosomal evolution in this genus.

The genomic RNA1 and RNA2 sequences of the tobacco rattle virus isolates found in Polish potato fields.

Four tobacco rattle virus (TRV) isolates were identified from tobacco bait seedlings planted in soil samples from Polish potato fields. Sequence analysis of the genomic RNA1 of the isolates revealed significant similarity to the isolates Ho and AL recently found in Germany. Multiple sequence alignments of the genomic RNA2 indicated that the two isolates from northern Poland (Deb57 and Slu24) are in a cluster with the isolates PSG and PLB found in the Netherlands. The remaining two isolates, from central Poland (11r21 and Mlo7), are in a distinct group with the unique isolate SYM found in England. The RNA2 sequences of the studied isolates range from 1998 nt to 2739 nt in length, and all carry deletions of the 2b and/or 2c genes. The isolate Mlo7 has an atypical RNA2 structure, having its cp gene located in its central region.

Karyotype analysis and chromosomal distribution of repetitive DNA sequences of cucumis metuliferus using fluorescence in situ hybridization.

Cucumis metuliferus (2n = 24) is a cultivated species of the Cucumis genus which is a potential genetic resource for Cucumis crops. Although some cytogenetic research has been reported, there is no study of karyotyping in this species. Here, we used 4′,6-diamidino-2-phenylindole and chromomycin A3 staining to identify 12 pairs of chromosomes in early-metaphase cells. Fluorescence in situ hybridization revealed the chromosomal distribution patterns of the 5S and 45S ribosomal DNA (rDNA) genes, telomeres, and 3 different satellite repeats. The 2 major signals of the 45S rDNA were located on the satellite of chromosome 11, and the 2 signals of the 5S rDNA and 2 minor signals of the 45S rDNA were located on chromosome 12. The telomere probes hybridized to the ends of all chromosomes. The 3 satellite DNAs were localized at the ends of chromosomes 1, 2, 4-10, and at the end of the short arm of chromosome 3. In summary, we reported the identification of all chromosomes of C. metuliferus. We also depicted the location of 5S and 45S rDNA, the telomere motif sequence, CmetSat1, CmetSatT2, and CmetmSat1 in an ideogram.

DArT Markers Effectively Target Gene Space in the Rye Genome

Large genome size and complexity hamper considerably the genomics research in relevant species. Rye (Secale cereale L.) has one of the largest genomes among cereal crops and repetitive sequences account for over 90% of its length. Diversity Arrays Technology is a high-throughput genotyping method, in which a preferential sampling of gene-rich regions is achieved through the use of methylation sensitive restriction enzymes. We obtained sequences of 6,177 rye DArT markers and following a redundancy analysis assembled them into 3,737 non-redundant sequences, which were then used in homology searches against five Pooideae sequence sets. In total 515 DArT sequences could be incorporated into publicly available rye genome zippers providing a starting point for the integration of DArT- and transcript-based genomics resources in rye. Using Blast2Go pipeline we attributed putative gene functions to 1101 (29.4%) of the non-redundant DArT marker sequences, including 132 sequences with putative disease resistance-related functions, which were found to be preferentially located in the 4RL and 6RL chromosomes. Comparative analysis based on the DArT sequences revealed obvious inconsistencies between two recently published high density consensus maps of rye. Furthermore we demonstrated that DArT marker sequences can be a source of SSR polymorphisms. Obtained data demonstrate that DArT markers effectively target gene space in the large, complex, and repetitive rye genome. Through the annotation of putative gene functions and the alignment of DArT sequences relative to reference genomes we obtained information, that will complement the results of the studies, where DArT genotyping was deployed, by simplifying the gene ontology and microcolinearity based identification of candidate genes.

Genomes correction and assembling: present methods and tools

Recent rapid development of next generation sequencing (NGS) technologies provided significant impact into genomics field of study enabling implementation of many de novo sequencing projects of new species which was previously confined by technological costs. Along with advancement of NGS there was need for adjustment in assembly programs. New algorithms must cope with massive amounts of data computation in reasonable time limits and processing power and hardware is also an important factor. In this paper, we address the issue of assembly pipeline for de novo genome assembly provided by programs presently available for scientist both as commercial and as open – source software. The implementation of four different approaches – Greedy, Overlap – Layout – Consensus (OLC), De Bruijn and Integrated resulting in variation of performance is the main focus of our discussion with additional insight into issue of short and long reads correction.

Single-nucleotide polymorphisms and reading frame shifts in RNA2 recombinant regions of tobacco rattle virus isolates Slu24 and Deb57

Two previously sequenced tobacco rattle virus (TRV) isolates, Slu24 and Deb57, from Polish potato fields have recombinant RNA2 species. The 3’-proximal region of the Slu24 RNA2 is derived from the 3’ terminus of its supporting RNA1, while that of the Deb57 RNA2 is derived from the 3’ terminus of the unrelated RNA1 from the isolate SYM or PpK20. Gene structure annotation revealed open reading frames encoding truncated 16-kDa proteins in the recombinant regions of the RNA2 of Deb57 and Slu24. Reading frame shifts, single nucleotide substitutions and deletions occurred during recombination, including shifts from a stop codon or replacements of an internal stop codon. In the recombinant region of the Deb57 RNA2, the first reading frameshift event starts from the AUG start codon of the truncated 16-kDa protein. The second frameshift event, caused by a single nucleotide deletion upstream of the stop codon, leads to the splitting of the stop codon into two amino acid codons and the continuation of translation. In addition, a U-to-A substitution changes a potential internal stop codon UAA, which is caused by recombination-related frame shifts, into the codon AAA, encoding lysine. The replacement of this internal stop codon with an amino acid codon prevented the premature translation termination of the truncated 16-kDa protein. These recombination-related reading frame shifts are the driving force underlying the major differences in the translated amino acids, consequently leading to their translation into distinct polypeptides. Conversely, single nucleotide substitutions in the recombinant regions of the RNA2 of Deb57 or Slu24 resulted in only minor changes in the translated amino acids.

Comparison of bioinformatics programs for analysis of single nucleotide variants

Changes in genomic sequence might influence the gene expression, protein function and, what is related to phenotype of the organism. The Next Generation Sequencing provides a big amount of data that could be used in predicting the single nucleotide variants between analyzed and reference genome. Herein we compare three tools for predicting the structural variants: Freebayes, GATK toolkit and DeepVariant. Predictions with usage of each program were made on cucumber lines and the results were compared. Our analysis indicates that in order to obtain more precise and reliable variant predictions it is worth to use more than one program for detecting polymorphisms and cross-check the results.