Aggregated gene co-expression networks predict transcription factor regulatory landscapes in grapevine.

Gene co-expression networks (GCNs) have not been extensively studied in non-model plants. However, the rapid accumulation of transcriptome datasets in certain species represents an opportunity to explore underutilized network aggregation approaches. In fact, aggregated GCNs (aggGCNs) highlight robust co-expression interactions and improve functional connectivity. We applied and evaluated two different aggregation methods on public grapevine RNA-Seq datasets from three different tissues (leaf, berry, and 'all organs'). Our results show that co-occurrence-based aggregation generally yielded the best-performing networks. We applied aggGCNs to study several transcription factor gene families, showing their capacity for detecting both already-described and novel regulatory relationships between R2R3-MYBs, bHLH/MYC, and multiple specialized metabolic pathways. Specifically, transcription factor gene- and pathway-centered network analyses successfully ascertained the previously established role of VviMYBPA1 in controlling the accumulation of proanthocyanidins while providing insights into its novel role as a regulator of p-coumaroyl-CoA biosynthesis as well as the shikimate and aromatic amino acid pathways. This network was validated using DNA affinity purification sequencing data, demonstrating that co-expression networks of transcriptional activators can serve as a proxy of gene regulatory networks. This study presents an open repository to reproduce networks in other crops and a GCN application within the Vitviz platform, a user-friendly tool for exploring co-expression relationships.

An improved reference of the grapevine genome reasserts the origin of the PN40024 highly homozygous genotype.

The genome sequence of the diploid and highly homozygous Vitis vinifera genotype PN40024 serves as the reference for many grapevine studies. Despite several improvements to the PN40024 genome assembly, its current version PN12X.v2 is quite fragmented and only represents the haploid state of the genome with mixed haplotypes. In fact, being nearly homozygous, this genome contains several heterozygous regions that are yet to be resolved. Taking the opportunity of improvements that long-read sequencing technologies offer to fully discriminate haplotype sequences, an improved version of the reference, called PN40024.v4, was generated. Through incorporating long genomic sequencing reads to the assembly, the continuity of the 12X.v2 scaffolds was highly increased with a total number decreasing from 2,059 to 640 and a reduction in N bases of 88%. Additionally, the full alternative haplotype sequence was built for the first time, the chromosome anchoring was improved and the number of unplaced scaffolds was reduced by half. To obtain a high-quality gene annotation that outperforms previous versions, a liftover approach was complemented with an optimized annotation workflow for Vitis. Integration of the gene reference catalogue and its manual curation have also assisted in improving the annotation, while defining the most reliable estimation of 35,230 genes to date. Finally, we demonstrated that PN40024 resulted from 9 selfings of cv. "Helfensteiner" (cross of cv. "Pinot noir" and "Schiava grossa") instead of a single "Pinot noir". These advances will help maintain the PN40024 genome as a gold-standard reference, also contributing toward the eventual elaboration of the grapevine pangenome.

Direct regulation of shikimate, early phenylpropanoid, and stilbenoid pathways by Subgroup 2 R2R3-MYBs in grapevine.

The stilbenoid pathway is responsible for the production of resveratrol in grapevine (Vitis vinifera L.). A few transcription factors (TFs) have been identified as regulators of this pathway but the extent of this control has not been deeply studied. Here we show how DNA affinity purification sequencing (DAP-Seq) allows for the genome-wide TF-binding site interrogation in grape. We obtained 5190 and 4443 binding events assigned to 4041 and 3626 genes for MYB14 and MYB15, respectively (approximately 40% of peaks located within −10 kb of transcription start sites). DAP-Seq of MYB14/MYB15 was combined with aggregate gene co-expression networks (GCNs) built from more than 1400 transcriptomic datasets from leaves, fruits, and flowers to narrow down bound genes to a set of high confidence targets. The analysis of MYB14, MYB15, and MYB13, a third uncharacterized member of Subgroup 2 (S2), showed that in addition to the few previously known stilbene synthase (STS) targets, these regulators bind to 30 of 47 STS family genes. Moreover, all three MYBs bind to several PAL, C4H, and 4CL genes, in addition to shikimate pathway genes, the WRKY03 stilbenoid co-regulator and resveratrol-modifying gene candidates among which ROMT2-3 were validated enzymatically. A high proportion of DAP-Seq bound genes were induced in the activated transcriptomes of transient MYB15-overexpressing grapevine leaves, validating our methodological approach for delimiting TF targets. Overall, Subgroup 2 R2R3-MYBs appear to play a key role in binding and directly regulating several primary and secondary metabolic steps leading to an increased flux towards stilbenoid production. The integration of DAP-Seq and reciprocal GCNs offers a rapid framework for gene function characterization using genome-wide approaches in the context of non-model plant species and stands up as a valid first approach for identifying gene regulatory networks of specialized metabolism.

Vitis OneGenE: A Causality-Based Approach to Generate Gene Networks in Vitis vinifera Sheds Light on the Laccase and Dirigent Gene Families.

The abundance of transcriptomic data and the development of causal inference methods have paved the way for gene network analyses in grapevine. Vitis OneGenE is a transcriptomic data mining tool that finds direct correlations between genes, thus producing association networks. As a proof of concept, the stilbene synthase gene regulatory network obtained with OneGenE has been compared with published co-expression analysis and experimental data, including cistrome data for MYB stilbenoid regulators. As a case study, the two secondary metabolism pathways of stilbenoids and lignin synthesis were explored. Several isoforms of laccase, peroxidase, and dirigent protein genes, putatively involved in the final oxidative oligomerization steps, were identified as specifically belonging to either one of these pathways. Manual curation of the predicted sequences exploiting the last available genome assembly, and the integration of phylogenetic and OneGenE analyses, identified a group of laccases exclusively present in grapevine and related to stilbenoids. Here we show how network analysis by OneGenE can accelerate knowledge discovery by suggesting new candidates for functional characterization and application in breeding programs.

The Grape Gene Reference Catalogue as a Standard Resource for Gene Selection and Genetic Improvement.

Effective crop improvement, whether through selective breeding or biotech strategies, is largely dependent on the cumulative knowledge of a species' pangenome and its containing genes. Acquiring this knowledge is specially challenging in grapevine, one of the oldest fruit crops grown worldwide, which is known to have more than 30,000 genes. Well-established research communities studying model organisms have created and maintained, through public and private funds, a diverse range of online tools and databases serving as repositories of genomes and gene function data. The lack of such resources for the non-model, but economically important, Vitis vinifera species has driven the need for a standardised collection of genes within the grapevine community. In an effort led by the Integrape COST Action CA17111, we have recently developed the first grape gene reference catalogue, where genes are ascribed to functional data, including their accession identifiers from different genome-annotation versions (https://integrape.eu/resources/genes-genomes/). We present and discuss this gene repository together with a validation-level scheme based on varied supporting evidence found in current literature. The catalogue structure and online submission form provided permits community curation. Finally, we present the Gene Cards tool, developed within the Vitis Visualization (VitViz) platform, to visualize the data collected in the catalogue and link gene function with tissue-specific expression derived from public transcriptomic data. This perspective article aims to present these resources to the community as well as highlight their potential use, in particular for plant-breeding applications.

Targeting female flight for genetic control of mosquitoes.

Aedes aegypti Act4 is a paralog of the Drosophila melanogaster indirect flight muscle actin gene Act88F. Act88F has been shown to be haploinsufficient for flight in both males and females (amorphic mutants are dominant). Whereas Act88F is expressed in indirect flight muscles of both males and females, expression of Act4 is substantially female-specific. We therefore used CRISPR/Cas9 and homology directed repair to examine the phenotype of Act4 mutants in two Culicine mosquitoes, Aedes aegypti and Culex quinquefasciatus. A screen for dominant female-flightless mutants in Cx. quinquefasciatus identified one such mutant associated with a six base pair deletion in the CxAct4 coding region. A similar screen in Ae. aegypti identified no dominant mutants. Disruption of the AeAct4 gene by homology-dependent insertion of a fluorescent protein marker cassette gave a recessive female-flightless phenotype in Ae. aegypti. Reproducing the six-base deletion from Cx. quinquefasciatus in Ae. aegypti using oligo-directed mutagenesis generated dominant female-flightless mutants and identified additional dominant female-flightless mutants with other in-frame insertions or deletions. Our data indicate that loss of function mutations in the AeAct4 gene are recessive but that short in-frame deletions produce dominant-negative versions of the AeAct4 protein that interfere with flight muscle function. This makes Act4 an interesting candidate for genetic control methods, particularly population-suppression gene drives targeting female viability/fertility.