Single nucleotide polymorphism (SNP) analysis reveals ancestry and genetic diversity of cultivated and wild grapevines in Croatia.

BACKGROUND: Croatia is a geographically small country with a remarkable diversity of cultivated and spontaneous grapevines. Local germplasm has been characterised by microsatellite markers, but a detailed analysis based on single nucleotide polymorphisms (SNPs) is still lacking. Here we characterize the genetic diversity of 149 accessions from three germplasm repositories and four natural sites using 516,101 SNPs to identify complete parent-offspring trios and their relations with spontaneous populations, offering a proof-of-concept for the use of reduced-representation genome sequencing in population genetics and genome-wide association studies (GWAS). RESULTS: Principal component analysis revealed a clear discontinuity between cultivated (V. vinifera subsp. sativa) and spontaneous grapevines, supporting the notion that the latter represent local populations of the wild progenitor (V. vinifera subsp. sylvestris). ADMIXTURE identified three ancestry components. Two sativa components are alternatively predominant in cultivars grown either in northern Adriatic Croatia and Continental Croatia or in Dalmatia (i.e. central and southern Adriatic Croatia). A sylvestris component, which is predominant in accessions from spontaneous populations, is a minor ancestry component in cultivated accessions. TREEMIX provided evidence of unidirectional migration from the vineyards to natural sites, suggesting that gene flow has gone preferentially from the introduced domesticated germplasm into local wild populations rather than vice versa. Identity-by-descent analysis indicated an extensive kinship network, including 14 complete parent-offspring trios, involving only cultivated accessions, six full-sibling relationships and invalidated a presumed pedigree of one of the most important varieties in Croatia, 'Plavac Mali'. Despite this strong population structure, significant association was found between 143 SNPs and berry skin colour and between 2 SNPs and leaf hairiness, across two previously known genomic regions. CONCLUSIONS: The clear genetic separation between Croatian cultivars and sylvestris ruled out the hypothesis that those cultivars originated from local domestication events. On the other hand, the evidence of a crop-to-wild gene flow signals the need for an urgent adoption of conservation strategies that preserve the residual genetic integrity of wild relatives. The use of this reduced-representation genome sequencing protocol in grapevine enables an accurate pedigree reconstruction and can be recommended for GWAS experiments.

Open chromatin in grapevine marks candidate CREs and with other chromatin features correlates with gene expression.

Vitis vinifera is an economically important crop and a useful model in which to study chromatin dynamics. In contrast to the small and relatively simple genome of Arabidopsis thaliana, grapevine contains a complex genome of 487 Mb that exhibits extensive colonization by transposable elements. We used Hi-C, ChIP-seq and ATAC-seq to measure how chromatin features correlate to the expression of 31 845 grapevine genes. ATAC-seq revealed the presence of more than 16 000 open chromatin regions, of which we characterize nearly 5000 as possible distal enhancer candidates that occur in intergenic space > 2 kb from the nearest transcription start site (TSS). A motif search identified more than 480 transcription factor (TF) binding sites in these regions, with those for TCP family proteins in greatest abundance. These open chromatin regions are typically within 15 kb from their nearest promoter, and a gene ontology analysis indicated that their nearest genes are significantly enriched for TF activity. The presence of a candidate cis-regulatory element (cCRE) > 2 kb upstream of the TSS, location in the active nuclear compartment as determined by Hi-C, and the enrichment of H3K4me3, H3K4me1 and H3K27ac at the gene are correlated with gene expression. Taken together, these results suggest that regions of intergenic open chromatin identified by ATAC-seq can be considered potential candidates for cis-regulatory regions in V. vinifera. Our findings enhance the characterization of a valuable agricultural crop, and help to clarify the understanding of unique plant biology.

ddRAD-seq reveals the genetic structure and detects signals of selection in Italian brown trout.

BACKGROUND: Brown trout is one of the most widespread fresh-water fish species in Europe. The evolutionary history of and phylogenetic relationships between brown trout populations are complex, and this is especially true for Italian populations, which are heavily influenced in different ways by stocking practices. The characterization of the genetic structure of Italian brown trout populations may give information on the risk of losing endemic Italian populations due to lack of genetic diversity or to admixture with stocking populations. The identification of signatures of selection, and the information deriving from dense genotyping data will help genotype-informed breeding programs. We used a ddRAD-seq approach to obtain more than 100,000 single nucleotide polymorphisms (SNPs), and to characterize the population structure and signatures of selection in 90 brown trout samples. RESULTS: Italian brown trout populations are genetically differentiated, although the stocking practices have introduced strong admixture in endemic Italian trout, especially with the Atlantic lineage. Most of the analysed populations showed high levels of kinship and inbreeding. We detected putative signatures of selection using different approaches, and investigated if the regions were enriched for functional categories. Several regions putatively under selection and characterized by a reduction in heterozygosity across all the studied populations are enriched for genes involved in the response to viral infections. CONCLUSIONS: Our results, which show evidence of admixture with the Atlantic lineage (commonly used for stocking), confirm the need for controlling stocking practices, in order to avoid the erosion of the endemic gene pool; given the apparently high levels of kinship and inbreeding in local populations, our results also show the need to take action for increasing gene diversity. In addition, we used the genetically-distinct lineages to detect signatures of selection and we identified putative signatures of selection in several regions associated with resistance to infectious diseases. These constitute candidate regions for the study of resistance to infections in wild and farmed trout.

Gene duplication and transposition of mobile elements drive evolution of the Rpv3 resistance locus in grapevine.

A wild grape haplotype (Rpv3-1) confers resistance to Plasmopara viticola. We mapped the causal factor for resistance to an interval containing a TIR-NB-LRR (TNL) gene pair that originated 1.6-2.6 million years ago by a tandem segmental duplication. Transient coexpression of the TNL pair in Vitis vinifera leaves activated pathogen-induced necrosis and reduced sporulation compared with control leaves. Even though transcripts of the TNL pair from the wild haplotype appear to be partially subject to nonsense-mediated mRNA decay, mature mRNA levels in a homozygous resistant genotype were individually higher than the mRNA trace levels observed for the orthologous single-copy TNL in sensitive genotypes. Allelic expression imbalance in a resistant heterozygote confirmed that cis-acting regulatory variation promotes expression in the wild haplotype. The movement of transposable elements had a major impact on the generation of haplotype diversity, altering the DNA context around similar TNL coding sequences and the GC-content in their proximal 5'-intergenic regions. The wild and domesticated haplotypes also diverged in conserved single-copy intergenic DNA, but the highest divergence was observed in intraspecific and not in interspecific comparisons. In this case, introgression breeding did not transgress the genetic boundaries of the domesticated species, because haplotypes present in modern varieties sometimes predate speciation events between wild and cultivated species.

A draft genome of sweet cherry (Prunus avium L.) reveals genome-wide and local effects of domestication.

Sweet cherry (Prunus avium L.) trees are both economically important fruit crops but also important components of natural forest ecosystems in Europe, Asia and Africa. Wild and domesticated trees currently coexist in the same geographic areas with important questions arising on their historical relationships. Little is known about the effects of the domestication process on the evolution of the sweet cherry genome. We assembled and annotated the genome of the cultivated variety "Big Star*" and assessed the genetic diversity among 97 sweet cherry accessions representing three different stages in the domestication and breeding process (wild trees, landraces and modern varieties). The genetic diversity analysis revealed significant genome-wide losses of variation among the three stages and supports a clear distinction between wild and domesticated trees, with only limited gene flow being detected between wild trees and domesticated landraces. We identified 11 domestication sweeps and five breeding sweeps covering, respectively, 11.0 and 2.4 Mb of the P. avium genome. A considerable fraction of the domestication sweeps overlaps with those detected in the related species, Prunus persica (peach), indicating that artificial selection during domestication may have acted independently on the same regions and genes in the two species. We detected 104 candidate genes in sweep regions involved in different processes, such as the determination of fruit texture, the regulation of flowering and fruit ripening and the resistance to pathogens. The signatures of selection identified will enable future evolutionary studies and provide a valuable resource for genetic improvement and conservation programs in sweet cherry.

Genetic, epigenetic and genomic effects on variation of gene expression among grape varieties.

The transcriptional regulatory structure of plant genomes is still relatively unexplored, and little is known about factors that influence expression variation in plants. We used a genetic system consisting of 10 heterozygous grape varieties with high consanguinity and high haplotypic diversity to: (i) identify regions of haplotype sharing through whole-genome resequencing and single-nucleotide polymorphism (SNP) genotyping; (ii) analyse gene expression through RNA-seq in four stages of berry development; and (iii) associate gene expression variation with genetic and epigenetic properties. We found that haplotype sharing in and around genes was positively correlated with similarity in expression and was negatively correlated with the fraction of differentially expressed genes. Genetic and epigenetic properties of the gene and the surrounding region showed significant effects on the extent of expression variation, with negative associations for the level of gene body methylation and mean expression level, and with positive associations for nucleotide diversity, structural diversity and ratio of non-synonymous to synonymous nucleotide diversity. We also observed a spatial dependency of covariation of gene expression among varieties. These results highlight relevant roles for cis-acting factors, selective constraints and epigenetic features of the gene, and the regional context in which the gene is located, in the determination of expression variation. OPEN RESEARCH BADGES: This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://www.ncbi.nlm.nih.gov/bioproject/PRJNA385116; https://www.ncbi.nlm.nih.gov/bioproject/PRJNA392287; https://www.ncbi.nlm.nih.gov/bioproject/PRJNA373967 (released upon publication); https://www.ncbi.nlm.nih.gov/bioproject/PRJNA490160 (released upon publication); https://www.ncbi.nlm.nih.gov/bioproject/PRJNA265039; https://www.ncbi.nlm.nih.gov/bioproject/PRJNA265040.

Grapevine field experiments reveal the contribution of genotype, the influence of environment and the effect of their interaction (G×E) on the berry transcriptome.

Changes in the performance of genotypes in different environments are defined as genotype × environment (G×E) interactions. In grapevine (Vitis vinifera), complex interactions between different genotypes and climate, soil and farming practices yield unique berry qualities. However, the molecular basis of this phenomenon remains unclear. To dissect the basis of grapevine G×E interactions we characterized berry transcriptome plasticity, the genome methylation landscape and within-genotype allelic diversity in two genotypes cultivated in three different environments over two vintages. We identified, through a novel data-mining pipeline, genes with expression profiles that were: unaffected by genotype or environment, genotype-dependent but unaffected by the environment, environmentally-dependent regardless of genotype, and G×E-related. The G×E-related genes showed different degrees of within-cultivar allelic diversity in the two genotypes and were enriched for stress responses, signal transduction and secondary metabolism categories. Our study unraveled the mutual relationships between genotypic and environmental variables during G×E interaction in a woody perennial species, providing a reference model to explore how cultivated fruit crops respond to diverse environments. Also, the pivotal role of vineyard location in determining the performance of different varieties, by enhancing berry quality traits, was unraveled.

Single primer enrichment technology as a tool for massive genotyping: a benchmark on black poplar and maize.

BACKGROUND AND AIMS: The advent of molecular breeding is advocated to improve the productivity and sustainability of second-generation bioenergy crops. Advanced molecular breeding in bioenergy crops relies on the ability to massively sample the genetic diversity. Genotyping-by-sequencing has become a widely adopted method for cost-effective genotyping. It basically requires no initial investment for design as compared with array-based platforms which have been shown to offer very robust assays. The latter, however, has the drawback of being limited to analyse only the genetic diversity accounted during selection of a set of polymorphisms and design of the assay. In contrast, genotyping-by-sequencing with random sampling of genomic loci via restriction enzymes or random priming has been shown to be fast and convenient but lacks the ability to target specific regions of the genome and to maintain high reproducibility across laboratories. METHODS: Here we present a first adoption of single-primer enrichment technology (SPET) which provides a highly efficient and scalable system to obtain targeted sequence-based large genotyping data sets, bridging the gaps between array-based systems and traditional sequencing-based protocols. To fully explore SPET performance, we conducted a benchmark study in ten Zea mays lines and a large-scale study of a natural black poplar population of 540 individuals with the aim of discovering polymorphisms associated with biomass-related traits. KEY RESULTS: Our results showed the ability of this technology to provide dense genotype information on a customized panel of selected polymorphisms, while yielding hundreds of thousands of untargeted variable sites. This provided an ideal resource for association analysis of natural populations harbouring unexplored allelic diversities and structure such as in black poplar. CONCLUSION: The improvement of sequencing throughput and the development of efficient library preparation methods has made it feasible to carry out targeted genotyping-by-sequencing experiments cost-competitively with either random complexity reduction systems or traditional array-based platforms, while maintaining the key advantages of both technologies.

A chromosome-scale assembly reveals chromosomal aberrations and exchanges generating genetic diversity in Coffea arabica germplasm.

In order to better understand the mechanisms generating genetic diversity in the recent allotetraploid species Coffea arabica, here we present a chromosome-level assembly obtained with long read technology. Two genomic compartments with different structural and functional properties are identified in the two homoeologous genomes. The resequencing data from a large set of accessions reveals low intraspecific diversity in the center of origin of the species. Across a limited number of genomic regions, diversity increases in some cultivated genotypes to levels similar to those observed within one of the progenitor species, Coffea canephora, presumably as a consequence of introgressions deriving from the so-called Timor hybrid. It also reveals that, in addition to few, early-occurring exchanges between homoeologous chromosomes, there are numerous recent chromosomal aberrations including aneuploidies, deletions, duplications and exchanges. These events are still polymorphic in the germplasm and could represent a fundamental source of genetic variation in such a lowly variable species.

Evaluation of sensitivity and specificity in RNA-Seq-based detection of grapevine viral pathogens.

Virus detection is a crucial step for the implementation of clean stock programs that preserve healthy crop species. Viral infections in grapevine, a vegetatively propagated perennial crop, cannot be eradicated from the vineyards by the application of agrochemicals and must be curtailed at the stage of nursery production during the propagation of planting material. Viral detection is routinely performed using enzyme-linked immunosorbent assays (ELISA) or Reverse Transcription-quantitative Polymerase Chain Reactions (RT-qPCR). High throughput sequencing (HTS) approaches have the potential to detect all viral pathogens in a plant specimen. However, to date, no published HTS-based study has used threshold selection based on ROC curves for discriminating positive from negative samples. To fill this gap, we assessed the specificity and sensitivity of different sequencing and bioinformatics approaches for nine common viruses, which were tested in the same specimens using ELISA and/or RT-qPCR. The normalized detection thresholds giving the best results were 19.28 Fragments Per Kilobase of transcript per Million mapped reads (FPKM) for alignment-based total RNA-Seq approaches, 386 Reads Per Million mapped reads (RPM) for metagenomics-based total RNA-Seq, 1572 FPKM for alignment-based small RNA-Seq analysis and 0.97 % of contigs for de novo analysis of small RNA-Seq data. Validation of the proposed thresholds using independent specimens collected over time from the same stocks and other specimens collected from nearby stocks that had derived from the same propagating material showed that HTS approaches are accurate, with RNA-Seq approaches showing better performance than small RNA-Seq.

Extent of wild-to-crop interspecific introgression in grapevine (Vitis vinifera) as a consequence of resistance breeding and implications for the crop species definition.

Over the past two centuries, introgression through repeated backcrossing has introduced disease resistance from wild grape species into the domesticated lineage Vitis vinifera subsp. sativa. Introgression lines are being cultivated over increasing vineyard surface areas, as their wines now rival in quality those obtained from preexisting varieties. There is, however, a lot of debate about whether and how wine laws defining commercial product categories, which are based on the classification of V. vinifera and interspecific hybrid grapes, should be revised to accommodate novel varieties that do not fit either category. Here, we developed a method of multilocus genotype analysis using short-read resequencing to identify haplotypic blocks of wild ancestry in introgression lines and quantify the physical length of chromosome segments free-of-introgression or with monoallelic and biallelic introgression. We used this genomic data to characterize species, hybrids and introgression lines and show that newly released resistant varieties contain 76.5-94.8% of V. vinifera DNA. We found that varietal wine ratings are not always commensurate with the percentage of V. vinifera ancestry and linkage drag of wild alleles around known resistance genes persists over at least 7.1-11.5 Mb, slowing down the recovery of the recurrent parental genome. This method also allowed us to identify the donor species of known resistance haplotypes, define the ancestry of wild genetic background in introgression lines with complex pedigrees, validate the ancestry of the historic varieties Concord and Norton, and unravel sample curation errors in public databases.

The genomes of 204 Vitis vinifera accessions reveal the origin of European wine grapes.

In order to elucidate the still controversial processes that originated European wine grapes from its wild progenitor, here we analyse 204 genomes of Vitis vinifera and show that all analyses support a single domestication event that occurred in Western Asia and was followed by numerous and pervasive introgressions from European wild populations. This admixture generated the so-called international wine grapes that have diffused from Alpine countries worldwide. Across Europe, marked differences in genomic diversity are observed in local varieties that are traditionally cultivated in different wine producing countries, with Italy and France showing the largest diversity. Three genomic regions of reduced genetic diversity are observed, presumably as a consequence of artificial selection. In the lowest diversity region, two candidate genes that gained berry-specific expression in domesticated varieties may contribute to the change in berry size and morphology that makes the fruit attractive for human consumption and adapted for winemaking.

First detection of an Italian human-to-cat outbreak of SARS-CoV-2 Alpha variant - lineage B.1.1.7.

The emergence of new SARS-CoV-2 variants and their rapid spread pose a threat to both human and animal health and may conceal unknown risks. This report describes an Italian human-to-cat outbreak of SARS-CoV-2 lineage B.1.1.7 (the Alpha variant) . On March 7th, 2021, approximately ten days after COVID-19 appeared in the family, the onset of respiratory signs in a cat by COVID-19-affected owners led to an in-depth diagnostic investigation, combining clinical and serological data with rt-qPCR-based virus detection and whole genome sequencing. The Alpha variant was confirmed first in the owners and a few days later in the cat that was then monitored weekly: the course was similar with one-week lag time in the cat. In addition, based on comparative analysis of genome sequences from our study and from 200 random Italian cases of Alpha variant, the familial cluster was confirmed. The temporal sequence along with the genomic data support a human-to-animal transmission. Such an event emphasizes the importance of studying the circulation and dynamics of SARS-CoV-2 variants in humans and animals to better understand and prevent potential spillover risks or unwarranted alerts involving our pet populations.

Two-omics data revealed commonalities and differences between Rpv12- and Rpv3-mediated resistance in grapevine.

Plasmopara viticola is the causal agent of grapevine downy mildew (DM). DM resistant varieties deploy effector-triggered immunity (ETI) to inhibit pathogen growth, which is activated by major resistance loci, the most common of which are Rpv3 and Rpv12. We previously showed that a quick metabolome response lies behind the ETI conferred by Rpv3 TIR-NB-LRR genes. Here we used a grape variety operating Rpv12-mediated ETI, which is conferred by an independent locus containing CC-NB-LRR genes, to investigate the defence response using GC/MS, UPLC, UHPLC and RNA-Seq analyses. Eighty-eight metabolites showed significantly different concentration and 432 genes showed differential expression between inoculated resistant leaves and controls. Most metabolite changes in sugars, fatty acids and phenols were similar in timing and direction to those observed in Rpv3-mediated ETI but some of them were stronger or more persistent. Activators, elicitors and signal transducers for the formation of reactive oxygen species were early observed in samples undergoing Rpv12-mediated ETI and were paralleled and followed by the upregulation of genes belonging to ontology categories associated with salicylic acid signalling, signal transduction, WRKY transcription factors and synthesis of PR-1, PR-2, PR-5 pathogenesis-related proteins.

A single polyploidization event at the origin of the tetraploid genome of Coffea arabica is responsible for the extremely low genetic variation in wild and cultivated germplasm.

The genome of the allotetraploid species Coffea arabica L. was sequenced to assemble independently the two component subgenomes (putatively deriving from C. canephora and C. eugenioides) and to perform a genome-wide analysis of the genetic diversity in cultivated coffee germplasm and in wild populations growing in the center of origin of the species. We assembled a total length of 1.536 Gbp, 444 Mb and 527 Mb of which were assigned to the canephora and eugenioides subgenomes, respectively, and predicted 46,562 gene models, 21,254 and 22,888 of which were assigned to the canephora and to the eugeniodes subgenome, respectively. Through a genome-wide SNP genotyping of 736 C. arabica accessions, we analyzed the genetic diversity in the species and its relationship with geographic distribution and historical records. We observed a weak population structure due to low-frequency derived alleles and highly negative values of Taijma's D, suggesting a recent and severe bottleneck, most likely resulting from a single event of polyploidization, not only for the cultivated germplasm but also for the entire species. This conclusion is strongly supported by forward simulations of mutation accumulation. However, PCA revealed a cline of genetic diversity reflecting a west-to-east geographical distribution from the center of origin in East Africa to the Arabian Peninsula. The extremely low levels of variation observed in the species, as a consequence of the polyploidization event, make the exploitation of diversity within the species for breeding purposes less interesting than in most crop species and stress the need for introgression of new variability from the diploid progenitors.

Genetic Mapping of the Incompatibility Locus in Olive and Development of a Linked Sequence-Tagged Site Marker.

The genetic control of self-incompatibility (SI) has been recently disclosed in olive. Inter-varietal crossing confirmed the presence of only two incompatibility groups (G1 and G2), suggesting a simple Mendelian inheritance of the trait. A double digest restriction associated DNA (ddRAD) sequencing of a biparental population segregating for incompatibility groups has been performed and high-density linkage maps were constructed in order to map the SI locus and identify gene candidates and linked markers. The progeny consisted of a full-sib family of 229 individuals derived from the cross 'Leccino' (G1) × 'Dolce Agogia' (G2) varieties, segregating 1:1 (G1:G2), in accordance with a diallelic self-incompatibility (DSI) model. A total of 16,743 single nucleotide polymorphisms was identified, 7,006 in the female parent 'Leccino' and 9,737 in the male parent 'Dolce Agogia.' Each parental map consisted of 23 linkage groups and showed an unusual large size (5,680 cM in 'Leccino' and 3,538 cM in 'Dolce Agogia'). Recombination was decreased across all linkage groups in pollen mother cells of 'Dolce Agogia,' the parent with higher heterozygosity, compared to megaspore mother cells of 'Leccino,' in a context of a species that showed exceptionally high recombination rates. A subset of 109 adult plants was assigned to either incompatibility group by a stigma test and the diallelic self-incompatibility (DSI) locus was mapped to an interval of 5.4 cM on linkage group 18. This region spanned a size of approximately 300 Kb in the olive genome assembly. We developed a sequence-tagged site marker in the DSI locus and identified five haplotypes in 57 cultivars with known incompatibility group assignment. A combination of two single-nucleotide polymorphisms (SNPs) was sufficient to predict G1 or G2 phenotypes in olive cultivars, enabling early marker-assisted selection of compatible genotypes and allowing for a rapid screening of inter-compatibility among cultivars in order to guarantee effective fertilization and increase olive production. The construction of high-density linkage maps has led to the development of the first functional marker in olive and provided positional candidate genes in the SI locus.

Genetic properties of the MAGIC maize population: a new platform for high definition QTL mapping in Zea mays.

BACKGROUND: Maize (Zea mays) is a globally produced crop with broad genetic and phenotypic variation. New tools that improve our understanding of the genetic basis of quantitative traits are needed to guide predictive crop breeding. We have produced the first balanced multi-parental population in maize, a tool that provides high diversity and dense recombination events to allow routine quantitative trait loci (QTL) mapping in maize. RESULTS: We produced 1,636 MAGIC maize recombinant inbred lines derived from eight genetically diverse founder lines. The characterization of 529 MAGIC maize lines shows that the population is a balanced, evenly differentiated mosaic of the eight founders, with mapping power and resolution strengthened by high minor allele frequencies and a fast decay of linkage disequilibrium. We show how MAGIC maize may find strong candidate genes by incorporating genome sequencing and transcriptomics data. We discuss three QTL for grain yield and three for flowering time, reporting candidate genes. Power simulations show that subsets of MAGIC maize might achieve high-power and high-definition QTL mapping. CONCLUSIONS: We demonstrate MAGIC maize's value in identifying the genetic bases of complex traits of agronomic relevance. The design of MAGIC maize allows the accumulation of sequencing and transcriptomics layers to guide the identification of candidate genes for a number of maize traits at different developmental stages. The characterization of the full MAGIC maize population will lead to higher power and definition in QTL mapping, and lay the basis for improved understanding of maize phenotypes, heterosis included. MAGIC maize is available to researchers.