Analysis of high-identity segmental duplications in the grapevine genome

Distinct structural variants and repeat landscape shape the genomes of the ancient grapes Aglianico and Falanghina

Mounting evidence recognizes structural variations (SVs) and repetitive DNA sequences as crucial players in shaping the existing grape phenotypic diversity at intra- and inter-species levels. To deepen our understanding on the abundance, diversity, and distribution of SVs and repetitive DNAs, including transposable elements (TEs) and tandemly repeated satellite DNA (satDNAs), we re-sequenced the genomes of the ancient grapes Aglianico and Falanghina. The analysis of large copy number variants (CNVs) detected candidate polymorphic genes that are involved in the enological features of these varieties. In a comparative analysis of Aglianico and Falanghina sequences with 21 publicly available genomes of cultivated grapes, we provided a genome-wide annotation of grape TEs at the lineage level. We disclosed that at least two main clusters of grape cultivars could be identified based on the TEs content. Multiple TEs families appeared either significantly enriched or depleted. In addition, in silico and cytological analyses provided evidence for a diverse chromosomal distribution of several satellite repeats between Aglianico, Falanghina, and other grapes. Overall, our data further improved our understanding of the intricate grape diversity held by two Italian traditional varieties, unveiling a pool of unique candidate genes never so far exploited in breeding for improved fruit quality.

Comparative Genomic Analysis of Asian Cultivated Rice and Its Wild Progenitor (Oryza rufipogon) Has Revealed Evolutionary Innovation of the Pentatricopeptide Repeat Gene Family through Gene Duplication

The pentatricopeptide repeat (PPR) gene family is one of the largest gene families in land plants. However, current knowledge about the evolution of the PPR gene family remains largely limited. In this study, we performed a comparative genomic analysis of the PPR gene family in O. sativa and its wild progenitor, O. rufipogon, and outlined a comprehensive landscape of gene duplications. Our findings suggest that the majority of PPR genes originated from dispersed duplications. Although segmental duplications have only expanded approximately 11.30% and 13.57% of the PPR gene families in the O. sativa and O. rufipogon genomes, we interestingly obtained evidence that segmental duplication promotes the structural diversity of PPR genes through incomplete gene duplications. In the O. sativa and O. rufipogon genomes, 10 (~33.33%) and 22 pairs of gene duplications (~45.83%) had non-PPR paralogous genes through incomplete gene duplication. Segmental duplications leading to incomplete gene duplications might result in the acquisition of domains, thus promoting functional innovation and structural diversification of PPR genes. This study offers a unique perspective on the evolution of PPR gene structures and underscores the potential role of segmental duplications in PPR gene structural diversity.

Vitis labrusca genome assembly reveals diversification between wild and cultivated grapevine genomes

Wild grapevines are important genetic resources in breeding programs to confer adaptive fitness traits and unique fruit characteristics, but the genetics underlying these traits, and their evolutionary origins, are largely unknown. To determine the factors that contributed to grapevine genome diversification, we performed comprehensive intragenomic and intergenomic analyses with three cultivated European (including the PN40024 reference genome) and two wild North American grapevine genomes, including our newly released Vitis labrusca genome. We found the heterozygosity of the cultivated grapevine genomes was twice as high as the wild grapevine genomes studied. Approximately 30% of V. labrusca and 48% of V. vinifera Chardonnay genes were heterozygous or hemizygous and a considerable number of collinear genes between Chardonnay and V. labrusca had different gene zygosity. Our study revealed evidence that supports gene gain-loss events in parental genomes resulted in the inheritance of hemizygous genes in the Chardonnay genome. Thousands of segmental duplications supplied source material for genome-specific genes, further driving diversification of the genomes studied. We found an enrichment of recently duplicated, adaptive genes in similar functional pathways, but differential retention of environment-specific adaptive genes within each genome. For example, large expansions of NLR genes were discovered in the two wild grapevine genomes studied. Our findings support variation in transposable elements contributed to unique traits in grapevines. Our work revealed gene zygosity, segmental duplications, gene gain-and-loss variations, and transposable element polymorphisms can be key driving forces for grapevine genome diversification.

Genome-wide analysis of the Populus trichocarpa laccase gene family and functional identification of PtrLAC23

INTRODUCTION

Biofuel is a kind of sustainable, renewable and environment friendly energy. Lignocellulose from the stems of woody plants is the main raw material for "second generation biofuels". Lignin content limits fermentation yield and is therefore a major obstacle in biofuel production. Plant laccase plays an important role in the final step of lignin formation, which provides a new strategy for us to obtain ideal biofuels by regulating the expression of laccase genes to directly gain the desired lignin content or change the composition of lignin.

METHODS

Multiple sequence alignment and phylogenetic analysis were used to classify PtrLAC genes; sequence features of PtrLACs were revealed by gene structure and motif composition analysis; gene duplication, interspecific collinearity and Ka/Ks analysis were conducted to identify ancient PtrLACs; expression levels of PtrLAC genes were measured by RNA-Seq data and qRT-PCR; domain analysis combine with cis-acting elements prediction together showed the potential function of PtrLACs. Furthermore, Alphafold2 was used to simulate laccase 3D structures, proLAC23::LAC23-eGFP transgenic Populus stem transects were applied to fluorescence observation.

RESULTS

A comprehensive analysis of the P. trichocarpa laccase gene (PtLAC) family was performed. Some ancient PtrLAC genes such as PtrLAC25, PtrLAC19 and PtrLAC41 were identified. Gene structure and distribution of conserved motifs clearly showed sequence characteristics of each PtrLAC. Combining published RNA-Seq data and qRT-PCR analysis, we revealed the expression pattern of PtrLAC gene family. Prediction results of cis-acting elements show that PtrLAC gene regulation was closely related to light. Through above analyses, we selected 5 laccases and used Alphafold2 to simulate protein 3D structures, results showed that PtrLAC23 may be closely related to the lignification. Fluorescence observation of proLAC23::LAC23-eGFP transgenic Populus stem transects and qRT-PCR results confirmed our hypothesis again.

DISCUSSION

In this study, we fully analyzed the Populus trichocarpa laccase gene family and identified key laccase genes related to lignification. These findings not only provide new insights into the characteristics and functions of Populus laccase, but also give a new understanding of the broad prospects of plant laccase in lignocellulosic biofuel production.

An efficient CRISPR-Cas9 enrichment sequencing strategy for characterizing complex and highly duplicated genomic regions. A case study in the Prunus salicina LG3-MYB10 genes cluster

BACKGROUND

Genome complexity is largely linked to diversification and crop innovation. Examples of regions with duplicated genes with relevant roles in agricultural traits are found in many crops. In both duplicated and non-duplicated genes, much of the variability in agronomic traits is caused by large as well as small and middle scale structural variants (SVs), which highlights the relevance of the identification and characterization of complex variability between genomes for plant breeding.

RESULTS

Here we improve and demonstrate the use of CRISPR-Cas9 enrichment combined with long-read sequencing technology to resolve the MYB10 region in the linkage group 3 (LG3) of Japanese plum (Prunus salicina). This region, which has a length from 90 to 271 kb according to the P. salicina genomes available, is associated with fruit color variability in Prunus species. We demonstrate the high complexity of this region, with homology levels between Japanese plum varieties comparable to those between Prunus species. We cleaved MYB10 genes in five plum varieties using the Cas9 enzyme guided by a pool of crRNAs. The barcoded fragments were then pooled and sequenced in a single MinION Oxford Nanopore Technologies (ONT) run, yielding 194 Mb of sequence. The enrichment was confirmed by aligning the long reads to the plum reference genomes, with a mean read on-target value of 4.5% and a depth per sample of 11.9x. From the alignment, 3261 SNPs and 287 SVs were called and phased. A de novo assembly was constructed for each variety, which also allowed detection, at the haplotype level, of the variability in this region.

CONCLUSIONS

CRISPR-Cas9 enrichment is a versatile and powerful tool for long-read targeted sequencing even on highly duplicated and/or polymorphic genomic regions, being especially useful when a reference genome is not available. Potential uses of this methodology as well as its limitations are further discussed.

Chromosome-scale genome assembly and population genomics provide insights into the adaptation, domestication, and flavonoid metabolism of Chinese plum

Globally, commercialized plum cultivars are mostly diploid Chinese plums (Prunus salicina Lindl.), also known as Japanese plums, and are one of the most abundant and variable fruit tree species. To advance Prunus genomic research, we present a chromosome-scale P. salicina genome assembly, constructed using an integrated strategy that combines Illumina, Oxford Nanopore, and high-throughput chromosome conformation capture (Hi-C) sequencing. The high-quality genome assembly consists of a 318.6-Mb sequence (contig N50 length of 2.3 Mb) with eight pseudo-chromosomes. The expansion of the P. salicina genome is led by recent segmental duplications and a long terminal repeat burst of approximately 0.2 Mya. This resulted in a significant expansion of gene families associated with flavonoid metabolism and plant resistance, which impacted fruit flavor and increased species adaptability. Population structure and domestication history suggest that Chinese plum may have originated from South China and provides a domestication route with accompanying genomic variations. Selection sweep and genetic diversity analysis enabled the identification of several critical genes associated with flowering time, stress tolerance, and flavonoid metabolism, demonstrating the essential roles of related pathways during domestication. Furthermore, we reconstructed and exploited flavonoid-anthocyanin metabolism using multi-omics analysis in Chinese plum and proposed a complete metabolic pathway. Collectively, our results will facilitate further candidate gene discovery for important agronomic traits in Chinese plum and provide insights into future functional genomic studies and DNA-informed breeding.

Seeing the Forest through the (Phylogenetic) Trees: Functional Characterisation of Grapevine Terpene Synthase (VviTPS) Paralogues and Orthologues

Gene families involved in specialised metabolism play a key role in a myriad of ecophysiological and biochemical functions. The Vitis vinifera sesquiterpene synthases represent the largest subfamily of grapevine terpene synthase (VviTPS) genes and are important volatile metabolites for wine flavour and aroma, as well as ecophysiological interactions. The functional characterisation of VviTPS genes is complicated by a reliance on a single reference genome that greatly underrepresents this large gene family, exacerbated by extensive duplications and paralogy. The recent release of multiple phased diploid grapevine genomes, as well as extensive whole-genome resequencing efforts, provide a wealth of new sequence information that can be utilised to overcome the limitations of the reference genome. A large cluster of sesquiterpene synthases, localised to chromosome 18, was explored by means of comparative sequence analyses using the publicly available grapevine reference genome, three PacBio phased diploid genomes and whole-genome resequencing data from multiple genotypes. Two genes, VviTPS04 and -10, were identified as putative paralogues and/or allelic variants. Subsequent gene isolation from multiple grapevine genotypes and characterisation by means of a heterologous in planta expression and volatile analysis resulted in the identification of genotype-specific structural variations and polymorphisms that impact the gene function. These results present novel insight into how grapevine domestication likely shaped the VviTPS landscape to result in genotype-specific functions.

Comprehensive Genome-Wide Exploration of C2H2 Zinc Finger Family in Grapevine (Vitis vinifera L.): Insights into the Roles in the Pollen Development Regulation

Some C2H2 zinc-finger proteins (ZFP) transcription factors are involved in the development of pollen in plants. In grapevine (Vitis vinifera L.), it has been suggested that abnormalities in pollen development lead to the phenomenon called parthenocarpy that occurs in some varieties of this cultivar. At present, a network involving several transcription factors types has been revealed and key roles have been assigned to members of the C2H2 zinc-finger proteins (ZFP) family in model plants. However, particularities of the regulatory mechanisms controlling pollen formation in grapevine remain unknown. In order to gain insight into the participation of ZFPs in grapevine gametophyte development, we performed a genome-wide identification and characterization of genes encoding ZFP (VviZFP family). A total of 98 genes were identified and renamed based on the gene distribution into grapevine genome. The analysis performed indicate significant changes throughout VviZFP genes evolution explained by high heterogeneity in sequence, length, number of ZF and presence of another conserved domains. Moreover, segmental duplication participated in the gene family expansion in grapevine. The VviZFPs were classified based on domain and phylogenetic analysis into three sets and different groups. Heat-map demonstrated differential and tissue-specific expression patterns of these genes and k-means clustering allowed to identify a group of putative orthologs to some ZFPs related to pollen development. In transgenic plants carrying the promVviZFP13::GUS and promVviZFP68::GUS constructs, GUS signals were detectable in the anther and mature pollen grains. Expression profiling of selected VviZFP genes showed differential expression pattern during flower development and provides a basis for deepening in the understanding of VviZFPs role on grapevine reproductive development.

Segmental duplications are hot spots of copy number variants affecting barley gene content

Copy number variants (CNVs) are pervasive in several animal and plant genomes and contribute to shaping genetic diversity. In barley, there is evidence that changes in gene copy number underlie important agronomic traits. The recently released reference sequence of barley represents a valuable genomic resource for unveiling the incidence of CNVs that affect gene content and for identifying sequence features associated with CNV formation. Using exome sequencing and read count data, we detected 16 605 deletions and duplications that affect barley gene content by surveying a diverse panel of 172 cultivars, 171 landraces, 22 wild relatives and other 32 uncategorized domesticated accessions. The quest for segmental duplications (SDs) in the reference sequence revealed many low-copy repeats, most of which overlap predicted coding sequences. Statistical analyses revealed that the incidence of CNVs increases significantly in SD-rich regions, indicating that these sequence elements act as hot spots for the formation of CNVs. The present study delivers a comprehensive genome-wide study of CNVs affecting barley gene content and implicates SDs in the molecular mechanisms that lead to the formation of this class of CNVs.

The Soybean Laccase Gene Family: Evolution and Possible Roles in Plant Defense and Stem Strength Selection

Laccase is a widely used industrial oxidase for food processing, dye synthesis, paper making, and pollution remediation. At present, laccases used by industries come mainly from fungi. Plants contain numerous genes encoding laccase enzymes that show properties which are distinct from that of the fungal laccases. These plant-specific laccases may have better potential for industrial purposes. The aim of this work was to conduct a genome-wide search for the soybean laccase genes and analyze their characteristics and specific functions. A total of 93 putative laccase genes (GmLac) were identified from the soybean genome. All 93 GmLac enzymes contain three typical Cu-oxidase domains, and they were classified into five groups based on phylogenetic analysis. Although adjacent members on the tree showed highly similar exon/intron organization and motif composition, there were differences among the members within a class for both conserved and differentiated functions. Based on the expression patterns, some members of laccase were expressed in specific tissues/organs, while some exhibited a constitutive expression pattern. Analysis of the transcriptome revealed that some laccase genes might be involved in providing resistance to oomycetes. Analysis of the selective pressures acting on the laccase gene family in the process of soybean domestication revealed that 10 genes could have been under artificial selection during the domestication process. Four of these genes may have contributed to the transition of the soft and thin stem of wild soybean species into strong, thick, and erect stems of the cultivated soybean species. Our study provides a foundation for future functional studies of the soybean laccase gene family.

Genome Wide Identification, Molecular Characterization, and Gene Expression Analyses of Grapevine NHX Antiporters Suggest Their Involvement in Growth, Ripening, Seed Dormancy, and Stress Response

Plant NHX antiporters are critical for cellular pH, Na⁺, and K⁺ homeostasis and salt tolerance. Even though their genomic and functional studies have been conducted in many species, the grapevine NHX family has not been described yet. Our work highlights the presence of six VvNHX genes whose phylogenetic analysis revealed their classification in two distinct groups: group I vacuolar (VvNHX1-5) and group II endosomal (VvNHX6). Several cis-acting regulatory elements related to tissue-specific expression, transcription factor binding, abiotic/biotic stresses response, and light regulation elements were identified in their promoter. Expression profile analyses of VvNHX genes showed variable transcription within organs and tissues with diverse patterns according to biochemical, environmental, and biotic treatments. All VvNHXs are involved in berry growth, except VvNHX5 that seems to be rather implicated in seed maturation. VvNHX4 would be more involved in floral development, while VvNHX2 and 3 display redundant roles. QPCR expression analyses of VvNHX1 showed its induction by NaCl and KNO₃ treatments, whereas VvNHX6 was induced by ABA application and strongly repressed by PEG treatment. VvNHX1 plays a crucial role in a bunch of grape developmental steps and adaptation responses through mechanisms of phyto-hormonal signaling. Overall, VvNHX family members could be valuable candidate genes for grapevine improvement.

Transcriptomic and genomic structural variation analyses on grape cultivars reveal new insights into the genotype-dependent responses to water stress

Grapevine (Vitis vinifera L.) is importantly cultivated worldwide for table grape and wine production. Its cultivation requires irrigation supply, especially in arid and semiarid areas. Water deficiency can affect berry and wine quality mostly depending on the extent of plant perceived stress, which is a cultivar-specific trait. We tested the physiological and molecular responses to water deficiency of two table grape cultivars, Italia and Autumn royal, and we highlighted their different adaptation. Microarray analyses revealed that Autumn royal reacts involving only 29 genes, related to plant stress response and ABA/hormone signal transduction, to modulate the response to water deficit. Instead, cultivar Italia orchestrates a very broad response (we found 1037 differentially expressed genes) that modifies the cell wall organization, carbohydrate metabolism, response to reactive oxygen species, hormones and osmotic stress. For the first time, we integrated transcriptomic data with cultivar-specific genomics and found that ABA-perception and -signalling are key factors mediating the varietal-specific behaviour of the early response to drought. We were thus able to isolate candidate genes for the genotype-dependent response to drought. These insights will allow the identification of reliable plant stress indicators and the definition of sustainable cultivar-specific protocols for water management.

Identification of jumonjiC domain containing gene family among the Oryza species and their expression analysis in FL478, a salt tolerant rice genotype

The jumonji (JMJ)-C domain containing proteins belong to histone demethylases family with the ability to demethylate the tri-methylated histone residues. They act as chromatin regulators to regulate many physiological functions in plants. The present study deals with the characterization of JMJ-C gene family members in wild as well as cultivated rice species and their expression analysis in salt tolerant rice genotype, FL478. The genome wide study identified 151 members belonging to JMJ-C gene family in 11 different Oryza species. We also studied their structure, genomic location, gene duplication events, phylogenetic relationship, in silico expression analysis and identified cis elements in their promoters. We also found a few JMJ-C gene family members in rice which underwent duplication before the whole genome duplication event of the rice. The qRT-PCR based expression profiling revealed that out of the total 15 rice JMJ-C members, two were highly expressed in the flag leaf stage of FL478 under salt treatment. These two candidate JMJ-C members were also found to render salinity tolerance when over-expressed in yeast cells. Thus, the present study helps in further structural as well as functional characterization of JMJ-C genes under salinity stress in Oryza species.

Revised Phylogeny of the Cellulose Synthase Gene Superfamily: Insights into Cell Wall Evolution

Cell walls are crucial for the integrity and function of all land plants and are of central importance in human health, livestock production, and as a source of renewable bioenergy. Many enzymes that mediate the biosynthesis of cell wall polysaccharides are encoded by members of the large cellulose synthase (CesA) gene superfamily. Here, we analyzed 29 sequenced genomes and 17 transcriptomes to revise the phylogeny of the CesA gene superfamily in angiosperms. Our results identify ancestral gene clusters that predate the monocot-eudicot divergence and reveal several novel evolutionary observations, including the expansion of the Poaceae-specific cellulose synthase-like CslF family to the graminids and restiids and the characterization of a previously unreported eudicot lineage, CslM, that forms a reciprocally monophyletic eudicot-monocot grouping with the CslJ clade. The CslM lineage is widely distributed in eudicots, and the CslJ clade, which was thought previously to be restricted to the Poales, is widely distributed in monocots. Our analyses show that some members of the CslJ lineage, but not the newly identified CslM genes, are capable of directing (1,3;1,4)-β-glucan biosynthesis, which, contrary to current dogma, is not restricted to Poaceae.

Analysis of Copy Number Variation in the Abp Gene Regions of Two House Mouse Subspecies Suggests Divergence during the Gene Family Expansions

The Androgen-binding protein (Abp) gene region of the mouse genome contains 64 genes, some encoding pheromones that influence assortative mating between mice from different subspecies. Using CNVnator and quantitative PCR, we explored copy number variation in this gene family in natural populations of Mus musculus domesticus (Mmd) and Mus musculus musculus (Mmm), two subspecies of house mice that form a narrow hybrid zone in Central Europe. We found that copy number variation in the center of the Abp gene region is very common in wild Mmd, primarily representing the presence/absence of the final duplications described for the mouse genome. Clustering of Mmd individuals based on this variation did not reflect their geographical origin, suggesting no population divergence in the Abp gene cluster. However, copy number variation patterns differ substantially between Mmd and other mouse taxa. Large blocks of Abp genes are absent in Mmm, Mus musculus castaneus and an outgroup, Mus spretus, although with differences in variation and breakpoint locations. Our analysis calls into question the reliance on a reference genome for interpreting the detailed organization of genes in taxa more distant from the Mmd reference genome. The polymorphic nature of the gene family expansion in all four taxa suggests that the number of Abp genes, especially in the central gene region, is not critical to the survival and reproduction of the mouse. However, Abp haplotypes of variable length may serve as a source of raw genetic material for new signals influencing reproductive communication and thus speciation of mice.

Sequence Polymorphisms and Structural Variations among Four Grapevine (Vitis vinifera L.) Cultivars Representing Sardinian Agriculture

The genetic diversity among grapevine (Vitis vinifera L.) cultivars that underlies differences in agronomic performance and wine quality reflects the accumulation of single nucleotide polymorphisms (SNPs) and small indels as well as larger genomic variations. A combination of high throughput sequencing and mapping against the grapevine reference genome allows the creation of comprehensive sequence variation maps. We used next generation sequencing and bioinformatics to generate an inventory of SNPs and small indels in four widely cultivated Sardinian grape cultivars (Bovale sardo, Cannonau, Carignano and Vermentino). More than 3,200,000 SNPs were identified with high statistical confidence. Some of the SNPs caused the appearance of premature stop codons and thus identified putative pseudogenes. The analysis of SNP distribution along chromosomes led to the identification of large genomic regions with uninterrupted series of homozygous SNPs. We used a digital comparative genomic hybridization approach to identify 6526 genomic regions with significant differences in copy number among the four cultivars compared to the reference sequence, including 81 regions shared between all four cultivars and 4953 specific to single cultivars (representing 1.2 and 75.9% of total copy number variation, respectively). Reads mapping at a distance that was not compatible with the insert size were used to identify a dataset of putative large deletions with cultivar Cannonau revealing the highest number. The analysis of genes mapping to these regions provided a list of candidates that may explain some of the phenotypic differences among the Bovale sardo, Cannonau, Carignano and Vermentino cultivars.

Genome-wide characterisation and expression profile of the grapevine ATL ubiquitin ligase family reveal biotic and abiotic stress-responsive and development-related members

The Arabidopsis Tóxicos en Levadura (ATL) protein family is a class of E3 ubiquitin ligases with a characteristic RING-H2 Zn-finger structure that mediates diverse physiological processes and stress responses in plants. We carried out a genome-wide survey of grapevine (Vitis vinifera L.) ATL genes and retrieved 96 sequences containing the canonical ATL RING-H2 domain. We analysed their genomic organisation, gene structure and evolution, protein domains and phylogenetic relationships. Clustering revealed several clades, as already reported in Arabidopsis thaliana and rice (Oryza sativa), with an expanded subgroup of grapevine-specific genes. Most of the grapevine ATL genes lacked introns and were scattered among the 19 chromosomes, with a high level of duplication retention. Expression profiling revealed that some ATL genes are expressed specifically during early or late development and may participate in the juvenile to mature plant transition, whereas others may play a role in pathogen and/or abiotic stress responses, making them key candidates for further functional analysis. Our data offer the first genome-wide overview and annotation of the grapevine ATL family, and provide a basis for investigating the roles of specific family members in grapevine physiology and stress responses, as well as potential biotechnological applications.

Inter-varietal structural variation in grapevine genomes

Grapevine (Vitis vinifera L.) is one of the world's most important crop plants, which is of large economic value for fruit and wine production. There is much interest in identifying genomic variations and their functional effects on inter-varietal, phenotypic differences. Using an approach developed for the analysis of human and mammalian genomes, which combines high-throughput sequencing, array comparative genomic hybridization, fluorescent in situ hybridization and quantitative PCR, we created an inter-varietal atlas of structural variations and single nucleotide variants (SNVs) for the grapevine genome analyzing four economically and genetically relevant table grapevine varieties. We found 4.8 million SNVs and detected 8% of the grapevine genome to be affected by genomic variations. We identified more than 700 copy number variation (CNV) regions and more than 2000 genes subjected to CNV as potential candidates for phenotypic differences between varieties.

Odintifier--A computational method for identifying insertions of organellar origin from modern and ancient high-throughput sequencing data based on haplotype phasing

Background

Cellular organelles with genomes of their own (e.g. plastids and mitochondria) can pass genetic sequences to other organellar genomes within the cell in many species across the eukaryote phylogeny. The extent of the occurrence of these organellar-derived inserted sequences (odins) is still unknown, but if not accounted for in genomic and phylogenetic studies, they can be a source of error. However, if correctly identified, these inserted sequences can be used for evolutionary and comparative genomic studies. Although such insertions can be detected using various laboratory and bioinformatic strategies, there is currently no straightforward way to apply them as a standard organellar genome assembly on next-generation sequencing data. Furthermore, most current methods for identification of such insertions are unsuitable for use on non-model organisms or ancient DNA datasets.

Results

We present a bioinformatic method that uses phasing algorithms to reconstruct both source and inserted organelle sequences. The method was tested in different shotgun and organellar-enriched DNA high-throughput sequencing (HTS) datasets from ancient and modern samples. Specifically, we used datasets from lions (Panthera leo ssp. and Panthera leo leo) to characterize insertions from mitochondrial origin, and from common grapevine (Vitis vinifera) and bugle (Ajuga reptans) to characterize insertions derived from plastid genomes. Comparison of the results against other available organelle genome assembly methods demonstrated that our new method provides an improvement in the sequence assembly.

Conclusion

Using datasets from a wide range of species and different levels of complexity we showed that our novel bioinformatic method based on phasing algorithms can be used to achieve the next two goals: i) reference-guided assembly of chloroplast/mitochondrial genomes from HTS data and ii) identification and simultaneous assembly of odins. This method represents the first application of haplotype phasing for automatic detection of odins and reference-based organellar genome assembly.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0682-1) contains supplementary material, which is available to authorized users.

Precision breeding of grapevine (Vitis vinifera L.) for improved traits

This review provides an overview of recent technological advancements that enable precision breeding to genetically improve elite cultivars of grapevine (Vitis vinifera L.). Precision breeding, previously termed "cisgenic" or "intragenic" genetic improvement, necessitates a better understanding and use of genomic resources now becoming accessible. Although it is now a relatively simple task to identify genetic elements and genes from numerous "omics" databases, the control of major agronomic and enological traits often involves the currently unknown participation of many genes and regulatory machineries. In addition, genetic evolution has left numerous vestigial genes and sequences without tangible functions. Thus, it is critical to functionally test each of these genetic entities to determine their real-world functionality or contribution to trait attributes. Toward this goal, several diverse techniques now are in place, including cell culture systems to allow efficient plant regeneration, advanced gene insertion techniques, and, very recently, resources for genomic analyses. Currently, these techniques are being used for high-throughput expression analysis of a wide range of grapevine-derived promoters and disease-related genes. It is envisioned that future research efforts will be extended to the study of promoters and genes functioning to enhance other important traits, such as fruit quality and vigor.

A UDP-Glucose:Monoterpenol Glucosyltransferase Adds to the Chemical Diversity of the Grapevine Metabolome

Terpenoids represent one of the major classes of natural products and serve different biological functions. In grape (Vitis vinifera), a large fraction of these compounds is present as nonvolatile terpene glycosides. We have extracted putative glycosyltransferase (GT) sequences from the grape genome database that show similarity to Arabidopsis (Arabidopsis thaliana) GTs whose encoded proteins glucosylate a diversity of terpenes. Spatial and temporal expression levels of the potential VvGT genes were determined in five different grapevine varieties. Heterologous expression and biochemical assays of candidate genes led to the identification of a UDP-glucose:monoterpenol β-d-glucosyltransferase (VvGT7). The VvGT7 gene was expressed in various tissues in accordance with monoterpenyl glucoside accumulation in grape cultivars. Twelve allelic VvGT7 genes were isolated from five cultivars, and their encoded proteins were biochemically analyzed. They varied in substrate preference and catalytic activity. Three amino acids, which corresponded to none of the determinants previously identified for other plant GTs, were found to be important for enzymatic catalysis. Site-specific mutagenesis along with the analysis of allelic proteins also revealed amino acids that impact catalytic activity and substrate tolerance. These results demonstrate that VvGT7 may contribute to the production of geranyl and neryl glucoside during grape ripening.

The role of retrotransposons in gene family expansions: insights from the mouse Abp gene family

Background

Retrotransposons have been suggested to provide a substrate for non-allelic homologous recombination (NAHR) and thereby promote gene family expansion. Their precise role, however, is controversial. Here we ask whether retrotransposons contributed to the recent expansions of the Androgen-binding protein (Abp) gene families that occurred independently in the mouse and rat genomes.

Results

Using dot plot analysis, we found that the most recent duplication in the Abp region of the mouse genome is flanked by L1Md_T elements. Analysis of the sequence of these elements revealed breakpoints that are the relicts of the recombination that caused the duplication, confirming that the duplication arose as a result of NAHR using L1 elements as substrates. L1 and ERVII retrotransposons are considerably denser in the Abp regions than in one Mb flanking regions, while other repeat types are depleted in the Abp regions compared to flanking regions. L1 retrotransposons preferentially accumulated in the Abp gene regions after lineage separation and roughly followed the pattern of Abp gene expansion. By contrast, the proportion of shared vs. lineage-specific ERVII repeats in the Abp region resembles the rest of the genome.

Conclusions

We confirmed the role of L1 repeats in Abp gene duplication with the identification of recombinant L1Md_T elements at the edges of the most recent mouse Abp gene duplication. High densities of L1 and ERVII repeats were found in the Abp gene region with abrupt transitions at the region boundaries, suggesting that their higher densities are tightly associated with Abp gene duplication. We observed that the major accumulation of L1 elements occurred after the split of the mouse and rat lineages and that there is a striking overlap between the timing of L1 accumulation and expansion of the Abp gene family in the mouse genome. Establishing a link between the accumulation of L1 elements and the expansion of the Abp gene family and identification of an NAHR-related breakpoint in the most recent duplication are the main contributions of our study.

Comparative genomic and transcriptomic analysis of tandemly and segmentally duplicated genes in rice

Tandem and segmental duplications significantly contribute to gene family expansion and genome evolution. Genome-wide identification of tandem and segmental genes has been analyzed before in several plant genomes. However, comparative studies in functional bias, expression divergence and their roles in species domestication are still lacking. We have carried out a genome-wide identification and comparative analysis of tandem and segmental genes in the rice genome. A total of 3,646 and 3,633 pairs of tandem and segmental genes, respectively, were identified in the genome. They made up around 30% of total annotated rice genes (excluding transposon-coding genes). Both tandem and segmental duplicates showed different physical locations and exhibited a biased subset of functions. These two types of duplicated genes were also under different functional constrains as shown by nonsynonymous substitutions per site (Ka) and synonymous substitutions per site (Ks) analysis. They are also differently regulated depending on the tissues and abiotic and biotic stresses based on transcriptomics data. The expression divergence might be related to promoter differentiation and DNA methylation status after tandem or segmental duplications. Both tandem and segmental duplications differ in their contribution to genetic novelty but evidence suggests that they play their role in species domestication and genome evolution.

In with the old, in with the new: the promiscuity of the duplication process engenders diverse pathways for novel gene creation

The gene duplication process has exhibited far greater promiscuity in the creation of paralogs with novel exon-intron structures than anticipated even by Ohno. In this paper I explore the history of the field, from the neo-Darwinian synthesis through Ohno's formulation of the canonical model for the evolution of gene duplicates and culminating in the present genomic era. I delineate the major tenets of Ohno's model and discuss its failure to encapsulate the full complexity of the duplication process as revealed in the era of genomics. I discuss the diverse classes of paralogs originating from both DNA- and RNA-mediated duplication events and their evolutionary potential for assuming radically altered functions, as well as the degree to which they can function unconstrained from the pressure of gene conversion. Lastly, I explore theoretical population-genetic considerations of how the effective population size (N(e)) of a species may influence the probability of emergence of genes with radically altered functions.

The grapevine expression atlas reveals a deep transcriptome shift driving the entire plant into a maturation program

We developed a genome-wide transcriptomic atlas of grapevine (Vitis vinifera) based on 54 samples representing green and woody tissues and organs at different developmental stages as well as specialized tissues such as pollen and senescent leaves. Together, these samples expressed ∼91% of the predicted grapevine genes. Pollen and senescent leaves had unique transcriptomes reflecting their specialized functions and physiological status. However, microarray and RNA-seq analysis grouped all the other samples into two major classes based on maturity rather than organ identity, namely, the vegetative/green and mature/woody categories. This division represents a fundamental transcriptomic reprogramming during the maturation process and was highlighted by three statistical approaches identifying the transcriptional relationships among samples (correlation analysis), putative biomarkers (O2PLS-DA approach), and sets of strongly and consistently expressed genes that define groups (topics) of similar samples (biclustering analysis). Gene coexpression analysis indicated that the mature/woody developmental program results from the reiterative coactivation of pathways that are largely inactive in vegetative/green tissues, often involving the coregulation of clusters of neighboring genes and global regulation based on codon preference. This global transcriptomic reprogramming during maturation has not been observed in herbaceous annual species and may be a defining characteristic of perennial woody plants.

Genome-wide analysis of the grapevine stilbene synthase multigenic family: genomic organization and expression profiles upon biotic and abiotic stresses

BACKGROUND

Plant stilbenes are a small group of phenylpropanoids, which have been detected in at least 72 unrelated plant species and accumulate in response to biotic and abiotic stresses such as infection, wounding, UV-C exposure and treatment with chemicals. Stilbenes are formed via the phenylalanine/polymalonate-route, the last step of which is catalyzed by the enzyme stilbene synthase (STS), a type III polyketide synthase (PKS). Stilbene synthases are closely related to chalcone synthases (CHS), the key enzymes of the flavonoid pathway, as illustrated by the fact that both enzymes share the same substrates. To date, STSs have been cloned from peanut, pine, sorghum and grapevine, the only stilbene-producing fruiting-plant for which the entire genome has been sequenced. Apart from sorghum, STS genes appear to exist as a family of closely related genes in these other plant species.

RESULTS

In this study a complete characterization of the STS multigenic family in grapevine has been performed, commencing with the identification, annotation and phylogenetic analysis of all members and integration of this information with a comprehensive set of gene expression analyses including healthy tissues at differential developmental stages and in leaves exposed to both biotic (downy mildew infection) and abiotic (wounding and UV-C exposure) stresses. At least thirty-three full length sequences encoding VvSTS genes were identified, which, based on predicted amino acid sequences, cluster in 3 principal groups designated A, B and C. The majority of VvSTS genes cluster in groups B and C and are located on chr16 whereas the few gene family members in group A are found on chr10. Microarray and mRNA-seq expression analyses revealed different patterns of transcript accumulation between the different groups of VvSTS family members and between VvSTSs and VvCHSs. Indeed, under certain conditions the transcriptional response of VvSTS and VvCHS genes appears to be diametrically opposed suggesting that flow of carbon between these two competing metabolic pathways is tightly regulated at the transcriptional level.

CONCLUSIONS

This study represents an overview of the expression pattern of each member of the STS gene family in grapevine under both constitutive and stress-induced conditions. The results strongly indicate the existence of a transcriptional subfunctionalization amongst VvSTSs and provide the foundation for further functional investigations about the role and evolution of this large gene family. Moreover, it represents the first study to clearly show the differential regulation of VvCHS and VvSTS genes, suggesting the involvement of transcription factors (TFs) in both the activation and repression of these genes.