Genomic variation, environmental adaptation and feralization in ramie, an ancient fiber crop.

Feralization is an important evolutionary process, but the mechanisms behind it remain poorly understood. Here, we use the ancient fiber crop, ramie (Boehmeria nivea (L.) Gaudich.) as a model to investigate genomic changes associated with both domestication and fertilization. We first produced a chromosome-scale de novo genome assembly of feral ramie and investigated structural variations between feral and domesticated ramie genomes. Next, 915 accessions from 20 countries were gathered, comprising cultivars, major landraces, feral populations and wild progenitor. Based on whole genome resequencing of these accessions, the most comprehensive ramie genomic variation map to date was constructed. Phylogenetic, demographic, and admixture signal detection analyses indicate that feral ramie is of exoferal or exo-endo origin, i.e., descended from hybridization between domesticated ramie and wild progenitor or ancient landraces. Feral ramie has greater genetic diversity than wild or domesticated ramie, and genomic regions affected by natural selection during feralization are different from those under selection during domestication. Ecological analyses showed that feral and domesticated ramie have similar ecological niches which are substantially different from the niche of the wild progenitor, and three environmental variables were associated with habitat-specific adaptation in feral ramie. Our findings advance our understanding of feralization, providing a scientific basis for the excavation of new crop germplasm resources and offering novel insights into the evolution of feralization in nature.

Genetic analyses of ancient tea trees provide insights into the breeding history and dissemination of Chinese Assam tea (Camellia sinensis var. assamica).

Chinese Assam tea (Camellia sinensis var. assamica) is an important tea crop with a long history of cultivation in Yunnan, China. Despite its potential value as a genetic resource, its genetic diversity and domestication/breeding history remain unclear. To address this issue, we genotyped 469 ancient tea plant trees representing 26 C. sinensis var. assamica populations, plus two of its wild relatives (six and three populations of C. taliensis and C. crassicolumna, respectively) using 16 nuclear microsatellite loci. Results showed that Chinese Assam tea has a relatively high, but comparatively lower gene diversity (HS = 0.638) than the wild relative C. crassicolumna (HS = 0.658). Clustering in STRUCTURE indicated that Chinese Assam tea and its two wild relatives formed distinct genetic groups, with considerable interspecific introgression. The Chinese Assam tea accessions clustered into three gene pools, corresponding well with their geographic distribution. However, NewHybrids analysis indicated that 68.48% of ancient Chinese Assam tea plants from Xishuangbanna were genetic intermediates between the Puer and Lincang gene pools. In addition, 10% of the ancient Chinese Assam tea individuals were found to be hybrids between Chinese Assam tea and C. taliensis. Our results suggest that Chinese Assam tea was domesticated separately in three gene pools (Puer, Lincang and Xishuangbanna) in the Mekong River valley and that the hybrids were subsequently selected during the domestication process. Although the domestication history of Chinese Assam tea in southwestern Yunnan remains complex, our results will help to identify valuable genetic resources that may be useful in future tea breeding programs.

Deep genome skimming reveals the hybrid origin of Pseudosasa gracilis (Poaceae: Bambusoideae).

Pseudosasa gracilis (Poaceae: Bambusoideae) is a temperate woody bamboo species endemic to South-central China with a narrow distribution. Previous phylogenetic studies revealed an unexpected, isolated phylogenetic position of Ps. gracilis. Here we conducted phylogenomic analysis by sampling populations of Ps. gracilis and its sympatric species Ps. nanunica and Sinosasa polytricha reflecting different genomic signals, by deep genome skimming. Integrating molecular evidence from chloroplast genes and genome-wide SNPs, we deciphered the phylogenetic relationships of Ps. gracilis. Both plastid and nuclear data indicate that Ps. gracilis is more closely related to Sinosasa, which is discordant with the taxonomic treatment. To further explore this molecular-morphological conflict, we screened 411 "perfect-copy" syntenic genes to reconstruct phylogenies using both the concatenation and coalescent methods. We observed extensive discordance between gene trees and the putative species tree. A significant hybridization event was detected based on 411 genes from the D subgenome, showing Ps. gracilis was a hybrid descendant between Sinosasa longiligulata and Ps. nanunica, with 63.56% and 36.44% inheritance probabilities of each parent. Moreover, introgression events were detected in the C subgenome between Ps. gracilis and S. polytricha in the same distribution region. Our findings suggest that sympatric hybridization and introgression play a crucial role in the origin of Ps. gracilis. By providing an empirical example of bamboo of hybrid origin using comprehensive analyses based on genomic data from different inheritance systems and morphological characters, our study represents a step forward in understanding of reticulate evolution of bamboos.

Niche convergence and biogeographic history shape elevational tree community assembly in a subtropical mountain forest.

Niche convergence or conservatism have been proposed as essential mechanisms underlying elevational plant community assembly in tropical mountain ecosystems. Subtropical mountains, compared to tropical mountains, are likely to be shaped by a mixing of different geographic affinities of species and remain somehow unclear. Here, we used 31 0.1-ha permanent plots distributed in subtropical forests on the eastern and western aspects of the Gaoligong Mountains, southwest China between 1498 m and 3204 m a.sl. to evaluate how niche-based and biogeographic processes shape tree community assembly along elevational gradients. We analyzed the elevational patterns of taxonomic, phylogenetic and functional diversity, as well as of individual traits, and assessed the relative importance of environmental effects on these diversity measures. We then classified tree species as being either tropical affiliated or temperate affiliated and estimated their contribution to the composition of biogeographic affinities. Species richness decreased with elevation, and species composition showed apparent turnover across the aspects and elevations. Most traits exhibited convergent patterns across the entire elevational gradient. Phylogenetic and functional diversity showed opposing patterns, with phylogenetic diversity increasing and functional diversity decreasing with elevation. Soil nutrients, especially phosphorus and nitrogen, appeared to be the main abiotic variables driving the elevational diversity patterns. Communities at lower elevations were occupied by tropical genera, while highlands contained species of tropical and temperate biogeographic affinities. Moreover, the high phylogenetic diversity at high elevations were likely due to differences in evolutionary history between temperate and tropical species. Our results highlight the importance of niche convergence of tropical species and the legacy of biogeographic history on the composition and structure of subtropical mountain forests. Furthermore, limited soil phosphorus caused traits divergence and the partitioning for different forms of phosphorus may explain the high biodiversity found in phosphorus-limited subtropical forests.

Unprecedented variation pattern of plastid genomes and the potential role in adaptive evolution in Poales.

BACKGROUND: The plastid is the photosynthetic organelle in plant cell, and the plastid genomes (plastomes) are generally conserved in evolution. As one of the most economically and ecologically important order of angiosperms, Poales was previously documented to exhibit great plastomic variation as an order of photoautotrophic plants. RESULTS: We acquired 93 plastomes, representing all the 16 families and 5 major clades of Poales to reveal the extent of their variation and evolutionary pattern. Extensive variation including the largest one in monocots with 225,293 bp in size, heterogeneous GC content, and a wide variety of gene duplication and loss were revealed. Moreover, rare occurrences of three inverted repeat (IR) copies in angiosperms and one IR loss were observed, accompanied by short IR (sIR) and small direct repeat (DR). Widespread structural heteroplasmy, diversified inversions, and unusual genomic rearrangements all appeared in Poales, occasionally within a single species. Extensive repeats in the plastomes were found to be positively correlated with the observed inversions and rearrangements. The variation all showed a "small-large-moderate" trend along the evolution of Poales, as well as for the sequence substitution rate. Finally, we found some positively selected genes, mainly in C4 lineages, while the closely related lineages of those experiencing gene loss tended to have undergone more relaxed purifying selection. CONCLUSIONS: The variation of plastomes in Poales may be related to its successful diversification into diverse habitats and multiple photosynthetic pathway transitions. Our order-scale analyses revealed unusual evolutionary scenarios for plastomes in the photoautotrophic order of Poales and provided new insights into the plastome evolution in angiosperms as a whole.

Variations and reduction of plastome are associated with the evolution of parasitism in Convolvulaceae.

Parasitic lifestyle can often relax the constraint on the plastome, leading to gene pseudogenization and loss, and resulting in diverse genomic structures and rampant genome degradation. Although several plastomes of parasitic Cuscuta have  been reported, the evolution of parasitism in the family Convolvulaceae which is linked to structural variations and reduction of plastome has not been well investigated. In this study, we assembled and collected 40 plastid genomes belonging to 23 species representing four subgenera of Cuscuta and ten species of autotrophic Convolvulaceae. Our findings revealed nine types of structural variations and six types of inverted repeat (IR) boundary variations in the plastome of Convolvulaceae spp. These structural variations were associated with the shift of parasitic lifestyle, and IR boundary shift, as well as the abundance of long repeats. Overall, the degradation of Cuscuta plastome proceeded gradually, with one clade exhibiting an accelerated degradation rate. We observed five stages of gene loss in Cuscuta, including NAD(P)H complex → PEP complex → Photosynthesis-related → Ribosomal protein subunits → ATP synthase complex. Based on our results, we speculated that the shift of parasitic lifestyle in early divergent time promoted relaxed selection on plastomes, leading to the accumulation of microvariations, which ultimately resulted in the plastome reduction. This study provides new evidence towards a better understanding of plastomic evolution, variation, and reduction in the genus Cuscuta.

Phylogenomic profiles of whole-genome duplications in Poaceae and landscape of differential duplicate retention and losses among major Poaceae lineages.

Poaceae members shared a whole-genome duplication called rho. However, little is known about the evolutionary pattern of the rho-derived duplicates among Poaceae lineages and implications in adaptive evolution. Here we present phylogenomic/phylotranscriptomic analyses of 363 grasses covering all 12 subfamilies and report nine previously unknown whole-genome duplications. Furthermore, duplications from a single whole-genome duplication were mapped to multiple nodes on the species phylogeny; a whole-genome duplication was likely shared by woody bamboos with possible gene flow from herbaceous bamboos; and recent paralogues of a tetraploid Oryza are implicated in tolerance of seawater submergence. Moreover, rho duplicates showing differential retention among subfamilies include those with functions in environmental adaptations or morphogenesis, including ACOT for aquatic environments (Oryzoideae), CK2ß for cold responses (Pooideae), SPIRAL1 for rapid cell elongation (Bambusoideae), and PAI1 for drought/cold responses (Panicoideae). This study presents a Poaceae whole-genome duplication profile with evidence for multiple evolutionary mechanisms that contribute to gene retention and losses.

Evolution of the FLOWERING LOCUS T-like genes in angiosperms: a core-Lamiales-specific diversification.

Plant life-history is determined by two transitions, the germination and the flowering times, in which the phosphatidylethanolamine-binding proteins (PEBP) FLOWERING LOCUS T (FT) and TERMINAL FLOWER1 (TFL1) play key regulatory roles. Compared to the highly conserved TFL1-likes, FT-like genes vary in copy numbers significantly in gymnosperms and monocots of the angiosperms, while sporadic duplications can be observed in eudicots. Here, via a systematic analysis of the PEBPs in angiosperms with a special focus on twelve representative species featuring high-quality genomes in the Lamiales order, we identified a successive lineage-specific but systematic expansion of FT-like genes in the families of core Lamiales. The first expansion event generated FT1-likes mainly via a core-Lamiales-specific whole-genome-duplication (cL-WGD), while on the other hand, a likely random duplication produced the FT2-likes in the lineages containing Scrophulariaceae and rest of the core Lamiales. Both FT1- and FT2-like genes were further amplified tandemly in some families. These expanded FT-likes featured highly diverged expression patterns and structural variation, indicating functional diversification. Intriguingly, some core Lamiales contained the relict MOTHER OF FT AND TFL1 like 2 (MFT2) that likely expanded in the common ancestor of angiosperms. Our data showcase the highly dynamic lineage-specific expansion of the FT-like genes, thus provide important and fresh evolutionary insights into the gene-regulatory-network underpinning flowering time diversity in Lamiales, and more generally, in angiosperms.

Dryopterisjinpingensis, a critically endangered diploid new species of Dryopteridaceae from Yunnan, China.

Dryopterisjinpingensis, a new species of diploid, sexually reproductive ferns of Dryopteridaceae from Yunnan, southwestern China, is described and illustrated. Morphologically, D.jinpingensis is similar to D.gaoligongensis but unique in elongated lanceolate laminae, sessile or subsessile pinna stalks, and overlapping membranous scales adnate to stipe base. Phylogenetic analyses based on both plastome and the nuclear AK1 gene sequences showed that D.jinpingensis is sister to D.gaoligongensis. A detailed taxonomic description with line drawings is provided, and its conservation status is evaluated to be critically endangered.

Genomes of Meniocus linifolius and Tetracme quadricornis reveal the ancestral karyotype and genomic features of core Brassicaceae.

Brassicaceae represents an important plant family from both a scientific and economic perspective. However, genomic features related to the early diversification of this family have not been fully characterized, especially upon the uplift of the Tibetan Plateau, which was followed by increasing aridity in the Asian interior, intensifying monsoons in Eastern Asia, and significantly fluctuating daily temperatures. Here, we reveal the genomic architecture that accompanied early Brassicaceae diversification by analyzing two high-quality chromosome-level genomes for Meniocus linifolius (Arabodae; clade D) and Tetracme quadricornis (Hesperodae; clade E), together with genomes representing all major Brassicaceae clades and the basal Aethionemeae. We reconstructed an ancestral core Brassicaceae karyotype (CBK) containing 9 pseudochromosomes with 65 conserved syntenic genomic blocks and identified 9702 conserved genes in Brassicaceae. We detected pervasive conflicting phylogenomic signals accompanied by widespread ancient hybridization events, which correlate well with the early divergence of core Brassicaceae. We identified a successive Brassicaceae-specific expansion of the class I TREHALOSE-6-PHOSPHATE SYNTHASE 1 (TPS1) gene family, which encodes enzymes with essential regulatory roles in flowering time and embryo development. The TPS1s were mainly randomly amplified, followed by expression divergence. Our results provide fresh insights into historical genomic features coupled with Brassicaceae evolution and offer a potential model for broad-scale studies of adaptive radiation under an ever-changing environment.

Genome assemblies of 11 bamboo species highlight diversification induced by dynamic subgenome dominance.

Polyploidy (genome duplication) is a pivotal force in evolution. However, the interactions between parental genomes in a polyploid nucleus, frequently involving subgenome dominance, are poorly understood. Here we showcase analyses of a bamboo system (Poaceae: Bambusoideae) comprising a series of lineages from diploid (herbaceous) to tetraploid and hexaploid (woody), with 11 chromosome-level de novo genome assemblies and 476 transcriptome samples. We find that woody bamboo subgenomes exhibit stunning karyotype stability, with parallel subgenome dominance in the two tetraploid clades and a gradual shift of dominance in the hexaploid clade. Allopolyploidization and subgenome dominance have shaped the evolution of tree-like lignified culms, rapid growth and synchronous flowering characteristic of woody bamboos as large grasses. Our work provides insights into genome dominance in a remarkable polyploid system, including its dependence on genomic context and its ability to switch which subgenomes are dominant over evolutionary time.

Positive selection and relaxed purifying selection contribute to rapid evolution of male-biased genes in a dioecious flowering plant.

Sex-biased genes offer insights into the evolution of sexual dimorphism. Sex-biased genes, especially those with male bias, show elevated evolutionary rates of protein sequences driven by positive selection and relaxed purifying selection in animals. Although rapid sequence evolution of sex-biased genes and evolutionary forces have been investigated in animals and brown algae, less is known about evolutionary forces in dioecious angiosperms. In this study, we separately compared the expression of sex-biased genes between female and male floral buds and between female and male flowers at anthesis in dioecious Trichosanthes pilosa (Cucurbitaceae). In floral buds, sex-biased gene expression was pervasive, and had significantly different roles in sexual dimorphism such as physiology. We observed higher rates of sequence evolution for male-biased genes in floral buds compared to female-biased and unbiased genes. Male-biased genes under positive selection were mainly associated with functions to abiotic stress and immune responses, suggesting that high evolutionary rates are driven by adaptive evolution. Additionally, relaxed purifying selection may contribute to accelerated evolution in male-biased genes generated by gene duplication. Our findings, for the first time in angiosperms, suggest evident rapid evolution of male-biased genes, advance our understanding of the patterns and forces driving the evolution of sexual dimorphism in dioecious plants.

Remarkable mitochondrial genome heterogeneity in Meniocus linifolius (Brassicaceae).

KEY MESSAGE: Detailed analyses of 16 genomes identified a remarkable acceleration of mutation rate, hence mitochondrial sequence and structural heterogeneity, in Meniocus linifolius (Brassicaceae). The powerhouse, mitochondria, in plants feature high levels of structural variation, while the encoded genes are normally conserved. However, the substitution rates and spectra of mitochondria DNA within the Brassicaceae, a family with substantial scientific and economic importance, have not been adequately deciphered. Here, by analyzing three newly assembled and 13 known mitochondrial genomes (mitogenomes), we report the highly variable genome structure and mutation rates in Brassicaceae. The genome sizes and GC contents are 196,604 bp and 46.83%, 288,122 bp and 44.79%, and 287,054 bp and 44.93%, for Meniocus linifolius (Mli), Crucihimalaya lasiocarpa (Cla), and Lepidium sativum (Lsa), respectively. In total, 29, 33, and 34 protein-coding genes (PCGs) and 14, 18, and 18 tRNAs are annotated for Mli, Cla, and Lsa, respectively, while all mitogenomes contain one complete circular molecule with three rRNAs and abundant RNA editing sites. The Mli mitogenome features four conformations likely mediated by the two pairs of long repeats, while at the same time seems to have an unusual evolutionary history due to higher GC content, loss of more genes and sequences, but having more repeats and plastid DNA insertions. Corroborating with these, an ambiguous phylogenetic position with long branch length and elevated synonymous substitution rate in nearly all PCGs are observed for Mli. Taken together, our results reveal a high level of mitogenome heterogeneity at the family level and provide valuable resources for further understanding the evolutionary pattern of organelle genomes in Brassicaceae.

Parallel evolution of morphological and genomic selfing syndromes accompany the breakdown of heterostyly.

Evolutionary transitions from outcrossing to selfing in flowering plants have convergent morphological and genomic signatures and can involve parallel evolution within related lineages. Adaptive evolution of morphological traits is often assumed to evolve faster than nonadaptive features of the genomic selfing syndrome. We investigated phenotypic and genomic changes associated with transitions from distyly to homostyly in the Primula oreodoxa complex. We determined whether the transition to selfing occurred more than once and investigated stages in the evolution of morphological and genomic selfing syndromes using 22 floral traits and both nuclear and plastid genomic data from 25 populations. Two independent transitions were detected representing an earlier and a more recently derived selfing lineage. The older lineage exhibited classic features of the morphological and genomic selfing syndrome. Although features of both selfing syndromes were less developed in the younger selfing lineage, they exhibited parallel development with the older selfing lineage. This finding contrasts with the prediction that some genomic changes should lag behind adaptive changes to morphological traits. Our findings highlight the value of comparative studies on the timing and extent of transitions from outcrossing to selfing between related lineages for investigating the tempo of morphological and molecular evolution.

Anthropogenic disturbance driving population decline of a dominant tree in East Asia evergreen broadleaved forests over the last 11,000 years.

Current biodiversity loss is generally considered to have been caused by anthropogenic disturbance, but it is unclear when anthropogenic activities began to affect biodiversity loss. One hypothesis suggests it began with the Industrial Revolution, whereas others propose that anthropogenic disturbance has been associated with biodiversity decline since the early Holocene. To test these hypotheses, we examined the unique vegetation of evergreen broadleaved forests (EBLFs) in East Asia, where humans have affected landscapes since the early Holocene. We adopted a genomic approach to infer the demographic history of a dominant tree (Litsea elongata) of EBLFs. We used Holocene temperature and anthropogenic disturbance factors to calculate the correlation between these variables and the historical effective population size of L. elongata with Spearman statistics and integrated the maximum-entropy niche model to determine the impact of climate change and anthropogenic disturbance on fluctuation in its effective population size. We identified 9 well-defined geographic clades for the populations of L. elongata. Based on the estimated historical population sizes of these clades, all the populations contracted, indicating persistent population decline over the last 11,000 years. Demographic history of L. elongata and human population change, change in cropland use, and change in irrigated rice area were significantly negatively correlated, whereas climate change in the Holocene was not correlated with demographic history. Our results support the early human impact hypothesis and provide comprehensive evidence that early anthropogenic disturbance may contribute to the current biodiversity crisis in East Asia.

Distinctive plastome evolution in carnivorous angiosperms.

BACKGROUND: Independent origins of carnivory in multiple angiosperm families are fabulous examples of convergent evolution using a diverse array of life forms and habitats. Previous studies have indicated that carnivorous plants have distinct evolutionary trajectories of plastid genome (plastome) compared to their non-carnivorous relatives, yet the extent and general characteristics remain elusive. RESULTS: We compared plastomes from 9 out of 13 carnivorous families and their non-carnivorous relatives to assess carnivory-associated evolutionary patterns. We identified inversions in all sampled Droseraceae species and four species of Utricularia, Pinguicula, Darlingtonia and Triphyophyllum. A few carnivores showed distinct shifts in inverted repeat boundaries and the overall repeat contents. Many ndh genes, along with some other genes, were independently lost in several carnivorous lineages. We detected significant substitution rate variations in most sampled carnivorous lineages. A significant overall substitution rate acceleration characterizes the two largest carnivorous lineages of Droseraceae and Lentibulariaceae. We also observe moderate substitution rates acceleration in many genes of Cephalotus follicularis, Roridula gorgonias, and Drosophyllum lusitanicum. However, only a few genes exhibit significant relaxed selection. CONCLUSION: Our results indicate that the carnivory of plants have different effects on plastome evolution across carnivorous lineages. The complex mechanism under carnivorous habitats may have resulted in distinctive plastome evolution with conserved plastome in the Brocchinia hechtioides to strongly reconfigured plastomes structures in Droseraceae. Organic carbon obtained from prey and the efficiency of utilizing prey-derived nutrients might constitute possible explanation.

Five new and noteworthy species of Epidendroideae (Orchidaceae) from southwestern China based on morphological and phylogenetic evidence.

Five new orchid species from southwestern China's Yunnan Province and the Tibetan Autonomous Region, Neottialihengiae, Neottiachawalongensis, Papilionanthemotuoensis, Gastrochiluslihengiae, and Gastrochilusbernhardtianus, are described and illustrated. To confirm their identities, and to resolve phylogenetic relationships, we sequenced the complete plastomes of these taxa with their congeneric species, adding new plastomes of three Neottia species, two Papilionanthe species and nine Gastrochilus species. Combined with published plastid sequences, our well-resolved phylogeny supported the alliance of N.lihengiae with the the N.grandiflora + N.pinetorum clade. Neottiachawalongensis is now sister to N.alternifolia, while P.motuoensis is closely related to P.subulata + P.teres. Conversely, phylogenetic analyses based on complete plastomes and plastid sequences showed inconsistent relationships among taxa in the genus Gastrochilus, but the two new species, G.lihengiae and G.bernhardtianus were supported by all datasets.

Dynamics of imprinted genes and their epigenetic mechanisms in castor bean seed with persistent endosperm.

Genomic imprinting refers to parent-of-origin-dependent gene expression and primarily occurs in the endosperm of flowering plants, but its functions and epigenetic mechanisms remain to be elucidated in eudicots. Castor bean, a eudicot with large and persistent endosperm, provides an excellent system for studying the imprinting. Here, we identified 131 imprinted genes in developing endosperms and endosperm at seed germination phase of castor bean, involving into the endosperm development, accumulation of storage compounds and specially seed germination. Our results showed that the transcriptional repression of maternal allele of DNA METHYLTRANSFERASE 1 (MET1) may be required for maternal genome demethylation in the endosperm. DNA methylation analysis showed that only a small fraction of imprinted genes was associated with allele-specific DNA methylation, and most of them were closely associated with constitutively unmethylated regions (UMRs), suggesting a limited role for DNA methylation in controlling genomic imprinting. Instead, histone modifications can be asymmetrically deposited in maternal and paternal genomes in a DNA methylation-independent manner to control expression of most imprinted genes. These results expanded our understanding of the occurrence and biological functions of imprinted genes and showed the evolutionary flexibility of the imprinting machinery and mechanisms in plants.