The telomere-to-telomere (T2T) genome of Peucedanum praeruptorum Dunn provides insights into the genome evolution and coumarin biosynthesis.

BACKGROUND: Traditional Chinese medicine has used Peucedanum praeruptorum Dunn (Apiaceae) for a long time. Various coumarins, including the significant constituents praeruptorin (A-E), are the active constituents in the dried roots of P. praeruptorum. Previous transcriptomic and metabolomic studies have attempted to elucidate the distribution and biosynthetic network of these medicinal-valuable compounds. However, the lack of a high-quality reference genome impedes an in-depth understanding of genetic traits and thus the development of better breeding strategies. RESULTS: A telomere-to-telomere (T2T) genome was assembled for P. praeruptorum by combining PacBio HiFi, ONT ultra-long, and Hi-C data. The final genome assembly was approximately 1.798 Gb, assigned to 11 chromosomes with genome completeness >98%. Comparative genomic analysis suggested that P. praeruptorum experienced 2 whole-genome duplication events. By the transcriptomic and metabolomic analysis of the coumarin metabolic pathway, we presented coumarins' spatial and temporal distribution and the expression patterns of critical genes for its biosynthesis. Notably, the COSY and cytochrome P450 genes showed tandem duplications on several chromosomes, which may be responsible for the high accumulation of coumarins. CONCLUSIONS: A T2T genome for P. praeruptorum was obtained, providing molecular insights into the chromosomal distribution of the coumarin biosynthetic genes. This high-quality genome is an essential resource for designing engineering strategies for improving the production of these valuable compounds.

The rise of baobab trees in Madagascar.

The baobab trees (genus Adansonia) have attracted tremendous attention because of their striking shape and distinctive relationships with fauna1. These spectacular trees have also influenced human culture, inspiring innumerable arts, folklore and traditions. Here we sequenced genomes of all eight extant baobab species and argue that Madagascar should be considered the centre of origin for the extant lineages, a key issue in their evolutionary history2,3. Integrated genomic and ecological analyses revealed the reticulate evolution of baobabs, which eventually led to the species diversity seen today. Past population dynamics of Malagasy baobabs may have been influenced by both interspecific competition and the geological history of the island, especially changes in local sea levels. We propose that further attention should be paid to the conservation status of Malagasy baobabs, especially of Adansonia suarezensis and Adansonia grandidieri, and that intensive monitoring of populations of Adansonia za is required, given its propensity for negatively impacting the critically endangered Adansonia perrieri.

Integrating genomic and multiomic data for Angelica sinensis provides insights into the evolution and biosynthesis of pharmaceutically bioactive compounds.

Angelica sinensis roots (Angelica roots) are rich in many bioactive compounds, including phthalides, coumarins, lignans, and terpenoids. However, the molecular bases for their biosynthesis are still poorly understood. Here, an improved chromosome-scale genome for A. sinensis var. Qinggui1 is reported, with a size of 2.16 Gb, contig N50 of 4.96 Mb and scaffold N50 of 198.27 Mb, covering 99.8% of the estimated genome. Additionally, by integrating genome sequencing, metabolomic profiling, and transcriptome analysis of normally growing and early-flowering Angelica roots that exhibit dramatically different metabolite profiles, the pathways and critical metabolic genes for the biosynthesis of these major bioactive components in Angelica roots have been deciphered. Multiomic analyses have also revealed the evolution and regulation of key metabolic genes for the biosynthesis of pharmaceutically bioactive components; in particular, TPSs for terpenoid volatiles, ACCs for malonyl CoA, PKSs for phthalide, and PTs for coumarin biosynthesis were expanded in the A. sinensis genome. These findings provide new insights into the biosynthesis of pharmaceutically important compounds in Angelica roots for exploration of synthetic biology and genetic improvement of herbal quality.

Multiple reinforcement measures of flysch landslide.

This work is mainly intended to investigate the flysch landslide reinforcement measures used in the Smokovac-Matesevo section of the North-South Expressway project in Montenegro. Bentley's Plaxis software is used for a numerical analysis of sliding surface parameters of flysch strata in the limit equilibrium state. This study analyzes the slope safety factor for rreinforcement measures such as rock bolts, retaining walls, anti-sliding piles, slope unloading and bolt anchoring and obtains an optimal combination of reinforcement application for the flysch landslide. The effects of seismic action on complex stress and the discontinuous stress boundary conditions arising from various reinforcement measures on landslide stability are also examined. The measures applied in this paper can be used as a reference for flysch landslide reinforcement or other similar slope engineering measures.

Research on construction deformation prediction and disaster warning of karst slope based on mutation theory.

In the study of deformation prediction and disaster warning during karst slope construction, the influencing factors and deformation law should be comprehensively considered. The layout of the deformation monitoring points of karst slope is affected by the thickness of karst overburden soil, dissolution and fragmentation degree, karst development degree, slope cracking degree, fault or weak interlayer and other factors. In this paper, the author aimed at the problem of construction deformation prediction and disaster warning of karst slope, proposed an improved model of cusp mutation by applying and optimizing the cusp mutation model, analysed the deformation trend and sudden change type of the slope, and obtained the critical control early warning value of slope deformation. Therefore, it is feasible to analyse the deformation and mutation characteristics of karstified slope by using a virtual reality-mutation model. In addition, based on the empirical formula of the slope sliding limit deformation rate and grey prediction model, the critical control warning value of slope deformation is obtained, which provides a basis to quantify the deformation index of risk evaluation. This method provides a new idea to predict karst slope construction deformation and catastrophic deformation warning and has a reference value for similar engineering examples.

Comparative genome anatomy reveals evolutionary insights into a unique amphitriploid fish.

Triploids are rare in nature because of difficulties in meiotic and gametogenic processes, especially in vertebrates. The Carassius complex of cyprinid teleosts contains sexual tetraploid crucian carp/goldfish (C. auratus) and unisexual hexaploid gibel carp/Prussian carp (C. gibelio) lineages, providing a valuable model for studying the evolution and maintenance mechanism of unisexual polyploids in vertebrates. Here we sequence the genomes of the two species and assemble their haplotypes, which contain two subgenomes (A and B), to the chromosome level. Sequencing coverage analysis reveals that C. gibelio is an amphitriploid (AAABBB) with two triploid sets of chromosomes; each set is derived from a different ancestor. Resequencing data from different strains of C. gibelio show that unisexual reproduction has been maintained for over 0.82 million years. Comparative genomics show intensive expansion and alterations of meiotic cell cycle-related genes and an oocyte-specific histone variant. Cytological assays indicate that C. gibelio produces unreduced oocytes by an alternative ameiotic pathway; however, sporadic homologous recombination and a high rate of gene conversion also exist in C. gibelio. These genomic changes might have facilitated purging deleterious mutations and maintaining genome stability in this unisexual amphitriploid fish. Overall, the current results provide novel insights into the evolutionary mechanisms of the reproductive success in unisexual polyploid vertebrates.

Deep sequencing of 1320 genes reveals the landscape of protein-truncating variants and their contribution to psoriasis in 19,973 Chinese individuals.

Protein-truncating variants (PTVs) have important impacts on phenotype diversity and disease. However, their population genetics characteristics in more globally diverse populations are not well defined. Here, we describe patterns of PTVs in 1320 genes sequenced in 10,539 healthy controls and 9434 patients with psoriasis, all of Han Chinese ancestry. We identify 8720 PTVs, of which 77% are novel, and estimate 88% of all PTVs are deleterious and subject to purifying selection. Furthermore, we show that individuals with psoriasis have a significantly higher burden of PTVs compared to controls (P = 0.02). Finally, we identified 18 PTVs in 14 genes with unusually high levels of population differentiation, consistent with the action of local adaptation. Our study provides insights into patterns and consequences of PTVs.

Evaluation of exosomal miRNAs as potential diagnostic biomarkers for acute myocardial infarction using next-generation sequencing.

BACKGROUND: Acute myocardial infarction (AMI) is one of the most common global causes of death. Although considerable progress has been made in AMI diagnosis, there remains an urgent need for novel diagnostic biomarkers for its prevention and treatment. Functional exosomal microRNAs (miRNAs) are recognized as potential biomarkers in many diseases. This study's objective was to identify specific plasma exosomal miRNAs with biomarker potential for early AMI detection. METHODS: Exosomes from the plasma of 26 coronary artery disease (CAD) patients, 55 AMI patients, and 37 healthy controls were isolated and characterized by transmission electron microcopy (TEM), western blotting, and nanoparticle tracking analysis (NTA). The miRNAs were purified from exosomes, and unique molecular identifier (UMI) small RNA sequencing was performed. The random forest (RF) model was trained to predict potential biomarkers. RESULTS: NTA demonstrated that nanoparticle concentration did not change after AMI, while nanoparticle size distribution significantly decreased. The CAD and AMI groups' miRNA expression profiles significantly differed from the healthy group's profile. The RF classifier could be used to distinguish the healthy group from the AMI group, but could not be used to distinguish the CAD group from the other groups, which caused a high classification error rate. Eighteen miRNAs were selected as biomarkers based on their RF classifier significance. The diagnostic accuracy of 18 miRNAs was evaluated using AUC values of 0.93, 0.87, and 0.75 to detect healthy controls, AMI, and CAD, respectively. CONCLUSIONS: Nanoparticle diameter and the 18 miRNAs may serve as simple and accessible fingerprints for early AMI diagnosis.

Whole genome sequencing of silver carp (Hypophthalmichthys molitrix) and bighead carp (Hypophthalmichthys nobilis) provide novel insights into their evolution and speciation.

The edible silver carp (Hypophthalmichthys molitrix) and bighead carp (H. nobilis), which are two of the "Four Domesticated Fish" of China, are cultivated intensively worldwide. Here, we constructed 837- and 845-Mb draft genome assemblies for the silver carp and the bighead carp, respectively, including 24,571 and 24,229 annotated protein-coding genes. Genetic maps, anchoring 71.7% and 83.8% of all scaffolds, were obtained for the silver and bighead carp, respectively. Phylogenetic analysis showed that the bighead carp formed a clade with the silver carp, with an estimated divergence time of 3.6 million years ago; the time of divergence between the silver carp and zebrafish was 50.7 million years ago. An East Asian cyprinid genome-specific chromosome fusion took place ~9.2 million years after this clade diverged from the clade containing the common carp and Sinocyclocheilus. KEGG and GO analyses indicated that the expanded gene families in the silver and bighead carp were associated with diseases, the immune system and environmental adaptations. Genomic regions differentiating the silver and bighead carp populations were detected based on the whole-genome sequences of 42 individuals. Genes associated with the divergent regions were associated with reproductive system development and the development of primary female sexual characteristics. Thus, our results provided a novel systematic genomic analysis of the East Asian cyprinids, as well as the evolution and speciation of the silver carp and bighead carp.

Resequencing of 296 cultivated and wild lotus accessions unravels its evolution and breeding history.

Lotus (family: Nelumbonaceae) are perennial aquatic plants that represent one of the most ancient basal dicots. In the present study, we resequenced 296 lotus accessions from various geographical locations and germplasms to explore their genomic diversity and population structure. This germplasm set consisted of four accessions of American wild lotus and 292 accessions of Asian lotus, which were divided into four subgroups: wild, rhizome, flower and seed. Total single nucleotide polymorphisms (SNPs) suggested that the wild lotus had the highest variant number (7 191 010). Population structure and genome diversity analysis indicated that the American wild lotus demonstrated a distant genetic relationship with the Asian lotus. Furthermore, the seed and rhizome lotus groups had not originated from a single source but rather had a more complex multisource origin. Besides that, the seed lotus showed higher genetic diversity, which might have been due to the gene flow from the flower lotus to seed lotus by artificial crossing, and the rhizome lotus showed a much lower genetic diversity than the other groups. The present study provides SNP markers for lotus genomic diversity analysis, which will be useful for guiding lotus breeding.

Genome sequence of Gossypium herbaceum and genome updates of Gossypium arboreum and Gossypium hirsutum provide insights into cotton A-genome evolution.

Upon assembling the first Gossypium herbaceum (A1) genome and substantially improving the existing Gossypium arboreum (A2) and Gossypium hirsutum ((AD)1) genomes, we showed that all existing A-genomes may have originated from a common ancestor, referred to here as A0, which was more phylogenetically related to A1 than A2. Further, allotetraploid formation was shown to have preceded the speciation of A1 and A2. Both A-genomes evolved independently, with no ancestor-progeny relationship. Gaussian probability density function analysis indicates that several long-terminal-repeat bursts that occurred from 5.7 million years ago to less than 0.61 million years ago contributed compellingly to A-genome size expansion, speciation and evolution. Abundant species-specific structural variations in genic regions changed the expression of many important genes, which may have led to fiber cell improvement in (AD)1. Our findings resolve existing controversial concepts surrounding A-genome origins and provide valuable genomic resources for cotton genetic improvement.

Correction to: Viral integration drives multifocal HCC during the occult HBV infection.

In the original publication of this article [1], Fig. 6 is wrong and the updated figure is shown below.

Global transcriptome analysis of different stages of preimplantation embryo development in river buffalo.

BACKGROUND: Water buffalo (Bubalus bubalis) are divided into river buffalo and swamp buffalo subspecies and are essential livestock for agriculture and the local economy. Studies on buffalo reproduction have primarily focused on optimal fertility and embryonic mortality. There is currently limited knowledge on buffalo embryonic development, especially during the preimplantation period. Assembly of the river buffalo genome offers a reference for omics studies and facilitates transcriptomic analysis of preimplantation embryo development (PED). METHODS: We revealed transcriptomic profile of four stages (2-cell, 8-cell, Morula and Blastocyst) of PED via RNA-seq (Illumina HiSeq4000). Each stage comprised three biological replicates. The data were analyzed according to the basic RNA-seq analysis process. Ingenuity analysis of cell lineage control, especially transcription factor (TF) regulatory networks, was also performed. RESULTS: A total of 21,519 expressed genes and 67,298 transcripts were predicted from approximately 81.94 Gb of raw data. Analysis of transcriptome-wide expression, gene coexpression networks, and differentially expressed genes (DEGs) allowed for the characterization of gene-specific expression levels and relationships for each stage. The expression patterns of TFs, such as POU5F1, TEAD4, CDX4 and GATAs, were elucidated across diverse time series; most TF expression levels were increased during the blastocyst stage, during which time cell differentiation is initiated. All of these TFs were involved in the composition of the regulatory networks that precisely specify cell fate. These findings offer a deeper understanding of PED at the transcriptional level in the river buffalo.

Genome assembly of a tropical maize inbred line provides insights into structural variation and crop improvement.

Maize is one of the most important crops globally, and it shows remarkable genetic diversity. Knowledge of this diversity could help in crop improvement; however, gold-standard genomes have been elucidated only for modern temperate varieties. Here, we present a high-quality reference genome (contig N50 of 15.78 megabases) of the maize small-kernel inbred line, which is derived from a tropical landrace. Using haplotype maps derived from B73, Mo17 and SK, we identified 80,614 polymorphic structural variants across 521 diverse lines. Approximately 22% of these variants could not be detected by traditional single-nucleotide-polymorphism-based approaches, and some of them could affect gene expression and trait performance. To illustrate the utility of the diverse SK line, we used it to perform map-based cloning of a major effect quantitative trait locus controlling kernel weight-a key trait selected during maize improvement. The underlying candidate gene ZmBARELY ANY MERISTEM1d provides a target for increasing crop yields.

Viral integration drives multifocal HCC during the occult HBV infection.

BACKGROUND & AIMS: Although the prognosis of patients with occult hepatitis B virus (HBV) infection (OBI) is usually benign, a small portion may undergo cirrhosis and subsequently hepatocellular carcinoma (HCC). We studied the mechanism of life-long Integration of virus DNA into OBI host's genome, of which may induce hepatocyte transformation. METHODS: We applied HBV capture sequencing on single cells from an OBI patient who, developed multiple HCC tumors and underwent liver resection in May 2013 at Tongji Hospital in China. Despite with the undetectable virus DNA in serum, we determined the pattern of viral integration in tumor cells and adjacent non-tumor cells and obtained the details of the viral arrangement in host genome, and furthermore the HBV integrated region in cancer genome. RESULTS: HBV captured sequencing of tissues and individual cells revealed that samples from multiple tumors shared two viral integration sites that could affect three host genes, including CSMD2 on chr1 and MED30/EXT1 on chr8. Whole genome sequencing further indicated one hybrid chromosome formed by HBV integrations between chr1 and chr8 that was shared by multiple tumors. Additional 50 poorly differentiated liver tumors and the paired adjacent non-tumors were evaluated and functional studies suggested up-regulated EXT1 expression promoted HCC growth. We further observed that the most somatic mutations within the tumor cell genome were common among the multiple tumors, suggesting that HBV associated, multifocal HCC is monoclonal in origin. CONCLUSION: Through analyzing the HBV integration sites in multifocal HCC, our data suggested that the tumor cells were monoclonal in origin and formed in the absence of active viral replication, whereas the affected host genes may subsequently contribute to carcinogenesis.

Genetic Study on Small Insertions and Deletions in Psoriasis Reveals a Role in Complex Human Diseases.

Genetic studies based on single-nucleotide polymorphisms have provided valuable insights into the genetic architecture of complex diseases. However, a large fraction of heritability for most of these diseases remains unexplained, and the impact of small insertions and deletions (InDels) has been neglected. We performed a comprehensive screen on the exome sequence data of 1,326 genes using the SOAP-PopIndel method for InDels in 32,043 Chinese Han individuals and identified 29 unreported InDels within 25 susceptibility genes associated with psoriasis. Specifically, we identified 12 common, 9 low-frequency, and 8 rare InDels that explained approximately 1.29% of the heritability of psoriasis. Further analyses identified KIAA0319, RELN, NCAPG, ABO, AADACL2, LMAN1, FLG, HERC5, CCDC66, LEKR1, AFF3, ABCG2, ANXA7, SYTL2,GIPR, METTL1, and FYCO1 as unreported genes for psoriasis. In addition, identified InDels were associated with the following reported genes: IFIH1, ERAP1, ERAP2, LNPEP, UBLCP1, and STAT3; unreported independent associations for exonic InDels were found within GJB2 and ZNF816A. Our study enriched the genetic basis and pathogenesis of psoriasis and highlighted the non-negligible impact of InDels on complex human diseases.

Population genomic data reveal genes related to important traits of quail.

Background: Japanese quail (Coturnix japonica), a recently domesticated poultry species, is important not only as an agricultural product, but also as a model bird species for genetic research. However, most of the biological questions concerning genomics, phylogenetics, and genetics of some important economic traits have not been answered. It is thus necessary to complete a high-quality genome sequence as well as a series of comparative genomics, evolution, and functional studies. Results: Here, we present a quail genome assembly spanning 1.04 Gb with 86.63% of sequences anchored to 30 chromosomes (28 autosomes and 2 sex chromosomes Z/W). Our genomic data have resolved the long-term debate of phylogeny among Perdicinae (Japanese quail), Meleagridinae (turkey), and Phasianinae (chicken). Comparative genomics and functional genomic data found that four candidate genes involved in early maturation had experienced positive selection, and one of them encodes follicle stimulating hormone beta (FSHß), which is correlated with different FSHß levels in quail and chicken. We re-sequenced 31 quails (10 wild, 11 egg-type, and 10 meat-type) and identified 18 and 26 candidate selective sweep regions in the egg-type and meat-type lines, respectively. That only one of them is shared between egg-type and meat-type lines suggests that they were subject to an independent selection. We also detected a haplotype on chromosome Z, which was closely linked with maroon/yellow plumage in quail using population resequencing and a genome-wide association study. This haplotype block will be useful for quail breeding programs. Conclusions: This study provided a high-quality quail reference genome, identified quail-specific genes, and resolved quail phylogeny. We have identified genes related to quail early maturation and a marker for plumage color, which is significant for quail breeding. These results will facilitate biological discovery in quails and help us elucidate the evolutionary processes within the Phasianidae family.

Chromosome-scale assembly of the Monopterus genome.

Background: The teleost fish Monopterus albus is emerging as a new model for biological studies due to its natural sex transition and small genome, in addition to its enormous economic and potential medical value. However, no genomic information for the Monopterus is currently available. Findings: Here, we sequenced and de novo assembled the genome of M. albus and report the de novochromosome assembly by FISH walking assisted by conserved synteny (Cafs). Using Cafs, 328 scaffolds were assembled into 12 chromosomes, which covered genomic sequences of 555 Mb, accounting for 81.3% of the sequences assembled in scaffolds (∼689 Mb). A total of 18 ,660 genes were mapped on the chromosomes and showed a nonrandom distribution along chromosomes. Conclusions: We report the first reference genome of the Monopterus and provide an efficient Cafs strategy for a de novo chromosome-level assembly of the Monopterus genome, which provides a valuable resource, not only for further studies in genetics, evolution, and development, particularly sex determination, but also for breed improvement of the species.

Draft genome sequence of ramie, Boehmeria nivea (L.) Gaudich.

Ramie, Boehmeria nivea (L.) Gaudich, family Urticaceae, is a plant native to eastern Asia, and one of the world's oldest fibre crops. It is also used as animal feed and for the phytoremediation of heavy metal-contaminated farmlands. Thus, the genome sequence of ramie was determined to explore the molecular basis of its fibre quality, protein content and phytoremediation. For further understanding ramie genome, different paired-end and mate-pair libraries were combined to generate 134.31 Gb of raw DNA sequences using the Illumina whole-genome shotgun sequencing approach. The highly heterozygous B. nivea genome was assembled using the Platanus Genome Assembler, which is an effective tool for the assembly of highly heterozygous genome sequences. The final length of the draft genome of this species was approximately 341.9 Mb (contig N50 = 22.62 kb, scaffold N50 = 1,126.36 kb). Based on ramie genome annotations, 30,237 protein-coding genes were predicted, and the repetitive element content was 46.3%. The completeness of the final assembly was evaluated by benchmarking universal single-copy orthologous genes (BUSCO); 90.5% of the 1,440 expected embryophytic genes were identified as complete, and 4.9% were identified as fragmented. Phylogenetic analysis based on single-copy gene families and one-to-one orthologous genes placed ramie with mulberry and cannabis, within the clade of urticalean rosids. Genome information of ramie will be a valuable resource for the conservation of endangered Boehmeria species and for future studies on the biogeography and characteristic evolution of members of Urticaceae.

The draft genome of blunt snout bream (Megalobrama amblycephala) reveals the development of intermuscular bone and adaptation to herbivorous diet.

The blunt snout bream Megalobrama amblycephala is the economically most important cyprinid fish species. As an herbivore, it can be grown by eco-friendly and resource-conserving aquaculture. However, the large number of intermuscular bones in the trunk musculature is adverse to fish meat processing and consumption. As a first towards optimizing this aquatic livestock, we present a 1.116-Gb draft genome of M. amblycephala, with 779.54 Mb anchored on 24 linkage groups. Integrating spatiotemporal transcriptome analyses, we show that intermuscular bone is formed in the more basal teleosts by intramembranous ossification and may be involved in muscle contractibility and coordinating cellular events. Comparative analysis revealed that olfactory receptor genes, especially of the beta type, underwent an extensive expansion in herbivorous cyprinids, whereas the gene for the umami receptor T1R1 was specifically lost in M. amblycephala. The composition of gut microflora, which contributes to the herbivorous adaptation of M. amblycephala, was found to be similar to that of other herbivores. As a valuable resource for the improvement of M. amblycephala livestock, the draft genome sequence offers new insights into the development of intermuscular bone and herbivorous adaptation.