Chromosomal-scale genomes of two Rosa species provide insights into genome evolution and ascorbate accumulation.

Rosa roxburghii and Rosa sterilis, two species belonging to the Rosaceae family, are widespread in the southwest of China. These species have gained recognition for their remarkable abundance of ascorbate in their fresh fruits, making them an ideal vitamin C resource. In this study, we generated two high-quality chromosome-scale genome assemblies for R. roxburghii and R. sterilis, with genome sizes of 504 and 981.2 Mb, respectively. Notably, we present a haplotype-resolved, chromosome-scale assembly for diploid R. sterilis. Our results indicated that R. sterilis originated from the hybridization of R. roxburghii and R. longicuspis. Genome analysis revealed the absence of recent whole-genome duplications in both species and identified a series of duplicated genes that possibly contributing to the accumulation of flavonoids. We identified two genes in the ascorbate synthesis pathway, GGP and GalLDH, that show signs of positive selection, along with high expression levels of GDP-d-mannose 3', 5'-epimerase (GME) and GDP-l-galactose phosphorylase (GGP) during fruit development. Furthermore, through co-expression network analysis, we identified key hub genes (MYB5 and bZIP) that likely regulate genes in the ascorbate synthesis pathway, promoting ascorbate biosynthesis. Additionally, we observed the expansion of terpene synthase genes in these two species and tissue expression patterns, suggesting their involvement in terpenoid biosynthesis. Our research provides valuable insights into genome evolution and the molecular basis of the high concentration of ascorbate in these two Rosa species.

Spatiotemporal transcriptomic atlas of rhizome formation in Oryza longistaminata.

Rhizomes are modified stems that grow underground and produce new individuals genetically identical to the mother plant. Recently, a breakthrough has been made in efforts to convert annual grains into perennial ones by utilizing wild rhizomatous species as donors, yet the developmental biology of this organ is rarely studied. Oryza longistaminata, a wild rice species featuring strong rhizomes, provides a valuable model for exploration of rhizome development. Here, we first assembled a double-haplotype genome of O. longistaminata, which displays a 48-fold improvement in contiguity compared to the previously published assembly. Furthermore, spatiotemporal transcriptomics was performed to obtain the expression profiles of different tissues in O. longistaminata rhizomes and tillers. Two spatially reciprocal cell clusters, the vascular bundle 2 cluster and the parenchyma 2 cluster, were determined to be the primary distinctions between the rhizomes and tillers. We also captured meristem initiation cells in the sunken area of parenchyma located at the base of internodes, which is the starting point for rhizome initiation. Trajectory analysis further indicated that the rhizome is regenerated through de novo generation. Collectively, these analyses revealed a spatiotemporal transcriptional transition underlying the rhizome initiation, providing a valuable resource for future perennial crop breeding.

Genetic Architecture and Selection of Chinese Cattle Revealed by Whole Genome Resequencing.

The bovine genetic resources in China are diverse, but their value and potential are yet to be discovered. To determine the genetic diversity and population structure of Chinese cattle, we analyzed the whole genomes of 46 cattle from six phenotypically and geographically representative Chinese cattle breeds, together with 18 Red Angus cattle genomes, 11 Japanese black cattle genomes and taurine and indicine genomes available from previous studies. Our results showed that Chinese cattle originated from hybridization between Bos taurus and Bos indicus. Moreover, we found that the level of genetic variation in Chinese cattle depends upon the degree of indicine content. We also discovered many potential selective sweep regions associated with domestication related to breed-specific characteristics, with selective sweep regions including genes associated with coat color (ERCC2, MC1R, ZBTB17, and MAP2K1), dairy traits (NCAPG, MAPK7, FST, ITFG1, SETMAR, PAG1, CSN3, and RPL37A), and meat production/quality traits (such as BBS2, R3HDM1, IGFBP2, IGFBP5, MYH9, MYH4, and MC5R). These findings substantially expand the catalogue of genetic variants in cattle and reveal new insights into the evolutionary history and domestication traits of Chinese cattle.

IDSSR: An Efficient Pipeline for Identifying Polymorphic Microsatellites from a Single Genome Sequence.

Simple sequence repeats (SSRs) are known as microsatellites, and consist of tandem 1-6-base motifs. They have become one of the most popular molecular markers, and are widely used in molecular ecology, conservation biology, molecular breeding, and many other fields. Previously reported methods identify monomorphic and polymorphic SSRs and determine the polymorphic SSRs via experimental validation, which is potentially time-consuming and costly. Herein, we present a new strategy named insertion/deletion (INDEL) SSR (IDSSR) to identify polymorphic SSRs by integrating SSRs with nucleotide insertions/deletions (INDEL) solely based on a single genome sequence and the sequenced pair-end reads. These INDEL indexes and polymorphic SSRs were identified, as well as the number of repeats, repeat motifs, chromosome location, annealing temperature, and primer sequences, enabling future experimental approaches to determine the correctness and polymorphism. Experimental validation with the giant panda demonstrated that our method has high reliability and stability. The efficient SSR pipeline would help researchers obtain high-quality genetic markers for plants and animals of interest, save labor, and reduce costly marker-screening experiments. IDSSR is freely available at https://github.com/Allsummerking/IDSSR.

The draft genome of Tibetan hulless barley reveals adaptive patterns to the high stressful Tibetan Plateau.

The Tibetan hulless barley (Hordeum vulgare L. var. nudum), also called "Qingke" in Chinese and "Ne" in Tibetan, is the staple food for Tibetans and an important livestock feed in the Tibetan Plateau. The diploid nature and adaptation to diverse environments of the highland give it unique resources for genetic research and crop improvement. Here we produced a 3.89-Gb draft assembly of Tibetan hulless barley with 36,151 predicted protein-coding genes. Comparative analyses revealed the divergence times and synteny between barley and other representative Poaceae genomes. The expansion of the gene family related to stress responses was found in Tibetan hulless barley. Resequencing of 10 barley accessions uncovered high levels of genetic variation in Tibetan wild barley and genetic divergence between Tibetan and non-Tibetan barley genomes. Selective sweep analyses demonstrate adaptive correlations of genes under selection with extensive environmental variables. Our results not only construct a genomic framework for crop improvement but also provide evolutionary insights of highland adaptation of Tibetan hulless barley.

A host plant genome (Zizania latifolia) after a century-long endophyte infection.

Despite the importance of host-microbe interactions in natural ecosystems, agriculture and medicine, the impact of long-term (especially decades or longer) microbial colonization on the dynamics of host genomes is not well understood. The vegetable crop 'Jiaobai' with enlarged edible stems was domesticated from wild Zizania latifolia (Oryzeae) approximately 2000 years ago as a result of persistent infection by a fungal endophyte, Ustilago esculenta. Asexual propagation via infected rhizomes is the only means of Jiaobai production, and the Z. latifolia-endophyte complex has been maintained continuously for two centuries. Here, genomic analysis revealed that cultivated Z. latifolia has a significantly smaller repertoire of immune receptors compared with wild Z. latifolia. There are widespread gene losses/mutations and expression changes in the plant-pathogen interaction pathway in Jiaobai. These results show that continuous long-standing endophyte association can have a major effect on the evolution of the structural and transcriptomic components of the host genome.

The chromosome-level genomes of the herbal magnoliids Warburgia ugandensis and Saururus chinensis.

Warburgia ugandensis and Saururus chinensis are two of the most important medicinal plants in magnoliids and are widely utilized in traditional Kenya and Chinese medicine, respectively. The absence of higher-quality reference genomes has hindered research on the medicinal compound biosynthesis mechanisms of these plants. We report the chromosome-level genome assemblies of W. ugandensis and S. chinensis, and generated 1.13 Gb and 0.53 Gb genomes from 74 and 27 scaffolds, respectively, using BGI-DIPSEQ, Nanopore, and Hi-C sequencing. The scaffold N50 lengths were 82.97 Mb and 48.53 Mb, and the assemblies were anchored to 14 and 11 chromosomes of W. ugandensis and S. chinensis, respectively. In total, 24,739 and 20,561 genes were annotated, and 98.5% and 98% of the BUSCO genes were fully represented, respectively. The chromosome-level genomes of W. ugandensis and S. chinensis will be valuable resources for understanding the genetics of these medicinal plants, studying the evolution of magnoliids and angiosperms and conserving plant genetic resources.

Balanophora genomes display massively convergent evolution with other extreme holoparasites and provide novel insights into parasite-host interactions.

Parasitic plants have evolved to be subtly or severely dependent on host plants to complete their life cycle. To provide new insights into the biology of parasitic plants in general, we assembled genomes for members of the sandalwood order Santalales, including a stem hemiparasite (Scurrula) and two highly modified root holoparasites (Balanophora) that possess chimaeric host-parasite tubers. Comprehensive genome comparisons reveal that hemiparasitic Scurrula has experienced a relatively minor degree of gene loss compared with autotrophic plants, consistent with its moderate degree of parasitism. Nonetheless, patterns of gene loss appear to be substantially divergent across distantly related lineages of hemiparasites. In contrast, Balanophora has experienced substantial gene loss for the same sets of genes as an independently evolved holoparasite lineage, the endoparasitic Sapria (Malpighiales), and the two holoparasite lineages experienced convergent contraction of large gene families through loss of paralogues. This unprecedented convergence supports the idea that despite their extreme and strikingly divergent life histories and morphology, the evolution of these and other holoparasitic lineages can be shaped by highly predictable modes of genome reduction. We observe substantial evidence of relaxed selection in retained genes for both hemi- and holoparasitic species. Transcriptome data also document unusual and novel interactions between Balanophora and host plants at the host-parasite tuber interface tissues, with evidence of mRNA exchange, substantial and active hormone exchange and immune responses in parasite and host.

Mechanisms of sex differentiation and sex reversal in hermaphrodite fish as revealed by the Epinephelus coioides genome.

Most grouper species are functional protogynous hermaphrodites, but the genetic basis and the molecular mechanisms underlying the regulation of this unique reproductive strategy remain enigmatic. In this study, we report a high-quality chromosome-level genome assembly of the representative orange-spotted grouper (Epinephelus coioides). No duplication or deletion of sex differentiation-related genes was found in the genome, suggesting that sex development in this grouper may be related to changes in regulatory sequences or environmental factors. Transcriptomic analyses showed that aromatase and retinoic acid are probably critical to promoting ovarian fate determination, and follicle-stimulating hormone triggers the female-to-male sex change. Socially controlled sex-change studies revealed that, in sex-changing fish, the brain's response to the social environment may be mediated by activation of the phototransduction cascade and the melatonin synthesis pathway. In summary, our genomic and experimental results provide novel insights into the molecular mechanisms of sex differentiation and sex change in the protogynous groupers.

Adaptations to a cold climate promoted social evolution in Asian colobine primates.

The biological mechanisms that underpin primate social evolution remain poorly understood. Asian colobines display a range of social organizations, which makes them good models for investigating social evolution. By integrating ecological, geological, fossil, behavioral, and genomic analyses, we found that colobine primates that inhabit colder environments tend to live in larger, more complex groups. Specifically, glacial periods during the past 6 million years promoted the selection of genes involved in cold-related energy metabolism and neurohormonal regulation. More-efficient dopamine and oxytocin pathways developed in odd-nosed monkeys, which may have favored the prolongation of maternal care and lactation, increasing infant survival in cold environments. These adaptive changes appear to have strengthened interindividual affiliation, increased male-male tolerance, and facilitated the stepwise aggregation from independent one-male groups to large multilevel societies.

Chloranthus genome provides insights into the early diversification of angiosperms.

Chloranthales remain the last major mesangiosperm lineage without a nuclear genome assembly. We therefore assemble a high-quality chromosome-level genome of Chloranthus spicatus to resolve enigmatic evolutionary relationships, as well as explore patterns of genome evolution among the major lineages of mesangiosperms (eudicots, monocots, magnoliids, Chloranthales, and Ceratophyllales). We find that synteny is highly conserved between genomic regions of Amborella, Vitis, and Chloranthus. We identify an ancient single whole-genome duplication (WGD) (κ) prior to the divergence of extant Chloranthales. Phylogenetic inference shows Chloranthales as sister to magnoliids. Furthermore, our analyses indicate that ancient hybridization may account for the incongruent phylogenetic placement of Chloranthales + magnoliids relative to monocots and eudicots in nuclear and chloroplast trees. Long genes and long introns are found to be prevalent in both Chloranthales and magnoliids compared to other angiosperms. Overall, our findings provide an improved context for understanding mesangiosperm relationships and evolution and contribute a valuable genomic resource for future investigations.

Chromosome-scale genomes provide new insights into subspecies divergence and evolutionary characteristics of the giant panda.

Extant giant pandas are divided into Sichuan and Qinling subspecies. The giant panda has many species-specific characteristics, including comparatively small organs for body size, small genitalia of male individuals, and low reproduction. Here, we report the most contiguous, high-quality chromosome-level genomes of two extant giant panda subspecies to date, with the first genome assembly of the Qinling subspecies. Compared with the previously assembled giant panda genomes based on short reads, our two assembled genomes increased contiguity over 200-fold at the contig level. Additional sequencing of 25 individuals dated the divergence of the Sichuan and Qinling subspecies into two distinct clusters from 10,000 to 12,000 years ago. Comparative genomic analyses identified the loss of regulatory elements in the dachshund family transcription factor 2 (DACH2) gene and specific changes in the synaptotagmin 6 (SYT6) gene, which may be responsible for the reduced fertility of the giant panda. Positive selection analysis between the two subspecies indicated that the reproduction-associated IQ motif containing D (IQCD) gene may at least partly explain the different reproduction rates of the two subspecies. Furthermore, several genes in the Hippo pathway exhibited signs of rapid evolution with giant panda-specific variants and divergent regulatory elements, which may contribute to the reduced inner organ sizes of the giant panda.

Predicting the state of cysteines based on sequence information.

A three-stage support vector machine (SVM) was constructed to predict the state of cysteines by fusing sequence information, evolution information and annotation information of protein sequences. The first and second stages were for predicting whether the protein sequences contain disulfide bonds and whether all of the cysteines are involved in disulfide bonds. In the last stage, one SVM was constructed for predicting which cysteines are involved in disulfide bonds, among all these cysteines in proteins. The three SVMs give a good performance and the overall prediction accuracy are 90.05%, 96.36% and 80.00%, respectively, which indicates that the features selected in this work are effective for predicting the state of cysteines. In addition, current methods only paid too much attention to the prediction performance and never showed us how much important the roles of these features played in the prediction. As a result a feature importance measurement designated as F-score function was used to evaluate these features. The result shows that among these protein descriptors; evolution information is the most important feature for representing the disulfide-containing proteins. The prediction software and data sets used in this article are freely available at http://cic.scu.edu.cn/bioinformatics/Predict_Cys.zip.

The draft genome of mandrill (Mandrillus sphinx): An Old World monkey.

Mandrill (Mandrillus sphinx) is a primate species, which belongs to the Old World monkey (Cercopithecidae) family. It is closely related to human, serving as a model for human health related research. However, the genetic studies on and genomic resources of mandrill are limited, especially in comparison to other primate species. Here we produced 284 Gb data, providing 96-fold coverage (considering the estimated genome size of 2.9 Gb), to construct a reference genome for the mandrill. The assembled draft genome was 2.79 Gb with contig N50 of 20.48 Kb and scaffold N50 of 3.56 Mb. We annotated the mandrill genome to find 43.83% repeat elements, as well as 21,906 protein-coding genes. The draft genome was of good quality with 98% gene annotation coverage by Benchmarking Universal Single-Copy Orthologs (BUSCO). Based on comparative genomic analyses of  the Major Histocompatibility Complex (MHC) of the immune system in mandrill and human, we found that 17 genes in the mandrill that have been associated with disease phenotypes in human such as Lung cancer, cranial volume and asthma, barbored amino acids changing mutations. Gene family analyses revealed expansion of several genes, and several genes associated with stress environmental adaptation and innate immunity responses exhibited signatures of positive selection. In summary, we established the first draft genome of  the mandrill of value for studies on evolution and human health.

Genomic analysis of the domestication and post-Spanish conquest evolution of the llama and alpaca.

BACKGROUND: Despite their regional economic importance and being increasingly reared globally, the origins and evolution of the llama and alpaca remain poorly understood. Here we report reference genomes for the llama, and for the guanaco and vicuña (their putative wild progenitors), compare these with the published alpaca genome, and resequence seven individuals of all four species to better understand domestication and introgression between the llama and alpaca. RESULTS: Phylogenomic analysis confirms that the llama was domesticated from the guanaco and the alpaca from the vicuña. Introgression was much higher in the alpaca genome (36%) than the llama (5%) and could be dated close to the time of the Spanish conquest, approximately 500 years ago. Introgression patterns are at their most variable on the X-chromosome of the alpaca, featuring 53 genes known to have deleterious X-linked phenotypes in humans. Strong genome-wide introgression signatures include olfactory receptor complexes into both species, hypertension resistance into alpaca, and fleece/fiber traits into llama. Genomic signatures of domestication in the llama include male reproductive traits, while in alpaca feature fleece characteristics, olfaction-related and hypoxia adaptation traits. Expression analysis of the introgressed region that is syntenic to human HSA4q21, a gene cluster previously associated with hypertension in humans under hypoxic conditions, shows a previously undocumented role for PRDM8 downregulation as a potential transcriptional regulation mechanism, analogous to that previously reported at high altitude for hypoxia-inducible factor 1α. CONCLUSIONS: The unprecedented introgression signatures within both domestic camelid genomes may reflect post-conquest changes in agriculture and the breakdown of traditional management practices.

In silico method for systematic analysis of feature importance in microRNA-mRNA interactions.

BACKGROUND: MicroRNA (miRNA), which is short non-coding RNA, plays a pivotal role in the regulation of many biological processes and affects the stability and/or translation of mRNA. Recently, machine learning algorithms were developed to predict potential miRNA targets. Most of these methods are robust but are not sensitive to redundant or irrelevant features. Despite their good performance, the relative importance of each feature is still unclear. With increasing experimental data becoming available, research interest has shifted from higher prediction performance to uncovering the mechanism of microRNA-mRNA interactions. RESULTS: Systematic analysis of sequence, structural and positional features was carried out for two different data sets. The dominant functional features were distinguished from uninformative features in single and hybrid feature sets. Models were developed using only statistically significant sequence, structural and positional features, resulting in area under the receiver operating curves (AUC) values of 0.919, 0.927 and 0.969 for one data set and of 0.926, 0.874 and 0.954 for another data set, respectively. Hybrid models were developed by combining various features and achieved AUC of 0.978 and 0.970 for two different data sets. Functional miRNA information is well reflected in these features, which are expected to be valuable in understanding the mechanism of microRNA-mRNA interactions and in designing experiments. CONCLUSIONS: Differing from previous approaches, this study focused on systematic analysis of all types of features. Statistically significant features were identified and used to construct models that yield similar accuracy to previous studies in a shorter computation time.

Genome Resequencing Reveals Congenital Causes of Embryo and Nestling Death in Crested Ibis (Nipponia nippon).

The crested ibis (Nipponia nippon) is endangered worldwide. Although a series of conservation measures have markedly increased the population size and distribution area of these birds, the high mortality of embryos and nestlings considerably decreases the survival potential of this bird species. High-throughput sequencing technology was utilized to compare whole genomes between ten samples from dead crested ibises (including six dead embryos and four dead nestlings aged 0-45 days) and 32 samples from living birds. The results indicated that the dead samples all shared the genetic background of a specific ancestral subpopulation. Furthermore, the dead individuals were less genetically diverse and suffered higher degrees of inbreeding compared with these measures in live birds. Several candidate genes (KLHL3, SETDB2, TNNT2, PKP1, AK1, and EXOSC3) associated with detrimental diseases were identified in the genomic regions that differed between the alive and dead samples, which are likely responsible for the death of embryos and nestlings. In addition, in these regions, we also found several genes involved in the protein catabolic process (UBE4A and LONP1), lipid metabolism (ACOT1), glycan biosynthesis and metabolism (HYAL1 and HYAL4), and the immune system (JAM2) that are likely to promote the normal development of embryos and nestlings. The aberrant conditions of these genes and biological processes may contribute to the death of embryos and nestlings. Our data identify congenital factors underlying the death of embryos and nestlings at the whole genome level, which may be useful toward informing more effective conservation efforts for this bird species.

Genome-wide Analysis of the MADS-Box Gene Family in Watermelon.

MADS-box genes comprise a family of transcription factors that function in the growth and development of plants. To obtain insights into their evolution in watermelon (Citrullus lanatus), we carried out a genome-wide analysis and identified 39 MADS-box genes. These genes were classified into MIKCc (25), MIKC*(3), Mα (5), Mß (3), and Mγ (3) clades according to their phylogenetic relationship with Arabidopsis thaliana and Cucumis sativus; moreover, these 25 genes in the MIKC clade could be classified into 13 subfamilies, and the Flowering Locus C (FLC) subfamily is absent in watermelon. Analysis of the conserved gene motifs showed similar motifs among clades. Continuing chromosomal localizations analysis indicated that MADS-box genes were distributed across 11 chromosomes in watermelon, and these genes were conditioned to be differentially expressed during plant growth and development. This research provides information that will aid further investigations into the evolution of the MADS-box gene family in plants.

Draft genome of the river water buffalo.

Water buffalo (Bubalus bubalis), a large-sized member of the Bovidae family, is considered as an important livestock species throughout Southeast Asia. In order to better understand the molecular basis of buffalo improvement and breeding, we sequenced and assembled the genome (2n=50) of a river buffalo species Bubalus bubalis from Bangladesh. Its genome size is 2.77 Gb, with a contig N50 of 25 kb and the scaffold N50 of 6.9 Mbp. Based on the assembled genome, we annotated 24,613 genes for future functional genomics studies. Phylogenetic tree analysis of cattle and water buffalo lineages showed that they diverged about 5.8-9.8 million years ago. Our findings provide an insight into the water buffalo genome which will contribute in further research on buffalo such as molecular breeding, understanding complex traits, conservation, and biodiversity.

Author Correction: Liriodendron genome sheds light on angiosperm phylogeny and species-pair differentiation.

In the Supplementary Information file originally published with this Letter, the authors mistakenly omitted Supplementary Table 14; this has now been amended.