Whole genome scanning of a Mediterranean basin hotspot collection provides new insights into olive tree biodiversity and biology.

Olive tree (Olea europaea L. subsp. europaea var. europaea) is one of the most important species of the Mediterranean region and one of the most ancient species domesticated. The availability of whole genome assemblies and annotations of olive tree cultivars and oleaster (O. europaea subsp. europaea var. sylvestris) has contributed to a better understanding of genetic and genomic differences between olive tree cultivars. However, compared to other plant species there is still a lack of genomic resources for olive tree populations that span the entire Mediterranean region. In the present study we developed the most complete genomic variation map and the most comprehensive catalog/resource of molecular variation to date for 89 olive tree genotypes originating from the entire Mediterranean basin, revealing the genetic diversity of this commercially significant crop tree and explaining the divergence/similarity among different variants. Additionally, the monumental ancient tree 'Throuba Naxos' was studied to characterize the potential origin or routes of olive tree domestication. Several candidate genes known to be associated with key agronomic traits, including olive oil quality and fruit yield, were uncovered by a selective sweep scan to be under selection pressure on all olive tree chromosomes. To further exploit the genomic and phenotypic resources obtained from the current work, genome-wide association analyses were performed for 23 morphological and two agronomic traits. Significant associations were detected for eight traits that provide valuable candidates for fruit tree breeding and for deeper understanding of olive tree biology.

Copy Number Variation on ABCC2-DNMBP Loci Affects the Diversity and Composition of the Fecal Microbiota in Pigs.

Genetic variation in the pig genome partially modulates the composition of porcine gut microbial communities. Previous studies have been focused on the association between single nucleotide polymorphisms (SNPs) and the gut microbiota, but little is known about the relationship between structural variants and fecal microbial traits. The main goal of this study was to explore the association between porcine genome copy number variants (CNVs) and the diversity and composition of pig fecal microbiota. For this purpose, we used whole-genome sequencing data to undertake a comprehensive identification of CNVs followed by a genome-wide association analysis between the estimated CNV status and the fecal bacterial diversity in a commercial Duroc pig population. A CNV predicted as gain (DUP) partially harboring ABCC2-DNMBP loci was associated with richness (P = 5.41 × 10-5, false discovery rate [FDR] = 0.022) and Shannon α-diversity (P = 1.42 × 10-4, FDR = 0.057). The in silico predicted gain of copies was validated by real-time quantitative PCR (qPCR), and its segregation, and positive association with the richness and Shannon α-diversity of the porcine fecal bacterial ecosystem was confirmed in an unrelated F1 (Duroc × Iberian) cross. Our results advise the relevance of considering the role of host-genome structural variants as potential modulators of microbial ecosystems and suggest the ABCC2-DNMBP CNV as a host-genetic factor for the modulation of the diversity and composition of the fecal microbiota in pigs. IMPORTANCE A better understanding of the environmental and host factors modulating gut microbiomes is a topic of greatest interest. Recent evidence suggests that genetic variation in the pig genome partially controls the composition of porcine gut microbiota. However, since previous studies have been focused on the association between single nucleotide polymorphisms and the fecal microbiota, little is known about the relationship between other sources of genetic variation, like the structural variants and microbial traits. Here, we identified, experimentally validated, and replicated in an independent population a positive link between the gain of copies of ABCC2-DNMBP loci and the diversity and composition of pig fecal microbiota. Our results advise the relevance of considering the role of host-genome structural variants as putative modulators of microbial ecosystems and open the possibility of implementing novel holobiont-based management strategies in breeding programs for the simultaneous improvement of microbial traits and host performance.

Development and Evaluation of an AxiomTM 60K SNP Array for Almond (Prunus dulcis).

A high-density single nucleotide polymorphism (SNP) array is essential to enable faster progress in plant breeding for new cultivar development. In this regard, we have developed an Axiom 60K almond SNP array by resequencing 81 almond accessions. For the validation of the array, a set of 210 accessions were genotyped and 82.8% of the SNPs were classified in the best recommended SNPs. The rate of missing data was between 0.4% and 2.7% for the almond accessions and less than 15.5% for the few peach and wild accessions, suggesting that this array can be used for peach and interspecific peach × almond genetic studies. The values of the two SNPs linked to the RMja (nematode resistance) and SK (bitterness) genes were consistent. We also genotyped 49 hybrids from an almond F2 progeny and could build a genetic map with a set of 1159 SNPs. Error rates, less than 1%, were evaluated by comparing replicates and by detection of departures from Mendelian inheritance in the F2 progeny. This almond array is commercially available and should be a cost-effective genotyping tool useful in the search for new genes and quantitative traits loci (QTL) involved in the control of agronomic traits.

A cryptic variation in a member of the Ovate Family Proteins is underlying the melon fruit shape QTL fsqs8.1.

KEY MESSAGE: The gene underlying the melon fruit shape QTL fsqs8.1 is a member of the Ovate Family Proteins. Variation in fruit morphology is caused by changes in gene expression likely due to a cryptic structural variation in this locus. Melon cultivars have a wide range of fruit morphologies. Quantitative trait loci (QTL) have been identified underlying such diversity. This research focuses on the fruit shape QTL fsqs8.1, previously detected in a cross between the accession PI 124112 (CALC, producing elongated fruit) and the cultivar 'Piel de Sapo' (PS, producing oval fruit). The CALC fsqs8.1 allele induced round fruit shape, being responsible for the transgressive segregation for this trait observed in that population. In fact, the introgression line CALC8-1, carrying the fsqs8.1 locus from CALC into the PS genetic background, produced perfect round fruit. Following a map-based cloning approach, we found that the gene underlying fsqs8.1 is a member of the Ovate Family Proteins (OFP), CmOFP13, likely a homologue of AtOFP1 and SlOFP20 from Arabidopsis thaliana and tomato, respectively. The induction of the round shape was due to the higher expression of the CALC allele at the early ovary development stage. The fsqs8.1 locus showed an important structural variation, being CmOFP13 surrounded by two deletions in the CALC genome. The deletions are present at very low frequency in melon germplasm. Deletions and single nucleotide polymorphisms in the fsqs8.1 locus could not be not associated with variation in fruit shape among different melon accessions, what indicates that other genetic factors should be involved to induce the CALC fsqs8.1 allele effects. Therefore, fsqs8.1 is an example of a cryptic variation that alters gene expression, likely due to structural variation, resulting in phenotypic changes in melon fruit morphology.

Resynthesis: Marker-Based Partial Reconstruction of Elite Genotypes in Clonally-Reproducing Plant Species.

We propose a method for marker-based selection of cultivars of clonally-reproducing plant species which keeps the basic genetic architecture of a top-performing cultivar (usually a partly heterozygous genotype), with the addition of some agronomically relevant differences (such as production time, product appearance or quality), providing added value to the product or cultivation process. The method is based on selecting a) two complementary nearly-inbred lines from successive selfing generations (ideally only F2 and F3) of large size, that may generate individuals with most of their genome identical to the original cultivar but being homozygous for either of the two component haplotypes in the rest, and b) individuals with such characteristics already occurring in the F2. Option a) allows for introgressing genes from other individuals in one or both of these nearly-inbred lines. Peach, a woody-perennial, clonally-reproduced species, was chosen as a model for a proof of concept of the Resynthesis process due to its biological characteristics: self-compatibility, compact and genetically well-known genome, low recombination rates and relatively short intergeneration time (3-4 years). From 416 F2 seedlings from cultivar Sweet Dream (SD), we obtained seven individuals with 76-94% identity with SD, and selected five pairs of complementary lines with average homozygosity of the two parents ≥0.70 such that crossing would produce some individuals highly similar to SD. The application of this scheme to other species with more complex genomes or biological features, including its generalization to F1 hybrids, is discussed.

Transposons played a major role in the diversification between the closely related almond and peach genomes: results from the almond genome sequence.

We sequenced the genome of the highly heterozygous almond Prunus dulcis cv. Texas combining short- and long-read sequencing. We obtained a genome assembly totaling 227.6 Mb of the estimated almond genome size of 238 Mb, of which 91% is anchored to eight pseudomolecules corresponding to its haploid chromosome complement, and annotated 27 969 protein-coding genes and 6747 non-coding transcripts. By phylogenomic comparison with the genomes of 16 additional close and distant species we estimated that almond and peach (Prunus persica) diverged around 5.88 million years ago. These two genomes are highly syntenic and show a high degree of sequence conservation (20 nucleotide substitutions per kb). However, they also exhibit a high number of presence/absence variants, many attributable to the movement of transposable elements (TEs). Transposable elements have generated an important number of presence/absence variants between almond and peach, and we show that the recent history of TE movement seems markedly different between them. Transposable elements may also be at the origin of important phenotypic differences between both species, and in particular for the sweet kernel phenotype, a key agronomic and domestication character for almond. Here we show that in sweet almond cultivars, highly methylated TE insertions surround a gene involved in the biosynthesis of amygdalin, whose reduced expression has been correlated with the sweet almond phenotype. Altogether, our results suggest a key role of TEs in the recent history and diversification of almond and its close relative peach.

Mapping Cucumber Vein Yellowing Virus Resistance in Cucumber (Cucumis sativus L.) by Using BSA-seq Analysis.

Cucumber vein yellowing virus (CVYV) causes severe yield losses in cucurbit crops across Mediterranean countries. The control of this virus is based on cultural practices to prevent the presence of its vector (Bemisia tabaci) and breeding for natural resistance, which requires the identification of the loci involved and the development of molecular markers for linkage analysis. In this work, we mapped a monogenic locus for resistance to CVYV in cucumber by using a Bulked Segregant Analysis (BSA) strategy coupled with whole-genome resequencing. We phenotyped 135 F3 families from a segregating population between a susceptible pickling cucumber and a resistant Long Dutch type cucumber for CVYV resistance. Phenotypic analysis determined the monogenic and incomplete dominance inheritance of the resistance. We named the locus CsCvy-1. For mapping this locus, 15 resistant and 15 susceptible homozygous F2 individuals were selected for whole genome resequencing. By using a customized bioinformatics pipeline, we identified a unique region in chromosome 5 associated to resistance to CVYV, explaining more than 80% of the variability. The resequencing data provided us with additional SNP markers to decrease the interval of CsCvy-1 to 625 kb, containing 24 annotated genes. Markers flanking CsCvy-1 in a 5.3 cM interval were developed for marker-assisted selection (MAS) in breeding programs and will be useful for the identification of the target gene in future studies.

The Evolutionary Consequences of Transposon-Related Pericentromer Expansion in Melon.

Transposable elements (TEs) are a major driver of plant genome evolution. A part from being a rich source of new genes and regulatory sequences, TEs can also affect plant genome evolution by modifying genome size and shaping chromosome structure. TEs tend to concentrate in heterochromatic pericentromeric regions and their proliferation may expand these regions. Here, we show that after the split of melon and cucumber, TEs have expanded the pericentromeric regions of melon chromosomes that, probably as a consequence, show a very low recombination frequency. In contrast, TEs have not proliferated to a high extent in cucumber, which has small TE-dense pericentromeric regions and shows a relatively constant recombination rate along chromosomes. These differences in chromosome structure also translate in differences in gene nucleotide diversity. Although gene nucleotide diversity is essentially constant along cucumber chromosomes, melon chromosomes show a bimodal pattern of genetic variability, with a gene-poor region where variability is negatively correlated with gene density. Interestingly, genes are not homogeneously distributed in melon, and the high variable low-recombining pericentromeric regions show a higher concentration of melon-specific genes whereas genes shared with cucumber and other plants are essentially found in gene-rich chromosomal arms. The results presented here suggest that melon pericentromeric regions may allow gene sequences to evolve more freely than in other chromosomal compartments which may allow new ORFs to arise and eventually be selected. These results show that TEs can drastically change the structure of chromosomes creating different chromosomal compartments imposing different constraints for gene evolution.

An improved assembly and annotation of the melon (Cucumis melo L.) reference genome.

We report an improved assembly (v3.6.1) of the melon (Cucumis melo L.) genome and a new genome annotation (v4.0). The optical mapping approach allowed correcting the order and the orientation of 21 previous scaffolds and permitted to correctly define the gap-size extension along the 12 pseudomolecules. A new comprehensive annotation was also built in order to update the previous annotation v3.5.1, released more than six years ago. Using an integrative annotation pipeline, based on exhaustive RNA-Seq collections and ad-hoc transposable element annotation, we identified 29,980 protein-coding loci. Compared to the previous version, the v4.0 annotation improved gene models in terms of completeness of gene structure, UTR regions definition, intron-exon junctions and reduction of fragmented genes. More than 8,000 new genes were identified, one third of them being well supported by RNA-Seq data. To make all the new resources easily exploitable and completely available for the scientific community, a redesigned Melonomics genomic platform was released at http://melonomics.net . The resources produced in this work considerably increase the reliability of the melon genome assembly and resolution of the gene models paving the way for further studies in melon and related species.

Genetic analysis of the wild strawberry (Fragaria vesca) volatile composition.

The volatile composition of wild strawberry (Fragaria vesca) fruit differs from that of the cultivated strawberry, having more intense and fruity aromas. Over the last few years, the diploid F. vesca has been recognized as a model species for genetic studies of cultivated strawberry (F. x ananassa), and here a previously developed F. vesca/F. bucharica Near Isogenic Line collection (NIL) was used to explore genetic variability of fruit quality traits. Analysis of fruit volatiles by GC-MS in our NIL collection revealed a complex and highly variable profile. One hundred compounds were unequivocally identified, including esters, aldehydes, ketones, alcohols, terpenoids, furans and lactones. Those in a subset, named key volatile compounds (KVCs), are likely contributors to the special aroma/flavour of wild strawberry. Genetic analysis revealed 50 major quantitative trait loci (QTL) including 14 QTL for KVCs, and one segregating as a dominant monogenetic trait for nerolidol. The most determinant regions affecting QTLs for KVCs, were mapped on LG5 and LG7. New candidate genes for the volatile QTL are proposed, based on differences in gene expression between NILs containing specific fragments of F. bucharica and the F. vesca recurrent genome. A high percentage of these candidate genes/alleles were colocalized within the boundaries of introgressed regions that contain QTLs, appearing to affect volatile metabolite accumulation acting in cis. A NIL collection is a good tool for the genetic dissection of volatile accumulation in wild strawberry fruit and a source of information for genes and alleles which may enhance aroma in cultivated strawberry.

A deletion affecting an LRR-RLK gene co-segregates with the fruit flat shape trait in peach.

In peach, the flat phenotype is caused by a partially dominant allele in heterozygosis (Ss), fruits from homozygous trees (SS) abort a few weeks after fruit setting. Previous research has identified a SSR marker (UDP98-412) highly associated with the trait, found suitable for marker assisted selection (MAS). Here we report a ∼10 Kb deletion affecting the gene PRUPE.6G281100, 400 Kb upstream of UDP98-412, co-segregating with the trait. This gene is a leucine-rich repeat receptor-like kinase (LRR-RLK) orthologous to the Brassinosteroid insensitive 1-associated receptor kinase 1 (BAK1) group. PCR markers suitable for MAS confirmed its strong association with the trait in a collection of 246 cultivars. They were used to evaluate the DNA from a round fruit derived from a somatic mutation of the flat variety 'UFO-4', revealing that the mutation affected the flat associated allele (S). Protein BLAST alignment identified significant hits with genes involved in different biological processes. Best protein hit occurred with AtRLP12, which may functionally complement CLAVATA2, a key regulator that controls the stem cell population size. RT-PCR analysis revealed the absence of transcription of the partially deleted allele. The data support PRUPE.6G281100 as a candidate gene for flat shape in peach.

Regulation of plant stem cell quiescence by a brassinosteroid signaling module.

The quiescent center (QC) maintains the activity of the surrounding stem cells within the root stem cell niche, yet specific molecular players sustaining the low rate of QC cell division remain poorly understood. Here, we identified a R2R3-MYB transcription factor, BRAVO (BRASSINOSTEROIDS AT VASCULAR AND ORGANIZING CENTER), acting as a cell-specific repressor of QC divisions in the primary root of Arabidopsis. Ectopic BRAVO expression restricts overall root growth and ceases root regeneration upon damage of the stem cells, demonstrating the role of BRAVO in counteracting Brassinosteroid (BR)-mediated cell division in the QC cells. Interestingly, BR-regulated transcription factor BES1 (BRI1-EMS SUPRESSOR 1) directly represses and physically interacts with BRAVO in vivo, creating a switch that modulates QC divisions at the root stem cell niche. Together, our results define a mechanism for BR-mediated regulation of stem cell quiescence in plants.

Plant organelle proteomics: collaborating for optimal cell function.

Organelle proteomics describes the study of proteins present in organelle at a particular instance during the whole period of their life cycle in a cell. Organelles are specialized membrane bound structures within a cell that function by interacting with cytosolic and luminal soluble proteins making the protein composition of each organelle dynamic. Depending on organism, the total number of organelles within a cell varies, indicating their evolution with respect to protein number and function. For example, one of the striking differences between plant and animal cells is the plastids in plants. Organelles have their own proteins, and few organelles like mitochondria and chloroplast have their own genome to synthesize proteins for specific function and also require nuclear-encoded proteins. Enormous work has been performed on animal organelle proteomics. However, plant organelle proteomics has seen limited work mainly due to: (i) inter-plant and inter-tissue complexity, (ii) difficulties in isolation of subcellular compartments, and (iii) their enrichment and purity. Despite these concerns, the field of organelle proteomics is growing in plants, such as Arabidopsis, rice and maize. The available data are beginning to help better understand organelles and their distinct and/or overlapping functions in different plant tissues, organs or cell types, and more importantly, how protein components of organelles behave during development and with surrounding environments. Studies on organelles have provided a few good reviews, but none of them are comprehensive. Here, we present a comprehensive review on plant organelle proteomics starting from the significance of organelle in cells, to organelle isolation, to protein identification and to biology and beyond. To put together such a systematic, in-depth review and to translate acquired knowledge in a proper and adequate form, we join minds to provide discussion and viewpoints on the collaborative nature of organelles in cell, their proper function and evolution.

A cyclin-dependent protein kinase, CDKC2, colocalizes with and modulates the distribution of spliceosomal components in Arabidopsis.

Cyclin-dependent kinases (CDKs) play key regulatory roles in diverse cellular functions, including cell-cycle progression, transcription and translation. In plants, CDKs have been classified into several groups, named A through to G, but the functions of most are poorly characterized. CDKCs are known to phosphorylate the C-terminal domain (CTD) of RNA polymerase II (RNAP II), and therefore the CDKC-cyclinT (CycT) complex may have a role similar to the animal CDK9-CycT complex of the positive transcription elongation factor b (P-TEFb). However, we found that the predicted structure of the Arabidopsis CDKC2 protein is more similar to the mammalian cdc2-related kinase, CRK7, than to CDK9. CRK7 is proposed to link transcription with splicing, and CDKC2 contains all the structural features of CRK7 that make the latter distinct from CDK9. Consistent with this, we show that GFP-CDKC2 fusion proteins co-localize with spliceosomal components, that the expression of CDKC2 modifies the location of these components, and that co-localization was dependent on the transcriptional status of the cells and on CDKC2-kinase activity. We propose, therefore, that the Arabidopsis CDKC2 combines the functions of both CRK7 and CDK9, and could also couple splicing with transcription.