The Chinese pine genome and methylome unveil key features of conifer evolution.

Conifers dominate the world's forest ecosystems and are the most widely planted tree species. Their giant and complex genomes present great challenges for assembling a complete reference genome for evolutionary and genomic studies. We present a 25.4-Gb chromosome-level assembly of Chinese pine (Pinus tabuliformis) and revealed that its genome size is mostly attributable to huge intergenic regions and long introns with high transposable element (TE) content. Large genes with long introns exhibited higher expressions levels. Despite a lack of recent whole-genome duplication, 91.2% of genes were duplicated through dispersed duplication, and expanded gene families are mainly related to stress responses, which may underpin conifers' adaptation, particularly in cold and/or arid conditions. The reproductive regulation network is distinct compared with angiosperms. Slow removal of TEs with high-level methylation may have contributed to genomic expansion. This study provides insights into conifer evolution and resources for advancing research on conifer adaptation and development.

Three amphioxus reference genomes reveal gene and chromosome evolution of chordates.

The slow-evolving invertebrate amphioxus has an irreplaceable role in advancing our understanding of the vertebrate origin and innovations. Here we resolve the nearly complete chromosomal genomes of three amphioxus species, one of which best recapitulates the 17 chordate ancestor linkage groups. We reconstruct the fusions, retention, or rearrangements between descendants of whole-genome duplications, which gave rise to the extant microchromosomes likely existed in the vertebrate ancestor. Similar to vertebrates, the amphioxus genome gradually establishes its three-dimensional chromatin architecture at the onset of zygotic activation and forms two topologically associated domains at the Hox gene cluster. We find that all three amphioxus species have ZW sex chromosomes with little sequence differentiation, and their putative sex-determining regions are nonhomologous to each other. Our results illuminate the unappreciated interspecific diversity and developmental dynamics of amphioxus genomes and provide high-quality references for understanding the mechanisms of chordate functional genome evolution.

Benchmarking variant callers in next-generation and third-generation sequencing analysis.

DNA variants represent an important source of genetic variations among individuals. Next- generation sequencing (NGS) is the most popular technology for genome-wide variant calling. Third-generation sequencing (TGS) has also recently been used in genetic studies. Although many variant callers are available, no single caller can call both types of variants on NGS or TGS data with high sensitivity and specificity. In this study, we systematically evaluated 11 variant callers on 12 NGS and TGS datasets. For germline variant calling, we tested DNAseq and DNAscope modes from Sentieon, HaplotypeCaller mode from GATK and WGS mode from DeepVariant. All the four callers had comparable performance on NGS data and 30× coverage of WGS data was recommended. For germline variant calling on TGS data, we tested DNAseq mode from Sentieon, HaplotypeCaller mode from GATK and PACBIO mode from DeepVariant. All the three callers had similar performance in SNP calling, while DeepVariant outperformed the others in InDel calling. TGS detected more variants than NGS, particularly in complex and repetitive regions. For somatic variant calling on NGS, we tested TNscope and TNseq modes from Sentieon, MuTect2 mode from GATK, NeuSomatic, VarScan2, and Strelka2. TNscope and Mutect2 outperformed the other callers. A higher proportion of tumor sample purity (from 10 to 20%) significantly increased the recall value of calling. Finally, computational costs of the callers were compared and Sentieon required the least computational cost. These results suggest that careful selection of a tool and parameters is needed for accurate SNP or InDel calling under different scenarios.

The reference genome and transcriptome of the limestone langur, Trachypithecus leucocephalus, reveal expansion of genes related to alkali tolerance.

BACKGROUND: Trachypithecus leucocephalus, the white-headed langur, is a critically endangered primate that is endemic to the karst mountains in the southern Guangxi province of China. Studying the genomic and transcriptomic mechanisms underlying its local adaptation could help explain its persistence within a highly specialized ecological niche. RESULTS: In this study, we used PacBio sequencing and optical assembly and Hi-C analysis to create a high-quality de novo assembly of the T. leucocephalus genome. Annotation and functional enrichment revealed many genes involved in metabolism, transport, and homeostasis, and almost all of the positively selected genes were related to mineral ion binding. The transcriptomes of 12 tissues from three T. leucocephalus individuals showed that the great majority of genes involved in mineral absorption and calcium signaling were expressed, and their gene families were significantly expanded. For example, FTH1 primarily functions in iron storage and had 20 expanded copies. CONCLUSIONS: These results increase our understanding of the evolution of alkali tolerance and other traits necessary for the persistence of T. leucocephalus within an ecologically unique limestone karst environment.

Individual and combined effects of ammonia-N and sulfide on the immune function and intestinal microbiota of Pacific white shrimp Litopenaeus vannamei.

In this study, we explored the individual and combined effects of ammonia-N and sulfide stress (1 mg/L sulfide and 15 mg/L ammonia-N) on the oxidation resistance, immune response and intestinal health of Litopenaeus vannamei during 72 h exposure. The total antioxidant capacity (T-AOC), malonaldehyde (MDA) and nitric oxide (NO) content, superoxide dismutase (SOD) and catalase activity (CAT), the immune-relative gene (caspase-3, hsp70 and IMD) expression in hepatopancreas and intestine of L.vannamei and the intestinal microbiota were measured. The result showed that MDA and NO contents in hepatopancreas of L. vannamei in all treatment groups increased and remain were at high levels at the end of the stress exposure. The L. vannamei employ antioxidant defense system by increasing the activities of T-AOC, SOD and CAT enzymes in hepatopancereas and intestine to reduce oxidant damage. More severe damages with combined ammonia-N and sulfide stress to antioxidant systems were observed. The gene expression results also demonstrated that antioxidant capacity of L. vannamei was severely impaired and the apoptosis cell was initiated under the ammonia-N and sulfide stress. In addition, the environmental stress also reshaped the intestinal microbial community structure of L. vannamei that a number of original genera decreased, such as Cellvibrio, Vibrio and Rheinheimera; some new genera increased or appeared, such as Photobacterium in all treatment groups, Arcobacter and Fusibacter in sulfide stress group. Therefore, the health of L. vannamei was severely impacted when exposed to the stress of ammonia nitrogen and sulfide and these two factors can have weak synergic effects.

Transcriptome Responses to Different Environments in Intertidal Zones in the Peanut Worm Sipunculus nudus.

The peanut worm (Sipunculus nudus) is an important intertidal species worldwide. Species living in the same aquaculture area might suffer different environmental impacts. To increase knowledge of the molecular mechanisms underlying the response to environmental fluctuations, we performed a transcriptome analysis of S. nudus from different intertidal zones using a combination of the SMRT platform and the Illumina sequencing platform. (1) A total of 105,259 unigenes were assembled, and 23,063 unigenes were perfectly annotated. The results of the PacBio Iso-Seq and IIIumina RNA-Seq enriched the genetic database of S. nudus. (2) A total of 830 DEGs were detected in S. nudus from the different groups. In particular, 33 DEGs had differential expression in the top nine KEGG pathways related to pathogens, protein synthesis, and cellular immune response and signaling. The results indicate that S. nudus from different zones experience different environmental stresses. (3) Several DEGs (HSPA1, NFKBIA, eEF1A, etc.) in pathways related to pathogens (influenza A, legionellosis, measles, and toxoplasmosis) had higher expression in groups M and L. HSPA1 was clearly enriched in most of the pathways, followed by NFKBIA. The results show that the peanut worms from the M and L tidal flats might have suffered more severe environmental conditions. (4) Some DEGs (MKP, MRAS, and HSPB1) were upregulated in peanut worms from the H tidal flat, and these DEGs were mainly involved in the MAPK signaling pathway. These results indicate that the MAPK pathway may play a vital role in the immune response of the peanut worm to the effects of different intertidal flats. This study provides a valuable starting point for further studies to elucidate the molecular basis of the response to different environmental stresses in S. nudus.

The role of transposon inverted repeats in balancing drought tolerance and yield-related traits in maize.

The genomic basis underlying the selection for environmental adaptation and yield-related traits in maize remains poorly understood. Here we carried out genome-wide profiling of the small RNA (sRNA) transcriptome (sRNAome) and transcriptome landscapes of a global maize diversity panel under dry and wet conditions and uncover dozens of environment-specific regulatory hotspots. Transgenic and molecular studies of Drought-Related Environment-specific Super eQTL Hotspot on chromosome 8 (DRESH8) and ZmMYBR38, a target of DRESH8-derived small interfering RNAs, revealed a transposable element-mediated inverted repeats (TE-IR)-derived sRNA- and gene-regulatory network that balances plant drought tolerance with yield-related traits. A genome-wide scan revealed that TE-IRs associate with drought response and yield-related traits that were positively selected and expanded during maize domestication. These results indicate that TE-IR-mediated posttranscriptional regulation is a key molecular mechanism underlying the tradeoff between crop environmental adaptation and yield-related traits, providing potential genomic targets for the breeding of crops with greater stress tolerance but uncompromised yield.

A high-quality chromosome-level genome of wild Rosa rugosa.

Rosa rugosa is an important shrub with economic, ecological, and pharmaceutical value. A high-quality chromosome-scale genome for R. rugosa sequences was assembled using PacBio and Hi-C technologies. The final assembly genome sequences size was about 407.1 Mb, the contig N50 size was 2.85 Mb, and the scaffold N50 size was 56.6 Mb. More than 98% of the assembled genome sequences were anchored to seven pseudochromosomes (402.9 Mb). The genome contained 37,512 protein-coding genes, with 37,016 genes (98.68%) that were functionally annotated, and 206.67 Mb (50.76%) of the assembled sequences are repetitive sequences. Phylogenetic analyses indicated that R. rugosa diverged from Rosa chinensis ∼6.6 million years ago, and no lineage-specific whole-genome duplication event occurred after divergence from R. chinensis. Chromosome synteny analysis demonstrated highly conserved synteny between R. rugosa and R. chinensis, between R. rugosa and Prunus persica as well. Comparative genome and transcriptome analysis revealed genes related to colour, scent, and environment adaptation. The chromosome-level reference genome provides important genomic resources for molecular-assisted breeding and horticultural comparative genomics research.

Restoring mammary gland structures and functions with autogenous cell therapy.

In somatic cell reprogramming, cells must escape the somatic cell-specific gene expression program to adopt other cell fates. Here, in vitro chemical induction with RepSox generated chemically induced mammary epithelial cells (CiMECs) with milk secreting functions from goat ear fibroblasts (GEFs). Transplanted CiMECs regenerated the normal mammary gland structure with milk-secreting functions in nude mice. Single-cell RNA sequencing revealed that during the reprogramming process, GEFs may sequentially undergo embryonic ectoderm (EE)-like and different MEC developmental states and finally achieve milk secreting functions, bypassing the pluripotent state. Mechanistically, Smad3 upregulation induced by transforming growth factor ß (TGFß) receptor 1 (TGFßR1) downregulation led to GEF reprogramming into CiMECs without other reprogramming factors. The TGFßR1-Smad3 regulatory effects will provide new insight into the TGFß signaling pathway regulation of somatic cell reprogramming. These findings suggest an innovative strategy for autogenous cell therapy for mammary gland defects and the production of transgenic mammary gland bioreactors.

Single cell RNA-seq reveals molecular pathways altered by 7, 12-dimethylbenz[a]anthracene treatment on pig oocytes.

Oocytes of better quality and developmental competence are highly demanded, which is affected by many intrinsic and external factors, including environmental pollutants. We have previously demonstrated that 7, 12-dimethylbenz [a]anthracene (DMBA) reduces the developmental competence of porcine oocytes, by desynchronizing nuclear and ooplasmic maturation. However, the underlying molecular mechanism remains obscure. Here we performed single cell RNA-seq to study the transcriptome changes in DMBA-treated porcine MII oocytes, and identified 19 protein-coding genes and 156 novel long non-coding RNAs (lncRNAs) with abundance to be significantly different (P < 0.05), which enriched in signaling pathways such as glycosphingolipid biosynthesis, nicotine addiction, basal transcription factors and nucleotide excision repair. RT-qPCR on oocyte pools confirmed ornithine aminotransferase (Oat) and serine/arginine-rich splicing factor 4 (Srsf4) to be significantly up- and down-regulated, respectively (P < 0.05). Treating porcine COCs with MAPK and PLC pathway inhibitors suppressed DMBA's effects on increasing PB1 extrusion rate. In addition, DMBA co-incubation with 250 µM vitamin C derivative (l-ascorbic acid 2-phosphate sesquimagnesium salt hydrate, AA2P) and 100 µM co-enzyme Q10 (CoQ10) could significantly reduce the DMBA-induced high ROS level, and partially alleviate the DMBA-induced high PB1 rate, whereas the cleavage and blastocyst rates of parthenotes derived from treated mature oocytes remained to be low. Collectively, our findings indicate that single cell RNA-seq can help reveal the dynamics of molecular signaling pathways for porcine oocytes treated by DMBA, and supplement of anti-oxidative reagents could not sufficiently rescue DMBA-induced defects of porcine oocytes.

Bioturbation of peanut worms Sipunculus nudus on the composition of prokaryotic communities in a tidal flat as revealed by 16S rRNA gene sequences.

To understand the impacts of peanut worms Sipunculus nudus on the prokaryotic community composition in a tidal flat, an onsite investigation was conducted in Suixi in the Beibu Gulf (109.82E, 21.35N) in the burrow sediments, non-burrow sediments and the sediments without peanut worm disturbance (control). The16S rRNA gene Illumina MiSeq sequencing was used to investigate the microbial communities and their response to bioturbation by S. nudus in a sandy tidal flat. A total of 18 bacteria phyla were detected, and Proteobacteria and Cyanobacteria constituted the majority of the prokaryotic community in the samples. The distribution of the relative abundances of genera showed that approximately 6.99%-17% of the reads in the samples were classified into 25 known genera. Sulfate-reducing bacteria (Desulfococcus and Desulfosarcina) were the most abundant taxa, followed by Thermodesulfovibrionaceae LCP-6, indicating that sulfate reduction is the main process in the sandy tidal flat. The abundances of Desulfococcus, LCP-6 and Cyanobacterium in the non-burrow sediment were greater than in the burrow sediment, suggesting that the anoxic condition is more suitable for Desulfococcus and LCP-6 when the activity of S. nudus is absent. The biomass of Cyanobacterium was decreased by the feeding bioturbation of S. nudus. Meanwhile, the relative abundance of the Bacteroidetes Luteimonas in the burrow sediments was significantly greater than in the non-burrow sediment, and there was a strong relationship between S. nudus bioturbation and increased in oxygen contents and oxidation-reduction potentials in the burrow sediment. The abundances of Desulfococcus and LCP-6 were greater in the middle layer (20-30 cm) than in the top layer in the non-burrow sediment. However, the middle and bottom layers (20-30, 30-40 cm) had higher abundances of these genera than did the upper layers (0-10, 10-20 cm) in the burrow sediments. The abundances of the Fusobacteria Propionigenium and the Spirochaetes Spirochaeta were greater in the middle and bottom layers (20-30 cm, 30-40 cm) than in the top layers (0-10, 10-20 cm) in the burrow sediment, but this phenomenon was not found in the non-burrow sediment. This study demonstrates that bioturbation by S. nudus plays an important role in reshaping the bacterial community composition in intertidal regions.

Genome assembly of Nannochloropsis oceanica provides evidence of host nucleus overthrow by the symbiont nucleus during speciation.

The species of the genus Nannochloropsis are unique in their maintenance of a nucleus-plastid continuum throughout their cell cycle, non-motility and asexual reproduction. These characteristics should have been endorsed in their gene assemblages (genomes). Here we show that N. oceanica has a genome of 29.3 Mb consisting of 32 pseudochromosomes and containing 7,330 protein-coding genes; and the host nucleus may have been overthrown by an ancient red alga symbiont nucleus during speciation through secondary endosymbiosis. In addition, N. oceanica has lost its flagella and abilities to undergo meiosis and sexual reproduction, and adopted a genome reduction strategy during speciation. We propose that N. oceanica emerged through the active fusion of a host protist and a photosynthesizing ancient red alga and the symbiont nucleus became dominant over the host nucleus while the chloroplast was wrapped by two layers of endoplasmic reticulum. Our findings evidenced an alternative speciation pathway of eukaryotes.

Comparison of hypoglycemic effects of ripened pu-erh tea and raw pu-erh tea in streptozotocin-induced diabetic rats.

Pu-erh tea is produced from the leaves of large-leaf tea species (Camellia sinensis var. assamica) in the Yunnan province of China and divided into ripened pu-erh tea (RIPT, with pile-fermentation) and raw pu-erh tea (RAPT) according to processing methods. RIPT extract showed more potent anti-diabetic effects on two-hour postprandial blood glucose (2h-PBG) and fasting blood glucose (FBG) than RAPT extract. UHPLC-Q-TOF/MS and UHPLC-PDA analyses found that 17 newly formed components and the increased components after fermentation, such as quinic acid, gallic acid, caffeine, puerin I and so on, might be the main contributors to the enhanced activities of RIPT. In addition, the probiotic role of RIPT to some beneficial gut bacteria, such as lactobacillus, Prevotellaceae NK3B31 group, Alloprevotella and Prevotella, was observed in our study. These results might provide a clue to anti-diabetic mechanism and active components of pu-erh tea, and use as a functional beverage worth to be further studied.

Alterations to the Lung Microbiome in Idiopathic Pulmonary Fibrosis Patients.

Lung microbiome ecosystem homeostasis in idiopathic pulmonary fibrosis (IPF) remains uncharacterized. The aims of this study were to identify unique microbial signatures of the lung microbiome and analyze microbial gene function in IPF patients. DNA isolated from BALF samples was obtained for high-throughput gene sequencing. Microbial metagenomic data were used for principal component analysis (PCA) and analyzed at different taxonomic levels. Shotgun metagenomic data were annotated using the KEGG database and were analyzed for functional and metabolic pathways. In this study, 17 IPF patients and 38 healthy subjects (smokers and non-smokers) were recruited. For the PCA, the first and the second principal component explained 16.3 and 13.4% of the overall variability, respectively. The ß diversity of microbiome was reduced in the IPF group. Signature of IPF's microbes was enriched of Streptococcus, Pseudobutyrivibrio, and Anaerorhabdus. The translocation of lung microbiome was shown that 32.84% of them were from oral. After analysis of gene function, ABC transporter systems, biofilm formation, and two-component regulatory system were enriched in IPF patients' microbiome. Here we shown the microbiology characteristics in IPF patients. The microbiome may participate in altering internal conditions and involving in generating antibiotic resistance in IPF patients.

Genomic methylation and transcriptomic profiling provides insights into heading depression in inbred Brassica rapa L. ssp. pekinensis.

Inbreeding depression is the reduction in fitness observed in inbred populations. In plants, it leads to disease, weaker resistance to adverse environmental conditions, inhibition of growth, and decrease of yield. To elucidate molecular mechanisms behind inbreeding depression, we compared global DNA methylation and transcriptome profiles of a normal and a highly inbred heading degenerated variety of the Chinese cabbage (Brassica rapa L. ssp. pekinensis). DNA methylation was reduced in inbred plants, suggesting a change in the epigenetic landscape. Transcriptome analysis by RNA-Seq revealed that genes in auxin-response and synthesis pathways were differentially expressed in the inbreeding depression lines. Interestingly, methylation levels of some of those genes were also changed. Furthermore, endogenous IAA content was decreased in inbred plants, in agreement with expression and methylation data. Chemical inhibition of auxin also replicated the degenerated phenotype in normal plants, while exogenous IAA application had no effect in inbred depression plants, suggesting a more complex mechanism. These data indicate DNA methylation-regulated auxin pathways play a role in establishing inbred depression phenotypes in plants. Our findings reveal new insights into inbreeding depression and leafy head development in Chinese cabbage.