Taxonomic structure and functional association of foxtail millet root microbiome

Diverse heterochromatin-associated proteins repress distinct classes of genes and repetitive elements

Heterochromatin, typically marked by histone H3 trimethylation at lysine 9 (H3K9me3) or lysine 27 (H3K27me3), represses different protein-coding genes in different cells, as well as repetitive elements. The basis for locus specificity is unclear. Previously, we identified 172 proteins that are embedded in sonication-resistant heterochromatin (srHC) harbouring H3K9me3. Here, we investigate in humans how 97 of the H3K9me3-srHC proteins repress heterochromatic genes. We reveal four groups of srHC proteins that each repress many common genes and repeat elements. Two groups repress H3K9me3-embedded genes with different extents of flanking srHC, one group is specific for srHC genes with H3K9me3 and H3K27me3, and one group is specific for genes with srHC as the primary feature. We find that the enhancer of rudimentary homologue (ERH) is conserved from Schizosaccharomyces pombe in repressing meiotic genes and, in humans, now represses other lineage-specific genes and repeat elements. The study greatly expands our understanding of H3K9me3-based gene repression in vertebrates.

Cas9 targeted enrichment of mobile elements using nanopore sequencing

Mobile element insertions (MEIs) are repetitive genomic sequences that contribute to genetic variation and can lead to genetic disorders. Targeted and whole-genome approaches using short-read sequencing have been developed to identify reference and non-reference MEIs; however, the read length hampers detection of these elements in complex genomic regions. Here, we pair Cas9-targeted nanopore sequencing with computational methodologies to capture active MEIs in human genomes. We demonstrate parallel enrichment for distinct classes of MEIs, averaging 44% of reads on-targeted signals and exhibiting a 13.4-54x enrichment over whole-genome approaches. We show an individual flow cell can recover most MEIs (97% L1Hs, 93% AluYb, 51% AluYa, 99% SVA_F, and 65% SVA_E). We identify seventeen non-reference MEIs in GM12878 overlooked by modern, long-read analysis pipelines, primarily in repetitive genomic regions. This work introduces the utility of nanopore sequencing for MEI enrichment and lays the foundation for rapid discovery of elusive, repetitive genetic elements.

A Comparative Analysis on the Structure and Function of the Panax notoginseng Rhizosphere Microbiome

Panax notoginseng, an important Chinese medicinal herb, can be mainly cultivated in two planting patterns, cropland planting (DT) and understory planting (LX). We speculate that the rhizosphere microbiome may vary in DT and LX and may play an important role in promoting the growth and health of P. notoginseng. In the present study, culture-independent Illumina HiSeq was employed to investigate the rhizosphere bacteria and fungi under DT and LX planting patterns. Predominant phyla include Proteobacteria, Acidobacteria, Actinobacteria, Gemmatimonadetes, and Ascomycota in the two planting patterns. DT has higher alpha diversity index than LX. The predominant LX-core genera include Bradyrhizobium, Streptomyces, and Actinomadura, and the predominant DT-core genera include Sphingomonas, Variovorax, and Novosphingobium. Total relative abundance of the disease-suppression phylum (Proteobacteria, Firmicutes, and Actinobacteria) and the potential plant growth-promoting rhizobacteria (PGPR) were both significantly higher in LX than in DT. We also identified over-presented microbial functional traits mediating plant-microbe and microbe-microbe interactions, nutrition acquisition, and plant growth promotion in P. notoginseng rhizosphere. Our findings provide a valuable reference for studying beneficial microbes and pathogens of P. notoginseng planted in DT and LX.

Enabling sustainable agriculture through understanding and enhancement of microbiomes

Harnessing plant-associated microbiomes offers an invaluable strategy to help agricultural production become more sustainable while also meeting growing demands for food, feed and fiber. A plethora of interconnected interactions among the host, environment and microbes, occurring both above and below ground, drive recognition, recruitment and colonization of plant-associated microbes, resulting in activation of downstream host responses and functionality. Dissecting these complex interactions by integrating multiomic approaches, high-throughput culturing, and computational and synthetic biology advances is providing deeper understanding of the structure and function of native microbial communities. Such insights are paving the way towards development of microbial products as well as microbiomes engineered with synthetic microbial communities capable of delivering agronomic solutions. While there is a growing market for microbial-based solutions to improve crop productivity, challenges with commercialization of these products remain. The continued translation of plant-associated microbiome knowledge into real-world scenarios will require concerted transdisciplinary research, cross-training of a next generation of scientists, and targeted educational efforts to prime growers and the general public for successful adoption of these innovative technologies.

Insights into the taxonomic and functional characterization of agricultural crop core rhizobiomes and their potential microbial drivers

While our understanding of plant-microbe interactions in the rhizosphere microbiome (rhizobiome) has increased, there is still limited information on which taxa and functions drive these rhizobiome interactions. Focusing on the core rhizobiome (members common to two or more microbial assemblages) of crops may reduce the number of targets for determining these interactions, as they are expected to have greater influence on soil nutrient cycling and plant growth than the rest of the rhizobiome. Here, we examined whether the characterization of a core rhizobiome on the basis of only taxonomic or functional traits rather than the combined analysis of taxonomic and functional traits provides a different assessment of the core rhizobiome of agricultural crops. Sequences of the bacterial 16S rRNA gene from six globally important crops were analyzed using two different approaches in order to identify and characterize the taxonomic and functional core rhizobiome. For all crops examined, we found significant differences in the taxonomic and functional composition between the core rhizobiomes, and different phyla, genera, and predicted microbial functions were dominant depending on the core rhizobiome type. Network analysis indicated potentially important taxa were present in both taxonomic and functional core rhizobiomes. A subset of genera and predicted functions were exclusively or predominately present in only one type of core rhizobiome while others were detected in both core rhizobiomes. These results highlight the necessity of including both taxonomy and function when assessing the core rhizobiome, as this will enhance our understanding of the relationships between microbial taxa and soil health, plant growth, and agricultural sustainability.

A practical guide to amplicon and metagenomic analysis of microbiome data

Advances in high-throughput sequencing (HTS) have fostered rapid developments in the field of microbiome research, and massive microbiome datasets are now being generated. However, the diversity of software tools and the complexity of analysis pipelines make it difficult to access this field. Here, we systematically summarize the advantages and limitations of microbiome methods. Then, we recommend specific pipelines for amplicon and metagenomic analyses, and describe commonly-used software and databases, to help researchers select the appropriate tools. Furthermore, we introduce statistical and visualization methods suitable for microbiome analysis, including alpha- and beta-diversity, taxonomic composition, difference comparisons, correlation, networks, machine learning, evolution, source tracing, and common visualization styles to help researchers make informed choices. Finally, a step-by-step reproducible analysis guide is introduced. We hope this review will allow researchers to carry out data analysis more effectively and to quickly select the appropriate tools in order to efficiently mine the biological significance behind the data.

Variants at the ASIP locus contribute to coat color darkening in Nellore cattle

BACKGROUND

Nellore cattle (Bos indicus) are well-known for their adaptation to warm and humid environments. Hair length and coat color may impact heat tolerance. The Nellore breed has been strongly selected for white coat, but bulls generally exhibit darker hair ranging from light grey to black on the head, neck, hump, and knees. Given the potential contribution of coat color variation to the adaptation of cattle populations to tropical and sub-tropical environments, our aim was to map positional and functional candidate genetic variants associated with darkness of hair coat (DHC) in Nellore bulls.

RESULTS

We performed a genome-wide association study (GWAS) for DHC using data from 432 Nellore bulls that were genotyped for more than 777 k single nucleotide polymorphism (SNP) markers. A single major association signal was detected in the vicinity of the agouti signaling protein gene (ASIP). The analysis of whole-genome sequence (WGS) data from 21 bulls revealed functional variants that are associated with DHC, including a structural rearrangement involving ASIP (ASIP-SV1). We further characterized this structural variant using Oxford Nanopore sequencing data from 13 Australian Brahman heifers, which share ancestry with Nellore cattle; we found that this variant originates from a 1155-bp deletion followed by an insertion of a transposable element of more than 150 bp that may impact the recruitment of ASIP non-coding exons.

CONCLUSIONS

Our results indicate that the variant ASIP sequence causes darker coat pigmentation on specific parts of the body, most likely through a decreased expression of ASIP and consequently an increased production of eumelanin.

Single cell transcriptome sequencing on the Nanopore platform with ScNapBar

The current ecosystem of single cell RNA-seq platforms is rapidly expanding, but robust solutions for single cell and single molecule full- length RNA sequencing are virtually absent. A high-throughput solution that covers all aspects is necessary to study the complex life of mRNA on the single cell level. The Nanopore platform offers long read sequencing and can be integrated with the popular single cell sequencing method on the 10x Chromium platform. However, the high error-rate of Nanopore reads poses a challenge in downstream processing (e.g. for cell barcode assignment). We propose a solution to this particular problem by using a hybrid sequencing approach on Nanopore and Illumina platforms. Our software ScNapBar enables cell barcode assignment with high accuracy, especially if sequencing satura- tion is low. ScNapBar uses unique molecular identifier (UMI) or Naıve Bayes probabilistic approaches in the barcode assignment, depending on the available Illumina sequencing depth. We have benchmarked the two approaches on simulated and real Nanopore datasets. We further applied ScNapBar to pools of cells with an active or a silenced non-sense mediated RNA decay pathway. Our Nanopore read assignment distinguishes the respective cell populations and reveals characteristic nonsense-mediated mRNA decay events depending on cell status.

Brockarchaeota, a novel archaeal phylum with unique and versatile carbon cycling pathways

Geothermal environments, such as hot springs and hydrothermal vents, are hotspots for carbon cycling and contain many poorly described microbial taxa. Here, we reconstructed 15 archaeal metagenome-assembled genomes (MAGs) from terrestrial hot spring sediments in China and deep-sea hydrothermal vent sediments in Guaymas Basin, Gulf of California. Phylogenetic analyses of these MAGs indicate that they form a distinct group within the TACK superphylum, and thus we propose their classification as a new phylum, 'Brockarchaeota', named after Thomas Brock for his seminal research in hot springs. Based on the MAG sequence information, we infer that some Brockarchaeota are uniquely capable of mediating non-methanogenic anaerobic methylotrophy, via the tetrahydrofolate methyl branch of the Wood-Ljungdahl pathway and reductive glycine pathway. The hydrothermal vent genotypes appear to be obligate fermenters of plant-derived polysaccharides that rely mostly on substrate-level phosphorylation, as they seem to lack most respiratory complexes. In contrast, hot spring lineages have alternate pathways to increase their ATP yield, including anaerobic methylotrophy of methanol and trimethylamine, and potentially use geothermally derived mercury, arsenic, or hydrogen. Their broad distribution and their apparent anaerobic metabolic versatility indicate that Brockarchaeota may occupy previously overlooked roles in anaerobic carbon cycling.

Analysis of the gene transcription patterns and DNA methylation characteristics of triploid sea cucumbers (Apostichopus japonicus)

Breeding of polyploid aquatic animals is still an important approach and research hotspot for realizing the economic benefits afforded by the improvement of aquatic animal germplasm. To better understand the molecular mechanisms of the growth of triploid sea cucumbers, we performed gene expression and genome-wide comparisons of DNA methylation using the body wall tissue of triploid sea cucumbers using RNA-seq and MethylRAD-seq technologies. We clarified the expression pattern of triploid sea cucumbers and found no dosage effect. DEGs were significantly enriched in the pathways of nucleic acid and protein synthesis, cell growth, cell division, and other pathways. Moreover, we characterized the methylation pattern changes and found 615 differentially methylated genes at CCGG sites and 447 differentially methylated genes at CCWGG sites. Integrative analysis identified 23 genes (such as Guf1, SGT, Col5a1, HAL, HPS1, etc.) that exhibited correlations between promoter methylation and expression. Altered DNA methylation and expression of various genes suggested their roles and potential functional interactions in the growth of triploid sea cucumbers. Our data provide new insights into the epigenetic and transcriptomic alterations of the body wall tissue of triploid sea cucumbers and preliminarily elucidate the molecular mechanism of their growth, which is of great significance for the breeding of fine varieties of sea cucumbers.

Research advances in and prospects of ornamental plant genomics

The term 'ornamental plant' refers to all plants with ornamental value, which generally have beautiful flowers or special plant architectures. China is rich in ornamental plant resources and known as the "mother of gardens". Genomics is the science of studying genomes and is useful for carrying out research on genome evolution, genomic variations, gene regulation, and important biological mechanisms based on detailed genome sequence information. Due to the diversity of ornamental plants and high sequencing costs, the progress of genome research on ornamental plants has been slow for a long time. With the emergence of new sequencing technologies and a reduction in costs since the whole-genome sequencing of the first ornamental plant (Prunus mume) was completed in 2012, whole-genome sequencing of more than 69 ornamental plants has been completed in <10 years. In this review, whole-genome sequencing and resequencing of ornamental plants will be discussed. We provide analysis with regard to basic data from whole-genome studies of important ornamental plants, the regulation of important ornamental traits, and application prospects.

Microbiome-wide association studies reveal correlations between the structure and metabolism of the rhizosphere microbiome and disease resistance in cassava

Cassava is one of the most important staple food crops in tropical regions. To date, an understanding of the relationship between microbial communities and disease resistance in cassava has remained elusive. In order to explore the relationship among microbiome and phenotypes for further targeted design of microbial community, 16S rRNA and ITS of microbiome of ten cassava varieties were analysed, and a distinctive microbial community in the rhizosphere showed significant interdependence with disease resistance. Shotgun metagenome sequencing was performed to elucidate the structure of microbiomes of cassava rhizosphere. Comprehensive microbiome studies were performed to assess the correlation between the rhizosphere microbiome and disease resistance. Subsequently, the metagenome of rhizosphere microbiome was annotated to obtain taxonomic information at species level and identify metabolic pathways that were significantly associated with cassava disease resistance. Notably, cassava disease resistance was significantly associated with Lactococcus sp., which specifically produces nisin. To definitively explain the role of nisin and underlying mechanism, analysis of nisin biosynthesis-associated genes together with in vitro and in vivo experiments highlighted the effect of nisin on inhibiting the growth of Xanthomonas axonopodis pv. manihotis (Xam) and activating immune response in cassava. The new insights between cassava rhizosphere microbiome especially Lactococcus sp. and disease resistance provide valuable information into further control of cassava disease.

The thermal dependence and molecular basis of physiological color change in Takydromus septentrionalis (Lacertidae)

One of the main functions of physiological color change is thermoregulation. This change occurs much more rapidly than morphological color change, but the underlying mechanism remains poorly understood. Here, we studied the thermal dependence and molecular basis of physiological color change in lizards using Takydromus septentrionalis (Lacertidae) as the model system. Body color was thermally sensitive, becoming increasingly light as body temperatures deviated from the level (∼30°C) preferred by this species. We identified 3389 differentially expressed genes (DEGs) between lizards at 24°C and 30°C, and 1,097 DEGs between lizards at 36°C and 30°C. Temperature affected the cAMP signal pathway, motor proteins, cytoskeleton, and the expression of genes related to melanocyte-stimulating hormone (MSH) and melanocyte-concentrating hormone (MCH). Our data suggest that the role of physiological color change in thermoregulation is achieved in T. septentrionalis by altering the arrangement of pigments and thus the amount of solar radiation absorbed and reflected. G protein-coupling system inhibits adenylate cyclase activity to transform ATP into cAMP and thereby causes rapid pigment aggregation. MCH deactivates the G proteins and thereby initiates pigment dispersion. This mechanism differs from that reported for teleost fish where MCH activates the G proteins and thereby causes pigment aggregation.

This article has an associated First Person interview with the first author of the paper.

Summary: G protein-coupling system inhibits adenylate cyclase activity to transform ATP into cAMP and thereby causes rapid pigment aggregation. MCH deactivates the G proteins and thereby initiates pigment dispersion.

Determination of Microbial Diversity and Community Composition in Unfermented and Fermented Washing Rice Water by High-Throughput Sequencing

Washing rice water (WRW) refers to the sewage produced by rice washing in China and other parts of Asia people's daily life. As in the WRW is rich a variety of nutrients, microorganisms are prone to multiply and pollute the environment. In this article, high-throughput sequencing is used to describe the microbial diversity in different fermentation time WRW. The results showed that the sequencing depth effectively covered the microbial species in the samples, and the bacterial community structure in the samples of WRW at different fermentation periods was rich in diversity. Preominant taxa included Proteobacteria (62%), Firmicutes (28%), approximately Cyanobacteria (10%) and Bacteroidetes (0.5%). The core WRW microbiome comprises Trabulsiella, Pseudomonas, Serratia, Lactobacillus, Erwinia, Enterobacter, Clostridium and Acinetobacter, some of which are potential beneficial microbes. The change of microbial community composition with the change of habitat was assessed. It was found that environmental factors had significant influence on the assembly structure of microbial community.

Genome sequences reveal global dispersal routes and suggest convergent genetic adaptations in seahorse evolution

Seahorses have a circum-global distribution in tropical to temperate coastal waters. Yet, seahorses show many adaptations for a sedentary, cryptic lifestyle: they require specific habitats, such as seagrass, kelp or coral reefs, lack pelvic and caudal fins, and give birth to directly developed offspring without pronounced pelagic larval stage, rendering long-range dispersal by conventional means inefficient. Here we investigate seahorses' worldwide dispersal and biogeographic patterns based on a de novo genome assembly of Hippocampus erectus as well as 358 re-sequenced genomes from 21 species. Seahorses evolved in the late Oligocene and subsequent circum-global colonization routes are identified and linked to changing dynamics in ocean currents and paleo-temporal seaway openings. Furthermore, the genetic basis of the recurring "bony spines" adaptive phenotype is linked to independent substitutions in a key developmental gene. Analyses thus suggest that rafting via ocean currents compensates for poor dispersal and rapid adaptation facilitates colonizing new habitats.

Evaluation of full-length nanopore 16S sequencing for detection of pathogens in microbial keratitis

BACKGROUND

Microbial keratitis is a leading cause of preventable blindness worldwide. Conventional sampling and culture techniques are time-consuming, with over 40% of cases being culture-negative. Nanopore sequencing technology is portable and capable of generating long sequencing reads in real-time. The aim of this study is to evaluate the potential of nanopore sequencing directly from clinical samples for the diagnosis of bacterial microbial keratitis.

METHODS

Using full-length 16S rRNA amplicon sequences from a defined mock microbial community, we evaluated and benchmarked our bioinformatics analysis pipeline for taxonomic assignment on three different 16S rRNA databases (NCBI 16S RefSeq, RDP and SILVA) with clustering at 97%, 99% and 100% similarities. Next, we optimised the sample collection using an ex vivo porcine model of microbial keratitis to compare DNA recovery rates of 12 different collection methods: 21-gauge needle, PTFE membrane (4 mm and 6 mm), Isohelix™ SK-2S, Sugi® Eyespear, Cotton, Rayon, Dryswab™, Hydraflock®, Albumin-coated, Purflock®, Purfoam and Polyester swabs. As a proof-of-concept study, we then used the sampling technique that provided the highest DNA recovery, along with the optimised bioinformatics pipeline, to prospectively collected samples from patients with suspected microbial keratitis. The resulting nanopore sequencing results were then compared to standard microbiology culture methods.

RESULTS

We found that applying alignment filtering to nanopore sequencing reads and aligning to the NCBI 16S RefSeq database at 100% similarity provided the most accurate bacterial taxa assignment. DNA concentration recovery rates differed significantly between the collection methods (p < 0.001), with the Sugi® Eyespear swab providing the highest mean rank of DNA concentration. Then, applying the optimised collection method and bioinformatics pipeline directly to samples from two patients with suspected microbial keratitis, sequencing results from Patient A were in agreement with culture results, whilst Patient B, with negative culture results and previous antibiotic use, showed agreement between nanopore and Illumina Miseq sequencing results.

CONCLUSION

We have optimised collection methods and demonstrated a novel workflow for identification of bacterial microbial keratitis using full-length 16S nanopore sequencing.

Investigation of ancestral alleles in the Bovinae subfamily

BACKGROUND

In evolutionary theory, divergence and speciation can arise from long periods of reproductive isolation, genetic mutation, selection and environmental adaptation. After divergence, alleles can either persist in their initial state (ancestral allele - AA), co-exist or be replaced by a mutated state (derived alleles -DA). In this study, we aligned whole genome sequences of individuals from the Bovinae subfamily to the cattle reference genome (ARS.UCD-1.2) for defining ancestral alleles necessary for selection signatures study.

RESULTS

Accommodating independent divergent of each lineage from the initial ancestral state, AA were defined based on fixed alleles on at least two groups of yak, bison and gayal-gaur-banteng resulting in ~ 32.4 million variants. Using non-overlapping scanning windows of 10 Kb, we counted the AA observed within taurine and zebu cattle. We focused on the extreme points, regions with top 0. 1% (high count) and regions without any occurrence of AA (null count). High count regions preserved gene functions from ancestral states that are still beneficial in the current condition, while null counts regions were linked to mutated ones. For both cattle, high count regions were associated with basal lipid metabolism, essential for survival of various environmental pressures. Mutated regions were associated to productive traits in taurine, i.e. higher metabolism, cell development and behaviors and in immune response domain for zebu.

CONCLUSIONS

Our findings suggest that retaining and losing AA in some regions are varied and made it species-specific with possibility of overlapping as it depends on the selective pressure they had to experience.

A robust unsupervised machine-learning method to quantify the morphological heterogeneity of cells and nuclei

Cell morphology encodes essential information on many underlying biological processes. It is commonly used by clinicians and researchers in the study, diagnosis, prognosis, and treatment of human diseases. Quantification of cell morphology has seen tremendous advances in recent years. However, effectively defining morphological shapes and evaluating the extent of morphological heterogeneity within cell populations remain challenging. Here we present a protocol and software for the analysis of cell and nuclear morphology from fluorescence or bright-field images using the VAMPIRE algorithm ( https://github.com/kukionfr/VAMPIRE_open ). This algorithm enables the profiling and classification of cells into shape modes based on equidistant points along cell and nuclear contours. Examining the distributions of cell morphologies across automatically identified shape modes provides an effective visualization scheme that relates cell shapes to cellular subtypes based on endogenous and exogenous cellular conditions. In addition, these shape mode distributions offer a direct and quantitative way to measure the extent of morphological heterogeneity within cell populations. This protocol is highly automated and fast, with the ability to quantify the morphologies from 2D projections of cells seeded both on 2D substrates or embedded within 3D microenvironments, such as hydrogels and tissues. The complete analysis pipeline can be completed within 60 minutes for a dataset of ~20,000 cells/2,400 images.

Developmental validation of the MGIEasy Signature Identification Library Prep Kit, an all-in-one multiplex system for forensic applications

Analyzing genetic markers in nuclear and mitochondrial genomes is helpful in various forensic applications, such as individual identifications and kinship analyses. However, most commercial kits detect these markers separately, which is time-consuming, laborious, and more error-prone (mislabelling, contamination, ...). The MGIEasy Signature Identification Library Prep Kit (hereinafter "MGIEasy identification system"; MGI Tech, Shenzhen, China) has been designed to provide a simple, fast, and robust way to detect appropriate markers in one multiplex PCR reaction: 52 autosomal STRs, 27 X-chromosomal STRs, 48 Y-chromosomal STRs, 145 identity-informative SNPs, 53 ancestry-informative SNPs, 29 phenotype-informative SNPs, and the hypervariable regions of mitochondrial DNA (mtDNA). Here, we validated the performance of MGIEasy identification system following the guidelines of the Scientific Working Group on DNA Analysis Methods (SWGDAM), assessing species specificity, sensitivity, mixture identification, stability under non-optimal conditions (degraded samples, inhibitor contamination, and various substrates), repeatability, and concordance. Libraries prepared using MGIEasy identification system were sequenced on a MGISEQ-2000 instrument (MGI Tech). MGIEasy-derived STR, SNP, and mtDNA genotypes were highly concordant with CE-based STR genotypes (99.79%), MiSeq FGx-based SNP genotypes (99.78%), and Sanger-based mtDNA genotypes (100%), respectively. This system was strongly human-specific, resistant to four common PCR inhibitors, and reliably amplified both low quantities of DNA (as low as 0.125 ng) and degraded DNA (~ 150 nt). Most of the unique alleles from the minor contributor were detected in 1:10 male-female and male-male mixtures; some minor Y-STR alleles were even detected in 1:1000 male-female mixtures. MGIEasy also successfully directly amplified markers from blood stains on FTA cards, filter papers, and swabs. Thus, our results demonstrated that MGIEasy identification system was suitable for use in forensic analyses due to its robust and reliable performance on samples of varying quality and quantity.

QTL mapping for foxtail millet plant height in multi-environment using an ultra-high density bin map

Using a fixed RIL population derived from a widely used foxtail millet backbone breeding line and an elite cultivar, we constructed a high-density bin map and identified six novel multi-environment effect QTLs and seven candidate genes for dwarf phenotype. Plant height is an important trait that determines tradeoffs between competition and resource allocation, which is crucial for yield potential. To improve the C4 model plant foxtail millet (Setaria italica) productivity, it is necessary to isolate plant height-related genes that contribute to ideal plant architecture in breeding. In the present study, we generated a foxtail millet population of 333 recombinant inbred lines (RILs) derived from a cross between a backbone line Ai 88 and an elite cultivar Liaogu 1. We evaluated plant height in 13 environmental conditions across 4 years, the mean plant height of the RIL population ranged from 89.5 to 149.9 cm. Using deep re-sequencing data, we constructed a high-density bin map with 3744 marker bins. Quantitative trait locus (QTL) mapping identified 26 QTLs significantly associated with plant height. Of these, 13 QTLs were repeatedly detected under multiple environments, including six novel QTLs that have not been reported before. Seita.1G242300, a gene encodes gibberellin 2-oxidase-8, which was detected in nine environments in a 1.54-Mb interval of qPH1.3, was considered as an important candidate gene. Moreover, other six genes involved in GA biosynthesis or signaling pathways, and fifteen genes encode F-box domain proteins which might function as E3 ligases, were also considered as candidate genes in different QTLs. These QTLs and candidate genes identified in this study will help to elucidate the genetic basis of foxtail millet plant height, and the linked markers will be useful for marker-assistant selection of varieties with ideal plant architecture and high yield potential.

Chromosome-level genome assembly of Ophiorrhiza pumila reveals the evolution of camptothecin biosynthesis

Plant genomes remain highly fragmented and are often characterized by hundreds to thousands of assembly gaps. Here, we report chromosome-level reference and phased genome assembly of Ophiorrhiza pumila, a camptothecin-producing medicinal plant, through an ordered multi-scaffolding and experimental validation approach. With 21 assembly gaps and a contig N50 of 18.49 Mb, Ophiorrhiza genome is one of the most complete plant genomes assembled to date. We also report 273 nitrogen-containing metabolites, including diverse monoterpene indole alkaloids (MIAs). A comparative genomics approach identifies strictosidine biogenesis as the origin of MIA evolution. The emergence of strictosidine biosynthesis-catalyzing enzymes precede downstream enzymes' evolution post γ whole-genome triplication, which occurred approximately 110 Mya in O. pumila, and before the whole-genome duplication in Camptotheca acuminata identified here. Combining comparative genome analysis, multi-omics analysis, and metabolic gene-cluster analysis, we propose a working model for MIA evolution, and a pangenome for MIA biosynthesis, which will help in establishing a sustainable supply of camptothecin.

Quantifying the influence of mutation detection on tumour subclonal reconstruction

Whole-genome sequencing can be used to estimate subclonal populations in tumours and this intra-tumoural heterogeneity is linked to clinical outcomes. Many algorithms have been developed for subclonal reconstruction, but their variabilities and consistencies are largely unknown. We evaluate sixteen pipelines for reconstructing the evolutionary histories of 293 localized prostate cancers from single samples, and eighteen pipelines for the reconstruction of 10 tumours with multi-region sampling. We show that predictions of subclonal architecture and timing of somatic mutations vary extensively across pipelines. Pipelines show consistent types of biases, with those incorporating SomaticSniper and Battenberg preferentially predicting homogenous cancer cell populations and those using MuTect tending to predict multiple populations of cancer cells. Subclonal reconstructions using multi-region sampling confirm that single-sample reconstructions systematically underestimate intra-tumoural heterogeneity, predicting on average fewer than half of the cancer cell populations identified by multi-region sequencing. Overall, these biases suggest caution in interpreting specific architectures and subclonal variants.

Whole genome variants across 57 pig breeds enable comprehensive identification of genetic signatures that underlie breed features

Background

A large number of pig breeds are distributed around the world, their features and characteristics vary among breeds, and they are valuable resources. Understanding the underlying genetic mechanisms that explain across-breed variation can help breeders develop improved pig breeds.

Results

In this study, we performed GWAS using a standard mixed linear model with three types of genome variants (SNP, InDel, and CNV) that were identified from public, whole-genome, sequencing data sets. We used 469 pigs of 57 breeds, and we identified and analyzed approximately 19 million SNPs, 1.8 million InDels, and 18,016 CNVs. We defined six biological phenotypes by the characteristics of breed features to identify the associated genome variants and candidate genes, which included coat color, ear shape, gradient zone, body weight, body length, and body height. A total of 37 candidate genes was identified, which included 27 that were reported previously (e.g., PLAG1 for body weight), but the other 10 were newly detected candidate genes (e.g., ADAMTS9 for coat color).

Conclusion

Our study indicated that using GWAS across a modest number of breeds with high density genome variants provided efficient mapping of complex traits.

Supplementary Information

Supplementary information accompanies this paper at 10.1186/s40104-020-00520-8.

Identification and validation of 174 COVID-19 vaccine candidate epitopes reveals low performance of common epitope prediction tools

The outbreak of SARS-CoV-2 (2019-nCoV) virus has highlighted the need for fast and efficacious vaccine development. Stimulation of a proper immune response that leads to protection is highly dependent on presentation of epitopes to circulating T-cells via the HLA complex. SARS-CoV-2 is a large RNA virus and testing of all of its overlapping peptides in vitro to deconvolute an immune response is not feasible. Therefore HLA-binding prediction tools are often used to narrow down the number of peptides to test. We tested NetMHC suite tools' predictions by using an in vitro peptide-MHC stability assay. We assessed 777 peptides that were predicted to be good binders across 11 MHC alleles in a complex-stability assay and tested a selection of 19 epitope-HLA-binding prediction tools against the assay. In this investigation of potential SARS-CoV-2 epitopes we found that current prediction tools vary in performance when assessing binding stability, and they are highly dependent on the MHC allele in question. Designing a COVID-19 vaccine where only a few epitope targets are included is therefore a very challenging task. Here, we present 174 SARS-CoV-2 epitopes with high prediction binding scores, validated to bind stably to 11 HLA alleles. Our findings may contribute to the design of an efficacious vaccine against COVID-19.

High quality genome assemblies of Mycoplasma bovis using a taxon-specific Bonito basecaller for MinION and Flongle long-read nanopore sequencing

BACKGROUND

Implementation of Third-Generation Sequencing approaches for Whole Genome Sequencing (WGS) all-in-one diagnostics in human and veterinary medicine, requires the rapid and accurate generation of consensus genomes. Over the last years, Oxford Nanopore Technologies (ONT) released various new devices (e.g. the Flongle R9.4.1 flow cell) and bioinformatics tools (e.g. the in 2019-released Bonito basecaller), allowing cheap and user-friendly cost-efficient introduction in various NGS workflows. While single read, overall consensus accuracies, and completeness of genome sequences has been improved dramatically, further improvements are required when working with non-frequently sequenced organisms like Mycoplasma bovis. As an important primary respiratory pathogen in cattle, rapid M. bovis diagnostics is crucial to allow timely and targeted disease control and prevention. Current complete diagnostics (including identification, strain typing, and antimicrobial resistance (AMR) detection) require combined culture-based and molecular approaches, of which the first can take 1-2 weeks. At present, cheap and quick long read all-in-one WGS approaches can only be implemented if increased accuracies and genome completeness can be obtained.

RESULTS

Here, a taxon-specific custom-trained Bonito v.0.1.3 basecalling model (custom-pg45) was implemented in various WGS assembly bioinformatics pipelines. Using MinION sequencing data, we showed improved consensus accuracies up to Q45.2 and Q46.7 for reference-based and Canu de novo assembled M. bovis genomes, respectively. Furthermore, the custom-pg45 model resulted in mean consensus accuracies of Q45.0 and genome completeness of 94.6% for nine M. bovis field strains. Improvements were also observed for the single-use Flongle sequencer (mean Q36.0 accuracies and 80.3% genome completeness).

CONCLUSIONS

These results implicate that taxon-specific basecalling of MinION and single-use Flongle Nanopore long reads are of great value to be implemented in rapid all-in-one WGS tools as evidenced for Mycoplasma bovis as an example.

Deciphering the Microbial Taxonomy and Functionality of Two Diverse Mangrove Ecosystems and Their Potential Abilities To Produce Bioactive Compounds

This study comprehensively described the taxonomy and functionality of mangrove microbiomes, including their capacity for secondary metabolite biosynthesis and their ability to resist antibiotics. The microbial taxonomic and functional characteristics differed between geographical locations, corresponding to the environmental condition of two diverse mangrove regions. A large number of microbial biosynthetic gene clusters encoding novel bioactivities were found, and this can serve as a valuable resource to guide novel bioactive compound discovery for potential clinical uses.

Mangroves, as important and special ecosystems, create unique ecological environments for examining the microbial gene capacity and potential for producing bioactive compounds. However, little is known about the biogeochemical implications of microbiomes in mangrove ecosystems, especially the variations between pristine and anthropogenic mangroves. To elucidate this, we investigated the microbial taxonomic and functional shifts of the mangrove microbiomes and their potential for bioactive compounds in two different coastal mangrove ecosystems in southern China. A gene catalogue, including 87 million unique genes, was constructed, based on deep shotgun metagenomic sequencing. Differentially enriched bacterial and archaeal taxa between pristine mangroves (Guangxi) and anthropogenic mangroves (Shenzhen) were found. The Nitrospira and ammonia-oxidizing archaea, specifically, were more abundant in Shenzhen mangroves, while sulfate-reducing bacteria and methanogens were more abundant in Guangxi mangroves. The results of functional analysis were consistent with the taxonomic results, indicating that the Shenzhen mangrove microbiome has a higher abundance of genes involved in nitrogen metabolism while the Guangxi mangrove microbiome has a higher capacity for sulfur metabolism and methanogenesis. Biosynthetic gene clusters were identified in the metagenome data and in hundreds of de novo reconstructed nonredundant microbial genomes, respectively. Notably, we found different biosynthetic potential in different taxa, and we identified three high quality and novel Acidobacteria genomes with a large number of BGCs. In total, 67,278 unique genes were annotated with antibiotic resistance, indicating the prevalence and persistence in multidrug-resistant genes in the mangrove microbiome.

IMPORTANCE This study comprehensively described the taxonomy and functionality of mangrove microbiomes, including their capacity for secondary metabolite biosynthesis and their ability to resist antibiotics. The microbial taxonomic and functional characteristics differed between geographical locations, corresponding to the environmental condition of two diverse mangrove regions. A large number of microbial biosynthetic gene clusters encoding novel bioactivities were found, and this can serve as a valuable resource to guide novel bioactive compound discovery for potential clinical uses.

Genomic and transcriptomic analysis unveils population evolution and development of pesticide resistance in fall armyworm Spodoptera frugiperda

The fall armyworm (FAW), Spodoptera frugiperda, is a destructive pest native to America and has recently become an invasive insect pest in China. Because of its rapid spread and great risks in China, understanding of FAW genetic background and pesticide resistance is urgent and essential to develop effective management strategies. Here, we assembled a chromosome-level genome of a male FAW (SFynMstLFR) and compared re-sequencing results of the populations from America, Africa, and China. Strain identification of 163 individuals collected from America, Africa and China showed that both C and R strains were found in the American populations, while only C strain was found in the Chinese and African populations. Moreover, population genomics analysis showed that populations from Africa and China have close relationship with significantly genetic differentiation from American populations. Taken together, FAWs invaded into China were most likely originated from Africa. Comparative genomics analysis displayed that the cytochrome p450 gene family is extremely expanded to 425 members in FAW, of which 283 genes are specific to FAW. Treatments of Chinese populations with twenty-three pesticides showed the variant patterns of transcriptome profiles, and several detoxification genes such as AOX, UGT and GST specially responded to the pesticides. These findings will be useful in developing effective strategies for management of FAW in China and other invaded areas.

Chromosome-level genome assembly of a parent species of widely cultivated azaleas

Azaleas (Ericaceae) comprise one of the most diverse ornamental plants, renowned for their cultural and economic importance. We present a chromosome-scale genome assembly for Rhododendron simsii, the primary ancestor of azalea cultivars. Genome analyses unveil the remnants of an ancient whole-genome duplication preceding the radiation of most Ericaceae, likely contributing to the genomic architecture of flowering time. Small-scale gene duplications contribute to the expansion of gene families involved in azalea pigment biosynthesis. We reconstruct entire metabolic pathways for anthocyanins and carotenoids and their potential regulatory networks by detailed analysis of time-ordered gene co-expression networks. MYB, bHLH, and WD40 transcription factors may collectively regulate anthocyanin accumulation in R. simsii, particularly at the initial stages of flower coloration, and with WRKY transcription factors controlling progressive flower coloring at later stages. This work provides a cornerstone for understanding the underlying genetics governing flower timing and coloration and could accelerate selective breeding in azalea.

Azaleas are one of the most diverse ornamental plants and have cultural and economic importance. Here, the authors report a chromosome-scale genome assembly for the primary ancestor of the azalea cultivar Rhododendro simsi and identify transcription factors that may function in flower coloration at different stages.

Analysis of muntjac deer genome and chromatin architecture reveals rapid karyotype evolution

Closely related muntjac deer show striking karyotype differences. Here we describe chromosome-scale genome assemblies for Chinese and Indian muntjacs, Muntiacus reevesi (2n = 46) and Muntiacus muntjak vaginalis (2n = 6/7), and analyze their evolution and architecture. The genomes show extensive collinearity with each other and with other deer and cattle. We identified numerous fusion events unique to and shared by muntjacs relative to the cervid ancestor, confirming many cytogenetic observations with genome sequence. One of these M. muntjak fusions reversed an earlier fission in the cervid lineage. Comparative Hi-C analysis showed that the chromosome fusions on the M. muntjak lineage altered long-range, three-dimensional chromosome organization relative to M. reevesi in interphase nuclei including A/B compartment structure. This reshaping of multi-megabase contacts occurred without notable change in local chromatin compaction, even near fusion sites. A few genes involved in chromosome maintenance show evidence for rapid evolution, possibly associated with the dramatic changes in karyotype.

A novel forensic panel of 186-plex SNPs and 123-plex STR loci based on massively parallel sequencing

The MiSeq® FGX Forensic system and the HID-Ion AmpliSeq Panel were previously developed for massively parallel sequencing (MPS) for forensic casework. Among the three major sequencing platforms, BGISEQ-500^TM, which is based on multiple PCRs, is still lacking in forensics. Here, a novel forensic panel was constructed to detect 186 single-nucleotide polymorphisms (SNPs) and 123 short tandem repeats (STRs) with MPS technology on the BGISEQ-500™ platform. First, the library preparation, sequencing process, and data analysis were performed, focusing on the average depth of coverage and heterozygote balance. We calculated the allelic frequencies and forensic parameters of STR and SNP loci in 73 unrelated Chinese Han individuals. In addition, performance was evaluated with accuracy, uniformity, sensitivity, PCR inhibitor, repeatability and reproducibility, mixtures, degraded samples, case-type samples, and pedigree analyses. The results showed that 100% accurate and concordant genotypes can be obtained, and the loci with an abundance in the interquartile range accounted for 92.90% of the total, suggesting reliable uniformity in this panel. We obtained a locus detection rate that was higher than 98.78% from 78 pg of input DNA, and the optimal amount was 1.25-10 ng. The maximum concentrations of hematin and humic acid were 200 and 100 μM, respectively (the ratios of detected loci were 96.52% and 92.41%), in this panel. As a mixture, compared with those of SNPs, minor-contributor alleles of STRs could be detected at higher levels. For the degraded sample, the ratio of detected loci was 98.41%, and most profiles from case-type samples were not significantly different in abundance in our studies. As a whole, this panel showed high-performance, reliable, robust, repeatable, and reproducible results, which are sufficient for paternity testing, individual identification, and use for potentially degraded samples in forensic science.

Hypoxic gene expression in chronic hepatitis B virus infected patients is not observed in state-of-the-art in vitro and mouse infection models

Hepatitis B virus (HBV) is the leading cause of hepatocellular carcinoma (HCC) worldwide. The prolyl hydroxylase domain (PHD)-hypoxia inducible factor (HIF) pathway is a key mammalian oxygen sensing pathway and is frequently perturbed by pathological states including infection and inflammation. We discovered a significant upregulation of hypoxia regulated gene transcripts in patients with chronic hepatitis B (CHB) in the absence of liver cirrhosis. We used state-of-the-art in vitro and in vivo HBV infection models to evaluate a role for HBV infection and the viral regulatory protein HBx to drive HIF-signalling. HBx had no significant impact on HIF expression or associated transcriptional activity under normoxic or hypoxic conditions. Furthermore, we found no evidence of hypoxia gene expression in HBV de novo infection, HBV infected human liver chimeric mice or transgenic mice with integrated HBV genome. Collectively, our data show clear evidence of hypoxia gene induction in CHB that is not recapitulated in existing models for acute HBV infection, suggesting a role for inflammatory mediators in promoting hypoxia gene expression.

Composition and function of rhizosphere microbiome of Panax notoginseng with discrepant yields

Background

Panax notoginseng is a highly valuable medicinal plant. Reduced P. notoginseng yield is a common and serious problem that arises in a continuous cropping system. Variation in the composition and function of soil microbial community is considered the primary cause of yield reduction.

Methods

This study used shotgun metagenomic sequencing approaches to describe the taxonomic and functional features of P. notoginseng rhizosphere microbiome and screen microbial taxa and functional traits related to yields.

Results

At the family and genus level, a total of 43 families and 45 genera (relative abundance > 0.1%) were obtained, and the correlation with the yield of P. notoginseng was further analyzed. Nitrosomonadaceae, Xanthomonadaceae, Mycobacterium and Arthrobacter that were enriched in soils with higher yields were positively correlated with P. notoginseng yields, thereby suggesting that they might increase yields. Negative correlation coefficients indicated that Xanthobacteraceae, Caulobacteraceae, Oxalobacteraceae, Chitinophagaceae, Sphingomonas, Hyphomicrobium, Variovorax and Phenylobacterium might be detrimental to P. notoginseng growth. A total of 85 functional traits were significantly (P < 0.05) correlated with P. notoginseng yields. Functional traits, likely steroid biosynthesis and MAPK signaling pathway were positively correlated with P. notoginseng yields. In contrast, functional traits, such as bacterial secretion system, ABC transporters, metabolism of xenobiotics by cytochrome P450 and drug metabolism–cytochrome P450, were negatively associated with yields.

Conclusions

This study describes an overview of the rhizosphere microbiome of P. notoginseng with discrepant yields and identifies the taxa and functional traits related to yields. Our results provide valuable information to guide the isolation and culture of potentially beneficial microorganisms and to utilize the power of the microbiome to increase plant yields in a continuous cropping system.

The scaffold-forming steps of plant alkaloid biosynthesis

Alkaloids from plants are characterised by structural diversity and bioactivity, and maintain a privileged position in both modern and traditional medicines. In recent years, there have been significant advances in elucidating the biosynthetic origins of plant alkaloids. In this review, I will describe the progress made in determining the metabolic origins of the so-called true alkaloids, specialised metabolites derived from amino acids containing a nitrogen heterocycle. By identifying key biosynthetic steps that feature in the majority of pathways, I highlight the key roles played by modifications to primary metabolism, iminium reactivity and spontaneous reactions in the molecular and evolutionary origins of these pathways.

Comparing and phylogenetic analysis chloroplast genome of three Achyranthes species

In this study, the chloroplast genome sequencing of the Achyranthes longifolia, Achyranthes bidentata and Achyranthes aspera were performed by Next-generation sequencing technology. The results revealed that there were a length of 151,520 bp (A. longifolia), 151,284 bp (A. bidentata), 151,486 bp (A. aspera), respectively. These chloroplast genome have a highly conserved structure with a pair of inverted repeat (IR) regions (25,150 bp; 25,145 bp; 25,150 bp), a large single copy (LSC) regions (83,732 bp; 83,933 bp; 83,966 bp) and a small single copy (SSC) regions (17,252 bp; 17,263 bp; 17,254 bp) in A. bidentate, A. aspera and A. longifolia. There were 127 genes were annotated, which including 8 rRNA genes, 37 tRNA genes and 82 functional genes. The phylogenetic analysis strongly revealed that Achyranthes is monophyletic, and A. bidentata was the closest relationship with A. aspera and A. longifolia. A. bidentata and A. longifolia were clustered together, the three Achyranthes species had the same origin, then the gunes of Achyranthes is the closest relative to Alternanthera, and that forms a group with Alternanthera philoxeroides. The research laid a foundation and provided relevant basis for the identification of germplasm resources in the future.

Whole-genome, transcriptome, and methylome analyses provide insights into the evolution of platycoside biosynthesis in Platycodon grandiflorus, a medicinal plant

Triterpenoid saponins (TSs) are common plant defense phytochemicals with potential pharmaceutical properties. Platycodon grandiflorus (Campanulaceae) has been traditionally used to treat bronchitis and asthma in East Asia. The oleanane-type TSs, platycosides, are a major component of the P. grandiflorus root extract. Recent studies show that platycosides exhibit anti-inflammatory, antiobesity, anticancer, antiviral, and antiallergy properties. However, the evolutionary history of platycoside biosynthesis genes remains unknown. In this study, we sequenced the genome of P. grandiflorus and investigated the genes involved in platycoside biosynthesis. The draft genome of P. grandiflorus is 680.1 Mb long and contains 40,017 protein-coding genes. Genomic analysis revealed that the CYP716 family genes play a major role in platycoside oxidation. The CYP716 gene family of P. grandiflorus was much larger than that of other Asterid species. Orthologous gene annotation also revealed the expansion of β-amyrin synthases (bASs) in P. grandiflorus, which was confirmed by tissue-specific gene expression. In these expanded gene families, we identified key genes showing preferential expression in roots and association with platycoside biosynthesis. In addition, whole-genome bisulfite sequencing showed that CYP716 and bAS genes are hypomethylated in P. grandiflorus, suggesting that epigenetic modification of these two gene families affects platycoside biosynthesis. Thus whole-genome, transcriptome, and methylome data of P. grandiflorus provide novel insights into the regulation of platycoside biosynthesis by CYP716 and bAS gene families.

COVID-19 pandemic reveals the peril of ignoring metadata standards

Efficient response to the pandemic through the mobilization of the larger scientific community is challenged by the limited reusability of the available primary genomic data. Here, the Genomic Standards Consortium board highlights the essential need for contextual genomic data FAIRness, for empowering key data-driven biological questions.

The Effect of Plant Geographical Location and Developmental Stage on Root-Associated Microbiomes of Gymnadenia conopsea

Gymnadenia conopsea (L.) R. Br. is an important perennial terrestrial photosynthetic orchid species whose microbiomes are considered to play an important role in helping its germination and growth. However, the assemblage of G. conopsea root-associated microbial communities is poorly understood. The compositions of fungal and bacterial communities from the roots and corresponding soil samples in G. conopsea across distinct biogeographical regions from two significantly different altitudes were characterized at the vegetative and reproductive growth stages. The geographical location, developmental stage and compartment were factors contributing to microbiome variation in G. conopsea. Predominant fungal taxa include Ascomycota, Basidiomycota, Mortierellomycota and Chytridiomycota, whereas Proteobacteria, Bacteroidetes, Acidobacteria, Actinobacteria, Verrucomicrobia, Chloroflexi, TM7 and Planctomycetes were predominant bacterial taxa. Using G. conopsea as a model, the structural and functional composition in G. conopsea root-associated microbiomes were comprehensive analyzed. Contrary to previous studies, biogeography was the main factor influencing the microbial community in this study. Besides, compartment and developmental stage should also be considered to analyze the variation of microbiota composition. Although the microbial composition varied greatly by location, the symbiotic microorganisms of G. conopsea still have certain specificity. This study gives an abundant information of G. conopsea root-associated microbiomes and provides new clues to better understanding the factors affecting the composition and diversity of fungal/bacterial communities associated with orchids. Our results also laying a foundation for harnessing the microbiome for sustainable G. conopsea cultivation. Moreover, these results might be generally applicable to other orchidaceae plants.

Tandem gene duplications drive divergent evolution of caffeine and crocin biosynthetic pathways in plants

Background

Plants have evolved a panoply of specialized metabolites that increase their environmental fitness. Two examples are caffeine, a purine psychotropic alkaloid, and crocins, a group of glycosylated apocarotenoid pigments. Both classes of compounds are found in a handful of distantly related plant genera (Coffea, Camellia, Paullinia, and Ilex for caffeine; Crocus, Buddleja, and Gardenia for crocins) wherein they presumably evolved through convergent evolution. The closely related Coffea and Gardenia genera belong to the Rubiaceae family and synthesize, respectively, caffeine and crocins in their fruits.

Results

Here, we report a chromosomal-level genome assembly of Gardenia jasminoides, a crocin-producing species, obtained using Oxford Nanopore sequencing and Hi-C technology. Through genomic and functional assays, we completely deciphered for the first time in any plant the dedicated pathway of crocin biosynthesis. Through comparative analyses with Coffea canephora and other eudicot genomes, we show that Coffea caffeine synthases and the first dedicated gene in the Gardenia crocin pathway, GjCCD4a, evolved through recent tandem gene duplications in the two different genera, respectively. In contrast, genes encoding later steps of the Gardenia crocin pathway, ALDH and UGT, evolved through more ancient gene duplications and were presumably recruited into the crocin biosynthetic pathway only after the evolution of the GjCCD4a gene.

Conclusions

This study shows duplication-based divergent evolution within the coffee family (Rubiaceae) of two characteristic secondary metabolic pathways, caffeine and crocin biosynthesis, from a common ancestor that possessed neither complete pathway. These findings provide significant insights on the role of tandem duplications in the evolution of plant specialized metabolism.

An integrated Asian human SNV and indel benchmark established using multiple sequencing methods

Sequencing technologies have been rapidly developed recently, leading to the breakthrough of sequencing-based clinical diagnosis, but accurate and complete genome variation benchmark would be required for further assessment of precision medicine applications. Despite the human cell line of NA12878 has been successfully developed to be a variation benchmark, population-specific variation benchmark is still lacking. Here, we established an Asian human variation benchmark by constructing and sequencing a stabilized cell line of a Chinese Han volunteer. By using seven different sequencing strategies, we obtained ~3.88 Tb clean data from different laboratories, hoping to reach the point of high sequencing depth and accurate variation detection. Through the combination of variations identified from different sequencing strategies and different analysis pipelines, we identified 3.35 million SNVs and 348.65 thousand indels, which were well supported by our sequencing data and passed our strict quality control, thus should be high confidence variation benchmark. Besides, we also detected 5,913 high-quality SNVs which had 969 sites were novel and  located in the high homologous regions supported by long-range information in both the co-barcoding single tube Long Fragment Read (stLFR) data and PacBio HiFi CCS data. Furthermore, by using the long reads data (stLFR and HiFi CCS), we were able to phase more than 99% heterozygous SNVs, which helps to improve the benchmark to be haplotype level. Our study provided comprehensive sequencing data as well as the integrated variation benchmark of an Asian derived cell line, which would be valuable for future sequencing-based clinical development.

Applying genomic resources to accelerate wheat biofortification

Wheat has low levels of the micronutrients iron and zinc in the grain, which contributes to 2 billion people suffering from micronutrient deficiency globally. While wheat flour is commonly fortified during processing, an attractive and more sustainable solution is biofortification, which could improve micronutrient content in the human diet, without the sustainability issues and costs associated with conventional fortification. Although many studies have used quantitative trait loci mapping and genome-wide association to identify genetic loci to improve micronutrient contents, recent developments in genomics offer an opportunity to accelerate marker discovery and use gene-focussed approaches to engineer improved micronutrient content in wheat. The recent publication of a high-quality wheat genome sequence, alongside gene expression atlases, variation datasets and sequenced mutant populations, provides a foundation to identify genetic loci and genes controlling micronutrient content in wheat. We discuss how novel genomic resources can identify candidate genes for biofortification, integrating knowledge from other cereal crops, and how these genes can be tested using gene editing, transgenic and TILLING approaches. Finally, we highlight key challenges remaining to develop wheat cultivars with high levels of iron and zinc.

Predominance of soil vs root effect in rhizosphere microbiota reassembly

Rhizosphere community assembly is simultaneously affected by both plants and bulk soils and is vital for plant health. However, it is still unclear how and to what extent disease-suppressive rhizosphere microbiota can be constructed from bulk soil, and the underlying agents involved in the process that render the rhizosphere suppressive against pathogenic microbes remain elusive. In this study, the evolutionary processes of the rhizosphere microbiome were explored based on transplanting plants previously growing in distinct disease-incidence soils to one disease-suppressive soil. Our results showed that distinct rhizoplane bacterial communities were assembled on account of the original bulk soil communities with different disease incidences. Furthermore, the bacterial communities in the transplanted rhizosphere were noticeably influenced by the second disease-suppressive microbial pool, rather than that of original formed rhizoplane microbiota and homogenous nontransplanted rhizosphere microbiome, contributing to a significant decrease in the pathogen population. In addition, Spearman's correlations between relative abundances of bacterial taxa and the abundance of Ralstonia solanacearum indicated Anoxybacillus, Flavobacterium, Permianibacter and Pseudomonas were predicted to be associated with disease-suppressive function formation. Altogether, our results showed that bulk soil played an important role in the process of assembling and reassembling the rhizosphere microbiome of plants.

Genome sequence and comparative analysis of reindeer (Rangifer tarandus) in northern Eurasia

Reindeer are semi-domesticated ruminants that have adapted to the challenging northern Eurasian environment characterized by long winters and marked annual fluctuations in daylight. We explored the genetic makeup behind their unique characteristics by de novo sequencing the genome of a male reindeer and conducted gene family analyses with nine other mammalian species. We performed a population genomics study of 23 additional reindeer representing both domestic and wild populations and several ecotypes from various geographic locations. We assembled 2.66 Gb (N50 scaffold of 5 Mb) of the estimated 2.92 Gb reindeer genome, comprising 27,332 genes. The results from the demographic history analysis suggested marked changes in the effective population size of reindeer during the Pleistocene period. We detected 160 reindeer-specific and expanded genes, of which zinc finger proteins (n = 42) and olfactory receptors (n = 13) were the most abundant. Comparative genome analyses revealed several genes that may have promoted the adaptation of reindeer, such as those involved in recombination and speciation (PRDM9), vitamin D metabolism (TRPV5, TRPV6), retinal development (PRDM1, OPN4B), circadian rhythm (GRIA1), immunity (CXCR1, CXCR2, CXCR4, IFNW1), tolerance to cold-triggered pain (SCN11A) and antler development (SILT2). The majority of these characteristic reindeer genes have been reported for the first time here. Moreover, our population genomics analysis suggested at least two independent reindeer domestication events with genetic lineages originating from different refugial regions after the Last Glacial Maximum. Taken together, our study has provided new insights into the domestication, evolution and adaptation of reindeer and has promoted novel genomic research of reindeer.

Single-molecule long-read sequencing of the full-length transcriptome of Rhododendron lapponicum L

Rhododendron lapponicum L. is a familiar ornamental plant worldwide with important ornamental and economic value. However, a full-length R. lapponicum transcriptome is still lacking. In the present study, we used the Pacific Biosciences single-molecule real-time sequencing technology to generate the R. lapponicum transcriptome. A total of 346,270 full-length non-chimeric reads were generated, from which we obtained 75,002 high-quality full-length transcripts. We identified 55,255 complete open reading frames, 7,140 alternative splicing events and 2,011 long non-coding RNAs. In gene annotation analyses, 71,155, 33,653, 30,359 and 31,749 transcripts were assigned to the Nr, GO, COG and KEGG databases, respectively. Additionally, 3,150 transcription factors were detected. KEGG pathway analysis showed that 96 transcripts were identified coding for the enzymes associated with anthocyanin synthesis. Furthermore, we identified 64,327 simple sequence repeats from 45,319 sequences, and 150 pairs of primers were randomly selected to develop SSR markers. This study provides a large number of full-length transcripts, which will facilitate the further study of the genetics of R. lapponicum.

Molecular evidence for adaptive evolution of olfactory-related genes in cervids

BACKGROUND

Cervids have evolved very successful means for survival and thriving to adapt to various climates and environments. One of these successful means might be the effective and efficient way of communication. To support this notion, cervids are well equipped with a variety of skin glands that distribute in different body regions. However, studies relevant to adaptive evolution in cervids, particularly on olfactory reception at the molecular level, have thus far not been reported.

OBJECTIVE

To provide valuable insights into molecular evidence for the adaptive evolution of olfactory-related gene in cervids.

METHODS

Based on recently sequenced genomes of cervids and closely-related-species, we performed comparative genomic analysis at genome level using bioinformatics tools.

RESULTS

Tree topology strongly supported that Bovidae was the sister group of Moschidae and both formed a branch that was then clustered with Cervidae. Expansion of heavy chain genes of the dynein family and 51 rapidly evolving genes could be associated with adaptation of cilia that serve as sensory organelles and act as cellular antennae. Based on the branch-site model test along the deer branch spanning 7-21 mammalian species, 14 deer olfactory receptor genes were found to be undergoing positive selection pressure and 89 positive selection sites (probability > 60%) had amino acid substitutions unique to deer.

CONCLUSION

This study, for the first time, provides significant molecular evidence for adaption of olfactory-related genes of cervids according to their olfactory behavior.

GenomeScope 2.0 and Smudgeplot for reference-free profiling of polyploid genomes

An important assessment prior to genome assembly and related analyses is genome profiling, where the k-mer frequencies within raw sequencing reads are analyzed to estimate major genome characteristics such as size, heterozygosity, and repetitiveness. Here we introduce GenomeScope 2.0 (https://github.com/tbenavi1/genomescope2.0), which applies combinatorial theory to establish a detailed mathematical model of how k-mer frequencies are distributed in heterozygous and polyploid genomes. We describe and evaluate a practical implementation of the polyploid-aware mixture model that quickly and accurately infers genome properties across thousands of simulated and several real datasets spanning a broad range of complexity. We also present a method called Smudgeplot (https://github.com/KamilSJaron/smudgeplot) to visualize and estimate the ploidy and genome structure of a genome by analyzing heterozygous k-mer pairs. We successfully apply the approach to systems of known variable ploidy levels in the Meloidogyne genus and the extreme case of octoploid Fragaria × ananassa.

Full-length transcript characterization of SF3B1 mutation in chronic lymphocytic leukemia reveals downregulation of retained introns

While splicing changes caused by somatic mutations in SF3B1 are known, identifying full-length isoform changes may better elucidate the functional consequences of these mutations. We report nanopore sequencing of full-length cDNA from CLL samples with and without SF3B1 mutation, as well as normal B cell samples, giving a total of 149 million pass reads. We present FLAIR (Full-Length Alternative Isoform analysis of RNA), a computational workflow to identify high-confidence transcripts, perform differential splicing event analysis, and differential isoform analysis. Using nanopore reads, we demonstrate differential 3' splice site changes associated with SF3B1 mutation, agreeing with previous studies. We also observe a strong downregulation of intron retention events associated with SF3B1 mutation. Full-length transcript analysis links multiple alternative splicing events together and allows for better estimates of the abundance of productive versus unproductive isoforms. Our work demonstrates the potential utility of nanopore sequencing for cancer and splicing research.

The genome evolution and domestication of tropical fruit mango

Background

Mango is one of the world’s most important tropical fruits. It belongs to the family Anacardiaceae, which includes several other economically important species, notably cashew, sumac and pistachio from other genera. Many species in this family produce family-specific urushiols and related phenols, which can induce contact dermatitis.

Results

We generate a chromosome-scale genome assembly of mango, providing a reference genome for the Anacardiaceae family. Our results indicate the occurrence of a recent whole-genome duplication (WGD) event in mango. Duplicated genes preferentially retained include photosynthetic, photorespiration, and lipid metabolic genes that may have provided adaptive advantages to sharp historical decreases in atmospheric carbon dioxide and global temperatures. A notable example of an extended gene family is the chalcone synthase (CHS) family of genes, and particular genes in this family show universally higher expression in peels than in flesh, likely for the biosynthesis of urushiols and related phenols. Genome resequencing reveals two distinct groups of mango varieties, with commercial varieties clustered with India germplasms and demonstrating allelic admixture, and indigenous varieties from Southeast Asia in the second group. Landraces indigenous in China formed distinct clades, and some showed admixture in genomes.

Conclusions

Analysis of chromosome-scale mango genome sequences reveals photosynthesis and lipid metabolism are preferentially retained after a recent WGD event, and expansion of CHS genes is likely associated with urushiol biosynthesis in mango. Genome resequencing clarifies two groups of mango varieties, discovers allelic admixture in commercial varieties, and shows distinct genetic background of landraces.

Microbial community analysis using high-throughput sequencing technology: a beginner's guide for microbiologists

Microbial communities present in diverse environments from deep seas to human body niches play significant roles in the complex ecosystem and human health. Characterizing their structural and functional diversities is indispensable, and many approaches, such as microscopic observation, DNA fingerprinting, and PCR-based marker gene analysis, have been successfully applied to identify microorganisms. Since the revolutionary improvement of DNA sequencing technologies, direct and high-throughput analysis of genomic DNA from a whole environmental community without prior cultivation has become the mainstream approach, overcoming the constraints of the classical approaches. Here, we first briefly review the history of environmental DNA analysis applications with a focus on profiling the taxonomic composition and functional potentials of microbial communities. To this end, we aim to introduce the shotgun metagenomic sequencing (SMS) approach, which is used for the untargeted ("shotgun") sequencing of all ("meta") microbial genomes ("genomic") present in a sample. SMS data analyses are performed in silico using various software programs; however, in silico analysis is typically regarded as a burden on wet-lab experimental microbiologists. Therefore, in this review, we present microbiologists who are unfamiliar with in silico analyses with a basic and practical SMS data analysis protocol. This protocol covers all the bioinformatics processes of the SMS analysis in terms of data preprocessing, taxonomic profiling, functional annotation, and visualization.

The sequence and de novo assembly of the wild yak genome

Vulnerable populations of wild yak (Bos mutus), the wild ancestral species of domestic yak, survive in extremely cold, harsh and oxygen-poor regions of the Qinghai-Tibetan Plateau (QTP) and adjacent high-altitude regions. In this study, we sequenced and assembled its genome de novo. In total, six different insert-size libraries were sequenced, and 662 Gb of clean data were generated. The assembled wild yak genome is 2.83 Gb in length, with an N50 contig size of 63.2 kb and a scaffold size of 16.3 Mb. BUSCO assessment indicated that 93.8% of the highly conserved mammal genes were completely present in the genome assembly. Annotation of the wild yak genome assembly identified 1.41 Gb (49.65%) of repetitive sequences and a total of 22,910 protein-coding genes, including 20,660 (90.18%) annotated with functional terms. This first construction of the wild yak genome provides a variable genetic resource that will facilitate further study of the genetic diversity of bovine species and accelerate yak breeding efforts.