Population genomics and epigenomics of Spirodela polyrhiza provide insights into the evolution of facultative asexuality.

Many plants are facultatively asexual, balancing short-term benefits with long-term costs of asexuality. During range expansion, natural selection likely influences the genetic controls of asexuality in these organisms. However, evidence of natural selection driving asexuality is limited, and the evolutionary consequences of asexuality on the genomic and epigenomic diversity remain controversial. We analyzed population genomes and epigenomes of Spirodela polyrhiza, (L.) Schleid., a facultatively asexual plant that flowers rarely, revealing remarkably low genomic diversity and DNA methylation levels. Within species, demographic history and the frequency of asexual reproduction jointly determined intra-specific variations of genomic diversity and DNA methylation levels. Genome-wide scans revealed that genes associated with stress adaptations, flowering and embryogenesis were under positive selection. These data are consistent with the hypothesize that natural selection can shape the evolution of asexuality during habitat expansions, which alters genomic and epigenomic diversity levels.

Will advancement in technologies bring fear and damage human employment? Evidence from China's manufacturing industry.

Advancement in technologies such as robotic industries and artificial intelligence bring fear among human being that jobs will be substituted by robots. Base on the panel data of 28 China's manufacturing industries, this research analyzed the impact of technical progress bias on employment. First, we calculate the technical progress bias index of 28 industries base on the stochastic frontier model with transcendental logarithm function found 16 industries were toward the skilled labor while the remaining 12 industries were toward the unskilled labor. Second, the empirical results show that technical progress bias has a positive impact on the total manufacturing employment and significant positive effect on the unskilled labor, while no significant impact on skilled labor employment. Third, the threshold effect test proves that if taking industry value-added per capita or R&D capital stock as threshold variable, the threshold about the impact exist, making the impact on skilled labor was insignificant.

Genetic Mechanism of Non-Targeted-Site Resistance to Diquat in Spirodela polyrhiza.

Understanding non-target-site resistance (NTSR) to herbicides represents a pressing challenge as NTSR is widespread in many weeds. Using giant duckweed (Spirodela polyrhiza) as a model, we systematically investigated genetic and molecular mechanisms of diquat resistance, which can only be achieved via NTSR. Quantifying the diquat resistance of 138 genotypes, we revealed an 8.5-fold difference in resistance levels between the most resistant and most susceptible genotypes. Further experiments suggested that diquat uptake and antioxidant-related processes jointly contributed to diquat resistance in S. polyrhiza. Using a genome-wide association approach, we identified several candidate genes, including a homolog of dienelactone hydrolase, that are associated with diquat resistance in S. polyrhiza. Together, these results provide new insights into the mechanisms and evolution of NTSR in plants.

A high-quality genome assembly of the waterlily aphid Rhopalosiphum nymphaeae.

Waterlily aphid, Rhopalosiphum nymphaeae (Linnaeus), is a host-alternating aphid known to feed on both terrestrial and aquatic hosts. It causes damage through direct herbivory and acting as a vector for plant viruses, impacting worldwide Prunus spp. fruits and aquatic plants. Interestingly, R. nymphaeae's ability to thrive in both aquatic and terrestrial conditions sets it apart from other aphids, offering a unique perspective on adaptation. We present the first high-quality R. nymphaeae genome assembly with a size of 324.4 Mb using PacBio long-read sequencing. The resulting assembly is highly contiguous with a contig N50 reached 12.7 Mb. The BUSCO evaluation suggested a 97.5% completeness. The R. nymphaeae genome consists of 16.9% repetitive elements and 16,834 predicted protein-coding genes. Phylogenetic analysis positioned R. nymphaeae within the Aphidini tribe, showing close relations to R. maidis and R. padi. The high-quality reference genome R. nymphaeae provides a unique resource for understanding genome evolution in aphids and paves the foundation for understanding host plant adaptation mechanisms and developing pest control strategies.

MAN2A1 predicts prognosis and progression through cancer-related pathways in colorectal cancer.

Background: Colorectal cancer (CRC) is a common malignant tumor leading to poor prognosis and high mortality. Mannosidase alpha class 2A member 1 (MAN2A1) turns to oncogene through fusing Fer tyrosine kinase (FER) and associates with multiple cancer occurrence. In order to determine whether MAN2A1 can promote tumorigenesis and metastasis in CRC, we conducted a series of studies. Methods: We obtained gene expression and clinical data of CRC from The Cancer Genome Atlas (TCGA) databases. RNA raw counts data was merged by Python. Batch processing of univariate Cox regression analysis was performed to preliminary identify the genes associated with prognosis. Differentially expressed genes (DEGs) between lymph node metastasis (LNM) patients and non-LNM patients were identified via edgR in Sangerbox tool. Protein-protein interactive (PPI) network was constructed using Search Tool for the Retrieval of Interacting Genes (STRING) and Cytoscape software. Kaplan-Meier (KM) survival of CRC patients was analyzed by Sangerbox tool. Clinicopathologic characteristics of CRC patients were analyzed by SPSS statistics software. Differences in RNA expression levels of genes were validated in our cohort by real-time polymerase chain reaction (RT-PCR). Analyses of Signaling pathways and gene ontology were explored by gene set enrichment analysis (GSEA). Results: We first obtained 4,455 genes associated with the prognosis of CRC, and 998 of these genes were also DEGs in CRC between metastatic CRC tissues and in situ tissues. Therein, MAN2A1 expression was downregulated in LNM CRC compared with CRC in situ, also downregulated in CRC compared with adjacent normal tissues, and high gene expression levels of MAN2A1 was associated with better survival. Conclusions: Our study suggested that MAN2A1 could be a potential biomarker significantly related to prognosis and LNM of CRC patients.

CCCTC-binding factor is an upstream regulator of the pluripotency factor Oct4 and functions in active transcription of linc1253 and linc1356 genes in pluripotent cells.

The embryonic stem cell- (ESC) specific transcription factor Oct4 is a well-known master regulator of pluripotency. The aim of this study was to identify upstream regulators of Oct4 and explore their functional link in regulating lincRNA expression in ESCs. By quantitative real-time PCR (RT-qPCR) analysis upon CCCTC-binding factor (CTCF) or Oct4 knockdown, here, we found that the chromatin insulator CTCF transcriptionally controls Oct4 gene expression in mouse ES cells. Furthermore, co-immunoprecipitation assays showed that CTCF physically interacts with Oct4. By analyzing CTCF and Oct4 ChIP-seq datasets in mouse ES cells and investigating their genomic occupancies, we demonstrated that CTCF and Oct4 share overlapping regulatory functions and are required for active transcription of long intergenic non-coding RNAs (lincRNAs) linc1253 and linc1356, which were reported to repress cellular lineage programs and maintain a pluripotent state. In summary, we propose an integrated model of transcriptional control mediated by CTCF, the master weaver of the genome, for the upstream regulation of Oct4-and ESC-associated genes. These results connect the chromatin insulator CTCF and the pluripotency factor Oct4 in the regulation of lincRNAs in pluripotent cells.

Transcriptomic Profiling Reveals Shared Signalling Networks Between Flower Development and Herbivory-Induced Responses in Tomato.

Most flowering plants must defend themselves against herbivores for survival and attract pollinators for reproduction. Although traits involved in plant defence and pollinator attraction are often localised in leaves and flowers, respectively, they will show a diffuse evolution if they share the same molecular machinery and regulatory networks. We performed RNA-sequencing to characterise and compare transcriptomic changes involved in herbivory-induced defences and flower development, in tomato leaves and flowers, respectively. We found that both the herbivory-induced responses and flower development involved alterations in jasmonic acid signalling, suppression of primary metabolism and reprogramming of secondary metabolism. We identified 411 genes that were involved in both processes, a number significantly higher than expected by chance. Genetic manipulation of key regulators of induced defences also led to the expression changes in the same genes in both leaves and flowers. Targeted metabolomic analysis showed that among closely related tomato species, jasmonic acid and α-tomatine are correlated in flower buds and herbivory-induced leaves. These findings suggest that herbivory-induced responses and flower development share a common molecular machinery and likely have coevolved in nature.

Genome Assembly and Analyses of the Macrofungus Macrocybe gigantea.

Macrocybe gigantea (M. gigantea) is a macrofungus genus that contains a big number of fairly fleshy gilled mushrooms with white spores. This macrofungus produces diverse bioactive compounds, antioxidants, and water-soluble polysaccharides. However, the genomic resources of this species remain unknown. Here, we assembled the genome of M. gigantea (41.23 Mb) into 336 scaffolds with a N50 size of 374,455 bp and compared it with the genomes of eleven other macrofungi. Comparative genomics study confirmed that M. gigantea belonged to the Macrocybe genus, a stand-alone genus different from the Tricholoma genus. In addition, we found that glycosyl hydrolase family 28 (GH28) in M. gigantea shared conserved motifs that were significantly different from their counterparts in Tricholoma. The genomic resource uncovered by this study will enhance our understanding of fungi biology, especially the differences in their growth rates and energy metabolism.

Improved de novo Assembly of the Achlorophyllous Orchid Gastrodia elata.

Achlorophyllous plants are full mycoheterotrophic plants with no chlorophyll and they obtain their nutrients from soil fungi. Gastrodia elata is a perennial, achlorophyllous orchid that displays distinctive evolutionary strategy of adaptation to the non-photosynthetic lifestyle. Here in this study, the genome of G. elata was assembled to 1.12 Gb with a contig N50 size of 110 kb and a scaffold N50 size of 1.64 Mb so that it helped unveil the genetic basics of those adaptive changes. Based on the genomic data, key genes related to photosynthesis, leaf development, and plastid division pathways were found to be lost or under relaxed selection during the course of evolution. Thus, the genome sequence of G. elata provides a good resource for future investigations of the evolution of orchids and other achlorophyllous plants.

Evolution of a Novel and Adaptive Floral Scent in Wild Tobacco.

Many plants emit diverse floral scents that mediate plant-environment interactions and attain reproductive success. However, how plants evolve novel and adaptive biosynthetic pathways for floral volatiles remains unclear. Here, we show that in the wild tobacco, Nicotiana attenuata, a dominant species-specific floral volatile (benzyl acetone, BA) that attracts pollinators and deters florivore is synthesized by phenylalanine ammonia-lyase 4 (NaPAL4), isoflavone reductase 3 (NaIFR3), and chalcone synthase 3 (NaCHAL3). Transient expression of NaFIR3 alone in N. attenuata leaves is sufficient and necessary for ectopic foliar BA emissions, and coexpressing NaIFR3 with NaPAL4 and NaCHAL3 increased the BA emission levels. Independent changes in transcription of NaPAL4 and NaCHAL3 contributed to intraspecific variations of floral BA emission. However, among species, the gain of expression of NaIFR3 resulted in the biosynthesis of BA, which was only found in N. attenuata. This study suggests that novel metabolic pathways associated with adaptation can arise via reconfigurations of gene expression.

The genome assembly of asparagus bean, Vigna unguiculata ssp. sesquipedialis.

Asparagus bean (Vigna. unguiculata ssp. sesquipedialis), known for its very long and tender green pods, is an important vegetable crop broadly grown in the developing Asian countries. In this study, we reported a 632.8 Mb assembly (549.81 Mb non-N size) of asparagus bean based on the whole genome shotgun sequencing strategy. We also generated a linkage map for asparagus bean, which helped anchor 94.42% of the scaffolds into 11 pseudo-chromosomes. A total of 42,609 protein-coding genes and 3,579 non-protein-coding genes were predicted from the assembly. Taken together, these genomic resources of asparagus bean will help develop a pan-genome of V. unguiculata and facilitate the investigation of economically valuable traits in this species, so that the cultivation of this plant would help combat the protein and energy malnutrition in the developing world.

CCCTC-Binding Factor Transcriptionally Targets Wdr5 to Mediate Somatic Cell Reprogramming.

It was previously reported that WD repeat domain 5 (Wdr5), a core member of the mammalian Trithorax complex, is a key regulator for the maintenance of embryonic stem (ES) cell pluripotency and somatic cell reprogramming. However, it remains unclear whether other factors are also involved in this process. Here, we show that CTCF is an upstream regulator of Wdr5: It physically associates with Wdr5 and further transcriptionally controls its expression by directly targeting the Wdr5 gene promoter. As a downstream effector, overexpression of Wdr5 can rescue ES cells from growth defects and decreased formation of induced pluripotent stem cells caused by CTCF knockdown. Our results reveal the functional relevance of CTCF and Wdr5, and they connect them to the re-establishment of pluripotency.

Genetic Variation of Goat Interferon Regulatory Factor 3 Gene and Its Implication in Goat Evolution.

The immune systems are fundamentally vital for evolution and survival of species; as such, selection patterns in innate immune loci are of special interest in molecular evolutionary research. The interferon regulatory factor (IRF) gene family control many different aspects of the innate and adaptive immune responses in vertebrates. Among these, IRF3 is known to take active part in very many biological processes. We assembled and evaluated 1356 base pairs of the IRF3 gene coding region in domesticated goats from Africa (Nigeria, Ethiopia and South Africa) and Asia (Iran and China) and the wild goat (Capra aegagrus). Five segregating sites with θ value of 0.0009 for this gene demonstrated a low diversity across the goats' populations. Fu and Li tests were significantly positive but Tajima's D test was significantly negative, suggesting its deviation from neutrality. Neighbor joining tree of IRF3 gene in domesticated goats, wild goat and sheep showed that all domesticated goats have a closer relationship than with the wild goat and sheep. Maximum likelihood tree of the gene showed that different domesticated goats share a common ancestor and suggest single origin. Four unique haplotypes were observed across all the sequences, of which, one was particularly common to African goats (MOCH-K14-0425, Poitou and WAD). In assessing the evolution mode of the gene, we found that the codon model dN/dS ratio for all goats was greater than one. Phylogenetic Analysis by Maximum Likelihood (PAML) gave a ω0 (dN/dS) value of 0.067 with LnL value of -6900.3 for the first Model (M1) while ω2 = 1.667 in model M2 with LnL value of -6900.3 with positive selection inferred in 3 codon sites. Mechanistic empirical combination (MEC) model for evaluating adaptive selection pressure on particular codons also confirmed adaptive selection pressure in three codons (207, 358 and 408) in IRF3 gene. Positive diversifying selection inferred with recent evolutionary changes in domesticated goat IRF3 led us to conclude that the gene evolution may have been influenced by domestication processes in goats.

Genome of Plant Maca (Lepidium meyenii) Illuminates Genomic Basis for High-Altitude Adaptation in the Central Andes.

Maca (Lepidium meyenii Walp, 2n = 8x = 64), belonging to the Brassicaceae family, is an economic plant cultivated in the central Andes sierra in Peru (4000-4500 m). Considering that the rapid uplift of the central Andes occurred 5-10 million years ago (Ma), an evolutionary question arises regarding how plants such as maca acquire high-altitude adaptation within a short geological period. Here, we report the high-quality genome assembly of maca, in which two closely spaced maca-specific whole-genome duplications (WGDs; ∼6.7 Ma) were identified. Comparative genomic analysis between maca and closely related Brassicaceae species revealed expansions of maca genes and gene families involved in abiotic stress response, hormone signaling pathway, and secondary metabolite biosynthesis via WGDs. The retention and subsequent functional divergence of many duplicated genes may account for the morphological and physiological changes (i.e., small leaf shape and self-fertility) in maca in a high-altitude environment. In addition, some duplicated maca genes were identified with functions in morphological adaptation (i.e., LEAF CURLING RESPONSIVENESS) and abiotic stress response (i.e., GLYCINE-RICH RNA-BINDING PROTEINS and DNA-DAMAGE-REPAIR/TOLERATION 2) under positive selection. Collectively, the maca genome provides useful information to understand the important roles of WGDs in the high-altitude adaptation of plants in the Andes.

Transcriptomic profile of tobacco in response to Alternaria longipes and Alternaria alternata infections.

Tobacco brown spot caused by Alternaria fungal species is one of the most damaging diseases, and results in significant yield losses. However, little is known about the systematic response of tobacco to this fungal infection. To fill this knowledge gap, de novo assemblies of tobacco leaf transcriptomes were obtained in cultivars V2 and NC89 after the inoculation of either Alternaria longipes (AL) or Alternaria alternata (AA) at three different time points. We studied the gene expression profile of each cultivar-pathogen combination, and identified eight differentially expressed genes shared among all combinations. Gene ontology enrichment analysis of the differentially expressed genes revealed key components during the fungal infection, which included regulation of gene expression (GO:0010468), regulation of RNA metabolic process (GO:0051252), tetrapyrrole binding (GO:0046906), and external encapsulating structure (GO:0030312). Further analyses of the continuously upregulated/downregulated genes and the resistance genes demonstrated that the gene expression profile upon fungal infection was contingent on the specific cultivar and pathogen. In conclusion, this study provides a solid foundation for the investigation of plant-pathogen interaction, and is of great importance for disease prevention and molecular breeding.

Single Base-Resolution Methylome of the Dizygotic Sheep.

Sheep is an important livestock in the world for meat, dairy and wool production. The third version of sheep reference genome has been recently assembled, but sheep DNA methylome has not been profiled yet. In this study, we report the comprehensive sheep methylome with 94.38% cytosine coverage at single base resolution by sequencing DNA samples from Longissimus dorsi of dizygotic Sunit sheep, which were bred in different habitats. We also compared methylomes between the twin sheep. DNA methylation status at genome-scale differentially methylated regions (DMRs), functional genomic regions and 248 DMR-containing genes were identified between the twin sheep. Gene ontology (GO) and KEGG annotations of these genes were performed to discover computationally predicted function. Lipid metabolism, sexual maturity and tumor-associated categories were observed to significantly enrich DMR-containing genes. These findings could be used to illustrate the relationship between phenotypic variations and gene methylation patterns.