Diverse and atypical manifestations of Q fever in a metropolitan city hospital: Emerging role of next-generation sequencing for laboratory diagnosis of Coxiella burnetii

Although Q fever has been widely reported in the rural areas of China, there is a paucity of data on the epidemiology and clinical characteristics of this disease in large metropolitan cities. In this study, we profile the epidemiology and clinical manifestations of Q fever from a tertiary hospital in Shenzhen, a Southern Chinese metropolitan city with a large immigrant population from other parts of China. A total of 14 patients were confirmed to have Q fever during a nine-year-and-six-month period, five of whom were retrospectively diagnosed during case review or incidentally picked up because of another research project on unexplained fever without localizing features. Some patients had the typical exposure histories and clinical features, while a few other patients had rare manifestations of Q fever, including one with heart failure and diffuse intracapillary proliferative glomerulonephritis, a patient presenting with a spontaneous bacterial peritonitis-like syndrome, and another one with concomitant Q fever and brucellosis. Using a combination of clinical manifestation, inflammatory marker levels, echocardiographic findings and serological or molecular test results, nine, three and two patients were diagnosed to have acute, chronic and convalescent Q fever, respectively. Seven, five and two patients were diagnosed to have Q fever by serological test, nested real-time PCR and next-generation sequencing respectively. Diverse and atypical manifestations are associated with Q fever. The incidence of Q fever is likely to be underestimated. Next-generation sequencing is becoming an important diagnostic modality for culture-negative infections, particularly those that the physicians fail to recognize clinically, such as Q fever.     

Identification of extracellular miRNA in human cerebrospinal fluid by next-generation sequencing

There has been a growing interest in using next-generation sequencing (NGS) to profile extracellular small RNAs from the blood and cerebrospinal fluid (CSF) of patients with neurological diseases, CNS tumors, or traumatic brain injury for biomarker discovery. Small sample volumes and samples with low RNA abundance create challenges for downstream small RNA sequencing assays. Plasma, serum, and CSF contain low amounts of total RNA, of which small RNAs make up a fraction. The purpose of this study was to maximize RNA isolation from RNA-limited samples and apply these methods to profile the miRNA in human CSF by small RNA deep sequencing. We systematically tested RNA isolation efficiency using ten commercially available kits and compared their performance on human plasma samples. We used RiboGreen to quantify total RNA yield and custom TaqMan assays to determine the efficiency of small RNA isolation for each of the kits. We significantly increased the recovery of small RNA by repeating the aqueous extraction during the phenol-chloroform purification in the top performing kits. We subsequently used the methods with the highest small RNA yield to purify RNA from CSF and serum samples from the same individual. We then prepared small RNA sequencing libraries using Illumina’s TruSeq sample preparation kit and sequenced the samples on the HiSeq 2000. Not surprisingly, we found that the miRNA expression profile of CSF is substantially different from that of serum. To our knowledge, this is the first time that the small RNA fraction from CSF has been profiled using next-generation sequencing.     

From genomics to functional markers in the era of next-generation sequencing

The availability of complete genome sequences, along with other genomic resources for Arabidopsis, rice, pigeon pea, soybean and other crops, has revolutionized our understanding of the genetic make-up of plants. Next-generation DNA sequencing (NGS) has facilitated single nucleotide polymorphism discovery in plants. Functionally-characterized sequences can be identified and functional markers (FMs) for important traits can be developed at an ever-increasing ease. FMs are derived from sequence polymorphisms found in allelic variants of a functional gene. Linkage disequilibrium-based association mapping and homologous recombinants have been developed for identification of “perfect” markers for their use in crop improvement practices. Compared with many other molecular markers, FMs derived from the functionally characterized sequence genes using NGS techniques and their use provide opportunities to develop high-yielding plant genotypes resistant to various stresses at a fast pace.     

Sporadic outbreaks of crimean-congo haemorrhagic fever in Uganda,July 2018-January 2019

IntroductionCrimean-Congo haemorrhagic fever (CCHF) is a tick-borne, zoonotic viral disease that causes haemorrhagic symptoms. Despite having eight confirmed outbreaks between 2013 and 2017, all within Uganda’s ‘cattle corridor’, no targeted tick control programs exist in Uganda to prevent disease. During a seven-month-period from July 2018-January 2019, the Ministry of Health confirmed multiple independent CCHF outbreaks. We investigated to identify risk factors and recommend interventions to prevent future outbreaks.MethodsWe defined a confirmed case as sudden onset of fever (≥37.5°C) with ≥4 of the following signs and symptoms: anorexia, vomiting, diarrhoea, headache, abdominal pain, joint pain, or sudden unexplained bleeding in a resident of the affected districts who tested positive for Crimean-Congo haemorrhagic fever virus (CCHFv) by RT-PCR from 1 July 2018–30 January 2019. We reviewed medical records and performed active case-finding. We conducted a case-control study and compared exposures of case-patients with age-, sex-, and sub-county-matched control-persons (1:4).ResultsWe identified 14 confirmed cases (64% males) with five deaths (case-fatality rate: 36%) from 11 districts in western and central region. Of these, eight (73%) case-patients resided in Uganda’s ‘cattle corridor’. One outbreak involved two case-patients and the remainder involved one. All case-patients had fever and 93% had unexplained bleeding. Case-patients were aged 6–36 years, with persons aged 20–44 years more affected (AR: 7.2/1,000,000) than persons ≤19 years (2.0/1,000,000), p = 0.015. Most (93%) case-patients had contact with livestock ≤2 weeks before symptom onset. Twelve (86%) lived <1 km from grazing fields compared with 27 (48%) controls (OR_M-H = 18, 95% CI = 3.2-∞) and 10 (71%) of 14 case-patients found ticks attached to their bodies ≤2 weeks before symptom onset, compared to 15 (27%) of 56 control-persons (OR_M-H = 9.3, 95%CI = 1.9–46).ConclusionsCCHF outbreaks occurred sporadically during 2018–2019, both within and outside ‘cattle corridor’ districts of Uganda. Most cases were associated with tick exposure. The Ministry of Health should partner with the Ministry of Agriculture, Animal Industry and Fisheries to develop joint nationwide tick control programs and strategies with shared responsibilities through a One Health approach.     

Random sampling in metagenomic sequencing leads to overestimated spatial scaling of microbial diversity

Revealing the spatial scaling patterns of microbial diversity is of special interest in microbial ecology. One critical question is whether the observed spatial turnover rate truly reflect the actual spatial patterns of extremely diverse microbial communities. Using simulated mock communities, this study suggested that the currently observed microbial spatial turnover rates were overestimated by random sampling processes associated with high-throughput metagenomic sequencing. The observed z values were largely contributed by accumulated microbial taxa due to cumulative number of samples. This is a crucial issue because microbial communities already have very low spatial turnover rate due to the small size and potential cosmopolitism nature of microorganisms. Further investigations suggested a linear relationship between the observed and expected z values, which can be applied to remove random sampling noises from the observed z values. Adjustment of z values for data sets from six American forests showed much lower spatial turnover rate than that before adjustment. However, the patterns of z values among these six forests remained unchanged. This study suggested that our current understanding of microbial taxa–area relationships could be inaccurate. Therefore, cautions and efforts should be made for more accurate estimation and interpretation of microbial spatial patterns.     

Use of next-generation DNA sequencing to analyze genetic variants in rheumatic disease

Next-generation DNA sequencing has revolutionized the field of genetics and genomics, providing researchers with the tools to efficiently identify novel rare and low frequency risk variants, which was not practical with previously available methodologies. These methods allow for the sequence capture of a specific locus or small genetic region all the way up to the entire six billion base pairs of the diploid human genome. Rheumatic diseases are a huge burden on the US population, affecting more than 46 million Americans. Those afflicted suffer from one or more of the more than 100 diseases characterized by inflammation and loss of function, mainly of the joints, tendons, ligaments, bones, and muscles. While genetics studies of many of these diseases (for example, systemic lupus erythematosus, rheumatoid arthritis, and inflammatory bowel disease) have had major successes in defining their genetic architecture, causal alleles and rare variants have still been elusive. This review describes the current high-throughput DNA sequencing methodologies commercially available and their application to rheumatic diseases in both case–control as well as family-based studies.     

Development of genomic simple sequence repeat markers in faba bean by next-generation sequencing

Faba bean (Vicia faba L.) is an important food legume crop with a huge genome. Development of genetic markers for faba bean is important to study diversity and for molecular breeding. In this study, we used Next Generation Sequencing (NGS) technology for the development of genomic simple sequence repeat (SSR) markers. A total of 14,027,500 sequence reads were obtained comprising 4,208 Mb. From these reads, 56,063 contigs were assembled (16,367 Mb) and 2138 SSRs were identified. Mono and dinucleotides were the most abundant, accounting for 57.5 % and 20.9 % of all SSR repeats, respectively. A total of 430 primer pairs were designed from contigs larger than 350 nucleotides and 50 primers pairs were tested for validation of SSR locus amplification. Nearly all (96 %) of the markers were found to produce clear amplicons and to be reproducible. Thirty-nine SSR markers were then applied to 46 faba bean accessions from worldwide origins, resulting in 161 alleles with 87.5 % polymorphism, and an average of 4.1 alleles per marker. Gene diversity (GD) of the markers ranged from 0 to 0.48 with an average of 0.27. Testing of the markers showed that they were useful in determining genetic relationships and population structure in faba bean accessions.     

Development of genomic simple sequence repeat markers in opium poppy by next-generation sequencing

Opium poppy (Papaver somniferum L.) is an important pharmaceutical crop with very few genetic marker resources. To expand these resources, we sequenced genomic DNA using pyrosequencing technology and examined the DNA sequences for simple sequence repeats (SSRs). A total of 1,244,412 sequence reads were obtained covering 474 Mb. Approximately half of the reads (52 %) were assembled into 166,724 contigs representing 105 Mb of the opium poppy genome. A total of 23,283 non-redundant SSRs were identified in 18,944 contigs (11.3 % of total contigs). Trinucleotide and tetranucleotide repeats were the most abundant SSR repeats, accounting for 49.0 and 27.9 % of all SSRs, respectively. The AAG/TTC repeat was the most abundant trinucleotide repeat, representing 19.7 % of trinucleotide repeats. Other SSR repeat types were AT-rich. A total of 23,126 primer pairs (98.7 % of total SSRs) were designed to amplify SSRs. Fifty-three genomic SSR markers were tested in 37 opium poppy accessions and seven Papaver species for determination of polymorphism and transferability. Intraspecific polymorphism information content (PIC) values of the genomic SSR markers were intermediate, with an average 0.17, while the interspecific average PIC value was slightly higher, 0.19. All markers showed at least 88 % transferability among related species. This study increases sequence coverage of the opium poppy genome by sevenfold and the number of opium poppy-specific SSR markers by sixfold. This is the first report of the development of genomic SSR markers in opium poppy, and the genomic SSR markers developed in this study will be useful in diversity, identification, mapping and breeding studies in opium poppy.     

The application of next-generation sequencing in the autozygosity mapping of human recessive diseases

Autozygosity, or the inheritance of two copies of an ancestral allele, has the potential to not only reveal phenotypes caused by biallelic mutations in autosomal recessive genes, but to also facilitate the mapping of such mutations by flagging the surrounding haplotypes as tractable runs of homozygosity (ROH), a process known as autozygosity mapping. Since SNPs replaced microsatellites as markers for the purpose of genomewide identification of ROH, autozygosity mapping of Mendelian genes has witnessed a significant acceleration. Historically, successful mapping traditionally required favorable family structure that permits the identification of an autozygous interval that is amenable to candidate gene selection and confirmation by Sanger sequencing. This requirement presented a major bottleneck that hindered the utilization of simplex cases and many multiplex families with autosomal recessive phenotypes. However, the advent of next-generation sequencing that enables massively parallel sequencing of DNA has largely bypassed this bottleneck and thus ushered in an era of unprecedented pace of Mendelian disease gene discovery. The ability to identify a single causal mutation among a massive number of variants that are uncovered by next-generation sequencing can be challenging, but applying autozygosity as a filter can greatly enhance the enrichment process and its throughput. This review will discuss the power of combining the best of both techniques in the mapping of recessive disease genes and offer some tips to troubleshoot potential limitations.     

Complete mitochondrial genome of Membranipora grandicella (Bryozoa: Cheilostomatida) determined with next-generation sequencing: The first representative of the suborder Malacostegina

Next-generation sequencing (NGS) has proven a valuable platform for fast and easy obtaining of large numbers of sequences at relatively low cost. In this study we use shot-gun sequencing method on Illumina HiSeq 2000, to obtain enough sequences for the assembly of the bryozoan Membranipora grandicella (Bryozoa: Cheilostomatida) mitochondrial genome, which is the first representative of the suborder Malacostegina. The complete mitochondrial genome is 15,861 bp in length, which is relatively larger than other studied bryozoans. The mitochondrial genome contains 13 protein-coding genes, 2 ribosomal RNAs and 20 transfer RNAs. To investigate the phylogenetic position and the inner relationships of the phylum Bryozoa, phylogenetic trees were constructed with amino acid sequences of 11 PCGs from 30 metazoans. Two superclades of protostomes, namely Lophotrochozoa and Ecdysozoa, are recovered as monophyletic with strong support in both ML and Bayesian analyses. Somewhat to surprise, Bryozoa appears as the sister group of Chaetognatha with moderate or high support. The relationship among five bryozoans is Tubulipora flabellaris + (M. grandicella + (Flustrellidra hispida + (Bugula neritina + Watersipora subtorquata))), which supports for the view that Cheilostomatida is not a natural, monophyletic clade. NGS proved to be a quick and easy method for sequencing a complete mitochondrial genome.     

下一代测序技术在干细胞转录调控研究中的应用

基因转录调控及其机制分析是后基因组时代生物学研究的重点之一。随着高通量测序技术的发展,人们可以从不同层面研究基因的转录调控行为,从转录组、转录因子结合,到染色质局部结构和整体空间构象,可系统分析转录调控的分子机制。干细胞分化过程的转录调控分析对研究再生医学和理解细胞癌变机制等具有重要意义。本文综述了下一代测序技术在干细胞转录调控研究中的应用,包括：(1)基于基因芯片或RNA测序的转录组分析;(2)基于染色体免疫共沉淀(chromatin immunoprecipitation, ChIP)测序的表观基因组和转录因子结合信息的分析;(3)基于DNase 酶切测序(DNase-Seq)的染色质开放性分析;(4)基于高通量染色质构象捕获(high-throughput chromosome conformation capture, Hi-C)技术的染色体远程相互作用分析。从基因表达谱、转录因子结合和基因组三维结构等层面展开介绍,重点关注了一些多能性转录因子(Oct4、Sox2和Nanog等)在维持干细胞干性和分化中的调控作用,以期为干细胞转录调控的研究提供借鉴和参考。