Chromosome-level reference genome of Tetrastigma hemsleyanum (Vitaceae) provides insights into genomic evolution and the biosynthesis of phenylpropanoids and flavonoids.

Here, we present a high-quality chromosome-scale genome assembly (2.19 Gb) and annotation of Tetrastigma hemsleyanum, a perennial herbaceous liana native to subtropical China with diverse medicinal applications. Approximately 73% of the genome was comprised of transposable elements (TEs), of which long terminal repeat retrotransposons (LTR-RTs) were a predominant group (69% of the genome). The genome size increase of T. hemsleyanum (relative to Vitis species) was mostly due to the proliferation of LTR-RTs. Of the different modes of gene duplication identified, transposed duplication (TRD) and dispersed duplication (DSD) were the predominant ones. Genes, particularly those involved in the phenylpropanoid-flavonoid (PF) pathway and those associated with therapeutic properties and environmental stress resistance, were significantly amplified through recent tandem duplications. We dated the divergence of two intraspecific lineages in Southwest (SW) versus Central-South-East (CSE) China to the late Miocene (approximately 5.2 million years ago). Of those, the former showed more upregulated genes and metabolites. Based on resequencing data of 38 individuals representing both lineages, we identified various candidate genes related to 'response to stimulus' and 'biosynthetic process', including ThFLS11, which is putatively involved in flavonoid accumulation. Overall, this study provides abundant genomic resources for future evolutionary, ecological, and functional genomics studies in T. hemsleyanum and related species.

What explains the high island endemicity of Philippine Rafflesia? A species distribution modeling analysis of three threatened parasitic plant species and their hosts.

PREMISE: Rafflesia are rare holoparasitic plants. In the Philippines, all but one species are found only on single islands. This study aimed to better understand the factors contributing to this distributional pattern. Specifically, we sought to determine whether narrow environmental tolerances of host and/or parasite species might explain their island endemicity. METHODS: We used Maxent species distribution modeling to identify areas with suitable habitat for R. lagascae, R. lobata, and R. speciosa and their Tetrastigma host species. These analyses were carried out for current climate conditions and two future climate change scenarios. RESULTS: Although species distribution models indicated suitable environmental conditions for the Tetrastigma host species in many parts of the Philippines, considerably fewer areas were inferred to have suitable conditions for the three Rafflesia species. Some of these areas are on islands from which they have not been reported. All three species will face significant threats as a result of climate change. CONCLUSIONS: Our results suggest that limited inter-island dispersibility and/or specific environmental requirements are likely responsible for the current pattern of island endemicity of the three Rafflesia species, rather than environmental requirements of their Tetrastigma host species.

Glycoside constituents from Cayratia geniculata.

Using various chromatographic separations, six glycoside derivatives (1-6), including one new ent-labdane glucoside named cayratioside (1), were isolated from the methanol extract of Cayratia geniculata stems and leaves. Their structures were elucidated by detailed analysis of the 1D, 2D NMR, and HRESIQTOF mass spectra. The inhibitory effect of 1-6 on LPS-induced NO production in RAW264.7 cells was also evaluated. Among isolated compounds, 1 exhibited moderate activity with an IC50 value of 59.65 ± 1.85 µM.

Herbarium specimens reveal century-long trait shifts in poison ivy due to anthropogenic CO2 emissions.

PREMISE: Previous experimental studies have shown that poison ivy (Toxicodendron radicans; Anacardicaceae) responds to elevated CO2 with increased leaf production, water-use efficiency, and toxicity (allergenic urushiol). However, long-term field data suggest no increase in poison ivy abundance over time. Using herbarium specimens, we examined whether poison ivy and other species shifted leaf traits under natural conditions with increasing atmospheric CO2 (pCO2 ) over the past century. METHODS: We measured stomatal density, leaf area, leaf N, leaf C:N, leaf carbon isotope discrimination (Δleaf ), and intrinsic water-use efficiency (iWUE) from 327 specimens collected from 1838 to 2020 across Pennsylvania. We compared poison ivy's responses to two evolutionarily related tree species, Toxicodendron vernix and Rhus typhina (Anacardiacae) and one ecological analog, Parthenocissus quinquefolia (Vitaceae), a common co-occurring liana. RESULTS: Stomatal density significantly decreased (P < 0.05) in poison ivy and the ecologically similar liana P. quinquefolia over the past century, but did not change in the related trees T. vernix and R. typhina. None of these species showed significant trends in changes in leaf N or C:N. Surprisingly, in poison ivy, but not the other species, Δleaf increased with increased pCO2 , corresponding to significant declines in iWUE over time. CONCLUSIONS: In contrast to the results of short-term experimental studies, iWUE decreased in poison ivy over the last century. Trait responses to pCO2 varied by species. Herbarium specimens suggest that realized long-term plant physiological responses to increased CO2 may not be reflected in short-term experimental growth studies, highlighting the value of collections.

A New Ceramide from Cissus Aralioides Baker (Vitaceae) and its Antimicrobial Activity.

Purification through repeated column chromatography over silica gel and Sephadex LH-20 of the ethanol extract of the stems of Cissus aralioides (Baker) Planch. resulted in the isolation of a new ceramide, aralioidamide A (1), along with five known compounds (2-6). Their structures were determined by the extensive analyses of their spectroscopic (1D and 2D NMR) and spectrometric data, and comparison with those reported in the literature. Aralioidamide A (1) displayed weak antibacterial activity (MIC=256 µg/mL) against Bacillus subtilis, Staphylococcus aureus and Shigella flexneri and was inactive (MIC>256 µg/mL) against the tested fungi.

Detecting signals of adaptive evolution in grape plastomes with a focus on the Cretaceous-Palaeogene (K/Pg) transition.

BACKGROUND AND AIMS: Although plastid genes are widely used in phylogenetic studies, signals of positive selection have been scarcely investigated in the grape family. The plastomes from 91 accessions of Vitaceae were examined to understand the extent to which positive selection is present and to identify which genes are involved. Moreover, the changes through time of genes under episodic positive selection were investigated and the hypothesis of an adaptive process following the Cretaceous-Palaeogene (K/Pg) transition about 66 million years ago was tested. METHODS: Different codon-substitution models were used to assess pervasive and episodic positive selection events on 70 candidate plastid genes. Divergence times between lineages were estimated and stochastic character mapping analysis was used to simulate variation over time of the genes found to be under episodic positive selection. KEY RESULTS: A total of 20 plastid genes (29 %) showed positive selection. Among them, 14 genes showed pervasive signatures of positive selection and nine genes showed episodic signatures of positive selection. In particular, four of the nine genes (psbK, rpl20, rpoB, rps11) exhibited a similar pattern showing an increase in the rate of variation close to the K/Pg transition. CONCLUSION: Multiple analyses have shown that the grape family has experienced ancient and recent positive selection events and that the targeted genes are involved in essential functions such as photosynthesis, self-replication and metabolism. Our results are consistent with the idea that the K/Pg transition has favoured an increased rate of change in some genes. Intense environmental perturbations have influenced the rapid diversification of certain lineages, and new mutations arising on some plastid genes may have been fixed by natural selection over the course of many generations.

Tannin phenotyping of the Vitaceae reveals a phylogenetic linkage of epigallocatechin in berries and leaves.

BACKGROUND AND AIMS: Condensed tannins, responsible for berry and wine astringency, may have been selected during grapevine domestication. This work examines the phylogenetic distribution of condensed tannins throughout the Vitaceae phylogenetic tree. METHODS: Green berries and mature leaves of representative true-to-type members of the Vitaceae were collected before 'véraison', freeze-dried and pulverized, and condensed tannins were measured following depolymerization by nucleophilic addition of 2-mercaptoethanol to the C4 of the flavan-3-ol units in an organic acidic medium. Reaction products were separated and quantified by ultrahigh pressure liquid chromatography/diode array detection/mass spectrometry. KEY RESULTS AND CONCLUSIONS: The original ability to incorporate epigallocatechin (EGC) into grapevine condensed tannins was lost independently in both the American and Eurasian/Asian branches of the Vitaceae, with exceptional cases of reversion to the ancestral EGC phenotype. This is particularly true in the genus Vitis, where we now find two radically distinct groups differing with respect to EGC content. While Vitis species from Asia are void of EGC, 50 % of the New World Vitis harbour EGC. Interestingly, the presence of EGC is tightly coupled with the degree of leaf margin serration. Noticeably, the rare Asian EGC-forming species are phylogenetically close to Vitis vinifera, the only remnant representative of Vitis in Eurasia. Both the wild ancestral V. vinifera subsp. sylvestris as well as the domesticated V. vinifera subsp. sativa can accumulate EGC and activate galloylation biosynthesis that compete for photoassimilates and reductive power.

Comparative transcriptome analyses between cultivated and wild grapes reveal conservation of expressed genes but extensive rewiring of co-expression networks.

KEY MESSAGE: The transcriptomes of wild and cultivated grapes consists of similar expressed genes but distinct wiring of co-expressed modules associated with environmental conditions. Grapevine is an important fruit crop worldwide, with high economic value and widespread distribution. Commercial production is based on Vitis vinifera, and, to a lesser extent, on hybrids with American grapes, such as V. labrusca. Wild grape relatives are important sources of resistance against biotic and abiotic factors; however, their global gene expression patterns remain poorly characterized. We associated genome-wide transcript profiling to phenotypic analyses to investigate the responses of cultivated and wild vines to vineyard conditions. The expressed genes in the Vitis reference transcriptome are largely shared by wild grapes, V. labrusca hybrids and vinifera cultivars. In contrast, significant differential regulation between wild and vinifera genotypes represents 80% of gene expression variation, regardless of the environment. In wild grapes, genes associated to regulatory processes are downregulated, whereas those involved in metabolic pathways are upregulated, in comparison to vinifera. Photosynthesis-related ontologies are overrepresented in the induced genes, in agreement with higher contents of chlorophyll in wild grapes. Co-regulated gene network analyses provide evidence of more complex transcriptome organization in vinifera. In wild grapes, genes involved in signaling pathways of stress-related hormones are overrepresented in modules associated with the environment. Consensus network analyses revealed high preservation within co-regulated gene modules between cultivated and wild grapes, but divergent relationships among the expression clusters. In conclusion, the distinct phenotypes of wild and cultivated grapes are underlain by differences in gene expression, but also by distinct higher-order organization of the transcriptome and contrasting association of co-expressed gene clusters with the environment.

Phylogenomic relationships and character evolution of the grape family (Vitaceae).

The grape family consists of 16 genera and ca. 950 species. It is best known for the economically important fruit crop - the grape Vitis vinifera. The deep phylogenetic relationships and character evolution of the grape family have attracted the attention of researchers in recent years. We herein reconstruct the phylogenomic relationships within Vitaceae using nuclear and plastid genes based on the Hyb-Seq approach and test the newly proposed classification system of the family. The five tribes of the grape family, including Ampelopsideae, Cayratieae, Cisseae, Parthenocisseae, and Viteae, are each robustly supported by both nuclear and chloroplast genomic data and the backbone relationships are congruent with previous reports. The cupular floral disc (raised above and free from ovary at the upper part) is an ancestral state of Vitaceae, with the inconspicuous floral disc as derived in the tribe Parthenocisseae, and the state of adnate to the ovary as derived in the tribe Viteae. The 5-merous floral pattern was inferred to be the ancestral in Vitaceae, with the 4-merous flowers evolved at least two times in the family. The compound dichasial cyme (cymose with two secondary axes) is ancestral in Vitaceae and the thyrse inflorescence (a combination of racemose and cymose branching) in tribe Viteae is derived. The ribbon-like trichome only evolved once in Vitaceae, as a synapomorphy for the tribe Viteae.

Vein-to-blade ratio is an allometric indicator of leaf size and plasticity.

PREMISE: As a leaf expands, its shape dynamically changes. Previously, we documented an allometric relationship between vein and blade area in grapevine leaves. Larger leaves have a smaller ratio of primary and secondary vein area relative to blade area compared to smaller leaves. We sought to use allometry as an indicator of leaf size and plasticity. METHODS: We measured the ratio of vein-to-blade area from the same 208 vines across four growing seasons (2013, 2015, 2016, and 2017). Matching leaves by vine and node, we analyzed the correlation between the size and shape of grapevine leaves as repeated measures with climate variables across years. RESULTS: The proportion of leaf area occupied by vein and blade exponentially decreased and increased, respectively, during leaf expansion making their ratio a stronger indicator of leaf size than area itself. Total precipitation and leaf wetness hours of the previous year but not the current showed strong negative correlations with vein-to-blade ratio, whereas maximum air temperature from the previous year was positively correlated. CONCLUSIONS: Our results demonstrate that vein-to-blade ratio is a strong allometric indicator of leaf size and plasticity in grapevines measured across years. Grapevine leaf primordia are initiated in buds the year before they emerge, and we found that total precipitation and maximum air temperature of the previous growing season exerted the largest statistically significant effects on leaf morphology. Vein-to-blade ratio is a promising allometric indicator of relationships between leaf morphology and climate, the robustness of which should be explored further.

Review of the Traditional Uses, Phytochemistry, and Pharmacological Activities of Rhoicissus Species (Vitaceae).

Species within the genus Rhoicissus (Vitaceae) are commonly used in South African traditional medicine. The current review discusses the occurrence, distribution, traditional uses, phytochemistry, and pharmacological properties of Rhoicissus species covering the period 1981-2020. The data reported were systematically collected, read, and analysed from scientific electronic databases including Scopus, Scifinder, Pubmed, and Google Scholar. Reported evidence indicates that species in this genus are used for the treatment of gastrointestinal complaints, sexually transmitted infections (STIs), and infertility, as well as to tone the uterus during pregnancy and to facilitate delivery. Pharmacological studies have further shown that members of the Rhoicissus genus display antidiabetic, uterotonic, ascaricidal, hepatoprotective, antioxidant, antimicrobial, anticancer, and anti-inflammatory properties. They are linked to the presence of bioactive compounds isolated from the genus. Hence, Rhoicissus species can potentially be an alternative therapeutic strategy to treat diseases and develop safer and more potent drugs to combat diseases. Plant species of this genus have valuable medicinal benefits due to their significant pharmacological potential. However, scientific investigation and information of the therapeutic potential of Rhoicissus remain limited as most of the species in the genus have not been fully exploited. Therefore, there is a need for further investigations to exploit the therapeutic potential of the genus Rhoicissus. Future studies should evaluate the phytochemical, pharmacological, and toxicological activities, as well as the mode of action, of Rhoicissus crude extracts and secondary compounds isolated from the species.

Dates and rates in grape's plastomes: evolution in slow motion.

The family Vitaceae includes the domesticated grapevine (Vitis vinifera), one of the most economically important crops in the world. Despite the importance of Vitaceae, there is still considerable controversy surrounding their phylogenetic relationships and evolutionary timescales. Moreover, variation in rates of molecular evolution among Vitaceae remains mostly unexplored. The present research aims to fill these knowledge gaps through the analysis of plastome sequences. Thirteen newly sequenced grape plastomes are presented and their phylogenetic relationships examined. Divergence times and absolute substitution rates are inferred under different molecular clocks by the analysis of 95 non-coding plastid regions and 43 representative accessions of the major lineages of Vitaceae. Furthermore, the phylogenetic informativeness of non-coding plastid regions is investigated. We find strong evidence in favor of the random local clock model and rate heterogeneity within Vitaceae. Substitution rates decelerate in Ampelocissus, Ampelopsis, Nekemias, Parthenocissus, Rhoicissus, and Vitis, with genus Vitis showing the lowest values up to a minimum of ~ 4.65 × 10-11 s/s/y. We suggest that liana-like species of Vitaceae evolve slower than erect growth habit plants and we invoke the "rate of mitosis hypothesis" to explain the observed pattern of the substitution rates. We identify a reduced set of 20 non-coding regions able to accurately reconstruct the phylogeny of Vitaceae and we provide a detailed description of all 152 non-coding regions identified in the plastomes of subg. Vitis. These polymorphic regions will find their applications in phylogenetics, phylogeography, and population genetics as well in grapes identification through DNA barcoding techniques.

Two new meroterpenoids from the aerial parts of Ampelopsis cantoniensis (Vitaceae).

From EtOAc-soluble fraction of the aerial parts of Ampelopsis cantoniensis (Vitaceae), two new meroterpenoids named ampechromonol A (1) and ampechromonol B (2), together with five known compounds (3-7), were isolated. Their structures were elucidated based on NMR spectroscopic analysis. The plausible biosynthesis pathway for the formation of two new meroterpenoids was proposed. This research is the first isolation of meroterpenoids from Ampelopsis genus. Compounds 1 and 2 showed weak cytotoxicity against MCF-7 breast cancer cells.

Inhibition potential of phenolic constituents from the aerial parts of Tetrastigma hemsleyanum against soluble epoxide hydrolase and nitric oxide synthase.

The aerial parts of Tetrastigma hemsleyanum (APTH) have been used as a functional tea in China. The purpose of the current study was to identify the bioactive constituents with inhibitory activity against soluble epoxide hydrolase (sEH) and inducible nitric oxide synthase (iNOS), which are jointly considered potential therapeutic targets for vascular system diseases. In the present study, 39 compounds (1-39) were isolated from the APTH. Among them, compounds 8, 10, 12, 16, 17, 19, and 32 displayed potential activities, with IC50 values ranging from 4.5 to 9.5 µM, respectively, and all in non-competitive inhibition mode. Compounds 5, 10, 12, 19, and 32 displayed potent iNOS inhibitory effects, with IC50 values ranging from 15.6 to 47.3 µM. The results obtained in this work contribute to a better understanding of the pharmacological activities of T. hemsleyanum and its potential application as a functional food.

High-throughput sequencing data clarify evolutionary relationships among North American Vitis species and improve identification in USDA Vitis germplasm collections.

PREMISE OF THE STUDY: Grapes are one of the most economically important berry crops worldwide, with the vast majority of production derived from the domesticated Eurasian species Vitis vinifera. Expansion of production into new areas, development of new cultivars, and concerns about adapting grapevines for changing climates necessitate the use of wild grapevine species in breeding programs. Diversity within Vitis has long been a topic of study; however, questions remain regarding relationships between species. Furthermore, the identity of some living accessions is unclear. METHODS: This study generated 11,020 single nucleotide polymorphism (SNP) markers for more than 300 accessions in the USDA-ARS grape germplasm repository using genotyping-by-sequencing. Resulting data sets were used to reconstruct evolutionary relationships among several North American and Eurasian Vitis species, and to suggest taxonomic labels for previously unidentified and misidentified germplasm accessions based on genetic distance. KEY RESULTS: Maximum likelihood analyses of SNP data support the monophyly of Vitis, subg. Vitis, a Eurasian subg. Vitis clade, and a North American subg. Vitis clade. Data delineate species groups within North America. In addition, analysis of genetic distance suggested taxonomic identities for 20 previously unidentified Vitis accessions and for 28 putatively misidentified accessions. CONCLUSIONS: This work advances understanding of Vitis evolutionary relationships and provides the foundation for ongoing germplasm enhancement. It supports conservation and breeding efforts by contributing to a growing genetic framework for identifying novel genetic variation and for incorporating new, unsampled populations into the germplasm repository system.

The Constituents of the Stems of Cissus assamica and Their Bioactivities.

Fifty-five compounds were isolated from the fresh stems of Cissus assamica, including 14 benzenoids, 11 triterpenes, nine steroids, five tocopherols, five chlorophylls, four flavonoids, two benzoquinones, two tannins, and three other compounds. Their structures were constructed by 1D and 2D nuclear magnetic resonance (NMR) and mass spectral data, and were also identified by a comparison of their spectral data with those reported in the literature. Among these isolates, 1,2-bis-(5--tocopheryl) ethane (51) was reported for the first time from natural sources. Some purified compounds were examined for their anti-inflammatory and anticancer bioactivities. The results indicated that betulinic acid (16) exhibited strong inhibition of superoxide anion generation with IC50 value of 0.2 ± 0.1 µM, while betulinic acid (16) and pheophytin-a (47) inhibited elastase release with IC50 value of 2.7 ± 0.3 and 5.3 ± 1.0 µM, respectively. In addition, betulinic acid (16) and epi-glut-5(6)-en-ol (18) exhibited potential cytotoxicity to non-small-cell lung carcinoma (NCI-H226) and colon cancer (HCT-116) cell lines with IC50 values in the range of 1.6 to 9.1 µM.

Anti-Dermatophyte and Anti-Malassezia Activity of Extracts Rich in Polymeric Flavan-3-ols Obtained from Vitis vinifera Seeds.

Several human skin diseases are associated with fungi as dermatophytes and Malassezia. Skin mycoses are increasing and new alternatives to conventional treatments with improved efficacy and/or safety profiles are desirable. For the first time, the anti-dermatophytes and the anti-Malassezia activities of Vitis vinifera seed extracts obtained from different table and wine cultivars have been evaluated. Geometric minimal inhibitory concentration ranged from 20 to 97 µg/mL for dermatophytes and from 32 to 161 µg/mL for Malassezia furfur. Dried grape seed extracts analyzed by HPLC/DAD/ESI/MS showed different quali-quantitative compositions in terms of monomeric and polymeric flavan-3-ols. The minimal inhibitory concentrations for Trichophyton mentagrophytes and for M. furfur were inversely correlated with the amount of the polymeric fraction (r = -0.7639 and r = -0.7228, respectively). Differently, the antifungal activity against T. mentagrophytes was not correlated to the content of flavan-3-ol monomers (r = 0.2920) and only weakly correlated for M. furfur (r = -0.53604). These results suggest that extracts rich in polymeric flavan-3-ols, recovered from V. vinifera seeds, could be used for the treatment of skin fungal infections. Copyright © 2016 John Wiley & Sons, Ltd.

Phylogeny of the Ampelocissus-Vitis clade in Vitaceae supports the New World origin of the grape genus.

The grapes and the close allies in Vitaceae are of great agronomic and economic importance. Our previous studies showed that the grape genus Vitis was closely related to three tropical genera, which formed the Ampelocissus-Vitis clade (including Vitis, Ampelocissus, Nothocissus and Pterisanthes). Yet the phylogenetic relationships of the four genera within this clade remain poorly resolved. Furthermore, the geographic origin of Vitis is still controversial, because the sampling of the close relatives of Vitis was too limited in the previous studies. This study reconstructs the phylogenetic relationships within the clade, and hypothesizes the origin of Vitis in a broader phylogenetic framework, using five plastid and two nuclear markers. The Ampelocissus-Vitis clade is supported to be composed of five main lineages. Vitis includes two described subgenera each as a monophyletic group. Ampelocissus is paraphyletic. The New World Ampelocissus does not form a clade and shows a complex phylogenetic relationship, with A. acapulcensis and A. javalensis forming a clade, and A. erdvendbergiana sister to Vitis. The majority of the Asian Ampelocissus species form a clade, within which Pterisanthes is nested. Pterisanthes is polyphyletic, suggesting that the lamellate inflorescence characteristic of the genus represents convergence. Nothocissus is sister to the clade of Asian Ampelocissus and Pterisanthes. The African Ampelocissus forms a clade with several Asian species. Based on the Bayesian dating and both the RASP and Lagrange analyses, Vitis is inferred to have originated in the New World during the late Eocene (39.4Ma, 95% HPD: 32.6-48.6Ma), then migrated to Eurasia in the late Eocene (37.3Ma, 95% HPD: 30.9-45.1Ma). The North Atlantic land bridges (NALB) are hypothesized to be the most plausible route for the Vitis migration from the New World to Eurasia, while intercontinental long distance dispersal (LDD) cannot be eliminated as a likely mechanism.

Cissus sicyoides: Pharmacological Mechanisms Involved in the Anti-Inflammatory and Antidiarrheal Activities.

The objective of this study was to evaluate the pharmacological mechanisms involved in anti-inflammatory and antidiarrheal actions of hydroalcoholic extract obtained from the leaves of Cissus sicyoides (HECS). The anti-inflammatory effect was evaluated by oral administration of HECS against acute model of edema induced by xylene, and the mechanisms of action were analysed by involvement of arachidonic acid (AA) and prostaglandin E2 (PGE2). The antidiarrheal effect of HECS was observed and we analyzed the motility and accumulation of intestinal fluid. We also analyzed the antidiarrheal mechanisms of action of HECS by evaluating the role of the opioid receptor, α2 adrenergic receptor, muscarinic receptor, nitric oxide (NO) and PGE2. The oral administration of HECS inhibited the edema induced by xylene and AA and was also able to significantly decrease the levels of PGE2. The extract also exhibited significant anti-diarrheal activity by reducing motility and intestinal fluid accumulation. This extract significantly reduced intestinal transit stimulated by muscarinic agonist and intestinal secretion induced by PGE2. Our data demonstrate that the mechanism of action involved in the anti-inflammatory effect of HECS is related to PGE2. The antidiarrheal effect of this extract may be mediated by inhibition of contraction by acting on the intestinal smooth muscle and/or intestinal transit.

Orthology inference in nonmodel organisms using transcriptomes and low-coverage genomes: improving accuracy and matrix occupancy for phylogenomics.

Orthology inference is central to phylogenomic analyses. Phylogenomic data sets commonly include transcriptomes and low-coverage genomes that are incomplete and contain errors and isoforms. These properties can severely violate the underlying assumptions of orthology inference with existing heuristics. We present a procedure that uses phylogenies for both homology and orthology assignment. The procedure first uses similarity scores to infer putative homologs that are then aligned, constructed into phylogenies, and pruned of spurious branches caused by deep paralogs, misassembly, frameshifts, or recombination. These final homologs are then used to identify orthologs. We explore four alternative tree-based orthology inference approaches, of which two are new. These accommodate gene and genome duplications as well as gene tree discordance. We demonstrate these methods in three published data sets including the grape family, Hymenoptera, and millipedes with divergence times ranging from approximately 100 to over 400 Ma. The procedure significantly increased the completeness and accuracy of the inferred homologs and orthologs. We also found that data sets that are more recently diverged and/or include more high-coverage genomes had more complete sets of orthologs. To explicitly evaluate sources of conflicting phylogenetic signals, we applied serial jackknife analyses of gene regions keeping each locus intact. The methods described here can scale to over 100 taxa. They have been implemented in python with independent scripts for each step, making it easy to modify or incorporate them into existing pipelines. All scripts are available from https://bitbucket.org/yangya/phylogenomic_dataset_construction.