A catalogue of recombination coldspots in interspecific tomato hybrids.

Increasing natural resistance and resilience in plants is key for ensuring food security within a changing climate. Breeders improve these traits by crossing cultivars with their wild relatives and introgressing specific alleles through meiotic recombination. However, some genomic regions are devoid of recombination especially in crosses between divergent genomes, limiting the combinations of desirable alleles. Here, we used pooled-pollen sequencing to build a map of recombinant and non-recombinant regions between tomato and five wild relatives commonly used for introgressive tomato breeding. We detected hybrid-specific recombination coldspots that underscore the role of structural variations in modifying recombination patterns and maintaining genetic linkage in interspecific crosses. Crossover regions and coldspots show strong association with specific TE superfamilies exhibiting differentially accessible chromatin between somatic and meiotic cells. About two-thirds of the genome are conserved coldspots, located mostly in the pericentromeres and enriched with retrotransposons. The coldspots also harbor genes associated with agronomic traits and stress resistance, revealing undesired consequences of linkage drag and possible barriers to breeding. We presented examples of linkage drag that can potentially be resolved by pairing tomato with other wild species. Overall, this catalogue will help breeders better understand crossover localization and make informed decisions on generating new tomato varieties.

Genome architecture and genetic diversity of allopolyploid okra (Abelmoschus esculentus).

The allopolyploid okra (Abelmoschus esculentus) unveiled telomeric repeats flanking distal gene-rich regions and short interstitial TTTAGGG telomeric repeats, possibly representing hallmarks of chromosomal speciation. Ribosomal RNA (rRNA) genes organize into 5S clusters, distinct from the 18S-5.8S-28S units, indicating an S-type rRNA gene arrangement. The assembly, in line with cytogenetic and cytometry observations, identifies 65 chromosomes and a 1.45 Gb genome size estimate in a haploid sibling. The lack of aberrant meiotic configurations implies limited to no recombination among sub-genomes. k-mer distribution analysis reveals 75% has a diploid nature and 15% heterozygosity. The configurations of Benchmarking Universal Single-Copy Ortholog (BUSCO), k-mer, and repeat clustering point to the presence of at least two sub-genomes one with 30 and the other with 35 chromosomes, indicating the allopolyploid nature of the okra genome. Over 130 000 putative genes, derived from mapped IsoSeq data and transcriptome data from public okra accessions, exhibit a low genetic diversity of one single nucleotide polymorphisms per 2.1 kbp. The genes are predominantly located at the distal chromosome ends, declining toward central scaffold domains. Long terminal repeat retrotransposons prevail in central domains, consistent with the observed pericentromeric heterochromatin and distal euchromatin. Disparities in paralogous gene counts suggest potential sub-genome differentiation implying possible sub-genome dominance. Amino acid query sequences of putative genes facilitated phenol biosynthesis pathway annotation. Comparison with manually curated reference KEGG pathways from related Malvaceae species reveals the genetic basis for putative enzyme coding genes that likely enable metabolic reactions involved in the biosynthesis of dietary and therapeutic compounds in okra.

Chasing breeding footprints through structural variations in Cucumis melo and wild relatives.

Cucumis melo (melon or muskmelon) is an important crop in the family of the Cucurbitaceae. Melon is cross pollinated and domesticated at several locations throughout the breeding history, resulting in highly diverse genetic structure in the germplasm. Yet, the relations among the groups and cultivars are still incomplete. We shed light on the melonbreeding history, analyzing structural variations ranging from 50 bp up to 100 kb, identified from whole genome sequences of 100 selected melon accessions and wild relatives. Phylogenetic trees based on SV types completely resolve cultivars and wild accessions into two monophyletic groups and clustering of cultivars largely correlates with their geographic origin. Taking into account morphology, we found six mis-categorized cultivars. Unique inversions are more often shared between cultivars, carrying advantageous genes and do not directly originate from wild species. Approximately 60% of the inversion breaks carry a long poly A/T motif, and following observations in other plant species, suggest that inversions in melon likely resulted from meiotic recombination events. We show that resistance genes in the linkage V region are expanded in the cultivar genomes compared to wild relatives. Furthermore, particular agronomic traits such as fruit ripening, fragrance, and stress response are specifically selected for in the melon subspecies. These results represent distinctive footprints of selective breeding that shaped today's melon. The sequences and genomic relations between land races, wild relatives, and cultivars will serve the community to identify genetic diversity, optimize experimental designs, and enhance crop development.

Meiotic recombination profiling of interspecific hybrid F1 tomato pollen by linked read sequencing.

Genome wide screening of pooled pollen samples from a single interspecific F1 hybrid obtained from a cross between tomato, Solanum lycopersicum and its wild relative, Solanum pimpinellifolium using linked read sequencing of the haploid nuclei, allowed profiling of the crossover (CO) and gene conversion (GC) landscape. We observed a striking overlap between cold regions of CO in the male gametes and our previously established F6 recombinant inbred lines (RILs) population. COs were overrepresented in non-coding regions in the gene promoter and 5'UTR regions of genes. Poly-A/T and AT rich motifs were found enriched in 1 kb promoter regions flanking the CO sites. Non-crossover associated allelic and ectopic GCs were detected in most chromosomes, confirming that besides CO, GC represents also a source for genetic diversity and genome plasticity in tomato. Furthermore, we identified processed break junctions pointing at the involvement of both homology directed and non-homology directed repair pathways, suggesting a recombination machinery in tomato that is more complex than currently anticipated.

Correcting palindromes in long reads after whole-genome amplification.

BACKGROUND: Next-generation sequencing requires sufficient DNA to be available. If limited, whole-genome amplification is applied to generate additional amounts of DNA. Such amplification often results in many chimeric DNA fragments, in particular artificial palindromic sequences, which limit the usefulness of long sequencing reads. RESULTS: Here, we present Pacasus, a tool for correcting such errors. Two datasets show that it markedly improves read mapping and de novo assembly, yielding results similar to these that would be obtained with non-amplified DNA. CONCLUSIONS: With Pacasus long-read technologies become available for sequencing targets with very small amounts of DNA, such as single cells or even single chromosomes.

A footprint of desiccation tolerance in the genome of Xerophyta viscosa.

Desiccation tolerance is common in seeds and various other organisms, but only a few angiosperm species possess vegetative desiccation tolerance. These 'resurrection species' may serve as ideal models for the ultimate design of crops with enhanced drought tolerance. To understand the molecular and genetic mechanisms enabling vegetative desiccation tolerance, we produced a high-quality whole-genome sequence for the resurrection plant Xerophyta viscosa and assessed transcriptome changes during its dehydration. Data revealed induction of transcripts typically associated with desiccation tolerance in seeds and involvement of orthologues of ABI3 and ABI5, both key regulators of seed maturation. Dehydration resulted in both increased, but predominantly reduced, transcript abundance of genomic 'clusters of desiccation-associated genes' (CoDAGs), reflecting the cessation of growth that allows for the expression of desiccation tolerance. Vegetative desiccation tolerance in X. viscosa was found to be uncoupled from drought-induced senescence. We provide strong support for the hypothesis that vegetative desiccation tolerance arose by redirection of genetic information from desiccation-tolerant seeds.

Botryococcus braunii strains compared for biomass productivity, hydrocarbon and carbohydrate content.

Botryococcus braunii can produce both long-chain hydrocarbons as well as carbohydrates in large quantities, and is therefore a promising industrial organism for the production of biopolymer building blocks. Many studies describe the use of different strains of Botryococcus braunii but differences in handling and cultivation conditions make the comparison between strains difficult. In this study, 16 B. braunii strains obtained from six culture collections were compared for their biomass productivity and hydrocarbon and carbohydrate content. Biomass productivity was highest for AC768 strain with 1.8gL-1day-1, while hydrocarbon production ranged from none to up to 42% per gram biomass dry weight, with Showa showing the highest hydrocarbon content followed by AC761. The total carbohydrate content varied from 20% to 76% per gram of the biomass dry weight, with CCALA777 as the highest producer. Glucose and galactose are the main monosaccharides in most strains and fucose content reached 463mgL-1 in CCALA778.

Expression of natural human ß1,4-GalT1 variants and of non-mammalian homologues in plants leads to differences in galactosylation of N-glycans.

ß1,4-Galactosylation of plant N-glycans is a prerequisite for commercial production of certain biopharmaceuticals in plants. Two different types of galactosylated N-glycans have initially been reported in plants as the result of expression of human ß1,4-galactosyltransferase 1 (GalT). Here we show that these differences are associated with differences at its N-terminus: the natural short variant of human GalT results in hybrid type N-glycans, whereas the long form generates bi-antennary complex type N-glycans. Furthermore, expression of non-mammalian, chicken and zebrafish GalT homologues with N-termini resembling the short human GalT N-terminus also induce hybrid type N-glycans. Providing both non-mammalian GalTs with a 13 amino acid N-terminal extension that distinguishes the two naturally occurring forms of human GalT, acted to increase the levels of bi-antennary galactosylated N-glycans when expressed in tobacco leaves. Replacement of the cytosolic tail and transmembrane domain of chicken and zebrafish GalTs with the corresponding region of rat α2,6-sialyltransferase yielded a gene whose expression enhanced the level of bi-antennary galactosylation even further.

Orientation of llama antibodies strongly increases sensitivity of biosensors.

Sensitivity of biosensors depends on the orientation of bio-receptors on the sensor surface. The objective of this study was to organize bio-receptors on surfaces in a way that their analyte binding site is exposed to the analyte solution. VHH proteins recognizing foot-and-mouth disease virus (FMDV) were used for making biosensors, and azides were introduced in the VHH to function as bioorthogonal reactive groups. The importance of the orientation of bio-receptors was addressed by comparing sensors with randomly oriented VHH (with multiple exposed azide groups) to sensors with uniformly oriented VHH (with only a single azide group). A surface plasmon resonance (SPR) chip exposing cyclooctyne was reacted to azide functionalized VHH domains, using click chemistry. Comparison between randomly and uniformly oriented bio-receptors showed up to 800-fold increase in biosensor sensitivity. This technique may increase the containment of infectious diseases such as FMDV as its strongly enhanced sensitivity may facilitate early diagnostics.

Identification of alg3 in the mushroom-forming fungus Schizophyllum commune and analysis of the Δalg3 knockout mutant.

Alg3 of Saccharomyces cerevisiae catalyzes the mannosyl transfer from Man-P-Dol to Man(5)GlcNAc(2)-PP-Dol resulting in the formation of Man(6)GlcNAc(2)-PP-Dol, which is then further processed to the final precursor oligosaccharide Glc(3)Man(9)GlcNAc(2) for N-glycosylation of proteins. Here, we identified the alg3 gene of the mushroom-forming fungus Schizophyllum commune by homology search. Its function was confirmed by the complementation of the Δalg3 strain of S. cerevisiae. Inactivation of alg3 in S. commune resulted in the production of predominantly Man(3)GlcNAc(2) protein-linked N-glycans. No impact on growth nor a developmental phenotype due to the deletion was observed. This provides a first step toward engineering of a homogeneous, humanized N-glycosylation pattern for the production of therapeutic glycoproteins in mushrooms.

Performance and long-term stability of the barley hordothionin gene in multiple transgenic apple lines.

Introduction of sustainable scab resistance in elite apple cultivars is of high importance for apple cultivation when aiming at reducing the use of chemical crop protectants. Genetic modification (GM) allows the rapid introduction of resistance genes directly into high quality apple cultivars. Resistance genes can be derived from apple itself but genetic modification also opens up the possibility to use other, non-host resistance genes. A prerequisite for application is the long-term performance and stability of the gene annex trait in the field. For this study, we produced and selected a series of transgenic apple lines of two cultivars, i.e. 'Elstar' and 'Gala' in which the barley hordothionin gene (hth) was introduced. After multiplication, the GM hth-lines, non-GM susceptible and resistant controls and GM non-hth controls were planted in a random block design in a field trial in 40 replicates. Scab resistance was monitored after artificial inoculation (first year) and after natural infection (subsequent years). After the trial period, the level of expression of the hth gene was checked by quantitative RT-PCR. Four of the six GM hth apple lines proved to be significantly less susceptible to apple scab and this trait was found to be stable for the entire 4-year period. Hth expression at the mRNA level was also stable.

Characterization of the single-chain Fv-Fc antibody MBP10 produced in Arabidopsis alg3 mutant seeds.

ER resident glycoproteins, including ectopically expressed recombinant glycoproteins, carry so-called high-mannose type N-glycans, which can be at different stages of processing. The presence of heterogeneous high-mannose type glycans on ER-retained therapeutic proteins is undesirable for specific therapeutic applications. Previously, we described an Arabidopsis alg3-2 glycosylation mutant in which aberrant Man(5)GlcNAc(2) mannose type N-glycans are transferred to proteins. Here we show that the alg3-2 mutation reduces the N-glycan heterogeneity on ER resident glycoproteins in seeds. We compared the properties of a scFv-Fc, with a KDEL ER retention tag (MBP10) that was expressed in seeds of wild type and alg3-2 plants. N-glycans on these antibodies from mutant seeds were predominantly of the intermediate Man(5)GlcNAc(2) compared to Man(8)GlcNAc(2) and Man(7)GlcNAc(2) isoforms on MBP10 from wild-type seeds. The presence of aberrant N-glycans on MBP10 did not seem to affect MBP10 dimerisation nor binding of MBP10 to its antigen. In alg3-2 the fraction of underglycosylated MBP10 protein forms was higher than in wild type. Interestingly, the expression of MBP10 resulted also in underglycosylation of other, endogenous glycoproteins.

Synthesis of Lewis X epitopes on plant N-glycans.

Glycoproteins from tobacco line xFxG1, in which expression of a hybrid beta-(1-->4)-galactosyltransferase (GalT) and a hybrid alpha-(1-->3)-fucosyltransferase IXa (FUT9a) is combined, contained an abundance of hybrid N-glycans with Lewis X (Le(X)) epitopes. A comparison with N-glycan profiles from plants expressing only the hybrid beta-(1-->4)-galactosyltransferase suggested that the fucosylation of the LacNAc residues in line xFxG1 protected galactosylated N-glycans from endogenous plant beta-galactosidase activity.

Solanum lycopersicum cv. Heinz 1706 chromosome 6: distribution and abundance of genes and retrotransposable elements.

We studied the physical and genetic organization of chromosome 6 of tomato (Solanum lycopersicum) cv. Heinz 1706 by combining bacterial artificial chromosome (BAC) sequence analysis, high-information-content fingerprinting, genetic analysis, and BAC-fluorescent in situ hybridization (FISH) mapping data. The chromosome positions of 81 anchored seed and extension BACs corresponded in most cases with the linear marker order on the high-density EXPEN 2000 linkage map. We assembled 25 BAC contigs and eight singleton BACs spanning 2.0 Mb of the short-arm euchromatin, 1.8 Mb of the pericentromeric heterochromatin and 6.9 Mb of the long-arm euchromatin. Sequence data were combined with their corresponding genetic and pachytene chromosome positions into an integrated map that covers approximately a third of the chromosome 6 euchromatin and a small part of the pericentromeric heterochromatin. We then compared physical length (Mb), genetic (cM) and chromosome distances (microm) for determining gap sizes between contigs, revealing relative hot and cold spots of recombination. Through sequence annotation we identified several clusters of functionally related genes and an uneven distribution of both gene and repeat sequences between heterochromatin and euchromatin domains. Although a greater number of the non-transposon genes were located in the euchromatin, the highly repetitive (22.4%) pericentromeric heterochromatin displayed an unexpectedly high gene content of one gene per 36.7 kb. Surprisingly, the short-arm euchromatin was relatively rich in repeats as well, with a repeat content of 13.4%, yet the ratio of Ty3/Gypsy and Ty1/Copia retrotransposable elements across the chromosome clearly distinguished euchromatin (2:3) from heterochromatin (3:2).

TOPAAS, a tomato and potato assembly assistance system for selection and finishing of bacterial artificial chromosomes.

We have developed the software package Tomato and Potato Assembly Assistance System (TOPAAS), which automates the assembly and scaffolding of contig sequences for low-coverage sequencing projects. The order of contigs predicted by TOPAAS is based on read pair information; alignments between genomic, expressed sequence tags, and bacterial artificial chromosome (BAC) end sequences; and annotated genes. The contig scaffold is used by TOPAAS for automated design of nonredundant sequence gap-flanking PCR primers. We show that TOPAAS builds reliable scaffolds for tomato (Solanum lycopersicum) and potato (Solanum tuberosum) BAC contigs that were assembled from shotgun sequences covering the target at 6- to 8-fold coverage. More than 90% of the gaps are closed by sequence PCR, based on the predicted ordering information. TOPAAS also assists the selection of large genomic insert clones from BAC libraries for walking. For this, tomato BACs are screened by automated BLAST analysis and in parallel, high-density nonselective amplified fragment length polymorphism fingerprinting is used for constructing a high-resolution BAC physical map. BLAST and amplified fragment length polymorphism analysis are then used together to determine the precise overlap. Assembly onto the seed BAC consensus confirms the BACs are properly selected for having an extremely short overlap and largest extending insert. This method will be particularly applicable where related or syntenic genomes are sequenced, as shown here for the Solanaceae, and potentially useful for the monocots Brassicaceae and Leguminosea.