Comparative performance of the BGISEQ-500 vs Illumina HiSeq2500 sequencing platforms for palaeogenomic sequencing

MTP, Carneiro M. Molecular parallelisms between pigmentation in the avian iris and the integument of ectothermic vertebrates

Birds exhibit striking variation in eye color that arises from interactions between specialized pigment cells named chromatophores. The types of chromatophores present in the avian iris are lacking from the integument of birds or mammals, but are remarkably similar to those found in the skin of ectothermic vertebrates. To investigate molecular mechanisms associated with eye coloration in birds, we took advantage of a Mendelian mutation found in domestic pigeons that alters the deposition of yellow pterin pigments in the iris. Using a combination of genome-wide association analysis and linkage information in pedigrees, we mapped variation in eye coloration in pigeons to a small genomic region of ~8.5kb. This interval contained a single gene, SLC2A11B, which has been previously implicated in skin pigmentation and chromatophore differentiation in fish. Loss of yellow pigmentation is likely caused by a point mutation that introduces a premature STOP codon and leads to lower expression of SLC2A11B through nonsense-mediated mRNA decay. There were no substantial changes in overall gene expression profiles between both iris types as well as in genes directly associated with pterin metabolism and/or chromatophore differentiation. Our findings demonstrate that SLC2A11B is required for the expression of pterin-based pigmentation in the avian iris. They further highlight common molecular mechanisms underlying the production of coloration in the iris of birds and skin of ectothermic vertebrates.

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.

Optimizing culture methods according to preoperative mNGS results can improve joint infection diagnosis

AIMS

Metagenomic next-generation sequencing (mNGS) is useful in the diagnosis of infectious disease. However, while it is highly sensitive at identifying bacteria, it does not provide information on the sensitivity of the organisms to antibiotics. The purpose of this study was to determine whether the results of mNGS can be used to guide optimization of culture methods to improve the sensitivity of culture from intraoperative samples.

METHODS

Between July 2014 and October 2019, patients with suspected joint infection (JI) from whom synovial fluid (SF) was obtained preoperatively were enrolled. Preoperative aspirated SF was analyzed by conventional microbial culture and mNGS. In addition to samples taken for conventional microbial culture, some samples were taken for intraoperative culture to optimize the culture method according to the preoperative mNGS results. The demographic characteristics, medical history, laboratory examination, mNGS, and culture results of the patients were recorded, and the possibility of the optimized culture methods improving diagnostic efficiency was evaluated.

RESULTS

A total of 56 cases were included in this study. There were 35 cases of JI and 21 cases of non-joint infection (NJI). The sensitivity, specificity, and accuracy of intraoperative microbial culture after optimization of the culture method were 94.29%, 76.19%, and 87.5%, respectively, while those of the conventional microbial culture method were 60%, 80.95%, and 67.86%, respectively.

CONCLUSION

Preoperative aspirated SF detected via mNGS can provide more aetiological information than preoperative culture, which can guide the optimization and improve the sensitivity of intraoperative culture. Cite this article: Bone Joint J 2021;103-B(1):39-45.

Genome diversity in Ukraine

Background

The main goal of this collaborative effort is to provide genome-wide data for the previously underrepresented population in Eastern Europe, and to provide cross-validation of the data from genome sequences and genotypes of the same individuals acquired by different technologies. We collected 97 genome-grade DNA samples from consented individuals representing major regions of Ukraine that were consented for public data release. BGISEQ-500 sequence data and genotypes by an Illumina GWAS chip were cross-validated on multiple samples and additionally referenced to 1 sample that has been resequenced by Illumina NovaSeq6000 S4 at high coverage.

Results

The genome data have been searched for genomic variation represented in this population, and a number of variants have been reported: large structural variants, indels, copy number variations, single-nucletide polymorphisms, and microsatellites. To our knowledge, this study provides the largest to-date survey of genetic variation in Ukraine, creating a public reference resource aiming to provide data for medical research in a large understudied population.

Conclusions

Our results indicate that the genetic diversity of the Ukrainian population is uniquely shaped by evolutionary and demographic forces and cannot be ignored in future genetic and biomedical studies. These data will contribute a wealth of new information bringing forth a wealth of novel, endemic and medically related alleles.

Genomes of Pleistocene Siberian Wolves Uncover Multiple Extinct Wolf Lineages

Extant Canis lupus genetic diversity can be grouped into three phylogenetically distinct clades: Eurasian and American wolves and domestic dogs.1 Genetic studies have suggested these groups trace their origins to a wolf population that expanded during the last glacial maximum (LGM)1-3 and replaced local wolf populations.4 Moreover, ancient genomes from the Yana basin and the Taimyr peninsula provided evidence of at least one extinct wolf lineage that dwelled in Siberia during the Pleistocene.35 Previous studies have suggested that Pleistocene Siberian canids can be classified into two groups based on cranial morphology. Wolves in the first group are most similar to present-day populations, although those in the second group possess intermediate features between dogs and wolves.67 However, whether this morphological classification represents distinct genetic groups remains unknown. To investigate this question and the relationships between Pleistocene canids, present-day wolves, and dogs, we resequenced the genomes of four Pleistocene canids from Northeast Siberia dated between >50 and 14 ka old, including samples from the two morphological categories. We found these specimens cluster with the two previously sequenced Pleistocene wolves, which are genetically more similar to Eurasian wolves. Our results show that, though the four specimens represent extinct wolf lineages, they do not form a monophyletic group. Instead, each Pleistocene Siberian canid branched off the lineage that gave rise to present-day wolves and dogs. Finally, our results suggest the two previously described morphological groups could represent independent lineages similarly related to present-day wolves and dogs.

The complete chloroplast genome of Cymbidium hookerianum (Orchidaceae): genome structure and basic analysis

Cymbidium hookerianum Rchb.f. is an ornamental orchid with large flowers and delicate aroma. Here, we reported the complete chloroplast genome sequence of C. hookerianum. The total chloroplast genome cycle was 155,447 bp. It displayed a typical structure including one large single-copy (LSC, 84,186 bp) region, one small single-copy (SSC, 17,839 bp) region, and two inverted repeat (IRs, 26,711 bp). 124 genes (78 CDSs, 38 tRNAs, 8 rRNAs) were encoded by the cp genome. The average GC content of this sequence was 36.8%. The phylogenetic analysis revealed C. hookerianum and C. changningense are sisters. The groundwork of chloroplast genome would provide available reference for molecular taxonomy and breeding.

An integrated morphological–molecular approach reveals new insights on the systematics of the octocoral Telestula humilis (Thomson, 1927) (Octocorallia:Alcyonacea:Clavulariidae)

Telestula humilis (Thomson, 1927) is a rare deep-sea stoloniferan octocoral distributed in the eastern Atlantic. Here we compared seven putative colonies of this species collected off Spain with the lectotype from the Oceanographic Museum of Monaco and found them to be identical morphologically. Phylogenetic analyses on both full mitogenomes and a concatenated alignment containing two mtDNA genes (mtMutS and Cox1) and nuclear 28S rRNA gene recovered Telestula humilis sister to Incrustatus and Inconstantia rather than to other species of Telestula. This therefore supports its taxonomic reassignment to Pseudotelestula gen. nov. as Pseudotelestula humilis comb. nov. The taxonomic reassignment is also supported by subtle differences observed between the morphology of the colony and the sclerome of Pseudotelestula humilis comb. nov. and the two sister genera. The occurrence of an intrusion tissue with sclerites in the basal part of the gastric cavity of the adult polyps is shared among Telestula and Pseudotelestula gen. nov. However, Pseudotelestula gen. nov. has sclerites arranged in a collaret and points below the tentacles, the sclerites of the calyx wall and the stolon are plump warty spindles, and the intrusion tissue has long sticks and spindles with cone-like spines.

Archaeological Central American maize genomes suggest ancient gene flow from South America

Maize is a global food staple with great economic and cultural importance. Archaeogenomic studies have revealed a process of protracted maize domestication and multiple waves of human-mediated dispersal in the Americas. Maize first arrived in South America as a partial domesticate, where the domestication syndrome became independently fixed and improved varieties developed away from the influence of wild gene flow. We demonstrate that hybrids of some of these improved varieties were likely reintroduced back to Central America. We hypothesize that this backflow of South American genetic material may have contributed to the development of a more productive staple, which was related to the growth and aggregation of human populations, and the formation of more complex social and political structures regionally.

Maize (Zea mays ssp. mays) domestication began in southwestern Mexico ∼9,000 calendar years before present (cal. BP) and humans dispersed this important grain to South America by at least 7,000 cal. BP as a partial domesticate. South America served as a secondary improvement center where the domestication syndrome became fixed and new lineages emerged in parallel with similar processes in Mesoamerica. Later, Indigenous cultivators carried a second major wave of maize southward from Mesoamerica, but it has been unclear until now whether the deeply divergent maize lineages underwent any subsequent gene flow between these regions. Here we report ancient maize genomes (2,300–1,900 cal. BP) from El Gigante rock shelter, Honduras, that are closely related to ancient and modern maize from South America. Our findings suggest that the second wave of maize brought into South America hybridized with long-established landraces from the first wave, and that some of the resulting newly admixed lineages were then reintroduced to Central America. Direct radiocarbon dates and cob morphological data from the rock shelter suggest that more productive maize varieties developed between 4,300 and 2,500 cal. BP. We hypothesize that the influx of maize from South America into Central America may have been an important source of genetic diversity as maize was becoming a staple grain in Central and Mesoamerica.

Efficient COI barcoding using high throughput single-end 400 bp sequencing

BACKGROUND

Over the last decade, the rapid development of high-throughput sequencing platforms has accelerated species description and assisted morphological classification through DNA barcoding. However, the current high-throughput DNA barcoding methods cannot obtain full-length barcode sequences due to read length limitations (e.g. a maximum read length of 300 bp for the Illumina's MiSeq system), or are hindered by a relatively high cost or low sequencing output (e.g. a maximum number of eight million reads per cell for the PacBio's SEQUEL II system).

RESULTS

Pooled cytochrome c oxidase subunit I (COI) barcodes from individual specimens were sequenced on the MGISEQ-2000 platform using the single-end 400 bp (SE400) module. We present a bioinformatic pipeline, HIFI-SE, that takes reads generated from the 5' and 3' ends of the COI barcode region and assembles them into full-length barcodes. HIFI-SE is written in Python and includes four function modules of filter, assign, assembly and taxonomy. We applied the HIFI-SE to a set of 845 samples (30 marine invertebrates, 815 insects) and delivered a total of 747 fully assembled COI barcodes as well as 70 Wolbachia and fungi symbionts. Compared to their corresponding Sanger sequences (72 sequences available), nearly all samples (71/72) were correctly and accurately assembled, including 46 samples that had a similarity score of 100% and 25 of ca. 99%.

CONCLUSIONS

The HIFI-SE pipeline represents an efficient way to produce standard full-length barcodes, while the reasonable cost and high sensitivity of our method can contribute considerably more DNA barcodes under the same budget. Our method thereby advances DNA-based species identification from diverse ecosystems and increases the number of relevant applications.

A Chromosome-Level Genome Assembly of the Anglerfish Lophius litulon

Anglerfishes are a highly diverse group of species with unique characteristics. Here, we report the first chromosome-level genome of a species in the order Lophiiformes, the yellow goosefish (Lophius litulon), obtained by whole genome shotgun sequencing and high-throughput chromatin conformation capture. Approximately 97.20% of the assembly spanning 709.23 Mb could be anchored to 23 chromosomes with a contig N50 of 164.91 kb. The BUSCO value was 95.4%, suggesting that the quality of the assembly was high. A comparative gene family analysis identified expanded and contracted gene families, and these may be associated with adaptation to the benthic environment and the lack of scales in the species. A majority of positively selected genes were related to metabolic processes, suggesting that digestive and metabolic system evolution expanded the diversity of yellow goosefish prey. Our study provides a valuable genetic resource for understanding the mechanisms underlying the unique features of the yellow goosefish and for investigating anglerfish evolution.

Cytogenetic and molecular genetic methods for chromosomal translocations detection with reference to the KMT2A/MLL gene

Acute leukemias (ALs) are often associated with chromosomal translocations, in particular, KMT2A/MLL gene rearrangements. Identification or confirmation of these translocations is carried out by a number of genetic and molecular methods, some of which are routinely used in clinical practice, while others are primarily used for research purposes. In the clinic, these methods serve to clarify diagnoses and monitor the course of disease and therapy. On the other hand, the identification of new translocations and the confirmation of known translocations are of key importance in the study of disease mechanisms and further molecular classification. There are multiple methods for the detection of rearrangements that differ in their principle of operation, the type of problem being solved, and the cost-result ratio. This review is intended to help researchers and clinicians studying AL and related chromosomal translocations to navigate this variety of methods. All methods considered in the review are grouped by their principle of action and include karyotyping, fluorescence in situ hybridization (FISH) with probes for whole chromosomes or individual loci, PCR and reverse transcription-based methods, and high-throughput sequencing. Another characteristic of the described methods is the type of problem being solved. This can be the discovery of new rearrangements, the determination of unknown partner genes participating in the rearrangement, or the confirmation of the proposed rearrangement between the two genes. We consider the specifics of the application, the basic principle of each method, and its pros and cons. To illustrate the application, examples of studying the rearrangements of the KMT2A/MLL gene, one of the genes that are often rearranged in AL, are mentioned.

Chloroplast characterizations of a Phalaenopsis native to China, Phalaenopsis mannii (Orchidaceae)

Phalaenopsis mannii, one of the native Phalaenopsis in China, is an important parent for breeding new varieties. However, its position has been unclear in Phalaenopsis. The obtained high-quality P. mannii chloroplast genome will provide useful information for phylogenetic and future breeding of Phalaenopsis. Herein, we reported a complete chloroplast genome of P. mannii from Yunnan, China. The sequencing data obtained from BGISEQ-500 platform were assembled. This sequence had a circular molecular length of 148,596 bp and contained a total of 127 genes with an average GC content of 36.7%. Phylogenetic analysis indicated that Phalaenopsis was monophyletic with strong support, in which the P.mannii was the sister-group of Phalaenopsis aphrodite subsp. formosas, Phalaenopsis ‘TinyStar’ and Phalaenopsis equestris.

Complete Genome Sequence of a Jumbo Bacteriophage, vB_pir03, against Vibrio harveyi

Vibrio harveyi is a persistent pathogen responsible for disease outbreaks in aquaculture. We have sequenced the genome of a jumbo Vibrio phage, vB_pir03, isolated in Greece. Here, we present the complete genome of vB_pir03, which consists of 286,284 bp and 336 open reading frames.

Advantages and Limitations of 16S rRNA Next-Generation Sequencing for Pathogen Identification in the Diagnostic Microbiology Laboratory: Perspectives from a Middle-Income Country

Bacterial culture and biochemical testing (CBtest) have been the cornerstone of pathogen identification in the diagnostic microbiology laboratory. With the advent of Sanger sequencing and later, next-generation sequencing, 16S rRNA next-generation sequencing (16SNGS) has been proposed to be a plausible platform for this purpose. Nevertheless, usage of the 16SNGS platform has both advantages and limitations. In addition, transition from the traditional methods of CBtest to 16SNGS requires procurement of costly equipment, timely and sustainable maintenance of these platforms, specific facility infrastructure and technical expertise. All these factors pose a challenge for middle-income countries, more so for countries in the lower middle-income range. In this review, we describe the basis for CBtest and 16SNGS, and discuss the limitations, challenges, advantages and future potential of using 16SNGS for bacterial pathogen identification in diagnostic microbiology laboratories of middle-income countries.

Multi-omic detection of Mycobacterium leprae in archaeological human dental calculus

Mineralized dental plaque (calculus) has proven to be an excellent source of ancient biomolecules. Here we present a Mycobacterium leprae genome (6.6-fold), the causative agent of leprosy, recovered via shotgun sequencing of sixteenth-century human dental calculus from an individual from Trondheim, Norway. When phylogenetically placed, this genome falls in branch 3I among the diversity of other contemporary ancient strains from Northern Europe. Moreover, ancient mycobacterial peptides were retrieved via mass spectrometry-based proteomics, further validating the presence of the pathogen. Mycobacterium leprae can readily be detected in the oral cavity and associated mucosal membranes, which likely contributed to it being incorporated into this individual's dental calculus. This individual showed some possible, but not definitive, evidence of skeletal lesions associated with early-stage leprosy. This study is the first known example of successful multi-omics retrieval of M. leprae from archaeological dental calculus. Furthermore, we offer new insights into dental calculus as an alternative sample source to bones or teeth for detecting and molecularly characterizing M. leprae in individuals from the archaeological record.

This article is part of the theme issue ‘Insights into health and disease from ancient biomolecules’.

The Application of Next-Generation Sequencing to Define Factors Related to Oral Cancer and Discover Novel Biomarkers

Despite the introduction of next-generation sequencing in the realm of DNA sequencing technology, it is not often used in the investigation of oral squamous cell carcinoma (OSCC). Oral cancer is one of the most frequently occurring malignancies in some parts of the world and has a high mortality rate. Patients with this malignancy are likely to have a poor prognosis and may suffer from severe facial deformity or mastication problems even after successful treatment. Therefore, a thorough understanding of this malignancy is essential to prevent and treat it. This review sought to highlight the contributions of next-generation sequencing (NGS) in unveiling the genetic alterations and differential expressions of miRNAs involved in OSCC progression. By applying an appropriate eligibility criterion, we selected relevant studies for review. Frequently identified mutations in genes such as TP53, NOTCH1, and PIK3CA are discussed. The findings of existing miRNAs (e.g., miR-21) as well as novel discoveries pertaining to OSCC are also covered. Lastly, we briefly mention the latest findings in targeted gene therapy and the potential use of miRNAs as biomarkers. Our goal is to encourage researchers to further adopt NGS in their studies and give an overview of the latest findings of OSCC treatment.

Molecular identification and phylogenetic analysis of a Callosciurus notatus complete mitogenome from Peninsular Malaysia

The mitogenome of a plantain squirrel, Callosciurus notatus, collected from Bukit Tarek Forest Reserve (Extension), Selangor, Malaysia was sequenced using BGISEQ-500RS technology. The 16,582 bp mitogenome consists of 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and 1 control region. A phylogenetic and BLASTn analysis against other available datasets showed that the mitogenome matched with 99.49% similarity to a previously published C. notatus mitogenome from Peninsular Malaysia. However, it also diverged by nearly 8% (92.24% match) from a second previously published mitogenome for the same species, sampled in East Kalimantan, Indonesia. This suggests a difference in landscape features between both localities might affect its genetic connectivity.

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.

Marked mitochondrial genetic variation in individuals and populations of the carcinogenic liver fluke Clonorchis sinensis

Clonorchiasis is a neglected tropical disease caused by the Chinese liver fluke, Clonorchis sinensis, and is often associated with a malignant form of bile duct cancer (cholangiocarcinoma). Although some aspects of the epidemiology of clonorchiasis are understood, little is known about the genetics of C. sinensis populations. Here, we conducted a comprehensive genetic exploration of C. sinensis from endemic geographic regions using complete mitochondrial protein gene sets. Genomic DNA samples from C. sinensis individuals (n = 183) collected from cats and dogs in China (provinces of Guangdong, Guangxi, Hunan, Heilongjiang and Jilin) as well as from rats infected with metacercariae from cyprinid fish from the Russian Far East (Primorsky Krai region) were deep sequenced using the BGISEQ-500 platform. Informatic analyses of mitochondrial protein gene data sets revealed marked genetic variation within C. sinensis; significant variation was identified within and among individual worms from distinct geographical locations. No clear affiliation with a particular location or host species was evident, suggesting a high rate of dispersal of the parasite across endemic regions. The present work provides a foundation for future biological, epidemiological and ecological studies using mitochondrial protein gene data sets, which could aid in elucidating associations between particular C. sinensis genotypes/haplotypes and the pathogenesis or severity of clonorchiasis and its complications (including cholangiocarcinoma) in humans.

Clonorchiasis is an important neglected tropical disease caused by the Chinese liver fluke, Clonorchis sinensis, which can induce malignant bile duct cancer (cholangiocarcinoma). Little precise information is available on the biology, epidemiology and population genetics of C. sinensis. For this reason, we explored here the genetic composition of C. sinensis populations in distinct endemic areas in China and Russia. Using a deep sequencing-informatic approach, we revealed marked mitochondrial genetic variation within and between individuals and populations of C. sinensis, with no particular affiliation with geographic or host origin. These molecular findings and the methodology established should underpin future genetic studies of C. sinensis causing human clonorchiasis and associated complications (cancer) as well as transmission patterns in endemic regions.

Pre-extinction Demographic Stability and Genomic Signatures of Adaptation in the Woolly Rhinoceros

Ancient DNA has significantly improved our understanding of the evolution and population history of extinct megafauna. However, few studies have used complete ancient genomes to examine species responses to climate change prior to extinction. The woolly rhinoceros (Coelodonta antiquitatis) was a cold-adapted megaherbivore widely distributed across northern Eurasia during the Late Pleistocene and became extinct approximately 14 thousand years before present (ka BP). While humans and climate change have been proposed as potential causes of extinction [1-3], knowledge is limited on how the woolly rhinoceros was impacted by human arrival and climatic fluctuations [2]. Here, we use one complete nuclear genome and 14 mitogenomes to investigate the demographic history of woolly rhinoceros leading up to its extinction. Unlike other northern megafauna, the effective population size of woolly rhinoceros likely increased at 29.7 ka BP and subsequently remained stable until close to the species' extinction. Analysis of the nuclear genome from a ∼18.5-ka-old specimen did not indicate any increased inbreeding or reduced genetic diversity, suggesting that the population size remained steady for more than 13 ka following the arrival of humans [4]. The population contraction leading to extinction of the woolly rhinoceros may have thus been sudden and mostly driven by rapid warming in the Bølling-Allerød interstadial. Furthermore, we identify woolly rhinoceros-specific adaptations to arctic climate, similar to those of the woolly mammoth. This study highlights how species respond differently to climatic fluctuations and further illustrates the potential of palaeogenomics to study the evolutionary history of extinct species.

Tagsteady: A metabarcoding library preparation protocol to avoid false assignment of sequences to samples

Metabarcoding of environmental DNA (eDNA) and DNA extracted from bulk specimen samples is a powerful tool in studies of biodiversity, diet and ecological interactions as its inherent labelling of amplicons allows sequencing of taxonomically informative genetic markers from many samples in parallel. However, the occurrence of so-called 'tag-jumps' can cause incorrect assignment of sequences to samples and artificially inflate diversity. Two steps during library preparation of pools of 5' nucleotide-tagged amplicons have been suggested to cause tag-jumps: (a) T4 DNA polymerase blunt-ending in the end-repair step and (b) postligation PCR amplification of amplicon libraries. The discovery of tag-jumps has led to recommendations to only carry out metabarcoding PCR amplifications with primers carrying twin-tags to ensure that tag-jumps cannot result in false assignments of sequences to samples. As this increases both cost and workload, a metabarcoding library preparation protocol which circumvents the two steps that causes tag-jumps is needed. Here, we demonstrate Tagsteady, a PCR-free metabarcoding Illumina library preparation protocol for pools of nucleotide-tagged amplicons that enables efficient and cost-effective generation of metabarcoding data with virtually no tag-jumps. We use pools of twin-tagged amplicons to investigate the effect of T4 DNA polymerase blunt-ending and postligation PCR on the occurrence of tag-jumps and demonstrate that both blunt-ending and postligation PCR, alone or together, can result in detrimental amounts of tag-jumps (here, up to ca. 49% of total sequences), while leaving both steps out (the Tagsteady protocol) results in amounts of sequences carrying new combinations of used tags (tag-jumps) comparable to background contamination.

Hi-C chromosome conformation capture sequencing of avian genomes using the BGISEQ-500 platform

BACKGROUND

Hi-C experiments couple DNA-DNA proximity with next-generation sequencing to yield an unbiased description of genome-wide interactions. Previous methods describing Hi-C experiments have focused on the industry-standard Illumina sequencing. With new next-generation sequencing platforms such as BGISEQ-500 becoming more widely available, protocol adaptations to fit platform-specific requirements are useful to give increased choice to researchers who routinely generate sequencing data.

RESULTS

We describe an in situ Hi-C protocol adapted to be compatible with the BGISEQ-500 high-throughput sequencing platform. Using zebra finch (Taeniopygia guttata) as a biological sample, we demonstrate how Hi-C libraries can be constructed to generate informative data using the BGISEQ-500 platform, following circularization and DNA nanoball generation. Our protocol is a modification of an Illumina-compatible method, based around blunt-end ligations in library construction, using un-barcoded, distally overhanging double-stranded adapters, followed by amplification using indexed primers. The resulting libraries are ready for circularization and subsequent sequencing on the BGISEQ series of platforms and yield data similar to what can be expected using Illumina-compatible approaches.

CONCLUSIONS

Our straightforward modification to an Illumina-compatible in situHi-C protocol enables data generation on the BGISEQ series of platforms, thus expanding the options available for researchers who wish to utilize the powerful Hi-C techniques in their research.

Complete Mitogenomes of Three Carangidae (Perciformes) Fishes: Genome Description and Phylogenetic Considerations

Carangidae are ecologically and economically important marine fish. The complete mitogenomes of three Carangidae species (Alectis indicus, Decapterus tabl, and Alepes djedaba) were sequenced, characterized, and compared with 29 other species of the family Carangidae in this study. The length of the three mitogenomes ranged from 16,530 to 16,610 bp, and the structures included 2 rRNA genes (12S rRNA and 16S rRNA), 1 control region (a non-coding region), 13 protein-coding genes, and 22 tRNA genes. Among the 22 tRNA genes, only tRNA-Ser (GCT) was not folded into a typical cloverleaf secondary structure and had no recognizable DHU stem. The full-length sequences and protein-coding genes (PCGs) of the mitogenomes of the three species all had obvious AT biases. The majority of the AT-skew and GC-skew values of the PCGs among the three species were negative, demonstrating bases T and C were more plentiful than A and G. Analyses of Ka/Ks and overall p-genetic distance demonstrated that ATP8 showed the highest evolutionary rate and COXI/COXII were the most conserved genes in the three species. The phylogenetic tree based on PCGs sequences of mitogenomes using maximum likelihood and Bayesian inference analyses showed that three clades were divided corresponding to the subfamilies Caranginae, Naucratinae, and Trachinotinae. The monophyly of each superfamily was generally well supported. The divergence time analyses showed that Carangidae evolved during three geological periods, the Cretaceous, Paleogene, and Neogene. A. indicus began to differentiate from other species about 27.20 million years ago (Mya) in the early Miocene, while D. tabl (21.25 Mya) and A. djedaba (14.67 Mya) differentiated in the middle Oligocene.

Comparative performance of the BGI and Illumina sequencing technology for single-cell RNA-sequencing

The libraries generated by high-throughput single cell RNA-sequencing (scRNA-seq) platforms such as the Chromium from 10× Genomics require considerable amounts of sequencing, typically due to the large number of cells. The ability to use these data to address biological questions is directly impacted by the quality of the sequence data. Here we have compared the performance of the Illumina NextSeq 500 and NovaSeq 6000 against the BGI MGISEQ-2000 platform using identical Single Cell 3′ libraries consisting of over 70 000 cells generated on the 10× Genomics Chromium platform. Our results demonstrate a highly comparable performance between the NovaSeq 6000 and MGISEQ-2000 in sequencing quality, and the detection of genes, cell barcodes, Unique Molecular Identifiers. The performance of the NextSeq 500 was also similarly comparable to the MGISEQ-2000 based on the same metrics. Data generated by both sequencing platforms yielded similar analytical outcomes for general single-cell analysis. The performance of the NextSeq 500 and MGISEQ-2000 were also comparable for the deconvolution of multiplexed cell pools via variant calling, and detection of guide RNA (gRNA) from a pooled CRISPR single-cell screen. Our study provides a benchmark for high-capacity sequencing platforms applied to high-throughput scRNA-seq libraries.

Differential Expression of Plasma Exo-miRNA in Neurodegenerative Diseases by Next-Generation Sequencing

Neurodegenerative diseases encompass a wide variety of pathological conditions caused by a loss of neurons in the central nervous system (CNS) and are severely debilitating. Exosome contains bio-signatures of great diagnostic and therapeutic value. There is proof that exosomal proteins can be biomarkers for Alzheimer’s disease (AD) and Parkinson’s disease (PD). MicroRNAs in exosome has potential to be an important source of biomarkers for neurodegenerative diseases. Here, we report exosomal microRNA performance of human plasma in neurodegenerative diseases by small RNA sequencing. A wide range of altered exo-miRNA expression levels were detected in both AD and PD patients. Down-regulated miRNAs in AD samples were enriched in ECM-receptor interaction pathway and both up-/down-regulated miRNAs in PD samples were enriched in fatty acid biosynthesis pathway. Compared to the control, 8 miRNAs were found to be significantly elevated/declined in AD and PD samples, of which 4 miRNAs were newly identified. Additionally, two exosome isolating methods were compared and the reproducibility of plasma exo-miRNA expression was confirmed, suggesting the feasibility of large-scale clinical application of this method. This study revealed exo-miRNA expression levels in neurodegenerative diseases, proposed new biomarkers and their potential functional pathway for AD and PD, confirmed the reproductivity of exo-miRNA profiles by using a different exosome isolating method, and compared the results with plasma miRNA expression. Therefore, this study also provides a precedent for identifying exosomal biomarkers of neurodegenerative diseases in plasma by high-throughput sequencing and it could extend the therapeutic repertoire of exosomal biomarkers.

Comparative analysis of novel MGISEQ-2000 sequencing platform vs Illumina HiSeq 2500 for whole-genome sequencing

The MGISEQ-2000 developed by MGI Tech Co. Ltd. (a subsidiary of the BGI Group) is a new competitor of such next-generation sequencing platforms as NovaSeq and HiSeq (Illumina). Its sequencing principle is based on the DNB and the cPAS technologies, which were also used in the previous version of the BGISEQ-500 device. However, the reagents for MGISEQ-2000 have been refined and the platform utilizes updated software. The cPAS method is an advanced technology based on the cPAL previously created by Complete Genomics. In this paper, the authors compare the results of the whole-genome sequencing of a DNA sample from a Russian female donor performed on MGISEQ-2000 and Illumina HiSeq 2500 (both PE150). Two platforms were compared in terms of sequencing quality, number of errors and performance. Additionally, we performed variant calling using four different software packages: Samtools mpileaup, Strelka2, Sentieon, and GATK. The accuracy of SNP detection was similar in the data generated by MGISEQ-2000 and HiSeq 2500, which was used as a reference. At the same time, a separate indel analysis of the overall error rate revealed similar FPR values and lower sensitivity. It may be concluded with confidence that the data generated by the analyzed sequencing systems is characterized by comparable magnitudes of error and that MGISEQ-2000 and HiSeq 2500 can be used interchangeably for similar tasks like whole genome sequencing.

Recent mitochondrial lineage extinction in the critically endangered Javan rhinoceros

The Javan rhinoceros (Rhinoceros sondaicus) is one of five extant rhinoceros species and among the rarest large mammals on Earth. Once widespread across Southeast Asia, it is now on the verge of extinction, with only one wild population remaining (estimated at ~60 individuals) on the island of Java, Indonesia. To assess the past genetic diversity of the female lineage of R. sondaicus, we generated mitochondrial genome data from eight museum specimens dating back to the 19th century, before the range of the Javan rhinoceros was dramatically reduced, for comparison against mitochondrial DNA (mtDNA) sequences of current R. sondaicus and other rhinoceros species. We succeeded in reconstructing five full and three partial ancient mitogenomes from the eight samples. We used BEAST to assess the phylogenetic relationship of the five extant rhinoceros species and the historical samples. The results show that the oldest and most diverse mtDNA lineages of R. sondaicus are found in historical samples, indicating a significant reduction of mtDNA diversity in modern Javan rhinos. We anticipate that the newly sequenced data will represent a useful resource for improving our understanding of evolutionary history of this species, should future studies be able to increase the available dataset. We hope this information may help in conservation efforts for this species.

Plasmodium vivax Malaria Viewed through the Lens of an Eradicated European Strain

The protozoan Plasmodium vivax is responsible for 42% of all cases of malaria outside Africa. The parasite is currently largely restricted to tropical and subtropical latitudes in Asia, Oceania, and the Americas. Though, it was historically present in most of Europe before being finally eradicated during the second half of the 20th century. The lack of genomic information on the extinct European lineage has prevented a clear understanding of historical population structuring and past migrations of P. vivax. We used medical microscope slides prepared in 1944 from malaria-affected patients from the Ebro Delta in Spain, one of the last footholds of malaria in Europe, to generate a genome of a European P. vivax strain. Population genetics and phylogenetic analyses placed this strain basal to a cluster including samples from the Americas. This genome allowed us to calibrate a genomic mutation rate for P. vivax, and to estimate the mean age of the last common ancestor between European and American strains to the 15th century. This date points to an introduction of the parasite during the European colonization of the Americas. In addition, we found that some known variants for resistance to antimalarial drugs, including Chloroquine and Sulfadoxine, were already present in this European strain, predating their use. Our results shed light on the evolution of an important human pathogen and illustrate the value of antique medical collections as a resource for retrieving genomic information on pathogens from the past.

The complete chloroplast genome of an ornamental orchid, Vanda coerulescens (Orchidaceae)

Though the chloroplast genomes of several Vanda species have been sequenced, there is little information about the complete chloroplast (cp) genome of Vanda coerulescens. Herein, we established the cp genome of V. coerulescens. The chloroplast genome circle was 149,410 bp in length, with the structure of an 85,954 bp large single-copy (LSC) region and a 11,526 bp small single-copy (SSC) region, which separated by two inverted repeat (IRs) regions of 25,965 bp. It encoded 130 genes, including 74 protein-coding genes, 38 tRNA genes and 8 rRNA genes. The overall GC-content of the whole plastome is 36.7%, whereas the corresponding values of the LSC, SSC, and IR regions ranged from 28.2% to 43.1%. In addition, the phylogenetic analysis base on 20 chloroplast genomes of Orchidaceae indicates that V. brunnea is closely related to V. coerulescens. This announcement of the complete V. coerulescens cp genome sequence could provide valuable information for further genetic modification and phylogenetic study in Vanda genus.

Chloroplast characterizations and phylogenetic position of an endangered orchid, Vanda coerulea (Orchidaceae)

Vanda coerulea possess a high ornamental value and medical effect against glaucoma and cataract. The whole complete chloroplast (cp) genome of V. coerulea and the phylogenetic position based on the cp sequences remain unclear. Herein, we report the complete chloroplast genome of V. coerulea. The chloroplast genome was 149,376 bp in length, including a large single-copy (LSC) region of 86,100 bp, a small single-copy (SSC) region of 11,702 bp, and two inverted repeat (IRs) regions of 25,787 bp. A total of 129 genes were characterized, including 74 protein-coding genes, 36 tRNA genes, and 8 rRNA genes. The overall GC content was 36.6%, and GC percentages range from 27.9% to 43.2% throughout LSC, IRs, and SSC regions. Phylogenetic analysis based on 20 chloroplast genomes of Orchidaceae indicated that V. coerulea is closely related to V. brunnea. Our study provides a valuable resource for the identification and distinction of Vanda genus, and will lay a foundation for further research and conservation measures of V. coerulea.

Evolutionary History, Genomic Adaptation to Toxic Diet, and Extinction of the Carolina Parakeet

As the only endemic neotropical parrot to have recently lived in the northern hemisphere, the Carolina parakeet (Conuropsis carolinensis) was an iconic North American bird. The last surviving specimen died in the Cincinnati Zoo in 1918 [1]. The cause of its extinction remains contentious: besides excessive mortality associated to habitat destruction and active hunting, their survival could have been negatively affected by its range having become increasingly patchy [2] or by the exposure to poultry pathogens [3, 4]. In addition, the Carolina parakeet showed a predilection for cockleburs, an herbaceous plant that contains a powerful toxin, carboxyatractyloside, or CAT [5], which did not seem to affect them but made the birds notoriously toxic to most predators [3]. To explore the demographic history of this bird, we generated the complete genomic sequence of a preserved specimen held in a private collection in Espinelves (Girona, Spain), as well as of a close extant relative, Aratinga solstitialis. We identified two non-synonymous genetic changes in two highly conserved proteins known to interact with CAT that could underlie a specific dietary adaptation to this toxin. Our genomic analyses did not reveal evidence of a dramatic past demographic decline in the Carolina parakeet; also, its genome did not exhibit the long runs of homozygosity that are signals of recent inbreeding and are typically found in endangered species. As such, our results suggest its extinction was an abrupt process and thus likely solely attributable to human causes.

Comparative Transcriptomics Analysis of Brassica napus L. during Seed Maturation Reveals Dynamic Changes in Gene Expression between Embryos and Seed Coats and Distinct Expression Profiles of Acyl-CoA-Binding Proteins for Lipid Accumulation

Production of vegetable oils is a vital agricultural resource and oilseed rape (Brassica napus) is the third most important oil crop globally. Although the regulation of lipid biosynthesis in oilseeds is still not fully defined, the acyl-CoA-binding proteins (ACBPs) have been reported to be involved in such metabolism, including oil accumulation, in several plant species. In this study, progressive changes in gene expression in embryos and seed coats at different stages of seed development were comprehensively investigated by transcriptomic analyses in B. napus, revealing dynamic changes in the expression of genes involved in lipid biosynthesis. We show that genes encoding BnACBP proteins show distinct changes in expression at different developmental stages of seed development and show markedly different expression between embryos and seed coats. Both isoforms of the ankyrin-repeat BnACBP2 increased during the oil accumulation period of embryo development. By contrast, the expression of the three most abundant isoforms of the small molecular mass BnACBP6 in embryos showed progressive reduction, despite having the highest overall expression level. In seed coats, BnACBP3, BnACBP4 and BnACBP5 expression remained constant during development, whereas the two major isoforms of BnACBP6 increased, contrasting with the data from embryos. We conclude that genes related to fatty acid and triacylglycerol biosynthesis showing dynamic expression changes may regulate the lipid distribution in embryos and seed coats of B. napus and that BnACBP2 and BnACBP6 are potentially important for oil accumulation.

Basic chloroplast characterizations of Prunus campanulata x kanzakura 'praecox', a warm-adapted cherry cultivar in South China (Rosaceae)

Flowering cherries are domesticated to cultivars by multiple hybridization from some original species, and there formed various adaptive abilities. Herein, we established the complete chloroplast genome of Prunus campanulata x kanzakura ‘Praecox’. The chloroplast genome circle (157,917 bp) presented a typical structure of one 85,928 bp LSC, one 19,117 bp SSC region, and two 26,436 bp IRs. It encoded 124 genes, and same conserved tRNA genes and rRNA genes with P. campanulata, 37 and 8, respectively. The overall GC content was 36.72%. Phylogenetic analysis confirmed the maternal inheritance relationship that P. campanulata x kanzakura ‘Praecox’ comes from Prunus campanulate with little chloroplast genome change, nested inside subgenus Cerasus. This announcement of chloroplast genome would provide helps in further genetic modification and distinction mark study in Prunus genus.

Chloroplast characterizations and phylogenetic location of a common ornamental cherry cultivar, Prunus campanulata 'Kanhizakura-plena' (Rosaceae)

Flowering cherries are well-grown in the world to develop a gorgeous landscape. Though several species are sequenced, there is huge mass of genome-level aberrances between wild species and highly domesticated cultivars. Herein, we established the complete chloroplast genome of Prunus campanulata ‘Kanhizakura-plena’. The chloroplast genome circle (157948 bp) was formed by an 85949 bp large single-copy (LSC) region, a 19127 bp small single-copy (SSC) region, and 2 inverted repeat (IRs) regions of 26436 bp. Independent annotation showed 124 genes were found and conserved tRNA genes and rRNA genes were 37 and 8, respectively. The overall GC content was 36.72%, same with the known species, P. campanulata. Phylogenetic tree confirmed the relationship that P. campanulata ‘Kanhizakura-plena’ is most closely related to P. campanulata, nested inside Prunus. This announcement of chloroplast genome helps genetic modification and phylogenetic study in Prunus genus with useful information. It is a valuable resource for further breeding.

The complete chloroplast genome of a widespread ornamental shrub in China, Magnolia figo (Magnoliaceae)

Plants in Michelia, presented by Magnolia figo DC, are wonderful resources in Magnoliaceae, covering a series of aromatic plants. Despite extensive studies in this family, the M. figo complete chloroplast genome and the taxonomical status based on the whole chloroplast sequences remain unclear. Herein, we report the complete chloroplast genome of M. figo. The chloroplast genome was 160,113 bp in length, with a large single-copy (LSC) region of 88,113 bp and a small single-copy (SSC) region of 18,797 bp, separated by two inverted repeat (protein-coding) regions of 26,602 bp. A total of 135 CDSs were found, including 129 genes, 85 protein-coding mRNAs, 36 tRNA genes, and eight rRNA genes. The overall GC content was 39.3%, and GC percentages range from 34.3% to 43.2% throughout LSC, IRs, and SSC regions. Phylogenetic analysis showed that M. figo is most closely to Michelia odora and displayed a relationship that three Michelia were nested inside Magnolia. This announcement of the complete M. figo cp genome sequence could provide valuable information for further breeding, cp genetic modification, and phylogenetic study in Magnoliaceae.

Ancient RNA from Late Pleistocene permafrost and historical canids shows tissue-specific transcriptome survival

While sequencing ancient DNA (aDNA) from archaeological material is now commonplace, very few attempts to sequence ancient transcriptomes have been made, even from typically stable deposition environments such as permafrost. This is presumably due to assumptions that RNA completely degrades relatively quickly, particularly when dealing with autolytic, nuclease-rich mammalian tissues. However, given the recent successes in sequencing ancient RNA (aRNA) from various sources including plants and animals, we suspect that these assumptions may be incorrect or exaggerated. To challenge the underlying dogma, we generated shotgun RNA data from sources that might normally be dismissed for such study. Here, we present aRNA data generated from two historical wolf skins, and permafrost-preserved liver tissue of a 14,300-year-old Pleistocene canid. Not only is the latter the oldest RNA ever to be sequenced, but it also shows evidence of biologically relevant tissue specificity and close similarity to equivalent data derived from modern-day control tissue. Other hallmarks of RNA sequencing (RNA-seq) data such as exon-exon junction presence and high endogenous ribosomal RNA (rRNA) content confirms our data’s authenticity. By performing independent technical library replicates using two high-throughput sequencing platforms, we show not only that aRNA can survive for extended periods in mammalian tissues but also that it has potential for tissue identification. aRNA also has possible further potential, such as identifying in vivo genome activity and adaptation, when sequenced using this technology.

Ancient DNA is known to survive in cold environments for tens of millennia, but it is assumed that ancient RNA could not persist in such a way due to its relative instability. However, this study shows that under permafrost conditions, ancient RNA can survive well enough to show tissue specificity even in mammalian soft tissues.

Impact of sequencing depth and technology on de novo RNA-Seq assembly

BACKGROUND

RNA-Seq data is inherently nonuniform for different transcripts because of differences in gene expression. This makes it challenging to decide how much data should be generated from each sample. How much should one spend to recover the less expressed transcripts? The sequencing technology used is another consideration, as there are inevitably always biases against certain sequences. To investigate these effects, we first looked at high-depth libraries from a set of well-annotated organisms to ascertain the impact of sequencing depth on de novo assembly. We then looked at libraries sequenced from the Universal Human Reference RNA (UHRR) to compare the performance of Illumina HiSeq and MGI DNBseq™ technologies.

RESULTS

On the issue of sequencing depth, the amount of exomic sequence assembled plateaued using data sets of approximately 2 to 8 Gbp. However, the amount of genomic sequence assembled did not plateau for many of the analyzed organisms. Most of the unannotated genomic sequences are single-exon transcripts whose biological significance will be questionable for some users. On the issue of sequencing technology, both of the analyzed platforms recovered a similar number of full-length transcripts. The missing "gap" regions in the HiSeq assemblies were often attributed to higher GC contents, but this may be an artefact of library preparation and not of sequencing technology.

CONCLUSIONS

Increasing sequencing depth beyond modest data sets of less than 10 Gbp recovers a plethora of single-exon transcripts undocumented in genome annotations. DNBseq™ is a viable alternative to HiSeq for de novo RNA-Seq assembly.

Improving Species Identification of Ancient Mammals Based on Next-Generation Sequencing Data

The taxonomical identification merely based on morphology is often difficult for ancient remains. Therefore, universal or specific PCR amplification followed by sequencing and BLAST (basic local alignment search tool) search has become the most frequently used genetic-based method for the species identification of biological samples, including ancient remains. However, it is challenging for these methods to process extremely ancient samples with severe DNA fragmentation and contamination. Here, we applied whole-genome sequencing data from 12 ancient samples with ages ranging from 2.7 to 700 kya to compare different mapping algorithms, and tested different reference databases, mapping similarities and query coverage to explore the best method and mapping parameters that can improve the accuracy of ancient mammal species identification. The selected method and parameters were tested using 152 ancient samples, and 150 of the samples were successfully identified. We further screened the BLAST-based mapping results according to the deamination characteristics of ancient DNA to improve the ability of ancient species identification. Our findings demonstrate a marked improvement to the normal procedures used for ancient species identification, which was achieved through defining the mapping and filtering guidelines to identify true ancient DNA sequences. The guidelines summarized in this study could be valuable in archaeology, paleontology, evolution, and forensic science. For the convenience of the scientific community, we wrote a software script with Perl, called AncSid, which is made available on GitHub.

Prospect and challenge of detecting dynamic gene copy number increases in stem cells by whole genome sequencing

Abstract

Gene amplification is an evolutionarily well-conserved and highly efficient mechanism to increase the amount of specific proteins. In humans, gene amplification is a hallmark of cancer and has recently been found during stem cell differentiation. Amplifications in stem cells are restricted to specific tissue areas and time windows, rendering their detection difficult. Here, we report on the performance of deep WGS sequencing (average 82-fold depth of coverage) on the BGISEQ with nanoball technology to detect amplifications in human mesenchymal and neural stem cells. As reference technology, we applied array-based comparative genomic hybridization (aCGH), fluorescence in situ hybridization (FISH), and qPCR. Using different in silico strategies for amplification detection, we analyzed the potential of WGS for amplification detection. Our results provide evidence that WGS accurately identifies changes of the copy number profiles in human stem cell differentiation. However, the identified changes are not in all cases consistent between WGS and aCGH. The results between WGS and the validation by qPCR were concordant in 83.3% of all tested 36 cases. In sum, both genome-wide techniques, aCGH and WGS, have unique advantages and specific challenges, calling for locus-specific confirmation by the low-throughput approaches qPCR or FISH.

Key messages

WGS allows for the identification of dynamic copy number changes in human stem cells.

Less stringent threshold setting is crucial for detection of copy number increase.

Broad knowledge of dynamic copy number is pivotal to estimate stem cell capabilities.

Electronic supplementary material

The online version of this article (10.1007/s00109-019-01792-y) contains supplementary material, which is available to authorized users.

Altered respiratory virome and serum cytokine profile associated with recurrent respiratory tract infections in children

Recurrent acute respiratory tract infections (ARTIs) affect a large population, yet the specific decisive factors are largely unknown. Here we study a population of 4407 children diagnosed with ARTI, comparing respiratory virome and serum cytokine profiles associated with multiple ARTIs and single ARTI during a six-year period. The relative abundance of Propionibacterium phages is significantly elevated in multiple ARTIs compared to single ARTI group. Serum levels of TIMP-1 and PDGF-BB are markedly increased in multiple ARTIs compared to single-ARTI and non-ARTI controls, making these two cytokines potential predictors for multiple ARTIs. The presence of Propionibacterium phages is associated with higher levels of TIMP-1 and PDGF-BB. Receiver operating characteristic (ROC) curve analyses show that the combination of TIMP-1, PDGF-BB and Propionibacterium phages could be a strong predictor for multiple ARTIs. These findings indicate that respiratory microbe homeostasis and specific cytokines are associated with the onset of multiple ARTIs over time.

Comparative analysis of sequencing technologies for single-cell transcriptomics

Single-cell RNA-seq technologies require library preparation prior to sequencing. Here, we present the first report to compare the cheaper BGISEQ-500 platform to the Illumina HiSeq platform for scRNA-seq. We generate a resource of 468 single cells and 1297 matched single cDNA samples, performing SMARTer and Smart-seq2 protocols on two cell lines with RNA spike-ins. We sequence these libraries on both platforms using single- and paired-end reads. The platforms have comparable sensitivity and accuracy in terms of quantification of gene expression, and low technical variability. Our study provides a standardized scRNA-seq resource to benchmark new scRNA-seq library preparation protocols and sequencing platforms.

Screening of the candidate genes related to low-temperature tolerance of Fenneropenaeus chinensis based on high-throughput transcriptome sequencing

In order to screen the candidate genes of Fenneropenaeus chinensis related to low-temperature tolerance, this research takes juvenile prawns of F. chinensis (P₄₀) in low temperature stress group (4°C) and normal temperature group (18°C) as experimental materials. The results showed that a total of 127,939 Unigenes with average length of 1,190 bp were obtained by assembly, of which 46% were annotated in the Nr database. A total of 1,698 differentially expressed genes were screened by differential gene expression analysis, of which 920 genes showed up-regulated expression and 778 genes showed down-regulated expression. Both GO and KEGG enrichment analysis revealed that differentially expressed genes were enriched in spliceosomes, ribosomes, bile secretion, ABC transport pathways, and cellular nitrogen compound synthesis. A further in-depth analysis obtained 8 genes that may be associated with low-temperature traits of F. chinensis. Five of them displayed up-regulated expression, including ATP-binding cassette protein C, acid ceramidase, glutathione transferase, C-type lectin and heat shock protein HSP70. The remaining three genes, γ-butyl betaine hydroxylase, β-hexosaminidase A and long chain fatty acid-CoA ligase displayed down-regulated expression. Eight differentially expressed genes were randomly selected and the real time RT-PCR verification showed that their expression levels were consistent with the sequencing results, demonstrating the accuracy of the sequencing results. The results of this study provide basic data for revealing the molecular mechanisms of F. chinensis in response to low temperature stress and the molecular assisted breeding of F. chinensis in low temperature.

Efficient and unique cobarcoding of second-generation sequencing reads from long DNA molecules enabling cost-effective and accurate sequencing, haplotyping, and de novo assembly

Here, we describe single-tube long fragment read (stLFR), a technology that enables sequencing of data from long DNA molecules using economical second-generation sequencing technology. It is based on adding the same barcode sequence to subfragments of the original long DNA molecule (DNA cobarcoding). To achieve this efficiently, stLFR uses the surface of microbeads to create millions of miniaturized barcoding reactions in a single tube. Using a combinatorial process, up to 3.6 billion unique barcode sequences were generated on beads, enabling practically nonredundant cobarcoding with 50 million barcodes per sample. Using stLFR, we demonstrate efficient unique cobarcoding of more than 8 million 20- to 300-kb genomic DNA fragments. Analysis of the human genome NA12878 with stLFR demonstrated high-quality variant calling and phase block lengths up to N50 34 Mb. We also demonstrate detection of complex structural variants and complete diploid de novo assembly of NA12878. These analyses were all performed using single stLFR libraries, and their construction did not significantly add to the time or cost of whole-genome sequencing (WGS) library preparation. stLFR represents an easily automatable solution that enables high-quality sequencing, phasing, SV detection, scaffolding, cost-effective diploid de novo genome assembly, and other long DNA sequencing applications.

Molecular digitization of a botanical garden: high-depth whole-genome sequencing of 689 vascular plant species from the Ruili Botanical Garden

BACKGROUND

Genome sequencing has been widely used in plant research to construct reference genomes and provide evolutionary insights. However, few plant species have had their whole genome sequenced, thus restraining the utility of these data. We collected 1,093 samples of vascular plant species growing in the Ruili Botanical Garden, located in southwest China. Of these, we sequenced 761 samples and collected voucher specimens stored in the Herbarium of China National GeneBank.

RESULTS

The 761 sequenced samples represented 689 vascular plant species from 137 families belonging to 49 orders. Of these, 257 samples were identified to the species level and 504 to the family level, using specimen and chloroplast sequences. In total, we generated 54 Tb of sequencing data, with an average sequencing depth of 60X per species, as estimated from genome sizes. A reference phylogeny was reconstructed with 78 chloroplast genes for molecular identification and other possible applications.

CONCLUSIONS

The large dataset of vascular plant genomes generated in this study, which includes both high-depth whole-genome sequencing data and associated voucher specimens, is valuable for plant genome research and other applications. This project also provides insight into the feasibility and technical requirements for "planetary-scale" projects such as the 10,000 Plant Genomes Project and the Earth BioGenome Project.

MobiSeq: De novo SNP discovery in model and non-model species through sequencing the flanking region of transposable elements

In recent years, the availability of reduced representation library (RRL) methods has catalysed an expansion of genome-scale studies to characterize both model and non-model organisms. Most of these methods rely on the use of restriction enzymes to obtain DNA sequences at a genome-wide level. These approaches have been widely used to sequence thousands of markers across individuals for many organisms at a reasonable cost, revolutionizing the field of population genomics. However, there are still some limitations associated with these methods, in particular the high molecular weight DNA required as starting material, the reduced number of common loci among investigated samples, and the short length of the sequenced site-associated DNA. Here, we present MobiSeq, a RRL protocol exploiting simple laboratory techniques, that generates genomic data based on PCR targeted enrichment of transposable elements and the sequencing of the associated flanking region. We validate its performance across 103 DNA extracts derived from three mammalian species: grey wolf (Canis lupus), red deer complex (Cervus sp.) and brown rat (Rattus norvegicus). MobiSeq enables the sequencing of hundreds of thousands loci across the genome and performs SNP discovery with relatively low rates of clonality. Given the ease and flexibility of MobiSeq protocol, the method has the potential to be implemented for marker discovery and population genomics across a wide range of organisms-enabling the exploration of diverse evolutionary and conservation questions.

Expanding the Clinical Spectrum of Osteogenesis Imperfecta Type V: 13 Additional Patients and Review

Osteogenesis imperfecta (OI) is an inherited connective tissue disorder characterized by bone fragility and is characterized by clinical and genetic heterogeneity. Previous studies showed that the same mutation (c.-14C> T) of the IFITM5 gene is responsible for autosomal dominant OI type V. However, the mutation has a variable expressivity. Clinical heterogeneity has been recognized in OI type V. In this study, we investigated 13 individuals with molecularly confirmed OI type V from seven Chinese families and explored the genotype-phenotype relationship. Increased callus formation is not observed in all individuals, and several novel clinical features were described: joint contractures (three individuals) and unexplained hip arthritis (six individuals). Significant clinical variability was observed even within families. Specific facial features were observed in six individuals from two families consistent with the facial features associated with OI type V reported so far in the literature. Interestingly, we report the process of hypertrophic callus formation in detail for the first time, and in five individuals with hyperplastic callus, increased erythrocyte sedimentation rate (ESR) and levels of C-reactive protein (C-RP) were measured, suggestive of inflammatory activation.