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.

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.

Genomic tools for durum wheat breeding: de novo assembly of Svevo transcriptome and SNP discovery in elite germplasm.

BACKGROUND: The tetraploid durum wheat (Triticum turgidum L. ssp. durum Desf. Husnot) is an important crop which provides the raw material for pasta production and a valuable source of genetic diversity for breeding hexaploid wheat (Triticum aestivum L.). Future breeding efforts to enhance yield potential and climate resilience will increasingly rely on genomics-based approaches to identify and select beneficial alleles. A deeper characterisation of the molecular and functional diversity of the durum wheat transcriptome will be instrumental to more effectively harness its genetic diversity. RESULTS: We report on the de novo transcriptome assembly of durum wheat cultivar 'Svevo'. The transcriptome of four tissues/organs (shoots and roots at the seedling stage, reproductive organs and developing grains) was assembled de novo, yielding 180,108 contigs, with a N50 length of 1121 bp and mean contig length of 883 bp. Alignment against the transcriptome of nine plant species identified 43% of transcripts with homology to at least one reference transcriptome. The functional annotation was completed by means of a combination of complementary software. The presence of differential expression between the A- and B-homoeolog copies of the durum wheat tetraploid genome was ascertained by phase reconstruction of polymorphic sites based on the T. urartu transcripts and inferring homoeolog-specific sequences. We observed greater expression divergence between A and B homoeologs in grains rather than in leaves and roots. The transcriptomes of 13 durum wheat cultivars spanning the breeding period from 1969 to 2005 were analysed for SNP diversity, leading to 95,358 non-rare, hemi-SNPs shared among two or more cultivars and 33,747 locus-specific (diploid inheritance) SNPs. CONCLUSIONS: Our study updates and expands the de novo transcriptome reference assembly available for durum wheat. Out of 180,108 assembled transcripts, 13,636 were specific to the Svevo cultivar as compared to the only other reference transcriptome available for durum, thus contributing to the identification of the tetraploid wheat pan-transcriptome. Additionally, the analysis of 13 historically relevant hallmark varieties produced a SNP dataset that could successfully validate the genotyping in tetraploid wheat and provide a valuable resource for genomics-assisted breeding of both tetraploid and hexaploid wheats.

The microbiota of hematophagous ectoparasites collected from migratory birds.

Arthropod vectors are responsible for the transmission of human pathogens worldwide. Several arthropod species are bird ectoparasites, however, no study to date has characterized their microbiota as a whole. We sampled hematophagous ectoparasites that feed on migratory birds and performed 16S rRNA gene metabarcoding to characterize their microbial community. A total of 194 ectoparasites were collected from 115 avian hosts and classified into three groups: a) Hippoboscidae diptera; b) ticks; c) other arthropods. Metabarcoding showed that endosymbionts were the most abundant genera of the microbial community, including Wolbachia for Hippoboscidae diptera, Candidatus Midichloria for ticks, Wolbachia and Arsenophonus for the other arthropod group. Genera including pathogenic species were: Rickettsia, Borrelia, Coxiella, Francisella, Bartonella, Anaplasma. Co-infection with Borrelia-Rickettsia and Anaplasma-Rickettsia was also observed. A global overview of the microbiota of ectoparasites sampled from migratory birds was obtained with the use of 16S rRNA gene metabarcoding. A novel finding is the first identification of Rickettsia in the common swift louse fly, Crataerina pallida. Given their possible interaction with pathogenic viruses and bacteria, the presence of endosymbionts in arthropods merits attention. Finally, molecular characterization of genera, including both pathogenic and symbiont species, plays a pivotal role in the design of targeted molecular diagnostics.

Regeneration of the entire human epidermis using transgenic stem cells.

Junctional epidermolysis bullosa (JEB) is a severe and often lethal genetic disease caused by mutations in genes encoding the basement membrane component laminin-332. Surviving patients with JEB develop chronic wounds to the skin and mucosa, which impair their quality of life and lead to skin cancer. Here we show that autologous transgenic keratinocyte cultures regenerated an entire, fully functional epidermis on a seven-year-old child suffering from a devastating, life-threatening form of JEB. The proviral integration pattern was maintained in vivo and epidermal renewal did not cause any clonal selection. Clonal tracing showed that the human epidermis is sustained not by equipotent progenitors, but by a limited number of long-lived stem cells, detected as holoclones, that can extensively self-renew in vitro and in vivo and produce progenitors that replenish terminally differentiated keratinocytes. This study provides a blueprint that can be applied to other stem cell-mediated combined ex vivo cell and gene therapies.

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.

Functional features of a single chromosome arm in wheat (1AL) determined from its structure.

Bread wheat (Triticum aestivum L.) is one of the most important crops globally and a high priority for genetic improvement, but its large and complex genome has been seen as intractable to whole genome sequencing. Isolation of individual wheat chromosome arms has facilitated large-scale sequence analyses. However, so far there is no such survey of sequences from the A genome of wheat. Greater understanding of an A chromosome could facilitate wheat improvement and future sequencing of the entire genome. We have constructed BAC library from the long arm of T. aestivum chromosome 1A (1AL) and obtained BAC end sequences from 7,470 clones encompassing the arm. We obtained 13,445 (89.99%) useful sequences with a cumulative length of 7.57 Mb, representing 1.43% of 1AL and about 0.14% of the entire A genome. The GC content of the sequences was 44.7%, and 90% of the chromosome was estimated to comprise repeat sequences, while just over 1% encoded expressed genes. From the sequence data, we identified a large number of sites suitable for development of molecular markers (362 SSR and 6,948 ISBP) which will have utility for mapping this chromosome and for marker assisted breeding. From 44 putative ISBP markers tested 23 (52.3%) were found to be useful. The BAC end sequence data also enabled the identification of genes and syntenic blocks specific to chromosome 1AL, suggesting regions of particular functional interest and targets for future research.

Large-scale detection of rare variants via pooled multiplexed next-generation sequencing: towards next-generation Ecotilling.

Common variants, such as those identified by genome-wide association scans, explain only a small proportion of trait variation. Growing evidence suggests that rare functional variants, which are usually missed by genome-wide association scans, play an important role in determining the phenotype. We used pooled multiplexed next-generation sequencing and a customized analysis workflow to detect mutations in five candidate genes for lignin biosynthesis in 768 pooled Populus nigra accessions. We identified a total of 36 non-synonymous single nucleotide polymorphisms, one of which causes a premature stop codon. The most common variant was estimated to be present in 672 of the 1536 tested chromosomes, while the rarest was estimated to occur only once in 1536 chromosomes. Comparison with individual Sanger sequencing in a selected sub-sample confirmed that variants are identified with high sensitivity and specificity, and that the variant frequency was estimated accurately. This proposed method for identification of rare polymorphisms allows accurate detection of variation in many individuals, and is cost-effective compared to individual sequencing.

The SSR-based molecular profile of 1005 grapevine (Vitis vinifera L.) accessions uncovers new synonymy and parentages, and reveals a large admixture amongst varieties of different geographic origin.

A collection of 1005 grapevine accessions was genotyped at 34 microsatellite loci (SSR) with the aim of analysing genetic diversity and exploring parentages. The comparison of molecular profiles revealed 200 groups of synonymy. The removal of perfect synonyms reduced the database to 745 unique genotypes, on which population genetic parameters were calculated. The analysis of kinship uncovered 74 complete pedigrees, with both parents identified. Many of these parentages were not previously known and are of considerable historical interest, e.g. Chenin blanc (Sauvignon × Traminer rot), Covè (Harslevelu selfed), Incrocio Manzoni 2-14 and 2-15 (Cabernet franc × Prosecco), Lagrein (Schiava gentile × Teroldego), Malvasia nera of Bolzano (Perera × Schiava gentile), Manzoni moscato (Raboso veronese × Moscato d'Amburgo), Moscato violetto (Moscato bianco × Duraguzza), Muscat of Alexandria (Muscat blanc à petit grain × Axina de tres bias) and others. Statistical robustness of unexpected pedigrees was reinforced with the analysis of an additional 7-30 SSRs. Grouping the accessions by profile resulted in a weak correlation with their geographical origin and/or current area of cultivation, revealing a large admixture of local varieties with those most widely cultivated, as a result of ancient commerce and population flow. The SSRs with tri- to penta-nucleotide repeats adopted for the present study showed a great capacity for discriminating amongst accessions, with probabilities of identity by chance as low as 1.45 × 10(-27) and 9.35 × 10(-12) for unrelated and full sib individuals, respectively. A database of allele frequencies and SSR profiles of 32 reference cultivars are provided.

Quick MLPA test for quantification of SMN1 and SMN2 copy numbers.

Spinal muscular atrophy (SMA) is an autosomal recessive disease caused in about 95% of SMA patients by homozygous deletion of the survival motor neuron 1 (SMN1) gene or its conversion to the highly homologous SMN2 gene. In the majority of cases, disease severity correlates inversely with increased SMN2 copy number. Because of the comparatively high incidence of healthy carriers and severity of the disease, detection of sequence alterations and quantification of SMN1 and SMN2 copy numbers are essential for exact diagnosis and genetic counselling. Several assays have been developed for this purpose. Multiplex ligation-dependent probe amplification (MLPA) is a versatile technique for relative quantification of different nucleic acid sequences in a single reaction. Here, we establish a quick MLPA-based assay for the detection of SMN1 and SMN2 copy numbers with high specificity and low complexity.

Physical mapping in highly heterozygous genomes: a physical contig map of the Pinot Noir grapevine cultivar.

BACKGROUND: Most of the grapevine (Vitis vinifera L.) cultivars grown today are those selected centuries ago, even though grapevine is one of the most important fruit crops in the world. Grapevine has therefore not benefited from the advances in modern plant breeding nor more recently from those in molecular genetics and genomics: genes controlling important agronomic traits are practically unknown. A physical map is essential to positionally clone such genes and instrumental in a genome sequencing project. RESULTS: We report on the first whole genome physical map of grapevine built using high information content fingerprinting of 49,104 BAC clones from the cultivar Pinot Noir. Pinot Noir, as most grape varieties, is highly heterozygous at the sequence level. This resulted in the two allelic haplotypes sometimes assembling into separate contigs that had to be accommodated in the map framework or in local expansions of contig maps. We performed computer simulations to assess the effects of increasing levels of sequence heterozygosity on BAC fingerprint assembly and showed that the experimental assembly results are in full agreement with the theoretical expectations, given the heterozygosity levels reported for grape. The map is anchored to a dense linkage map consisting of 994 markers. 436 contigs are anchored to the genetic map, covering 342 of the 475 Mb that make up the grape haploid genome. CONCLUSIONS: We have developed a resource that makes it possible to access the grapevine genome, opening the way to a new era both in grape genetics and breeding and in wine making. The effects of heterozygosity on the assembly have been analyzed and characterized by using several complementary approaches which could be easily transferred to the study of other genomes which present the same features.

Analysis of transposons and repeat composition of the sunflower (Helianthus annuus L.) genome.

A sample-sequencing strategy combined with slot-blot hybridization and FISH was used to study the composition of the repetitive component of the sunflower genome. One thousand six hundred thirty-eight sequences for a total of 954,517 bp were analyzed. The fraction of sequences that can be classified as repetitive using computational and hybridization approaches amounts to 62% in total. Almost two thirds remain as yet uncharacterized in nature. Of those characterized, most belong to the gypsy superfamily of LTR-retrotransposons. Unlike in other species, where single families can account for large fractions of the genome, it appears that no transposon family has been amplified to very high levels in sunflower. All other known classes of transposable elements were also found. One family of unknown nature (contig 61) was the most repeated in the sunflower genome. The evolution of the repetitive component in the Helianthus genus and in other Asteraceae was studied by comparative analysis of the hybridization of total genomic DNAs from these species to the sunflower small-insert library and compared to gene-based phylogeny. Very little similarity is observed between Helianthus species and two related Asteraceae species outside of the genus. Most repetitive elements are similar in annual and perennial Helianthus species indicating that sequence amplification largely predates such divergence. Gypsy-like elements are more represented in the annuals than in the perennials, while copia-like elements are similarly represented, attesting a different amplification history of the two superfamilies of LTR-retrotransposons in the Helianthus genus.

Multilocus patterns of nucleotide diversity, linkage disequilibrium and demographic history of Norway spruce [Picea abies (L.) Karst].

DNA polymorphism at 22 loci was studied in an average of 47 Norway spruce [Picea abies (L.) Karst.] haplotypes sampled in seven populations representative of the natural range. The overall nucleotide variation was limited, being lower than that observed in most plant species so far studied. Linkage disequilibrium was also restricted and did not extend beyond a few hundred base pairs. All populations, with the exception of the Romanian population, could be divided into two main domains, a Baltico-Nordic and an Alpine one. Mean Tajima's D and Fay and Wu's H across loci were both negative, indicating the presence of an excess of both rare and high-frequency-derived variants compared to the expected frequency spectrum in a standard neutral model. Multilocus neutrality tests based on D and H led to the rejection of the standard neutral model and exponential growth in the whole population as well as in the two main domains. On the other hand, in all three cases the data are compatible with a severe bottleneck occurring some hundreds of thousands of years ago. Hence, demographic departures from equilibrium expectations and population structure will have to be accounted for when detecting selection at candidate genes and in association mapping studies, respectively.