Demographically-Based Evaluation of Genomic Regions under Selection in Domestic Dogs.

Controlling for background demographic effects is important for accurately identifying loci that have recently undergone positive selection. To date, the effects of demography have not yet been explicitly considered when identifying loci under selection during dog domestication. To investigate positive selection on the dog lineage early in the domestication, we examined patterns of polymorphism in six canid genomes that were previously used to infer a demographic model of dog domestication. Using an inferred demographic model, we computed false discovery rates (FDR) and identified 349 outlier regions consistent with positive selection at a low FDR. The signals in the top 100 regions were frequently centered on candidate genes related to brain function and behavior, including LHFPL3, CADM2, GRIK3, SH3GL2, MBP, PDE7B, NTAN1, and GLRA1. These regions contained significant enrichments in behavioral ontology categories. The 3rd top hit, CCRN4L, plays a major role in lipid metabolism, that is supported by additional metabolism related candidates revealed in our scan, including SCP2D1 and PDXC1. Comparing our method to an empirical outlier approach that does not directly account for demography, we found only modest overlaps between the two methods, with 60% of empirical outliers having no overlap with our demography-based outlier detection approach. Demography-aware approaches have lower-rates of false discovery. Our top candidates for selection, in addition to expanding the set of neurobehavioral candidate genes, include genes related to lipid metabolism, suggesting a dietary target of selection that was important during the period when proto-dogs hunted and fed alongside hunter-gatherers.

Genomic and transcriptomic plasticity in treatment-naive ovarian cancer.

Intra-tumor heterogeneity is a hallmark of many cancers and may lead to therapy resistance or interfere with personalized treatment strategies. Here, we combined topographic mapping of somatic breakpoints and transcriptional profiling to probe intra-tumor heterogeneity of treatment-naïve stage IIIC/IV epithelial ovarian cancer. We observed that most substantial differences in genomic rearrangement landscapes occurred between metastases in the omentum and peritoneum versus tumor sites in the ovaries. Several cancer genes such as NF1, CDKN2A, and FANCD2 were affected by lesion-specific breakpoints. Furthermore, the intra-tumor variability involved different mutational hallmarks including lesion-specific kataegis (local mutation shower coinciding with genomic breakpoints), rearrangement classes, and coding mutations. In one extreme case, we identified two independent TP53 mutations in ovary tumors and omentum/peritoneum metastases, respectively. Examination of gene expression dynamics revealed up-regulation of key cancer pathways including WNT, integrin, chemokine, and Hedgehog signaling in only subsets of tumor samples from the same patient. Finally, we took advantage of the multilevel tumor analysis to understand the effects of genomic breakpoints on qualitative and quantitative gene expression changes. We show that intra-tumor gene expression differences are caused by site-specific genomic alterations, including formation of in-frame fusion genes. These data highlight the plasticity of ovarian cancer genomes, which may contribute to their strong capacity to adapt to changing environmental conditions and give rise to the high rate of recurrent disease following standard treatment regimes.

Somatic mutations found in the healthy blood compartment of a 115-yr-old woman demonstrate oligoclonal hematopoiesis.

The somatic mutation burden in healthy white blood cells (WBCs) is not well known. Based on deep whole-genome sequencing, we estimate that approximately 450 somatic mutations accumulated in the nonrepetitive genome within the healthy blood compartment of a 115-yr-old woman. The detected mutations appear to have been harmless passenger mutations: They were enriched in noncoding, AT-rich regions that are not evolutionarily conserved, and they were depleted for genomic elements where mutations might have favorable or adverse effects on cellular fitness, such as regions with actively transcribed genes. The distribution of variant allele frequencies of these mutations suggests that the majority of the peripheral white blood cells were offspring of two related hematopoietic stem cell (HSC) clones. Moreover, telomere lengths of the WBCs were significantly shorter than telomere lengths from other tissues. Together, this suggests that the finite lifespan of HSCs, rather than somatic mutation effects, may lead to hematopoietic clonal evolution at extreme ages.

Genome sequencing highlights the dynamic early history of dogs.

To identify genetic changes underlying dog domestication and reconstruct their early evolutionary history, we generated high-quality genome sequences from three gray wolves, one from each of the three putative centers of dog domestication, two basal dog lineages (Basenji and Dingo) and a golden jackal as an outgroup. Analysis of these sequences supports a demographic model in which dogs and wolves diverged through a dynamic process involving population bottlenecks in both lineages and post-divergence gene flow. In dogs, the domestication bottleneck involved at least a 16-fold reduction in population size, a much more severe bottleneck than estimated previously. A sharp bottleneck in wolves occurred soon after their divergence from dogs, implying that the pool of diversity from which dogs arose was substantially larger than represented by modern wolf populations. We narrow the plausible range for the date of initial dog domestication to an interval spanning 11-16 thousand years ago, predating the rise of agriculture. In light of this finding, we expand upon previous work regarding the increase in copy number of the amylase gene (AMY2B) in dogs, which is believed to have aided digestion of starch in agricultural refuse. We find standing variation for amylase copy number variation in wolves and little or no copy number increase in the Dingo and Husky lineages. In conjunction with the estimated timing of dog origins, these results provide additional support to archaeological finds, suggesting the earliest dogs arose alongside hunter-gathers rather than agriculturists. Regarding the geographic origin of dogs, we find that, surprisingly, none of the extant wolf lineages from putative domestication centers is more closely related to dogs, and, instead, the sampled wolves form a sister monophyletic clade. This result, in combination with dog-wolf admixture during the process of domestication, suggests that a re-evaluation of past hypotheses regarding dog origins is necessary.

Novel drug design for Chagas disease via targeting Trypanosoma cruzi tubulin: Homology modeling and binding pocket prediction on Trypanosoma cruzi tubulin polymerization inhibition by naphthoquinone derivatives.

Chagas disease, also called American trypanosomiasis, is a parasitic disease caused by Trypanosoma cruzi (T. cruzi). Recent findings have underscored the abundance of the causative organism, (T. cruzi), especially in the southern tier states of the US and the risk burden for the rural farming communities there. Due to a lack of safe and effective drugs, there is an urgent need for novel therapeutic options for treating Chagas disease. We report here our first scientific effort to pursue a novel drug design for treating Chagas disease via the targeting of T. cruzi tubulin. First, the anti T. cruzi tubulin activities of five naphthoquinone derivatives were determined and correlated to their anti-trypanosomal activities. The correlation between the ligand activities against the T. cruzi organism and their tubulin inhibitory activities was very strong with a Pearson's r value of 0.88 (P value <0.05), indicating that this class of compounds could inhibit the activity of the trypanosome organism via T. cruzi tubulin polymerization inhibition. Subsequent molecular modeling studies were carried out to understand the mechanisms of the anti-tubulin activities, wherein, the homology model of T. cruzi tubulin dimer was generated and the putative binding site of naphthoquinone derivatives was predicted. The correlation coefficient for ligand anti-tubulin activities and their binding energies at the putative pocket was found to be r=0.79, a high correlation efficiency that was not replicated in contiguous candidate pockets. The homology model of T. cruzi tubulin and the identification of its putative binding site lay a solid ground for further structure based drug design, including molecular docking and pharmacophore analysis. This study presents a new opportunity for designing potent and selective drugs for Chagas disease.

A small-cell lung cancer genome with complex signatures of tobacco exposure.

Cancer is driven by mutation. Worldwide, tobacco smoking is the principal lifestyle exposure that causes cancer, exerting carcinogenicity through >60 chemicals that bind and mutate DNA. Using massively parallel sequencing technology, we sequenced a small-cell lung cancer cell line, NCI-H209, to explore the mutational burden associated with tobacco smoking. A total of 22,910 somatic substitutions were identified, including 134 in coding exons. Multiple mutation signatures testify to the cocktail of carcinogens in tobacco smoke and their proclivities for particular bases and surrounding sequence context. Effects of transcription-coupled repair and a second, more general, expression-linked repair pathway were evident. We identified a tandem duplication that duplicates exons 3-8 of CHD7 in frame, and another two lines carrying PVT1-CHD7 fusion genes, indicating that CHD7 may be recurrently rearranged in this disease. These findings illustrate the potential for next-generation sequencing to provide unprecedented insights into mutational processes, cellular repair pathways and gene networks associated with cancer.

Genetic Characteristics of Polycythemia Vera and Essential Thrombocythemia in Korean Patients.

BACKGROUND: Despite recent advances in the investigation of myeloproliferative neoplasms (MPN), the impact of genetic heterogeneity on its molecular pathogenesis has not been fully elucidated. Thus, in this study, we aim to characterize the genetic complexity in Korean patients with polycythemia vera (PV) and essential thrombocythemia (ET). METHODS: We conducted association studies using 84 single-nucleotide polymorphisms (SNPs) in 229 patients (96 with PV and 133 with ET) and 170 controls. Further, whole-genome sequencing was performed in six patients (two with JAK2 V617F and four with wild-type JAK2), and putative somatic mutations were validated in a further 69 ET patients. Clinical and laboratory characteristics were also analyzed. RESULTS: Several germline SNPs and the 46 haplotype were significantly associated with PV and ET. Three somatic mutations in MPDZ, IQCH, and CALR genes were selected and validated. The frequency of the CALR mutation was 58.0% (40/69) in ET patients, who did not carry JAK2/MPL mutations. Moreover, compared with JAK2 V617F-positive patients, those with CALR mutations showed lower hemoglobin and hematocrit levels (P = 0.004 and P = 0.002, respectively), higher platelet counts (P =0.008), and a lower frequency of cytoreductive therapy (P = 0.014). CONCLUSION: This study was the first comprehensive investigation of the genetic characteristics of Korean patients with PV and ET. We found that somatic mutations and the 46 haplotype contribute to PV and ET pathogenesis in Korean patients.

Molecular assessment of differences in the duodenal microbiome in subjects with irritable bowel syndrome.

OBJECTIVE: Breath testing and duodenal culture studies suggest that a significant proportion of irritable bowel syndrome (IBS) patients have small intestinal bacterial overgrowth. In this study, we extended these data through 16S rDNA amplicon sequencing and quantitative PCR (qPCR) analyses of duodenal aspirates from a large cohort of IBS, non-IBS and control subjects. MATERIALS AND METHODS: Consecutive subjects presenting for esophagogastroduodenoscopy only and healthy controls were recruited. Exclusion criteria included recent antibiotic or probiotic use. Following extensive medical work-up, patients were evaluated for symptoms of IBS. DNAs were isolated from duodenal aspirates obtained during endoscopy. Microbial populations in a subset of IBS subjects and controls were compared by 16S profiling. Duodenal microbes were then quantitated in the entire cohort by qPCR and the results compared with quantitative live culture data. RESULTS: A total of 258 subjects were recruited (21 healthy, 163 non-healthy non-IBS, and 74 IBS). 16S profiling in five IBS and five control subjects revealed significantly lower microbial diversity in the duodenum in IBS, with significant alterations in 12 genera (false discovery rate < 0.15), including overrepresentation of Escherichia/Shigella (p = 0.005) and Aeromonas (p = 0.051) and underrepresentation of Acinetobacter (p = 0.024), Citrobacter (p = 0.031) and Microvirgula (p = 0.036). qPCR in all 258 subjects confirmed greater levels of Escherichia coli in IBS and also revealed increases in Klebsiella spp, which correlated strongly with quantitative culture data. CONCLUSIONS: 16S rDNA sequencing confirms microbial overgrowth in the small bowel in IBS, with a concomitant reduction in diversity. qPCR supports alterations in specific microbial populations in IBS.

Comparative genomics of 274 Vibrio cholerae genomes reveals mobile functions structuring three niche dimensions.

BACKGROUND: Vibrio cholerae is a globally dispersed pathogen that has evolved with humans for centuries, but also includes non-pathogenic environmental strains. Here, we identify the genomic variability underlying this remarkable persistence across the three major niche dimensions space, time, and habitat. RESULTS: Taking an innovative approach of genome-wide association applicable to microbial genomes (GWAS-M), we classify 274 complete V. cholerae genomes by niche, including 39 newly sequenced for this study with the Ion Torrent DNA-sequencing platform. Niche metadata were collected for each strain and analyzed together with comprehensive annotations of genetic and genomic attributes, including point mutations (single-nucleotide polymorphisms, SNPs), protein families, functions and prophages. CONCLUSIONS: Our analysis revealed that genomic variations, in particular mobile functions including phages, prophages, transposable elements, and plasmids underlie the metadata structuring in each of the three niche dimensions. This underscores the role of phages and mobile elements as the most rapidly evolving elements in bacterial genomes, creating local endemicity (space), leading to temporal divergence (time), and allowing the invasion of new habitats. Together, we take a data-driven approach for comparative functional genomics that exploits high-volume genome sequencing and annotation, in conjunction with novel statistical and machine learning analyses to identify connections between genotype and phenotype on a genome-wide scale.

A comprehensively molecular haplotype-resolved genome of a European individual.

Independent determination of both haplotype sequences of an individual genome is essential to relate genetic variation to genome function, phenotype, and disease. To address the importance of phase, we have generated the most complete haplotype-resolved genome to date, "Max Planck One" (MP1), by fosmid pool-based next generation sequencing. Virtually all SNPs (>99%) and 80,000 indels were phased into haploid sequences of up to 6.3 Mb (N50 ~1 Mb). The completeness of phasing allowed determination of the concrete molecular haplotype pairs for the vast majority of genes (81%) including potential regulatory sequences, of which >90% were found to be constituted by two different molecular forms. A subset of 159 genes with potentially severe mutations in either cis or trans configurations exemplified in particular the role of phase for gene function, disease, and clinical interpretation of personal genomes (e.g., BRCA1). Extended genomic regions harboring manifold combinations of physically and/or functionally related genes and regulatory elements were resolved into their underlying "haploid landscapes," which may define the functional genome. Moreover, the majority of genes and functional sequences were found to contain individual or rare SNPs, which cannot be phased from population data alone, emphasizing the importance of molecular phasing for characterizing a genome in its molecular individuality. Our work provides the foundation to understand that the distinction of molecular haplotypes is essential to resolve the (inherently individual) biology of genes, genomes, and disease, establishing a reference point for "phase-sensitive" personal genomics. MP1's annotated haploid genomes are available as a public resource.

Improving mammalian genome scaffolding using large insert mate-pair next-generation sequencing.

BACKGROUND: Paired-tag sequencing approaches are commonly used for the analysis of genome structure. However, mammalian genomes have a complex organization with a variety of repetitive elements that complicate comprehensive genome-wide analyses. RESULTS: Here, we systematically assessed the utility of paired-end and mate-pair (MP) next-generation sequencing libraries with insert sizes ranging from 170 bp to 25 kb, for genome coverage and for improving scaffolding of a mammalian genome (Rattus norvegicus). Despite a lower library complexity, large insert MP libraries (20 or 25 kb) provided very high physical genome coverage and were found to efficiently span repeat elements in the genome. Medium-sized (5, 8 or 15 kb) MP libraries were much more efficient for genome structure analysis than the more commonly used shorter insert paired-end and 3 kb MP libraries. Furthermore, the combination of medium- and large insert libraries resulted in a 3-fold increase in N50 in scaffolding processes. Finally, we show that our data can be used to evaluate and improve contig order and orientation in the current rat reference genome assembly. CONCLUSIONS: We conclude that applying combinations of mate-pair libraries with insert sizes that match the distributions of repetitive elements improves contig scaffolding and can contribute to the finishing of draft genomes.

Chromatin and sequence features that define the fine and gross structure of genomic methylation patterns.

Abnormalities of genomic methylation patterns are lethal or cause disease, but the cues that normally designate CpG dinucleotides for methylation are poorly understood. We have developed a new method of methylation profiling that has single-CpG resolution and can address the methylation status of repeated sequences. We have used this method to determine the methylation status of >275 million CpG sites in human and mouse DNA from breast and brain tissues. Methylation density at most sequences was found to increase linearly with CpG density and to fall sharply at very high CpG densities, but transposons remained densely methylated even at higher CpG densities. The presence of histone H2A.Z and histone H3 di- or trimethylated at lysine 4 correlated strongly with unmethylated DNA and occurred primarily at promoter regions. We conclude that methylation is the default state of most CpG dinucleotides in the mammalian genome and that a combination of local dinucleotide frequencies, the interaction of repeated sequences, and the presence or absence of histone variants or modifications shields a population of CpG sites (most of which are in and around promoters) from DNA methyltransferases that lack intrinsic sequence specificity.

A comprehensive assay for CFTR mutational analysis using next-generation sequencing.

BACKGROUND: Cystic fibrosis is a life-threatening genetic disorder that has been associated with mutations in the CFTR [cystic fibrosis transmembrane conductance regulator (ATP-binding cassette sub-family C, member 7)] gene. Hundreds of CFTR mutations have been detected to date. Current CFTR genotyping assays target a subset of these mutations, particularly a mutation panel recommended by the American College of Medical Genetics for carrier screening of the general population. Fast sequencing of the entire coding sequence in a scalable manner could expand the detection of CFTR mutations and facilitate management of costs and turnaround times in the clinical laboratory. METHODS: We describe a proof-of-concept CFTR assay that uses PCR target enrichment and next-generation sequencing on the Ion Torrent Personal Genome Machine™ (PGM™) platform. RESULTS: The scalability of the assay was demonstrated, with an average mean depth of coverage ranging from 500× to 3500×, depending on the number of multiplexed patient samples and the Ion Torrent chip used. In a blinded study of 79 previously genotyped patient DNA samples and cell lines, our assay detected most of the mutations, including single-nucleotide variants, small insertions and deletions, and large copy-number variants. The reproducibility was 100% for detecting mutations in independent runs. Our assay demonstrated high specificity, with only 2 false-positive calls (at 2184delA) found in 2 samples caused by a sequencing error in a homopolymer stretch of sequence. The detection rate for variants of unknown significance was very low in the targeted region. CONCLUSIONS: With continued optimization and system refinements, PGM sequencing promises to be a powerful, rapid, and scalable means of clinical diagnostic sequencing.

Cell-free DNA reveals distinct pathology of multisystem inflammatory syndrome in children.

Multisystem inflammatory syndrome in children (MIS-C) is a rare but life-threatening hyperinflammatory condition induced by infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) that causes pediatric COVID-19 (pCOVID-19). The relationship of the systemic tissue injury to the pathophysiology of MIS-C is poorly defined. We leveraged the high sensitivity of epigenomics analyses of plasma cell-free DNA (cfDNA) and plasma cytokine measurements to identify the spectrum of tissue injury and glean mechanistic insights. Compared with pediatric healthy controls (pHCs) and patients with pCOVID-19, patients with MIS-C had higher levels of cfDNA primarily derived from innate immune cells, megakaryocyte-erythroid precursor cells, and nonhematopoietic tissues such as hepatocytes, cardiac myocytes, and kidney cells. Nonhematopoietic tissue cfDNA levels demonstrated significant interindividual variability, consistent with the heterogenous clinical presentation of MIS-C. In contrast, adaptive immune cell-derived cfDNA levels were comparable in MIS-C and pCOVID-19 patients. Indeed, the cfDNA of innate immune cells in patients with MIS-C correlated with the levels of innate immune inflammatory cytokines and nonhematopoietic tissue-derived cfDNA, suggesting a primarily innate immunity-mediated response to account for the multisystem pathology. These data provide insight into the pathogenesis of MIS-C and support the value of cfDNA as a sensitive biomarker to map tissue injury in MIS-C and likely other multiorgan inflammatory conditions.

Sequence and structural variation in a human genome uncovered by short-read, massively parallel ligation sequencing using two-base encoding.

We describe the genome sequencing of an anonymous individual of African origin using a novel ligation-based sequencing assay that enables a unique form of error correction that improves the raw accuracy of the aligned reads to >99.9%, allowing us to accurately call SNPs with as few as two reads per allele. We collected several billion mate-paired reads yielding approximately 18x haploid coverage of aligned sequence and close to 300x clone coverage. Over 98% of the reference genome is covered with at least one uniquely placed read, and 99.65% is spanned by at least one uniquely placed mate-paired clone. We identify over 3.8 million SNPs, 19% of which are novel. Mate-paired data are used to physically resolve haplotype phases of nearly two-thirds of the genotypes obtained and produce phased segments of up to 215 kb. We detect 226,529 intra-read indels, 5590 indels between mate-paired reads, 91 inversions, and four gene fusions. We use a novel approach for detecting indels between mate-paired reads that are smaller than the standard deviation of the insert size of the library and discover deletions in common with those detected with our intra-read approach. Dozens of mutations previously described in OMIM and hundreds of nonsynonymous single-nucleotide and structural variants in genes previously implicated in disease are identified in this individual. There is more genetic variation in the human genome still to be uncovered, and we provide guidance for future surveys in populations and cancer biopsies.

mRNA-Seq whole-transcriptome analysis of a single cell.

Next-generation sequencing technology is a powerful tool for transcriptome analysis. However, under certain conditions, only a small amount of material is available, which requires more sensitive techniques that can preferably be used at the single-cell level. Here we describe a single-cell digital gene expression profiling assay. Using our mRNA-Seq assay with only a single mouse blastomere, we detected the expression of 75% (5,270) more genes than microarray techniques and identified 1,753 previously unknown splice junctions called by at least 5 reads. Moreover, 8-19% of the genes with multiple known transcript isoforms expressed at least two isoforms in the same blastomere or oocyte, which unambiguously demonstrated the complexity of the transcript variants at whole-genome scale in individual cells. Finally, for Dicer1(-/-) and Ago2(-/-) (Eif2c2(-/-)) oocytes, we found that 1,696 and 1,553 genes, respectively, were abnormally upregulated compared to wild-type controls, with 619 genes in common.

Human U6 promoter drives stronger shRNA activity than its schistosome orthologue in Schistosoma mansoni and human fibrosarcoma cells.

Blood flukes or schistosomes are the causative agents of human schistosomiasis, one of the major neglected tropical diseases. Draft genome sequences have been reported for schistosomes, but functional genomics tools are needed to investigate the role and essentiality of the newly reported genes. Vector based RNA interference can contribute to functional genomics analysis for schistosomes. Using mRNA encoding reporter firefly luciferase as a model target, we compared the performance of a schistosome and a human promoter from the U6 gene in driving shRNA in human fibrosarcoma cells and in cultured schistosomes. Further, both a retroviral [Murine leukemia virus (MLV)] and plasmid (piggyBac, pXL-Bac II) vector were utilized. The schistosome U6 gene promoter was 270 bp in length, the human U6 gene promoter was 264 bp; they shared 41% identity. Following transduction of both HT1080 fibrosarcoma cells and schistosomules of Schistosoma mansoni with pseudotyped MLV virions, stronger knockdown of luciferase activity was seen with the virions encoding the human U6 promoter driven shRNA than the schistosome U6 promoter. A similar trend was seen after transfection of HT1080 cells and schistosomules with the pXL-Bac-II constructs-stronger knockdown of luciferase activity was seen with constructs encoding the human compared to schistosome U6 promoter. The findings indicate that a human U6 gene promoter drives stronger shRNA activity than its schistosome orthologue, not only in a human cancer cell line but also in larval schistosomes. This RNA polymerase III promoter represents a potentially valuable component for vector based RNA interference studies in schistosomes and related platyhelminth parasites.

Trichomoniasis, a new look at a common but neglected STI in African descendance population in the United States and the Black Diaspora. A review of its incidence, research prioritization, and the resulting health disparities.

Structural racism plays a significant role in limited access to higher education, financial resources, employment opportunities, and high-quality healthcare for African Americans. The lack of healthcare equity and infrastructure has directly contributed to overall poor healthcare outcomes for the Black community. Studies have shown that adverse health outcomes such as sexually transmitted diseases (STDs) are more prevalent in African Americans, regardless of their socioeconomic factors and lifestyles. For example, trichomoniasis, transmitted sexually by its etiological agent, Trichomonas vaginalis (T. vaginalis), predisposes those infected to co-infections with other STDs, such as human immunodeficiency virus (HIV), herpes, and other related infections. Our review showcases the impact of trichomoniasis on the health of the Black community with an emphasis on African American women. A critical examination of the socio-demographic history of Black people in the United States (US) is vital to illustrate the origin of past and current racial health disparities. Further, we expand the complex and nuanced conversation on the intersectionality of racism, health equity, and innovative epidemiological and biomedical research strategies needed to eradicate this global public health threat.

Leishmaniasis: Recent epidemiological studies in the Middle East.

Leishmaniasis, one of the most neglected tropical diseases (NTDs), is the third most important vector-borne disease worldwide. This disease has a global impact and severity of the infection and is greatest in the Middle East. The agent of infection is a protozoan parasite of the genus, Leishmania, and is generally transmitted by blood-sucking female sandflies. In humans, there are three clinical forms of infection: (1) cutaneous (CL), (2) mucocutaneous (ML), and (3) visceral leishmaniasis (VL). This review aims to discuss the current epidemiological status of leishmaniasis in Saudi Arabia, Iraq, Syria, and Yemen with a consideration of treatment options. The elevated risk of leishmaniasis is influenced by the transmission of the disease across endemic countries into neighboring non-infected regions.

Stem cell transcriptome profiling via massive-scale mRNA sequencing.

We developed a massive-scale RNA sequencing protocol, short quantitative random RNA libraries or SQRL, to survey the complexity, dynamics and sequence content of transcriptomes in a near-complete fashion. This method generates directional, random-primed, linear cDNA libraries that are optimized for next-generation short-tag sequencing. We surveyed the poly(A)(+) transcriptomes of undifferentiated mouse embryonic stem cells (ESCs) and embryoid bodies (EBs) at an unprecedented depth (10 Gb), using the Applied Biosystems SOLiD technology. These libraries capture the genomic landscape of expression, state-specific expression, single-nucleotide polymorphisms (SNPs), the transcriptional activity of repeat elements, and both known and new alternative splicing events. We investigated the impact of transcriptional complexity on current models of key signaling pathways controlling ESC pluripotency and differentiation, highlighting how SQRL can be used to characterize transcriptome content and dynamics in a quantitative and reproducible manner, and suggesting that our understanding of transcriptional complexity is far from complete.