Conservation genetics of endangered Trillium govanianum Wall. ex D. Don - A pharmaceutically prized medicinal plant from the Himalaya and implications for species recovery.

Understanding the genetic diversity and population structure of pharmaceutically important endangered plant species is crucial for their conservation and sustainable use. Despite the continuous population decline in Trillium govanianum Wall. ex D. Don, a highly prized medicinal plant endemic to the Himalaya, information regarding its conservation genetics has been lacking. Here, we employed a conservation genetics approach to investigate how drastically declining populations in natural habitats impact population genetic diversity and structure of this endangered species across the Kashmir Himalaya. We used Start codon targeted (SCoT) and Simple sequence repeat (SSR) markers to assess the intra- and inter-population genetic variation in seven sites across the study region. Based on these markers, we found a very low genetic diversity in T. govanianum populations. Very low levels of observed heterozygosity (Ho = 0.000) and that expected (He = 0.064) in the populations indicate high heterozygote deficiency and high levels of inbreeding depression (FIS = 1.000). A high genetic differentiation was observed among the populations for both SCoT (Gst = 0.719) and SSR (Fst = 0.707) markers. Both the markers showed low gene flow, SCoT (Nm = 0.195) and SSR (Nm = 0.119), depicting high among-population variation than within-population variation. Analysis of molecular variance also indicated a higher genetic variation between the populations than within populations. We also observed a significant positive correlation between genetic divergence and geographical distance, indicating that genetic differentiation in T. govanianum follows a pattern of isolation by distance. Bayesian structure and cluster analysis grouped the populations according to their geographical proximity. Further, redundancy analysis (RDA) revealed the presence of one polymorphic locus for each marker with high discriminatory power. Overall, our findings reveal a very low genetic diversity, high levels of inbreeding, and high genetic differentiation among the populations; likely resulting from habitat fragmentation, population isolation, bottleneck effect, low gene flow, and predominantly asexual reproduction currently operative in the species. Finally, based on the insights gained, we discuss the potential implications of our findings in guiding species recovery and habitat rehabilitation of T. govanianum in the Himalaya with conservation lessons for elsewhere in the world.

Spatial transcriptome analysis provides insights of key gene(s) involved in steroidal saponin biosynthesis in medicinally important herb Trillium govanianum.

Trillium govanianum, an endangered medicinal herb native to the Himalaya, is less studied at the molecular level due to the non-availability of genomic resources. To facilitate the basic understanding of the key genes and regulatory mechanism of pharmaceutically important biosynthesis pathways, first spatial transcriptome sequencing of T. govanianum was performed. 151,622,376 (~11.5 Gb) high quality reads obtained using paired-end Illumina sequencing were de novo assembled into 69,174 transcripts. Functional annotation with multiple public databases identified array of genes involved in steroidal saponin biosynthesis and other secondary metabolite pathways including brassinosteroid, carotenoid, diterpenoid, flavonoid, phenylpropanoid, steroid and terpenoid backbone biosynthesis, and important TF families (bHLH, MYB related, NAC, FAR1, bZIP, B3 and WRKY). Differentially expressed large number of transcripts, together with CYPs and UGTs suggests involvement of these candidates in tissue specific expression. Combined transcriptome and expression analysis revealed that leaf and fruit tissues are the main site of steroidal saponin biosynthesis. In conclusion, comprehensive genomic dataset created in the current study will serve as a resource for identification of potential candidates for genetic manipulation of targeted bioactive metabolites and also contribute for development of functionally relevant molecular marker resource to expedite molecular breeding and conservation efforts in T. govanianum.

Analysis of the floral MADS-box genes from monocotyledonous Trilliaceae species indicates the involvement of SEPALLATA3-like genes in sepal-petal differentiation.

The evolution of greenish sepals from petaloid outer tepals has occurred repeatedly in various lineages of non-grass monocots. Studies in distinct monocot species showed that the evolution of sepals could be explained by the ABC model; for example, the defect of B-class function in the outermost whorl was linked to the evolution of sepals. Here, floral MADS-box genes from three sepal-bearing monocotyledonous Trilliaceae species, Trillium camschatcense, Paris verticillata, and Kinugasa japonica were examined. Unexpectedly, expression of not only A- but also B-class genes was detected in the sepals of all three species. Although the E-class gene is generally expressed across all floral whorls, no expression was detected in sepals in the three species examined here. Overexpression of the E-class SEPALLATA3-like gene from T. camschatcense (TcamSEP) in Arabidopsis thaliana produced phenotypes identical to those reported for orthologs in other monocots. Additionally, yeast hybrid experiments indicated that TcamSEP could form a higher-order complex with an endogenous heterodimer of B-class APETALA3/DEFICIENS-like (TcamDEF) and PISTILLATA/GLOBOSA-like (TcamGLO) proteins. These results suggest a conserved role for Trilliaceae SEPALLATA3-like genes in functionalization of the B-class genes, and that a lack of SEPALLATA3-like gene expression in the outermost whorl may be related to the formation of greenish sepals.

Evolution in stage-structured populations.

For many organisms, stage is a better predictor of demographic rates than age. Yet no general theoretical framework exists for understanding or predicting evolution in stage-structured populations. Here, we provide a general modeling approach that can be used to predict evolution and demography of stage-structured populations. This advances our ability to understand evolution in stage-structured populations to a level previously available only for populations structured by age. We use this framework to provide the first rigorous proof that Lande's theorem, which relates adaptive evolution to population growth, applies to stage-classified populations, assuming only normality and that evolution is slow relative to population dynamics. We extend this theorem to allow for different means or variances among stages. Our next major result is the formulation of Price's theorem, a fundamental law of evolution, for stage-structured populations. In addition, we use data from Trillium grandiflorum to demonstrate how our models can be applied to a real-world population and thereby show their practical potential to generate accurate projections of evolutionary and population dynamics. Finally, we use our framework to compare rates of evolution in age- versus stage-structured populations, which shows how our methods can yield biological insights about evolution in stage-structured populations.

A heterogeneity test for fine-scale genetic structure.

For organisms with limited vagility and/or occupying patchy habitats, we often encounter nonrandom patterns of genetic affinity over relatively small spatial scales, labelled fine-scale genetic structure. Both the extent and decay rate of that pattern can be expected to depend on numerous interesting demographic, ecological, historical, and mating system factors, and it would be useful to be able to compare different situations. There is, however, no heterogeneity test currently available for fine-scale genetic structure that would provide us with any guidance on whether the differences we encounter are statistically credible. Here, we develop a general nonparametric heterogeneity test, elaborating on standard autocorrelation methods for pairs of individuals. We first develop a 'pooled within-population' correlogram, where the distance classes (lags) can be defined as functions of distance. Using that pooled correlogram as our null-hypothesis reference frame, we then develop a heterogeneity test of the autocorrelations among different populations, lag-by-lag. From these single-lag tests, we construct an analogous test of heterogeneity for multilag correlograms. We illustrate with a pair of biological examples, one involving the Australian bush rat, the other involving toadshade trillium. The Australian bush rat has limited vagility, and sometimes occupies patchy habitat. We show that the autocorrelation pattern diverges somewhat between continuous and patchy habitat types. For toadshade trillium, clonal replication in Piedmont populations substantially increases autocorrelation for short lags, but clonal replication is less pronounced in mountain populations. Removal of clonal replicates reduces the autocorrelation for short lags and reverses the sign of the difference between mountain and Piedmont correlograms.

Evolutionary dynamics as a component of stage-structured matrix models: an example using Trillium grandiflorum.

Evolution by natural selection improves fitness and may therefore influence population trajectories. Demographic matrix models are often employed in conservation studies to project population dynamics, but such analyses have not incorporated evolutionary dynamics. We project evolutionarily informed population trajectories for a population of the perennial plant Trillium grandiflorum, which is declining due to high levels of herbivory by white-tailed deer. Individuals with later flowering times are less often consumed, so there is selection on this trait. We first incorporated selection analyses into a deterministic matrix model in three ways (corresponding to different methods that have been used for analyzing evolution in structured populations). Because it is not clear which of these methods works best for stage-structured models, we compared each with a more realistic, individual-based model. Deterministic models using fitness averaged over the phenotypic distribution gave trajectories that were similar to those of the individual-based model, whereas the deterministic model using fitness at the mean phenotype gave a much faster rate of evolution than that which was observed. This illustrates that subtle differences in the way in which one splices evolution into demographic models can have a large effect on expected outcomes. This study demonstrates that, by combining demographic and selection analyses, one can gauge the potential relevance of evolution to population dynamics and persistence.

Fine-scale spatial genetic structure within continuous and fragmented populations of Trillium camschatcense.

Spatial genetic structure (SGS) within populations was analyzed for the ling-lived understory perennial herb Trillium camschatcense using allozyme loci. We used Sp statistics to compare SGS between 2 life-history stages, juveniles (J) and reproductives (R), as well as between 2 populations, continuous and fragmented, with different habitat conditions. In the continuous population, significant SGS was detected in both stages but the extent was greatly reduced with the progress of the stage (J, Sp = 0.0475; R, Sp = 0.0053). We inferred that limited seed dispersal and subsequent random loss of individuals from the family patches are responsible for the J and R stage structures, respectively. The fragmented population differed in the patterns of SGS; significant structure was detected in the R stage, but not in the J stage (J, Sp = 0.0021; R, Sp = 0.0165) despite significant positive inbreeding coefficients (J, F(IS) = 0.251). The observed differences in the J-stage structures between populations may be explained by habitat fragmentation effects because reduced recruitment in the fragmented population prevents the development of maternal sibling cohort. Such comparative analysis between populations and life-history stages can be useful to understand the different underlying causes of SGS.

Pollen-mediated gene dispersal within continuous and fragmented populations of a forest understorey species, Trillium cuneatum.

Pollen movement plays a critical role in the distribution of genetic variation within and among plant populations. Direct measures of pollen movement in the large, continuous populations that characterize many herbaceous plant species are often technically difficult and biologically unreliable. Here, we studied contemporary pollen movement in four large populations of Trillium cuneatum. Three populations, located in the Georgia Piedmont, are exposed to strong anthropogenic disturbances, while the fourth population, located in the Southern Appalachian Mountains, is relatively undisturbed. Using the recently developed TwoGener analysis, we extracted estimates of the effective number of pollen donors (N(ep)), effective mating neighbourhood size (A(ep)) and the average distance of pollen movement (delta) for each population. We extended the TwoGener method by developing inference on the paternal gametic contribution to the embryo in situations where offspring genotypes are inferred from seeds and elaiosomes of species with bisporic megagametogenesis. Our estimates indicate that maternal plants do not sample pollen randomly from a global pool; rather, pollen movement in all four populations is highly restricted. Although the effective number of pollen donors per maternal plant is low (1.22-1.66) and pollen movement is highly localized in all populations, N(ep) in the disturbed Piedmont populations is higher and there is more pollen movement than in the mountains. The distance pollen moves is greater in disturbed sites and fragmented populations, possibly due to edge effects in Trillium habitats.

Genetic diversity and genetic structure of an endangered species, Trillium tschonoskii.

The genetic diversity and genetic structure of Trillium tschonoskii (Maxim) were investigated using amplified fragment length polymorphism markers. Eight primer combinations were carried out on 105 different individuals sampled from seven populations. Of the 619 discernible DNA fragments generated, 169 (27.3%) were polymorphic. The percentage of polymorphic bands within populations ranged from 4.52 to 10.50. Genetic diversity (H(E)) within populations ranged from 0.0130 to 0.0379, averaging 0.0536 at the species level. Genetic differentiation among populations was detected based on Nei's genetic diversity analysis (53.03%) and analysis of molecular variance (AMOVA) (52.43%). AMOVA indicated significant genetic differentiation among populations (52.43% of the variance) and within populations (47.57% of the variance) (p < 0.0002). Gene flow was low (0.4429) among populations. Species breeding system and limited gene flow among populations are plausible reasons for the high genetic differentiation observed for this species. We propose an appropriate strategy for conserving the genetic resources of T. tschonoskii in China.

Distribution of genetic diversity among disjunct populations of the rare forest understory herb, Trillium reliquum.

We assessed genetic diversity and its distribution in the rare southeastern US forest understory species, Trillium reliquum. In all, 21 loci were polymorphic (PS=95.5%) and the mean number of alleles per polymorphic locus was 3.05. However, genetic diversity was relatively low (Hes=0.120) considering the level of polymorphism observed for this outcrossing species. A relatively high portion of the genetic diversity (29.7%) was distributed among populations. There was no relationship between population size and genetic diversity, and we did not detect significant inbreeding. These results are best explained by the apparent self-incompatibility of this species, its longevity and clonal reproduction. To address questions regarding the history of T. reliquum's rarity, we compared results for T. reliquum with that of its more common and partially sympatric congener, T. cuneatum. Despite shared life history traits and history of land use, we observed significant genetic differences between the two species. Although T. cuneatum contains slightly lower polymorphism (Ps=85%), we detected significantly higher genetic diversity (Hes=0.217); most of its genetic diversity is contained within its populations (GST=0.092). Our results suggest that not only is there little gene flow among extant T. reliquum populations, but that rarity and population isolation in this species is of ancient origins, rather than due to more recent anthropogenic fragmentation following European colonization. The Chattahoochee River was identified as a major barrier to gene exchange.

Characterization of TrcMADS1 gene of Trillium camtschatcense (Trilliaceae) reveals functional evolution of the SOC1/TM3-like gene family.

Plant MADS-box genes encode transcriptional regulators that are critical for a number of developmental processes, such as the establishment of floral organ identity, flowering time, and fruit development. It appears that the MADS-box gene family has undergone considerable gene duplication and divergence within various angiosperm lineages. SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1)/Tomato MADS-box gene 3 (TM3)-like genes are members of the MADS-box gene family and have undergone repeated duplication events. Here, we isolated and characterized the SOC1/TM3-like gene TrcMADS1 from Trillium camtschatcense (Trilliaceae) to infer the ancestral function of SOC1/TM3-like genes. The alignment of SOC1/TM3-like genes revealed the presence of a highly conserved region in the C-terminal of predicted protein sequences, designated the SOC1 motif. Phylogenetic analysis indicated that TrcMADS1 is at the basal position of the SOC1/TM3-like gene family. The TrcMADS1 mRNA was detected in both vegetative and reproductive organs by RT-PCR. Our results suggest that duplicated copies of SOC1/TM3-like gene evolved to become variously functionally specialized.

[Chromosome maps of trilliaceae: II. A study of the genome composition in polyploid species of the genus Trillium by fluorescence nucleotide base-specific staining of heterochromatic chromosome regions]. / Karty khromosom rastenii semeistva Trilliaceae (II): izuchenie genomnogo sostava poliploidnykh vidov roda Trillium pri pomoshchi fluorestsentnogo nukleotid-spetsificheskogo okrashivaniia geterokhromatinovykh raionov khromosom.

Chromosome banding with nucleotide base-specific fluorochromes chromomycin A3 (CMA) and Hoechst 33258 (H33258) was used to study the karyotypes and to construct cytological maps for diploid Trillium camschatcense (2n = 10), tetraploid T. tschonoskii (2n = 20), hexaploid T. rhombifolium (2n = 30), and a triploid T. camschatcense x T. tschonoskii hybrid (T. x hagae, 2n = 15). With H33258, species- and genome-specific patterns with numerous AT-rich heterochromatin bands were obtained for each of the four forms; CMA revealed a few small, mostly telomeric GC-rich bands. In T. tschonoskii, the two subgenomes were similar to each other and differed from the T. camschatcense genome; on this evidence, the species was considered to be a segmental allotetraploid. In T. x hagae, one T. camschatcense and both T. tschonoskii subgenomes were identified. The subgenomes of T. rhombifolium only partly corresponded to the T. camschatcense and T. tschonoskii genomes, in contrast to the morphologically identical Japanese species T. hagae. This was assumed to indicate that allohexaploids T. rhombifolium and T. hagae originated independently at different times; i.e., their origin is polyphyletic. Based on the chromosome maps, a new nomenclature was proposed for the Trillium genomes examined: K1K1 for T. camschatcense, T1T1T2T2 for T. tschonoskii, T1T1T2T2 for T. x hagae, and K1RK1RT1RT1RT2RT2R for T. rhombifolium.