Network analysis of the relationship between non-suicidal self-injury, depression, and childhood trauma in adolescents.

BACKGROUND: Non-suicidal self-injury seriously harm the physical and mental health of adolescents. The aim of the current study was to explore the relationship between non-suicide self-injury, depression, and childhood trauma from the perspective of symptoms in adolescents. METHODS: A cross-sectional survey was conducted in four junior high middle schools and collected 2640 valid questionnaires. There were 1329 male students and 1311 female students. The age of the participants ranged from 11 to 17 years old, with a mean age of 13.3 (± 0.94) years. Non-suicidal self-injury (NSSI), depressive symptoms, and childhood trauma were assessed using the Adolescent Self-Harm Scale, the Childhood Depression Scale, and the Childhood Trauma Questionnaire, respectively. A network analysis was performed. RESULTS: In the network, NSSI, depressive symptoms, and childhood trauma were closely related. Negative self-esteem in the depressive symptoms and emotional abuse in childhood were the most central nodes. Negative self-esteem and negative mood were directly connected to NSSI, other nodes of depressive symptoms appeared to be indirectly connected to NSSI through these two nodes. Emotional abuse was the only node in childhood trauma categories directly connected to NSSI. Nodes of other categories of childhood trauma (physical neglect, physical abuse, emotional neglect, and sexual abuse) were indirectly connected to NSSI through emotional abuse. CONCLUSIONS: NSSI, depression, and childhood trauma of teenagers were closely related. Individuals who have suffered emotional abuse in childhood were more likely to have depressive symptoms and NSSI. Improving negative self-esteem and negative emotions and reducing emotional abuse may be beneficial in alleviating depression and reducing NSSI in adolescents.

Diversity and community assembly mechanism of soil ectomycorrhizal fungi in urban parks of Baotou City, China.

Ectomycorrhizal (EM) fungi play an important role in forest ecosystems. However, little is known about the mechanisms driving diversity and community composition of soil EM fungi in urban forest parks which are intensively affected by anthropogenic activities. In this study, we investigated the EM fungal community using Illumina high-throughput sequencing with soil samples collected from three typical forest parks, including Olympic Park, Laodong Park, and Aerding Botanical Garden of Baotou City. The results showed that soil EM fungi richness index followed a pattern of Laodong Park (146.43±25.17) > Aerding Botanical Garden (102.71±15.31) > Olympic Park (68.86±6.83). Russula, Geopora, Inocybe, Tomentella, Hebeloma, Sebacina, Amanita, Rhizopogon, Amphinema, and Lactarius were the dominant genera in the three parks. EM fungal community composition was significantly different among the three parks. Results of linear discriminant analysis effect size (LEfSe) indicated that all parks had biomarker EM fungi that exhibiting significantly different abundance. The normalized stochasticity ratio (NST) and the inferring community assembly mechanisms by phylogenetic-bin-based null model analysis (iCAMP) showed that both stochastic and deterministic processes determined soil EM fungal communities in the three urban parks, with a dominant role of the stochastic process. Drift and dispersal limitation in the stochastic process and homogeneous selection in the deterministic process were the dominant ecological processes of soil EM fungal community assembly in the three urban parks.

Comparative chloroplast genomics provides insights into the genealogical relationships of endangered Tetraena mongolica and the chloroplast genome evolution of related Zygophyllaceae species.

A comprehensive understanding of genetic background for rare species will provide an important theoretical basis for the future species management, monitoring and conservation. Tetraena mongolica is restrictedly distributed in the western Ordos plateau of China and has been listed as a national protected plant. We generated 13 chloroplast (cp) genomes of T. mongolica (size range of 106,062-106,230 bp) and conducted a series of comparative analyses of six Zygophyllaceae cp genomes. T. mongolica cp genome exhibited a quadripartite structure with drastically reduced inverted repeats (IRs, 4,315 bp) and undergone the loss of a suit of ndh genes and a copy of rRNAs. Furthermore, all the T. mongolica populations were divided into two genetic groups based on complete cp phylogenomics. In addition, notably variable genome size, gene order and structural changes had been observed among the six Zygophyllaceae cp genomes. Overall, our findings provide insights into the cp genome evolution mode and intraspecific relationships of T. mongolica, and provide a molecular basis for scientific conservation of this endangered plant.

De novo genome assembly of the medicinal plant Gentiana macrophylla provides insights into the genomic evolution and biosynthesis of iridoids.

Gentiana macrophylla is a perennial herb in the Gentianaceae family, whose dried roots are used in traditional Chinese medicine. Here, we assembled a chromosome-level genome of G. macrophylla using a combination of Nanopore, Illumina, and Hi-C scaffolding approaches. The final genome size was ~1.79 Gb (contig N50 = 720.804 kb), and 98.89% of the genome sequences were anchored on 13 pseudochromosomes (scaffold N50 = 122.73 Mb). The genome contained 55,337 protein-coding genes, and 73.47% of the assemblies were repetitive sequences. Genome evolution analysis indicated that G. macrophylla underwent two rounds of whole-genome duplication after the core eudicot γ genome triplication event. We further identified candidate genes related to the biosynthesis of iridoids, and the corresponding gene families mostly expanded in G. macrophylla. In addition, we found that root-specific genes are enriched in pathways involved in defense responses, which may greatly improve the biological adaptability of G. macrophylla. Phylogenomic analyses showed a sister relationship of asterids and rosids, and all Gentianales species formed a monophyletic group. Our study contributes to the understanding of genome evolution and active component biosynthesis in G. macrophylla and provides important genomic resource for the genetic improvement and breeding of G. macrophylla.

Characterization of the complete chloroplast genome of Hedysarum polybotrys var. alaschanicum (Fabaceae) and its phylogeny.

Hedysarum polybotrys var. alaschanicum is an important medicinal plant and is widely used in traditional Chinese medicine. The complete chloroplast genome of H. polybotrys var. alaschanicum was assembled from Illumina pair-end sequence reads. The whole chloroplast genome is 122,933 bp in length and encodes a total of 110 genes, including 76 protein-coding genes, 30 tRNA genes and 4 rRNA genes. The overall GC content of the cp genome is 35.3%. A maximum likelihood (ML) phylogenetic analysis revealed that H. polybotrys var. alaschanicum was close to Hedysarum semenovii.

Phylogenetic relationships in Chinese oaks (Fagaceae, Quercus): Evidence from plastid genome using low-coverage whole genome sequencing.

China is a second center of oak diversity but with less intensively systematic studies. Here, with 49 species representing all four sections in China, we firstly gave insight into the comprehensive phylogenetic relationships of Chinese oaks based on 54 complete plastid genomes. Our results recovered a robust phylogenetic framework and provided strong support for most nodes. The phylogenetic tree supported Quercus section Ilex as not monophyletic, in which Quercus section Cyclobalanopsis and Quercus section Cerris were nested. Most likely, incomplete lineage sorting and/or introgression among ancestral lineages in these three sections resulted in this complex pattern. The current distribution, diversification and molecular differentiation of Q. sect. Ilex in China are likely consequences of local adaptation to the geographic and paleoclimatic changes, which were driven by the uplift of Tibetan Plateau, the Hengduan Mountains and the Himalayas.

Characteristic and phylogenetic analyses of chloroplast genome for Ephedra monosperma (Ephedraceae), an important medicinal species.

Ephedra monosperma is an important medicinal plant of Ephedra (Ephedraceae). The complete chloroplast genome of E. monosperma was assembled from Illumina pair-end sequence reads. The whole chloroplast (cp) genome is 109,548 bp in length and presents a quadripartite structure consisting of two copies of inverted repeat (IR) regions (20,398) separated by a large single copy (LSC) region (60,674 bp) and a small single copy (SSC) region (8078 bp). The cp genome of E. monosperma encodes a total of 118 genes, including 73 protein-coding genes, 37 tRNA genes and 8 rRNA genes. The overall GC content of E. monosperma cp genome is 36.6%. A maximum likelihood (ML) phylogenetic analysis revealed that E. monosperma was close to Ephedra equisetina. The ML tree also showed Ephedraceae appeared more closely related to Gnetaceae than to the other families in Gymnospermae.

Plastid genomes of Elaeagnus mollis: comparative and phylogenetic analyses.

Plastomes, which are maternally inherited and show a moderate rate of evolution, play a critical role in phylogenetic reconstruction and assignment of plant species. However, little is known about the sequence divergence and molecular evolutionary patterns of plastid genomes in Elaeagnus mollis, a plant of great economic, medicinal, edible and ecological values. The plastid genome of E. mollis is 152,224-bp long and has 47 repeat sequences, including tandem (17), dispersed (12), and palindromic (18) types of repeat variations. Here, we reported six divergence hotspots (atpH-atpI, petN-psbM, trnT-psbD, trnP-psaJ, rpl32-trnL and ycf1) that could potentially be used as molecular genetic markers for population genetics and phylogenetic studies of E. mollis. A comparison of plastid genomes in the order Rosales showed that the trnH gene was duplicated only in Elaeagnaceae; therefore, it is an important marker in Elaeagnaceae. Phylogenetic analyses based on whole plastid genome sequences in 33 species revealed that Rosales is divided into two strongly supported clades and that the families Elaeagnaceae and Barbeyaceae are closely related.

Evaluating Population Genetic Structure and Demographic History of Quercus spinosa (Fagaceae) Based on Specific Length Amplified Fragment Sequencing.

Effectively identifying the genetic structure and related factors of a species can facilitate understanding the evolutionary history of the species. Phylogeographic patterns and genetic data are essential in investigating the species historical processes and diversification that response to environmental, climatic and geological influences. In this study, Specific Length Amplified Fragment Sequencing (SLAF-seq) data and ecological niche models (ENMs) are combined to identify the genetic structure and demographic modeling of Quercus spinosa, and evaluate the impacts of historical range shifts, climatic variation, and landscape factors on this species. The population topology and genetic divergence of the Cenozoic were inferred by a site frequency spectrum based composite-likelihood approach which is a novel strategy for maximizing the utility of linked SLAF markers. The overall genetic structure using model-based and model-free clustering methods was consistently identified as two geographically distinct genetic clusters. A deep divergence between two natural lineages (i.e., a western Himalaya-Hengduan Mountains lineage and an eastern Qin-ling Mountains lineage) was observed. The demographic modeling and Niche reconstruction indicated that the two groups were diverged in the late Miocene and then presented as two distinct genetic lineages. With the Quaternary glacial climate fluctuation, two groups had continuous asymmetrical secondary contact and gene exchange in the Sichuan Basin during the last glacial maximum. Besides, a significant relationship between genetic distance and geography in all individuals was identified by the Mantel test. Overall, this study 1) contributes to a better understanding of the role played by Quaternary climatic fluctuation in the present-day distributions of Q. spinosa; 2) provides a comprehensive view of the genome-wide variation of sclerophyllous forests in ecological adaptive evolution; 3) indicates that dispersal limitation and ecological divergence contribute to the genome-wide differentiation of Q. spinosa, which supports a hypothesis that complex geography and climatic changes strongly influence the evolutionary origin and history of the species.

Comparative Plastid Genomes of Primula Species: Sequence Divergence and Phylogenetic Relationships.

Compared to traditional DNA markers, genome-scale datasets can provide mass information to effectively address historically difficult phylogenies. Primula is the largest genus in the family Primulaceae, with members distributed mainly throughout temperate and arctic areas of the Northern Hemisphere. The phylogenetic relationships among Primula taxa still maintain unresolved, mainly due to intra- and interspecific morphological variation, which was caused by frequent hybridization and introgression. In this study, we sequenced and assembled four complete plastid genomes (Primula handeliana, Primula woodwardii, Primula knuthiana, and Androsace laxa) by Illumina paired-end sequencing. A total of 10 Primula species (including 7 published plastid genomes) were analyzed to investigate the plastid genome sequence divergence and their inferences for the phylogeny of Primula. The 10 Primula plastid genomes were similar in terms of their gene content and order, GC content, and codon usage, but slightly different in the number of the repeat. Moderate sequence divergence was observed among Primula plastid genomes. Phylogenetic analysis strongly supported that Primula was monophyletic and more closely related to Androsace in the Primulaceae family. The phylogenetic relationships among the 10 Primula species showed that the placement of P. knuthiana-P. veris clade was uncertain in the phylogenetic tree. This study indicated that plastid genome data were highly effective to investigate the phylogeny.

Plastid Genome Comparative and Phylogenetic Analyses of the Key Genera in Fagaceae: Highlighting the Effect of Codon Composition Bias in Phylogenetic Inference.

Fagaceae is one of the largest and economically important taxa within Fagales. Considering the incongruence among inferences from plastid and nuclear genes in the previous Fagaceae phylogeny studies, we assess the performance of plastid phylogenomics in this complex family. We sequenced and assembled four complete plastid genomes (Fagus engleriana, Quercus spinosa, Quercus aquifolioides, and Quercus glauca) using reference-guided assembly approach. All of the other 12 published plastid genomes in Fagaceae were retrieved for genomic analyses (including repeats, sequence divergence and codon usage) and phylogenetic inference. The genomic analyses reveal that plastid genomes in Fagaceae are conserved. Comparing the phylogenetic relationships of the key genera in Fagaceae inferred from different codon positions and gene function datasets, we found that the first two codon sites dataset recovered nearly all relationships and received high support. Thus, the result suggested that codon composition bias had great influence on Fagaceae phylogenetic inference. Our study not only provides basic understanding of Fagaceae plastid genomes, but also illuminates the effectiveness of plastid phylogenomics in resolving relationships of this intractable family.

The Complete Chloroplast Genome Sequences of Six Rehmannia Species.

Rehmannia is a non-parasitic genus in Orobanchaceae including six species mainly distributed in central and north China. Its phylogenetic position and infrageneric relationships remain uncertain due to potential hybridization and polyploidization. In this study, we sequenced and compared the complete chloroplast genomes of six Rehmannia species using Illumina sequencing technology to elucidate the interspecific variations. Rehmannia plastomes exhibited typical quadripartite and circular structures with good synteny of gene order. The complete genomes ranged from 153,622 bp to 154,055 bp in length, including 133 genes encoding 88 proteins, 37 tRNAs, and 8 rRNAs. Three genes (rpoA, rpoC2, accD) have potentially experienced positive selection. Plastome size variation of Rehmannia was mainly ascribed to the expansion and contraction of the border regions between the inverted repeat (IR) region and the single-copy (SC) regions. Despite of the conserved structure in Rehmannia plastomes, sequence variations provide useful phylogenetic information. Phylogenetic trees of 23 Lamiales species reconstructed with the complete plastomes suggested that Rehmannia was monophyletic and sister to the clade of Lindenbergia and the parasitic taxa in Orobanchaceae. The interspecific relationships within Rehmannia were completely different with the previous studies. In future, population phylogenomic works based on plastomes are urgently needed to clarify the evolutionary history of Rehmannia.

The complete chloroplast genome of Dipteronia sinensis (Aceraceae), an endangered endemic species to China.

The complete chloroplast sequence of D. sinensis was reported in this study. The total length was 157 080 bp containing a pair of 26 766 bp inverted repeat regions (IRa and IRb), which were separated by a small single copy regions and a large single copy regions (SSC and LSC) of 18 093 and 85 455 bp, respectively. A total of 138 functional genes were annotated, which included 90 protein-coding genes, 40 tRNAs, and eight rRNA genes. The overall GC content of the complete chloroplast genome is 37.8% and in LSC, SSC and IR regions were 35.9%, 32.1%, and 42.7%, respectively. The maximum likelihood (ML) phylogenetic analysis revealed that D. sinensis was closely related to A. buergerianum subsp. Ningpoense in Sapindales order.

The complete chloroplast genome of Quercus baronii (Quercus L.).

Quercus baronii is an endemic species to China. The total length of the Q. baronii chloroplast genome is 161 072 bp. It exhibits typical quadripartite structure, composed of two copies of inverted repeat (IR) regions (25 843 bp), large single copy (LSC) region (90 341 bp), and small single copy (SSC) region (19 045 bp). The genome contains 134 genes, including 86 protein-coding genes, 40 tRNA genes, and 8 rRNA genes. Phylogenetic analysis revealed its relationship with reported Quercus L. species.

Comparative Analysis of the Complete Chloroplast Genomes of Five Quercus Species.

Quercus is considered economically and ecologically one of the most important genera in the Northern Hemisphere. Oaks are taxonomically perplexing because of shared interspecific morphological traits and intraspecific morphological variation, which are mainly attributed to hybridization. Universal plastid markers cannot provide a sufficient number of variable sites to explore the phylogeny of this genus, and chloroplast genome-scale data have proven to be useful in resolving intractable phylogenetic relationships. In this study, the complete chloroplast genomes of four Quercus species were sequenced, and one published chloroplast genome of Quercus baronii was retrieved for comparative analyses. The five chloroplast genomes ranged from 161,072 bp (Q. baronii) to 161,237 bp (Q. dolicholepis) in length, and their gene organization and order, and GC content, were similar to those of other Fagaceae species. We analyzed nucleotide substitutions, indels, and repeats in the chloroplast genomes, and found 19 relatively highly variable regions that will potentially provide plastid markers for further taxonomic and phylogenetic studies within Quercus. We observed that four genes (ndhA, ndhK, petA, and ycf1) were subject to positive selection. The phylogenetic relationships of the Quercus species inferred from the chloroplast genomes obtained moderate-to-high support, indicating that chloroplast genome data may be useful in resolving relationships in this genus.