Flavonoid and chromone-rich extract from Euscaphis Konishii Hayata leaf attenuated alcoholic liver injury in mice.

ETHNOPHARMACOLOGICAL RELEVANCE: Euscaphis konishii Hayata is a traditional medicinal plant in China, and its leaves are usually used to make dishes for hepatic or gastrointestinal issues by Chinese She nationality. Pharmacological analysis showed that E. konishii leaves contain high levels of flavonoids and chromones with favorable anti-hepatoma effect. AIM OF THE STUDY: The extract from E. konishii leaves was detected to evaluate its chemical composition, and the alcoholic liver injury mice model was adopted to elucidate its hepatoprotective effects. MATERIALS AND METHODS: The total leaf extract from E. konishii was separated by polyamide column to get the flavonoid and chromone-rich extract (FCE). Single compounds from FCE was purified by gel and Sephadex LH-20 chromatography and analyzed by nuclear magnetic resonance (NMR). The chemical component of FCE was confirmed and quantified by HPLC-MS. The OH·, O2-, DPPH and ABTS + free radical assays were adopted to estimate the antioxidant activity of FCE in vitro. The alcohol-fed model mice were established to assess the hepatoprotective capacity of FCE in vivo, through biochemical determination, histopathological analysis, mitochondrial function measurement, quantitative Real-Time Polymerase Chain Reaction (qRT-PCR) detection and Western blot determination. RESULTS: 8 flavonoids and 2 chromones were recognized in the FCEextract by both NMR and HPLC-MS. FCE represented strong free radicals scavenging activity in vitro. With oral administration, FCE (50, 100 and 200 mg/kg) dose-dependently decreased the serum levels of alanine aminotransferase (ALT), alkaline phosphatase (ALP) and aspartate aminotransferase (AST) in alcohol-fed mice. FCE gradually reduced the malondialdehyde (MDA) content, increased the activity of superoxide dismutase (SOD) and glutathione peroxidase (GSH-Px) in the alcohol-treated liver tissues. FCE also alleviated the hepatic inflammation, inhibited the hepatocyte apoptosis and lessened the alcohol-induced histological alteration and lipid accumulation in the liver tissues. FCE administration inhibited the overexpression of endoplasmic reticulum (ER) chaperones signaling and unfolded protein response (UPR) pathways to defense the ER-induced apoptosis. Pretreatment with FCE also restored the mitochondrial membrane potentials andadenosine triphosphate (ATP) levels, which in turn suppressed the Cytochrome C release and mitochondria-induced apoptosis. CONCLUSIONS: FCE conferred great protection against alcoholic liver injury, which might be associated with its viability through suppressing reactive oxygen species (ROS) stress and hepatocyte apoptosis.

Genomes of leafy and leafless Platanthera orchids illuminate the evolution of mycoheterotrophy.

To improve our understanding of the origin and evolution of mycoheterotrophic plants, we here present the chromosome-scale genome assemblies of two sibling orchid species: partially mycoheterotrophic Platanthera zijinensis and holomycoheterotrophic Platanthera guangdongensis. Comparative analysis shows that mycoheterotrophy is associated with increased substitution rates and gene loss, and the deletion of most photoreceptor genes and auxin transporter genes might be linked to the unique phenotypes of fully mycoheterotrophic orchids. Conversely, trehalase genes that catalyse the conversion of trehalose into glucose have expanded in most sequenced orchids, in line with the fact that the germination of orchid non-endosperm seeds needs carbohydrates from fungi during the protocorm stage. We further show that the mature plant of P. guangdongensis, different from photosynthetic orchids, keeps expressing trehalase genes to hijack trehalose from fungi. Therefore, we propose that mycoheterotrophy in mature orchids is a continuation of the protocorm stage by sustaining the expression of trehalase genes. Our results shed light on the molecular mechanism underlying initial, partial and full mycoheterotrophy.

The Euscaphis japonica genome and the evolution of malvids.

Malvids is one of the largest clades of rosids, includes 58 families and exhibits remarkable morphological and ecological diversity. Here, we report a high-quality chromosome-level genome assembly for Euscaphis japonica, an early-diverging species within malvids. Genome-based phylogenetic analysis suggests that the unstable phylogenetic position of E. japonica may result from incomplete lineage sorting and hybridization event during the diversification of the ancestral population of malvids. Euscaphis japonica experienced two polyploidization events: the ancient whole genome triplication event shared with most eudicots (commonly known as the γ event) and a more recent whole genome duplication event, unique to E. japonica. By resequencing 101 samples from 11 populations, we speculate that the temperature has led to the differentiation of the evergreen and deciduous of E. japonica and the completely different population histories of these two groups. In total, 1012 candidate positively selected genes in the evergreen were detected, some of which are involved in flower and fruit development. We found that reddening and dehiscence of the E. japonica pericarp and long fruit-hanging time promoted the reproduction of E. japonica populations, and revealed the expression patterns of genes related to fruit reddening, dehiscence and abscission. The key genes involved in pentacyclic triterpene synthesis in E. japonica were identified, and different expression patterns of these genes may contribute to pentacyclic triterpene diversification. Our work sheds light on the evolution of E. japonica and malvids, particularly on the diversification of E. japonica and the genetic basis for their fruit dehiscence and abscission.

Genetic diversity and population structure of Euscaphis japonica, a monotypic species.

BACKGROUND: Understanding plant genetic diversity is important for effective conservation and utilization of genetic resources. Euscaphis japonica (Thunb.) Dippel, is a monotypic species with high phenotypic diversity, narrow distribution, and small population size. In this study, we estimated the genetic diversity and population structure of E. japonica using nine natural populations and inter-simple sequence repeat (ISSR) markers. Our results could provide a theoretical reference for future conservation and utilization of E. japonica. RESULTS: We obtained a total of 122 DNA bands, of which 121 (99.18%) were polymorphic. The average number of effective alleles (Ne = 1.4975), Nei's gene diversity index (H = 0.3016), and Shannon's information index (I = 0.4630) revealed that E. japonica possessed a high level of genetic diversity. We observed that E. japonica consisted of both deciduous and evergreen populations. UPGMA tree showed that the evergreen and deciduous E. japonica form a sister group. There is little genetic differentiation among geographic populations based on STRUCTURE analysis. The Dice's similarity coefficient between the deciduous and evergreen populations was low, and the Fst value was high, indicating that these two types of groups have high degree of differentiation. CONCLUSION: Rich genetic diversity has been found in E. japonica, deciduous E. japonica and evergreen E. japonica populations, and genetic variation mainly exists within the population. The low-frequency gene exchange between deciduous and evergreen populations may be the result of the differentiation of deciduous and evergreen populations. We suggest that in-situ protection, seed collection, and vegetative propagation could be the methods for maintenance and conservation of E. japonica populations.

[Chemical constituents from the leaves of Cinnamomum camphora var. linaloolifera and their anti-inflammatory activities].

Thirteen compounds were isolated and purified from the leaves of Cinnamomum camphora by the macroporous resin,silica gel,and Sephadex LH-20 column chromatographies. Those compounds were further identified by IR,UV,MS,and NMR techniques:( 2 S)-1-( 3″,4″-methylenedioxy phenyl)-3-( 2',6'-dimethoxy-4'-hydroxyphenyl)-propan-2-ol( 1),( 2 R,3 R)-5,7-dimethoxy-3',4'-methylenedioxy flavanol( 2),9-hydroxysesamin( 3),sesamin( 4),piperitol( 5),kobusin( 6),(-)-aptosimon( 7),acuminatolide( 8),1ß,11-dihydroxy-5-eudesmene( 9),lasiodiplodin( 10),vanillin( 11),p-hydroxybenzaldehyde( 12),and p-hydroxybenzoic acid ethyl ester( 13). Compound 1 was a novel compound,and compounds 2,6,7,9 and 10 were isolated from Cinnamomum plants for the first time. Compounds 4,7 and 10 were found to possess good inhibitory effect on IL-6 production in LPS-induced BV2 cells at a concentration of 20 µmol·L-1 in the in vitro bioassay,with inhibition rates of 51. 26% ± 4. 13%,67. 82% ± 3. 77% and85. 81%±1. 19%,respectively.

Wolfberry genomes and the evolution of Lycium (Solanaceae).

Wolfberry Lycium, an economically important genus of the Solanaceae family, contains approximately 80 species and shows a fragmented distribution pattern among the Northern and Southern Hemispheres. Although several herbaceous species of Solanaceae have been subjected to genome sequencing, thus far, no genome sequences of woody representatives have been available. Here, we sequenced the genomes of 13 perennial woody species of Lycium, with a focus on Lycium barbarum. Integration with other genomes provides clear evidence supporting a whole-genome triplication (WGT) event shared by all hitherto sequenced solanaceous plants, which occurred shortly after the divergence of Solanaceae and Convolvulaceae. We identified new gene families and gene family expansions and contractions that first appeared in Solanaceae. Based on the identification of self-incompatibility related-gene families, we inferred that hybridization hotspots are enriched for genes that might be functioning in gametophytic self-incompatibility pathways in wolfberry. Extremely low expression of LOCULE NUBER (LC) and COLORLESS NON-RIPENING (CNR) orthologous genes during Lycium fruit development and ripening processes suggests functional diversification of these two genes between Lycium and tomato. The existence of additional flowering locus C-like MADS-box genes might correlate with the perennial flowering cycle of Lycium. Differential gene expression involved in the lignin biosynthetic pathway between Lycium and tomato likely illustrates woody and herbaceous differentiation. We also provide evidence that Lycium migrated from Africa into Asia, and subsequently from Asia into North America. Our results provide functional insights into Solanaceae origins, evolution and diversification.

[A new lignan from Euscaphis konishii].

The chemical constituents from the extract of the twigs of Euscaphis konishii with anti-hepatoma activity were investigated, twelve compounds by repeated chromatography with silica gel, Sephadex LH-20 and preparative-HPLC. The structures of the chemical components were elucidated by spectroscopy methods, as konilignan(1),(7R, 8S)-dihydrodehydrodico-niferylalcohol-9-O-ß-D-glucopyranoside(2),illiciumlignan B(3),threo-1-(4-hydroxy-3-methoxyphenyl)-2-[4-(3-hydroxypropyl)-2-methoxyphenoxy]-1,3-panediol(4),erythro-1-(4-hydroxy-3-methoxyphenyl)-2-[4-(3-hydroxypropyl)-2-methoxyphenoxy]-1,3-panediol(5), matairesinol(6), wikstromol(7), isolariciresinol(8),(+)-lyoniresinol(9), 4-ketopinoresinol(10), syringaresin(11), and vladinol D(12). Among them, compound 1 is a new lignan. Compounds 10 and 12 had moderate inhibitory activity on HepG2 cells, with IC_(50) values of 107.12 µmol·L~(-1) and 183.56 µmol·L~(-1), respectively.

Orchid Bsister gene PeMADS28 displays conserved function in ovule integument development.

The ovules and egg cells are well developed to be fertilized at anthesis in many flowering plants. However, ovule development is triggered by pollination in most orchids. In this study, we characterized the function of a Bsister gene, named PeMADS28, isolated from Phalaenopsis equestris, the genome-sequenced orchid. Spatial and temporal expression analysis showed PeMADS28 predominantly expressed in ovules between 32 and 48 days after pollination, which synchronizes with integument development. Subcellular localization and protein-protein interaction analyses revealed that PeMADS28 could form a homodimer as well as heterodimers with D-class and E-class MADS-box proteins. In addition, ectopic expression of PeMADS28 in Arabidopsis thaliana induced small curled rosette leaves, short silique length and few seeds, similar to that with overexpression of other species' Bsister genes in Arabidopsis. Furthermore, complementation test revealed that PeMADS28 could rescue the phenotype of the ABS/TT16 mutant. Together, these results indicate the conserved function of Bsister PeMADS28 associated with ovule integument development in orchid.

Genome-Wide Identification and Co-Expression Analysis of ARF and IAA Family Genes in Euscaphis konishii: Potential Regulators of Triterpenoids and Anthocyanin Biosynthesis.

Euscaphis konishii is an evergreen plant that is widely planted as an industrial crop in Southern China. It produces red fruits with abundant secondary metabolites, giving E. konishii high medicinal and ornamental value. Auxin signaling mediated by members of the AUXIN RESPONSE FACTOR (ARF) and auxin/indole-3-acetic acid (Aux/IAA) protein families plays important roles during plant growth and development. Aux/IAA and ARF genes have been described in many plants but have not yet been described in E. konishii. In this study, we identified 34 EkIAA and 29 EkARF proteins encoded by the E. konishii genome through database searching using HMMER. We also performed a bioinformatic characterization of EkIAA and EkARF genes, including their phylogenetic relationships, gene structures, chromosomal distribution, and cis-element analysis, as well as conserved motifs in the proteins. Our results suggest that EkIAA and EkARF genes have been relatively conserved over evolutionary history. Furthermore, we conducted expression and co-expression analyses of EkIAA and EkARF genes in leaves, branches, and fruits, which identified a subset of seven EkARF genes as potential regulators of triterpenoids and anthocyanin biosynthesis. RT-qPCR, yeast one-hybrid, and transient expression analyses showed that EkARF5.1 can directly interact with auxin response elements and regulate downstream gene expression. Our results may pave the way to elucidating the function of EkIAA and EkARF gene families in E. konishii, laying a foundation for further research on high-yielding industrial products and E. konishii breeding.

The Phoebe genome sheds light on the evolution of magnoliids.

Lauraceae includes the genus Phoebe, and the family is linked to the evolution of magnoliids. We sequenced the genome of Phoebe bournei Nanmu. The assembled genome size was 989.19 Mb, with a contig N50 value of 2.05 Mb. A total of 28,198 protein-coding genes were annotated in P. bournei. Whole-genome duplication (WGD) analysis showed that Lauraceae has experienced two WGD events; the older WGD event occurred just before the divergence of Lauraceae and Magnoliales, and the more recent WGD was shared by all lineages of Lauraceae. The phylogenetic tree showed that magnoliids form a sister clade to monocots and eudicots. We also identified 63 MADS-box genes, including AGL12-like genes that may be related to the regulation of P. bournei roots and FIN219-like genes encoding GH3 proteins, which are involved in photomorphogenesis. SAUR50-like genes involved in light signal-mediated pedicel or stem development were also identified. Four ATMYB46- and three PtrEPSP-homologous genes related to lignin biosynthesis were identified. These genes may be associated with the formation of straight trunks in P. bournei. Overall, the P. bournei reference genome provides insight into the origin, evolution, and diversification of Phoebe and other magnoliids.

[Spatial variability and agglomeration of soil salinity in Minjiang estuary wetland, Southeast China]. / 闽江河口湿地土壤盐分的空间分异与集聚特征.

Understanding the spatial variability and agglomeration of soil salinity is of great significance for the sustainable development of estuarine wetland. Landsat 8 OLI remote sensing image, digital elevation mode and soil surface samples of Minjiang estuary wetland of Fuzhou were used as the data sources. The correlation analysis and principal component analysis were combined to select significant environmental variables and to reduce their dimensions. We analyzed the spatial variability of soil salinity with support vector regression ordinary kriging model (SVROK) and regression kri-ging model (RK), and quantified spatial agglomeration of soil salinity by the spatial autocorrelation analysis. The results showed that three principal components (PCs) extracted by the principal component analysis could explain at least 85% of the total variance in the original dataset and reflected the comprehensive information of vegetation cover, soil properties and topography. Both soil salinity and its residuals were affected by structural factors and random factors. The SVROK model based on principal component (PCs) as input variables can more accurately reflect the spatial variability of soil salinity, with a trend of "higher in the north and lower in the south". The Moran's I of soil salinity was more than 0.5, with significant positive spatial autocorrelation and a higher spatial aggregation degree, displaying the spatial agglomeration characteristics of "high value agglomeration, high value widespread, high value surrounded by low value".

Chromosome-scale assembly of the Kandelia obovata genome.

The mangrove Kandelia obovata (Rhizophoraceae) is an important coastal shelterbelt and landscape tree distributed in tropical and subtropical areas across East Asia and Southeast Asia. Herein, a chromosome-level reference genome of K. obovata based on PacBio, Illumina, and Hi-C data is reported. The high-quality assembled genome size is 177.99 Mb, with a contig N50 value of 5.74 Mb. A large number of contracted gene families and a small number of expanded gene families, as well as a small number of repeated sequences, may account for the small K. obovata genome. We found that K. obovata experienced two whole-genome polyploidization events: one whole-genome duplication shared with other Rhizophoreae and one shared with most eudicots (γ event). We confidently annotated 19,138 protein-coding genes in K. obovata and identified the MADS-box gene class and the RPW8 gene class, which might be related to flowering and resistance to powdery mildew in K. obovata and Rhizophora apiculata, respectively. The reference K. obovata genome described here will be very useful for further molecular elucidation of various traits, the breeding of this coastal shelterbelt species, and evolutionary studies with related taxa.

Comprehensive transcriptome analysis of reference genes for fruit development of Euscaphis konishii.

BACKGROUND: Quantitativereal-time reverse transcriptase polymerase chain reaction is the common method to quantify relative gene expression. Normalizating using reliable genes is critical in correctly interpreting expression data from qRT-PCR. Euscaphis konishii is a medicinal plant with a long history in China, which has various chemical compounds in fruit. However, there is no report describing the selection of reference genes in fruit development of Euscaphis konishii. METHODS: We selected eight candidate reference genes based on RNA-seq database analysis, and ranked expression stability using statistical algorithms GeNorm, NormFinder, BestKeeper and ReFinder. Finally, The nine genes related to the anthocyanin synthesis pathway of Euscaphis konishii were used to verify the suitability of reference gene. RESULTS: The results showed that the stability of EkUBC23, EkCYP38 and EkGAPDH2 was better, and the low expression reference genes (EkUBC23 and EkCYP38) were favourable for quantifying low expression target genes, while the high expression reference gene (EkGAPDH2) was beneficial for quantifying high expression genes. In this study, we present the suitable reference genes for fruit development of Euscaphis konishii based on transcriptome data, our study will contribute to further studies in molecular biology and gene function on Euscaphis konishii and other closely related species.

Total phenolic extract of Euscaphis konishii hayata Pericarp attenuates carbon tetrachloride (CCl4)-induced liver fibrosis in mice.

Liver fibrosis is a crucial pathological process involved in the hepatogenic morbidity and mortality. The pericarp of Euscaphis konishii Hayata is usually used in the cooking soup to improve liver function in Southern China, and high level of phenolic compounds has been found in the E. konishii pericarp. The total phenolic compounds extracted from E. konishii pericarp (TPEP) was obtained by polyamide column chromatograph, and 9 phenolic compounds of TPEP were identified through LC/MS and NMR. TPEP exhibited strong free radicals scavenging activity in vitro, and the chronic CCl4-induced liver fibrosis mice were established to elucidate the hepatoprotective mechanism of TPEP in vivo. TPEP treatment (50, 100 and 200 mg/kg) ameliorated the oxidative stress, immune dysfunction, inflammatory response and hepatic fibrosis induced by CCl4 introduction, alleviated the histopathological alteration and hepatocyte apoptosis in the liver tissue. Pretreatment with TPEP suppressed the activation of NADPH oxidase (NOX) signaling to attenuate oxidative stress in the liver tissue. TPEP administration inhibited the translocation of NF-κB into the nucleus to prevent the expression of downstream proinflammatory cytokines. TPEP treatment downregulated the activation of TGF-ß/Smad pathway, and facilitated the degradation of extracellular matrix through enhancing matrix metalloproteinases (MMPs) activity and decreasing the expression of matrix metalloproteinase inhibitors (TIMPs). In conclusion, TPEP inhibited CCl4-induced hepatic fibrosis through its antioxidant, anti-inflammatory and anti-fibrotic activities.

The complete chloroplast genome sequence of Tapiscia sinensis (Staphyleaceae).

Tapiscia sinensis, belong to Tapisciaceae, is endangered tree endemic to China. Here, we provide the complete plastid genomic data of T. sinensis with the aim of providing data for future conservation efforts research and revealing its phylogenetic position. The complete chloroplast sequence is 161,093 bp, including a large single-copy (LSC) region of 87,782 bp, a small single-copy (SSC) region of 18,517 bp, a pair of invert repeats (IR) regions of 27,387 bp. Plastid genome contains 131 genes, 85 protein-coding genes, 38 tRNA genes, and eight rRNA genes. Phylogenetic analysis base on 19 plastid genomes indicates that T. sinensis located Malvids branch, and is more closely related to the species of the order Sapindales than those of the order Malvales.

The complete chloroplast genome sequence of Turpinia montana (Staphyleaceae).

The Rosids are characterized by remarkable morphological and ecological diversity. Here, we provide the completed plasmid genome of Turpinia montana. The complete chloroplast size of T. montana is 160,111 bp, including a large single-copy (LSC) region of 89,631 bp, a small single-copy (SSC) region of 18,247 bp, a pair of invert repeats (IR) regions of 26,120 bp. Plastid genome contains 131 genes, 86 protein-coding genes, 37 tRNA genes, and eight rRNA genes. Phylogenetic analysis base on 23 plastid genomes indicates that T. montanas is clustered with the plants of the Euscaphis japonica and Staphylea bumalda.

Multivariate analysis reveals phenotypic diversity of Euscaphis japonica population.

Fruit traits affect population genetic diversity by affecting seed protection and dispersal strategies, thereby comprising important components of phenotypic variation. Understanding of the phenotypic variation is an indispensable first step for developing breeding strategies. However, little information is known about the genetic variation in E. japonica-a monotypic species with abundant phenotypes that is mainly distributed in southern China. In this study, we evaluated the phenotypic diversity of 67 E. japonica using 23 phenotypic traits. Our results showed that the Shannon-Wiener (I) index of qualitative traits ranged from 0.55 to 1.26, and the color traits had a relatively high I. The average coefficient of variation of compound leaf traits (14.74%) was higher than that of fruit and seed traits (12.77% and 11.47%, respectively). Principal component analysis also showed that compound leaf and fruit traits were important components of total variation. Furthermore, correlation analysis revealed a significant difference in elevation and fruit color, irregular ribs, leaf margin and texture. The F value within populations was smaller than among populations, indicating the variation in phenotypic traits among populations was much greater than within populations. Dehua and Zunyi populations had the highest coefficients of variation, whereas Wenzhou population had the smallest-which may be attributed to habitat destruction. According to Q-type clustering, 67 samples clustered into four groups, with those having similar phenotypes clustering into the same group. In general, leaf and fruit traits had abundant phenotypic diversity, representing the main sources of phenotypic variation. Combined with clustering results and field surveys, this study suggests that the phenotypes of E. japonica are classified into two main categories: The deciduous E. japonica present at high altitudes; and the evergreen E. japonica present at low altitudes. Excavating E. japonica variations provides a theoretical reference for its classification and diversity, and is of great significance for planning genetic resources and establishing conservation strategies.

Protective Effect of the Total Triterpenes of Euscaphis konishii Hayata Pericarp on Bacillus Calmette-Guérin Plus Lipopolysaccharide-Induced Liver Injury.

BACKGROUND: Liver injury has been recognized as a primary cause of hepatic morbidity and mortality. Euscaphis konishii Hayata, also called Euscaphis fukienensis Hsu, is usually used as a detumescent and analgesic agent to improve liver function in South China, but its mechanism of action and chemical composition are unclear. OBJECTIVE: The main aim of the study was to investigate the constituent and potential hepatoprotective mechanism of the total triterpenes of E. konishii pericarp (TTEP). METHODS: The constituent of TTEP was analyzed by a series of silica gel column to get single compounds and then identified by NMR and MS. In vitro assays were conducted to test the free radical scavenging activity of TTEP. The BCG/LPS-induced immunological livery injury mice model was established to clarify the hepatoprotective effect of TTEP in vivo. RESULTS: 8 pentacyclic triterpene acids were separated and identified by NMR and MS. TTEP treatment (50, 100, and 200 mg/Kg) improved the immune function of the BCG/LPS-infected mice, dose-dependently alleviated the BCG/LPS-induced inflammation and oxidative stress, and ameliorated the hepatocyte apoptosis in the liver tissue. CONCLUSION: The pericarp of E. konishii may be further considered as a potent natural food for liver disease treatment.

The complete chloroplast genome sequence of Euscaphis japonica (Staphyleaceae).

Euscaphis japonica is not only an ideal ornamental plant, but also a traditional medicinal plant, which is an extremely valuable species to study. We determined the complete chloroplast genome sequence for E. japonica using Illumina sequencing technology. The complete chloroplast sequence is 160,467 bp in length, including a large single-copy (LSC) region of 88,716 bp, a small single-copy (SSC) region of 18,614 bp, and a pair of invert repeats (IR) regions of 26,568 bp. Plastid genome contains 142 genes, 46 protein-coding genes, 39 tRNA genes, and eight rRNA genes. Phylogenetic analysis base on 14 chloroplast genomes indicates that E. japonica forms an isolated clade and sisters to Glycosmis-Gossypium clade with strong support.