Trichomes and unique gene expression confer insect herbivory resistance in Vitis labrusca grapevines.

BACKGROUND: Grapevine (Vitis) is one of the world's most valuable fruit crops, but insect herbivory can decrease yields. Understanding insect herbivory resistance is critical to mitigating these losses. Vitis labrusca, a wild North American grapevine species, has been leveraged in breeding programs to generate hybrid grapevines with enhanced abiotic and biotic stress resistance, rendering it a valuable genetic resource for sustainable viticulture. This study assessed the resistance of V. labrusca acc. 'GREM4' and Vitis vinifera cv. 'PN40024' grapevines to Popillia japonica (Japanese beetle) herbivory and identified morphological and genetic adaptations underlying this putative resistance. RESULTS: 'GREM4' displayed greater resistance to beetle herbivory compared to 'PN40024' in both choice and no-choice herbivory assays spanning periods of 30 min to 19 h. 'GREM4' had significantly higher average leaf trichome densities than 'PN40024' and beetles preferred to feed on the side of leaves with fewer trichomes. When leaves from each species that specifically did not differ in trichome densities were fed on by beetles, significantly less leaf area was damaged in 'GREM4' (3.29mm2) compared to 'PN40024' (9.80mm2), suggesting additional factors beyond trichomes contributed to insect herbivory resistance in 'GREM4'. Comparative transcriptomic analyses revealed 'GREM4' exhibited greater constitutive (0 h) expression of defense response and secondary metabolite biosynthesis genes compared to 'PN40024', indicative of heightened constitutive defenses. Upon herbivory, 'GREM4' displayed a greater number of differentially expressed genes (690) compared to 'PN40024' (502), suggesting a broader response. Genes up-regulated in 'GREM4' were enriched in terpene biosynthesis, flavonoid biosynthesis, phytohormone signaling, and disease defense-related functions, likely contributing to heighted insect herbivory defense, while genes differentially expressed in 'PN40024' under herbivory were enriched in xyloglucan, cell wall formation, and calcium ion binding. The majority of genes implicated in insect herbivory defense were orthologs with specific expression patterns in 'GREM4' and 'PN40024', but some paralogous and genome-specific genes also likely contributed to conferring resistance. CONCLUSIONS: Our findings suggest that 'GREM4' insect herbivory resistance was attributed to a combination of factors, including trichomes and unique constitutive and inducible expression of genes implicated in terpene, flavonoid, and phenylpropanoid biosynthesis, as well as pathogen defense.

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.

Evolution of weedy giant ragweed (Ambrosia trifida): Multiple origins and gene expression variability facilitates weediness.

Agricultural weeds may originate from wild populations, but the origination patterns and genetics underlying this transition remain largely unknown. Analysis of weedy-wild paired populations from independent locations may provide evidence to identify key genetic variation contributing to this adaptive shift. We performed genetic variation and expression analyses on transcriptome data from 67 giant ragweed samples collected from different locations in Ohio, Iowa, and Minnesota and found geographically separated weedy populations likely originated independently from their adjacent wild populations, but subsequent spreading of weedy populations also occurred locally. By using eight closely related weedy-wild paired populations, we identified thousands of unique transcripts in weedy populations that reflect shared or specific functions corresponding, respectively, to both convergently evolved and population-specific weediness processes. In addition, differential expression of specific groups of genes was detected between weedy and wild giant ragweed populations using gene expression diversity and gene co-expression network analyses. Our study suggests an integrated route of weedy giant ragweed origination, consisting of independent origination combined with the subsequent spreading of certain weedy populations, and provides several lines of evidence to support the hypothesis that gene expression variability plays a key role in the evolution of weedy species.

Rapid evolution of protein diversity by de novo origination in Oryza.

New protein-coding genes that arise de novo from non-coding DNA sequences contribute to protein diversity. However, de novo gene origination is challenging to study as it requires high-quality reference genomes for closely related species, evidence for ancestral non-coding sequences, and transcription and translation of the new genes. High-quality genomes of 13 closely related Oryza species provide unprecedented opportunities to understand de novo origination events. Here, we identify a large number of young de novo genes with discernible recent ancestral non-coding sequences and evidence of translation. Using pipelines examining the synteny relationship between genomes and reciprocal-best whole-genome alignments, we detected at least 175 de novo open reading frames in the focal species O. sativa subspecies japonica, which were all detected in RNA sequencing-based transcriptomes. Mass spectrometry-based targeted proteomics and ribosomal profiling show translational evidence for 57% of the de novo genes. In recent divergence of Oryza, an average of 51.5 de novo genes per million years were generated and retained. We observed evolutionary patterns in which excess indels and early transcription were favoured in origination with a stepwise formation of gene structure. These data reveal that de novo genes contribute to the rapid evolution of protein diversity under positive selection.

Publisher Correction: Genomes of 13 domesticated and wild rice relatives highlight genetic conservation, turnover and innovation across the genus Oryza.

This article was not made open access when initially published online, which was corrected before print publication. In addition, ORCID links were missing for 12 authors and have been added to the HTML and PDF versions of the article.

Genomes of 13 domesticated and wild rice relatives highlight genetic conservation, turnover and innovation across the genus Oryza.

The genus Oryza is a model system for the study of molecular evolution over time scales ranging from a few thousand to 15 million years. Using 13 reference genomes spanning the Oryza species tree, we show that despite few large-scale chromosomal rearrangements rapid species diversification is mirrored by lineage-specific emergence and turnover of many novel elements, including transposons, and potential new coding and noncoding genes. Our study resolves controversial areas of the Oryza phylogeny, showing a complex history of introgression among different chromosomes in the young 'AA' subclade containing the two domesticated species. This study highlights the prevalence of functionally coupled disease resistance genes and identifies many new haplotypes of potential use for future crop protection. Finally, this study marks a milestone in modern rice research with the release of a complete long-read assembly of IR 8 'Miracle Rice', which relieved famine and drove the Green Revolution in Asia 50 years ago.

New genes drive the evolution of gene interaction networks in the human and mouse genomes.

BACKGROUND: The origin of new genes with novel functions creates genetic and phenotypic diversity in organisms. To acquire functional roles, new genes must integrate into ancestral gene-gene interaction (GGI) networks. The mechanisms by which new genes are integrated into ancestral networks, and their evolutionary significance, are yet to be characterized. Herein, we present a study investigating the rates and patterns of new gene-driven evolution of GGI networks in the human and mouse genomes. RESULTS: We examine the network topological and functional evolution of new genes that originated at various stages in the human and mouse lineages by constructing and analyzing three different GGI datasets. We find a large number of new genes integrated into GGI networks throughout vertebrate evolution. These genes experienced a gradual integration process into GGI networks, starting on the network periphery and gradually becoming highly connected hubs, and acquiring pleiotropic and essential functions. We identify a few human lineage-specific hub genes that have evolved brain development-related functions. Finally, we explore the possible underlying mechanisms driving the GGI network evolution and the observed patterns of new gene integration process. CONCLUSIONS: Our results unveil a remarkable network topological integration process of new genes: over 5000 new genes were integrated into the ancestral GGI networks of human and mouse; new genes gradually acquire increasing number of gene partners; some human-specific genes evolved into hub structure with critical phenotypic effects. Our data cast new conceptual insights into the evolution of genetic networks.

Origin and domestication of papaya Yh chromosome.

Sex in papaya is controlled by a pair of nascent sex chromosomes. Females are XX, and two slightly different Y chromosomes distinguish males (XY) and hermaphrodites (XY(h)). The hermaphrodite-specific region of the Y(h) chromosome (HSY) and its X chromosome counterpart were sequenced and analyzed previously. We now report the sequence of the entire male-specific region of the Y (MSY). We used a BAC-by-BAC approach to sequence the MSY and resequence the Y regions of 24 wild males and the Y(h) regions of 12 cultivated hermaphrodites. The MSY and HSY regions have highly similar gene content and structure, and only 0.4% sequence divergence. The MSY sequences from wild males include three distinct haplotypes, associated with the populations' geographic locations, but gene flow is detected for other genomic regions. The Y(h) sequence is highly similar to one Y haplotype (MSY3) found only in wild dioecious populations from the north Pacific region of Costa Rica. The low MSY3-Y(h) divergence supports the hypothesis that hermaphrodite papaya is a product of human domestication. We estimate that Y(h) arose only ∼ 4000 yr ago, well after crop plant domestication in Mesoamerica >6200 yr ago but coinciding with the rise of the Maya civilization. The Y(h) chromosome has lower nucleotide diversity than the Y, or the genome regions that are not fully sex-linked, consistent with a domestication bottleneck. The identification of the ancestral MSY3 haplotype will expedite investigation of the mutation leading to the domestication of the hermaphrodite Y(h) chromosome. In turn, this mutation should identify the gene that was affected by the carpel-suppressing mutation that was involved in the evolution of males.

Evolution of gene structural complexity: an alternative-splicing-based model accounts for intron-containing retrogenes.

The structure of eukaryotic genes evolves extensively by intron loss or gain. Previous studies have revealed two models for gene structure evolution through the loss of introns: RNA-based gene conversion, dubbed the Fink model and retroposition model. However, retrogenes that experienced both intron loss and intron-retaining events have been ignored; evolutionary processes responsible for the variation in complex exon-intron structure were unknown. We detected hundreds of retroduplication-derived genes in human (Homo sapiens), fly (Drosophila melanogaster), rice (Oryza sativa), and Arabidopsis (Arabidopsis thaliana) and categorized them either as duplicated genes that have all introns lost or as duplicated genes that have at least lost one and retained one intron compared with the parental copy (intron-retaining [IR] type). Our new model attributes intron retention alternative splicing to the generation of these IR-type gene pairs. We presented 25 parental genes that have an intron retention isoform and have retained introns in the same locations in the IR-type duplicate genes, which directly support our hypothesis. Our alternative-splicing-based model in conjunction with the retroposition and Fink models can explain the IR-type gene observed. We discovered a greater percentage of IR-type genes in plants than in animals, which may be due to the abundance of intron retention cases in plants. Given the prevalence of intron retention in plants, this new model gives a support that plant genomes have very complex gene structures.

Genome of the long-living sacred lotus (Nelumbo nucifera Gaertn.).

BACKGROUND: Sacred lotus is a basal eudicot with agricultural, medicinal, cultural and religious importance. It was domesticated in Asia about 7,000 years ago, and cultivated for its rhizomes and seeds as a food crop. It is particularly noted for its 1,300-year seed longevity and exceptional water repellency, known as the lotus effect. The latter property is due to the nanoscopic closely packed protuberances of its self-cleaning leaf surface, which have been adapted for the manufacture of a self-cleaning industrial paint, Lotusan. RESULTS: The genome of the China Antique variety of the sacred lotus was sequenced with Illumina and 454 technologies, at respective depths of 101× and 5.2×. The final assembly has a contig N50 of 38.8 kbp and a scaffold N50 of 3.4 Mbp, and covers 86.5% of the estimated 929 Mbp total genome size. The genome notably lacks the paleo-triplication observed in other eudicots, but reveals a lineage-specific duplication. The genome has evidence of slow evolution, with a 30% slower nucleotide mutation rate than observed in grape. Comparisons of the available sequenced genomes suggest a minimum gene set for vascular plants of 4,223 genes. Strikingly, the sacred lotus has 16 COG2132 multi-copper oxidase family proteins with root-specific expression; these are involved in root meristem phosphate starvation, reflecting adaptation to limited nutrient availability in an aquatic environment. CONCLUSIONS: The slow nucleotide substitution rate makes the sacred lotus a better resource than the current standard, grape, for reconstructing the pan-eudicot genome, and should therefore accelerate comparative analysis between eudicots and monocots.

Rapid divergence and expansion of the X chromosome in papaya.

X chromosomes have long been thought to conserve the structure and gene content of the ancestral autosome from which the sex chromosomes evolved. We compared the recently evolved papaya sex chromosomes with a homologous autosome of a close relative, the monoecious Vasconcellea monoica, to infer changes since recombination stopped between the papaya sex chromosomes. We sequenced 12 V. monoica bacterial artificial chromosomes, 11 corresponding to the papaya X-specific region, and 1 to a papaya autosomal region. The combined V. monoica X-orthologous sequences are much shorter (1.10 Mb) than the corresponding papaya region (2.56 Mb). Given that the V. monoica genome is 41% larger than that of papaya, this finding suggests considerable expansion of the papaya X; expansion is supported by a higher repetitive sequence content of the X compared with the papaya autosomal sequence. The alignable regions include 27 transcript-encoding sequences, only 6 of which are functional X/V. monoica gene pairs. Sequence divergence from the V. monoica orthologs is almost identical for papaya X and Y alleles; the Carica-Vasconcellea split therefore occurred before the papaya sex chromosomes stopped recombining, making V. monoica a suitable outgroup for inferring changes in papaya sex chromosomes. The papaya X and the hermaphrodite-specific region of the Y(h) chromosome and V. monoica have all gained and lost genes, including a surprising amount of changes in the X.

Sequencing papaya X and Yh chromosomes reveals molecular basis of incipient sex chromosome evolution.

Sex determination in papaya is controlled by a recently evolved XY chromosome pair, with two slightly different Y chromosomes controlling the development of males (Y) and hermaphrodites (Y(h)). To study the events of early sex chromosome evolution, we sequenced the hermaphrodite-specific region of the Y(h) chromosome (HSY) and its X counterpart, yielding an 8.1-megabase (Mb) HSY pseudomolecule, and a 3.5-Mb sequence for the corresponding X region. The HSY is larger than the X region, mostly due to retrotransposon insertions. The papaya HSY differs from the X region by two large-scale inversions, the first of which likely caused the recombination suppression between the X and Y(h) chromosomes, followed by numerous additional chromosomal rearrangements. Altogether, including the X and/or HSY regions, 124 transcription units were annotated, including 50 functional pairs present in both the X and HSY. Ten HSY genes had functional homologs elsewhere in the papaya autosomal regions, suggesting movement of genes onto the HSY, whereas the X region had none. Sequence divergence between 70 transcripts shared by the X and HSY revealed two evolutionary strata in the X chromosome, corresponding to the two inversions on the HSY, the older of which evolved about 7.0 million years ago. Gene content differences between the HSY and X are greatest in the older stratum, whereas the gene content and order of the collinear regions are identical. Our findings support theoretical models of early sex chromosome evolution.

Construction of physical maps for the sex-specific regions of papaya sex chromosomes.

BACKGROUND: Papaya is a major fruit crop in tropical and subtropical regions worldwide. It is trioecious with three sex forms: male, female, and hermaphrodite. Sex determination is controlled by a pair of nascent sex chromosomes with two slightly different Y chromosomes, Y for male and Yh for hermaphrodite. The sex chromosome genotypes are XY (male), XYh (hermaphrodite), and XX (female). The papaya hermaphrodite-specific Yh chromosome region (HSY) is pericentromeric and heterochromatic. Physical mapping of HSY and its X counterpart is essential for sequencing these regions and uncovering the early events of sex chromosome evolution and to identify the sex determination genes for crop improvement. RESULTS: A reiterate chromosome walking strategy was applied to construct the two physical maps with three bacterial artificial chromosome (BAC) libraries. The HSY physical map consists of 68 overlapped BACs on the minimum tiling path, and covers all four HSY-specific Knobs. One gap remained in the region of Knob 1, the only knob structure shared between HSY and X, due to the lack of HSY-specific sequences. This gap was filled on the physical map of the HSY corresponding region in the X chromosome. The X physical map consists of 44 BACs on the minimum tiling path with one gap remaining in the middle, due to the nature of highly repetitive sequences. This gap was filled on the HSY physical map. The borders of the non-recombining HSY were defined genetically by fine mapping using 1460 F2 individuals. The genetically defined HSY spanned approximately 8.5 Mb, whereas its X counterpart extended about 5.4 Mb including a 900 Kb region containing the Knob 1 shared by the HSY and X. The 8.5 Mb HSY corresponds to 4.5 Mb of its X counterpart, showing 4 Mb (89%) DNA sequence expansion. CONCLUSION: The 89% increase of DNA sequence in HSY indicates rapid expansion of the Yh chromosome after genetic recombination was suppressed 2-3 million years ago. The genetically defined borders coincide with the common BACs on the minimum tiling paths of HSY and X. The minimum tiling paths of HSY and its X counterpart are being used for sequencing these X and Yh-specific regions.

The sex-specific region of sex chromosomes in animals and plants.

Our understanding of the evolution of sex chromosomes has increased greatly in recent years due to a number of molecular evolutionary investigations in divergent sex chromosome systems, and these findings are reshaping theories of sex chromosome evolution. In particular, the dynamics of the sex-determining region (SDR) have been demonstrated by recent findings in ancient and incipient sex chromosomes. Radical changes in genomic structure and gene content in the male specific region of the Y chromosome between human and chimpanzee indicated rapid evolution in the past 6 million years, defying the notion that the pace of evolution in the SDR was fast at early stages but slowed down overtime. The chicken Z and the human X chromosomes appeared to have acquired testis-expressed genes and expanded in intergenic regions. Transposable elements greatly contributed to SDR expansion and aided the trafficking of genes in the SDR and its X or Z counterpart through retrotransposition. Dosage compensation is not a destined consequence of sex chromosomes as once thought. Most X-linked microRNA genes escape silencing and are expressed in testis. Collectively, these findings are challenging many of our preconceived ideas of the evolutionary trajectory and fates of sex chromosomes.

Multiple developmental processes underlie sex differentiation in angiosperms.

The production of unisexual flowers has evolved numerous times in dioecious and monoecious plant taxa. Based on repeated evolutionary origins, a great variety of developmental and genetic mechanisms underlying unisexual flower development is predicted. Here, we comprehensively review the modes of development of unisexual flowers, test potential correlations with sexual system, and end with a synthesis of the genetics and hormonal regulation of plant sex determination. We find that the stage of organ abortion in male and female flowers is temporally correlated within species and also confirm that the arrest of development does not tend to occur preferentially at a particular stage, or via a common process.

Construction of papaya male and female BAC libraries and application in physical mapping of the sex chromosomes.

Papaya is a major fruit crop in the tropics and has recently evolved sex chromosomes. Towards sequencing the papaya sex chromosomes, two bacterial artificial chromosome (BAC) libraries were constructed from papaya male and female genomic DNA. The female BAC library was constructed using restriction enzyme BstY I and consists of 36,864 clones with an average insert size of 104 kb, providing 10.3x genome equivalents. The male BAC library was constructed using restriction enzyme EcoR I and consists of 55,296 clones with an average insert size of 101 kb, providing 15.0x genome equivalents. The male BAC library was used in constructing the physical map of the male-specific region of the male Y chromosome (MSY) and in filling gaps and extending the physical map of the hermaphrodite-specific region of the Y(h) chromosome (HSY) and the X chromosome physical map. The female BAC library was used to extend the X physical map gap. The MSY, HSY, and X physical maps offer a unique opportunity to study chromosomal rearrangements, Y chromosome degeneration, and dosage compensation of the papaya nascent sex chromosomes.

The draft genome of the transgenic tropical fruit tree papaya (Carica papaya Linnaeus).

Papaya, a fruit crop cultivated in tropical and subtropical regions, is known for its nutritional benefits and medicinal applications. Here we report a 3x draft genome sequence of 'SunUp' papaya, the first commercial virus-resistant transgenic fruit tree to be sequenced. The papaya genome is three times the size of the Arabidopsis genome, but contains fewer genes, including significantly fewer disease-resistance gene analogues. Comparison of the five sequenced genomes suggests a minimal angiosperm gene set of 13,311. A lack of recent genome duplication, atypical of other angiosperm genomes sequenced so far, may account for the smaller papaya gene number in most functional groups. Nonetheless, striking amplifications in gene number within particular functional groups suggest roles in the evolution of tree-like habit, deposition and remobilization of starch reserves, attraction of seed dispersal agents, and adaptation to tropical daylengths. Transgenesis at three locations is closely associated with chloroplast insertions into the nuclear genome, and with topoisomerase I recognition sites. Papaya offers numerous advantages as a system for fruit-tree functional genomics, and this draft genome sequence provides the foundation for revealing the basis of Carica's distinguishing morpho-physiological, medicinal and nutritional properties.