Genomic and transcriptomic analysis of the streptomycin-dependent Mycobacterium tuberculosis strain 18b.

The ability of Mycobacterium tuberculosis to establish a latent infection (LTBI) in humans confounds the treatment of tuberculosis. Consequently, there is a need to discover new therapeutic agents that can kill M. tuberculosis both during active disease and LTBI. The streptomycin-dependent strain of M. tuberculosis, 18b, provides a useful tool for this purpose since upon removal of streptomycin (STR) it enters a non-replicating state that mimics latency both in vitro and in animal models. The 4.41 Mb genome sequence of M. tuberculosis 18b was determined and this revealed the strain to belong to clade 3 of the ancient ancestral lineage of the Beijing family. STR-dependence was attributable to insertion of a single cytosine in the 530 loop of the 16S rRNA and to a single amino acid insertion in the N-terminal domain of initiation factor 3. RNA-seq was used to understand the genetic programme activated upon STR-withdrawal and hence to gain insight into LTBI. This revealed reconfiguration of gene expression and metabolic pathways showing strong similarities between non-replicating 18b and M. tuberculosis residing within macrophages, and with the core stationary phase and microaerophilic responses. The findings of this investigation confirm the validity of 18b as a model for LTBI, and provide insight into both the evolution of tubercle bacilli and the functioning of the ribosome.

Fragmentation of contaminant and endogenous DNA in ancient samples determined by shotgun sequencing; prospects for human palaeogenomics

Despite the successful retrieval of genomes from past remains, the prospects for human palaeogenomics remain unclear because of the difficulty of distinguishing contaminant from endogenous DNA sequences. Previous sequence data generated on high-throughput sequencing platforms indicate that fragmentation of ancient DNA sequences is a characteristic trait primarily arising due to depurination processes that create abasic sites leading to DNA breaks.

Hybridization Capture Using RAD Probes (hyRAD), a New Tool for Performing Genomic Analyses on Collection Specimens.

In the recent years, many protocols aimed at reproducibly sequencing reduced-genome subsets in non-model organisms have been published. Among them, RAD-sequencing is one of the most widely used. It relies on digesting DNA with specific restriction enzymes and performing size selection on the resulting fragments. Despite its acknowledged utility, this method is of limited use with degraded DNA samples, such as those isolated from museum specimens, as these samples are less likely to harbor fragments long enough to comprise two restriction sites making possible ligation of the adapter sequences (in the case of double-digest RAD) or performing size selection of the resulting fragments (in the case of single-digest RAD). Here, we address these limitations by presenting a novel method called hybridization RAD (hyRAD). In this approach, biotinylated RAD fragments, covering a random fraction of the genome, are used as baits for capturing homologous fragments from genomic shotgun sequencing libraries. This simple and cost-effective approach allows sequencing of orthologous loci even from highly degraded DNA samples, opening new avenues of research in the field of museum genomics. Not relying on the restriction site presence, it improves among-sample loci coverage. In a trial study, hyRAD allowed us to obtain a large set of orthologous loci from fresh and museum samples from a non-model butterfly species, with a high proportion of single nucleotide polymorphisms present in all eight analyzed specimens, including 58-year-old museum samples. The utility of the method was further validated using 49 museum and fresh samples of a Palearctic grasshopper species for which the spatial genetic structure was previously assessed using mtDNA amplicons. The application of the method is eventually discussed in a wider context. As it does not rely on the restriction site presence, it is therefore not sensitive to among-sample loci polymorphisms in the restriction sites that usually causes loci dropout. This should enable the application of hyRAD to analyses at broader evolutionary scales.

Identification of molecular subtypes and gene expression patterns of breast cancer analysing RNA-seq data

Breast cancer is the most common diagnosed cancer and the leading cause of cancer death among females worldwide. It is considered a highly heterogeneous disease and it must be classified into more homogeneous groups. Hence, the purpose of this study was to classify breast tumors based on variations in gene expression patterns derived from RNA sequencing by using different class discovery methods. 42 breast tumors paired-samples were sequenced by Illumine Genome Analyzer and the data was analyzed and prepared by TopHat2 and htseq-count. As reported previously, breast cancer could be grouped into five main groups known as basal epithelial-like group, HER2 group, normal breast-like group and two Luminal groups with a distinctive expression profile. Classifying breast tumor samples by using PAM50 method, the most common subtype was Luminal B and was significantly associated with ESR1 and ERBB2 high expression. Luminal A subtype had ESR1 and SLC39A6 significant high expression, whereas HER2 subtype had a high expression of ERBB2 and CNNE1 genes and low luminal epithelial gene expression. Basal-like and normal-like subtypes were associated with low expression of ESR1, PgR and HER2, and had significant high expression of cytokeratins 5 and 17. Our results were similar compared with TGCA breast cancer data results and with known studies related with breast cancer classification. Classifying breast tumors could add significant prognostic and predictive information to standard parameters, and moreover, identify marker genes for each subtype to find a better therapy for patients with breast cancer.

Genome-wide analysis of adaptive molecular evolution in the carnivorous plant Utricularia gibba

The genome of the bladderwort Utricularia gibba provides an unparalleled opportunity to uncover the adaptive landscape of an aquatic carnivorous plant with unique phenotypic features such as absence of roots, development of water-filled suction bladders, and a highly ramified branching pattern. Despite its tiny size, the U. gibba genome accommodates approximately as many genes as other plant genomes. To examine the relationship between the compactness of its genome and gene turnover, we compared the U. gibba genome with that of four other eudicot species, defining a total of 17,324 gene families (orthogroups). These families were further classified as either 1) lineage-specific expanded/contracted or 2) stable in size. The U. gibba-expanded families are generically related to three main phenotypic features: 1) trap physiology, 2) key plant morphogenetic/developmental pathways, and 3) response to environmental stimuli, including adaptations to life in aquatic environments. Further scans for signatures of protein functional specialization permitted identification of seven candidate genes with amino acid changes putatively fixed by positive Darwinian selection in the U. gibba lineage. The Arabidopsis orthologs of these genes (AXR, UMAMIT41, IGS, TAR2, SOL1, DEG9, and DEG10) are involved in diverse plant biological functions potentially relevant for U. gibba phenotypic diversification, including 1) auxin metabolism and signal transduction, 2) flowering induction and floral meristem transition, 3) root development, and 4) peptidases. Taken together, our results suggest numerous candidate genes and gene families as interesting targets for further experimental confirmation of their functional and adaptive roles in the U. gibba's unique lifestyle and highly specialized body plan.

Prototypical Arthropod Gene Content and Genome Organisation in the Centipede Strigamia maritima

Myriapods (e.g., centipedes and millipedes) display a simple homonomous body plan relative to other arthropods. All members of the class are terrestrial, but they attained terrestriality independently of insects. Myriapoda is the only arthropod class not represented by a sequenced genome. We present an analysis of the genome of the centipede Strigamia maritima. It retains a compact genome that has undergone less gene loss and shuffling than previously sequenced arthropods, and many orthologues of genes conserved from the bilaterian ancestor that have been lost in insects. Our analysis locates many genes in conserved macro-synteny contexts, and many small-scale examples of gene clustering. We describe several examples where S. maritima shows different solutions from insects to similar problems. The insect olfactory receptor gene family is absent from S. maritima, and olfaction in air is likely effected by expansion of other receptor gene families. For some genes S. maritima has evolved paralogues to generate coding sequence diversity, where insects use alternate splicing. This is most striking for the Dscam gene, which in Drosophila generates more than 100,000 alternate splice forms, but in S. maritima is encoded by over 100 paralogues. We see an intriguing linkage between the absence of any known photosensory proteins in a blind organism and the additional absence of canonical circadian clock genes. The phylogenetic position of myriapods allows us to identify where in arthropod phylogeny several particular molecular mechanisms and traits emerged. For example, we conclude that juvenile hormone signalling evolved with the emergence of the exoskeleton in the arthropods and that RR-1 containing cuticle proteins evolved in the lineage leading to Mandibulata. We also identify when various gene expansions and losses occurred. The genome of S. maritima offers us a unique glimpse into the ancestral arthropod genome, while also displaying many adaptations to its specific life history.

Novel Proteins Involved in Dynamics of The Photosynthetic Apparatus

During the past few years, a considerable number of research articles have been published relating to the structure and function of the major photosynthetic protein complexes, photosystem (PS) I, PSII, cytochrome (Cyt) b6f, and adenosine triphosphate (ATP) synthase. Sequencing of the Arabidopsis thaliana (Arabidopsis) genome together with several high-quality proteomics studies has, however, revealed that the thylakoid membrane network of plant chloroplasts still contains a number of functionally unknown proteins. These proteins may have a role as auxiliary proteins guiding the assembly, maintenance, and turnover of the thylakoid protein complexes, or they may be as yet unknown subunits of the photosynthetic complexes. Novel subunits are most likely to be found in the NAD(P)H dehydrogenase (NDH) complex, the structure and function of which have remained obscure in the absence of detailed crystallographic data, thus making this thylakoid protein complex a particularly interesting target of investigation. In this thesis, several novel thylakoid-associated proteins were identified by proteomics-based methods. The major goal of characterization of the stroma thylakoid associated polysome-nascent chain complexes was to determine the proteins that guide the dynamic life cycle of PSII. In addition, a large protein complex of ≥ 1,000 kDa, residing in the stroma thylakoid, was characterized in greater depth and it was found to be a supercomplex composed of the PSI and NDH complexes. A set of newly identified proteins from Arabidopsis thylakoids was subjected to detailed characterization using the reverse genetics approach and extensive biochemical and biophysical analysis. The role of the novel proteins, either as auxiliary proteins or subunits of the photosynthetic protein complexes, was revealed. Two novel thylakoid lumen proteins, TLP18.3 and AtCYP38, function as auxiliary proteins assisting specific steps of the assembly/repair of PSII. The role of the 10-kDa thylakoid lumen protein PsbR is related to the optimization of oxygen evolution of PSII by assisting the assembly of the PsbP protein. Two integral thylakoid membrane proteins, NDH45 and NDH48, are novel subunits of the chloroplast NDH complex. Finally, the thylakoid lumen immunophilin AtCYP20-2 is suggested to interact with the NDH complex, instead of PSII as was hypothesized earlier.

Analysis of the nucleotide sequence of the coat protein and 3'-untranslated region of two Brazilian Potato virus Y isolates

Two Brazilian Potato virus Y (PVY) isolates were biologically characterized as necrotic (PVY-NBR) and common (PVY-OBR) based upon symptoms on test plants. Additional characterization was performed by sequencing a cDNA corresponding to the 3' terminal region of the viral genome. The sequence consisted of 195 nucleotides (nt) coding part of the nuclear inclusion body b (NIb) gene, 804 nt of the coat protein (CP) gene, and 328 nt (PVY-OBR) or 326 nt (PVY-NBR) of the 3'-untranslated region (UTR). Translation of the sequence resulted in one single open reading frame with part of the NIb and a CP of 267 amino acids. The two isolates shared 95.1% similarity in the CP amino acid sequence. The CP and the 3'-UTR sequence of the Brazilian isolates were compared to those of other PVY isolates previously reported and unrooted phylogenetic trees were constructed. The trees revealed a separation of two distinct clusters, one comprising most of the common strains and the other comprising the necrotic strains. PVY-OBR was clustered in the common group and PVY-NBR in the necrotic one.

PCR amplification and sequence analyses of reverse transcriptase-like genes in Crinipellis perniciosa isolates

Reverse transcriptase (RT) sequence analysis is an important technique used to detect the presence of transposable elements in a genome. Putative RT sequences were analyzed in the genome of the pathogenic fungus C. perniciosa, the causal agent of witches' broom disease of cocoa. A 394 bp fragment was amplified from genomic DNA of different isolates of C. perniciosa belonging to C-, L-, and S-biotypes and collected from various geographical areas. The cleavage of PCR products with restriction enzymes and the sequencing of various RT fragments indicated the presence of several sequences showing transition events (G:C to A:T). Southern blot analysis revealed high copy numbers of RT signals, forming different patterns among C-, S-, and L-biotype isolates. Sequence comparisons of the predicted RT peptide indicate a close relationship with the RT protein from thegypsy family of LTR-retrotransposons. The possible role of these retrotransposons in generating genetic variability in the homothallic C. perniciosa is discussed.

Characterization of novel genes predominantly expressed in the epididymis – from genome analyses to genetically modified mice

In mammals, post-testicular sperm maturation taking place in the epididymis is required for the spermatozoa to acquire the abilities required to fertilize the egg in vivo. The epididymal epithelial cells secrete proteins and other small molecules into the lumen, where they interact with the spermatozoa and enable necessary maturational changes. In this study different in silico, in vitro and in vivo approaches were utilized in order to find novel genes responsible for the function of the epididymis and post-testicular sperm maturation in the mouse. Available online genomic databases were analyzed to identify genes potentially expressed in the epididymis, gene expression profiling was performed by studying their expression in different mouse tissues, and significance of certain genes to fertility was assessed by generating genetically modified mouse models. A recently discovered Pate (prostate and testis expression) gene family was found to be predominantly expressed in the epididymis. It represents one of the largest known gene families expressed in the epididymis, and the members code for proteins potentially involved in defense against microorganisms. Through genetically modified mouse models CRISP4 (cysteine-rich secretory protein 4) was identified to regulate sperm acrosome reaction, and BMYC to inhibit the expression of the Myc proto-oncogene in the developing testis. A mouse line expressing iCre recombinase specifically in the epididymis was also generated. This model can be used to generate conditional, epididymis-specific knock-out models, and will be a valuable tool in fertility studies.

Application of response surface methodology to study the biological removal of nitrogen from effluent of cattle slaughterhouse in a sequencing batch reactor

The aim of this study was to evaluate the efficiency of a sequencing batch reactor (SBR) on biological removal of nitrogen from cattle slaughterhouse wastewater by nitrification/denitrification processes. The effects of initial concentration of ammoniacal nitrogen were investigated at 100; 150 and 200 mg L-1 and air flow rate at 0.125; 0.375 and 0.625 L min¹ Lreactor-1 on the nitrogen compounds removal, by a Central Composite Rotational Design (CCRD) configuration. There were variations from 9.2 to 94.9%, 4.0 to 19.6% and 20.8 to 92.0% in the conversion of ammoniacal nitrogen to nitrate and nitrite concentration and removal of total nitrogen, respectively. The increase of air flow rate and decrease of the initial concentration of ammoniacal nitrogen resulted in higher efficiencies of total nitrogen removal, as well as the conversion of ammoniacal nitrogen to nitrate. During the pre-established intervals of this study, the removal and conversion efficiencies of nitrogen compounds above 85% were achieved in air flow rate variations from 0.375 to 0.725 L min-1 Lreactor-1 and initial concentration of ammoniacal nitrogen from 80 to 200 mg L-1. On denitrification process, we obtained efficiencies from 91.5 to 96.9% on the removal of nitrite/nitrate and from 78.3 to 87.9% on the removal of organic matter.

Effect of cycle time and airflow in biological nitrogen removal from poultry slaughterhouse wastewater using sequencing batch reactor

This study aimed to evaluate the influence of airflow (0.25, 0.50 and 0.75 L.L-1.min-1) and cycle time (10.45 h, 14.25 h and 17.35 h) on a sequencing batch reactor (SBR) performance in promoting nitrification and denitrification of poultry slaughterhouse wastewater. The operational stages included feeding, aerobic and anoxic reactions, sedimentation and discharge. SBR was operated in a laboratory scale with a working volume of 4 L, keeping 25% of biomass retained inside the reactor as inoculum for the next batch. In the anoxic stage, C: N ratio was maintained between 5 and 6 by adding cassava starch wastewater. A factorial design (22) with five repetitions was designed at the central point to evaluate the influence of cycle time and airflow on total inorganic nitrogen removal (N-NH4++N-NO2-+N-NO3-) and in the whole process (nitrification and denitrification). The highest total inorganic nitrogen removal (93.3%) was observed for airflow of 0.25 L.L-1.min‑1 and a cycle time of 14.25 h. At the end of the experiment, the sludge inside the reactor was characterized by fluorescent in situ hybridization (FISH), indicating the presence of ammonia and nitrite oxidizing bacteria.

Molecular diagnostics of rare inherited SYNDROMES with a view on diagnostic test development

CHARGE syndrome, Sotos syndrome and 3p deletion syndrome are examples of rare inherited syndromes that have been recognized for decades but for which the molecular diagnostics only have been made possible by recent advances in genomic research. Despite these advances, development of diagnostic tests for rare syndromes has been hindered by diagnostic laboratories having limited funds for test development, and their prioritization of tests for which a (relatively) high demand can be expected. In this study, the molecular diagnostic tests for CHARGE syndrome and Sotos syndrome were developed, resulting in their successful translation into routine diagnostic testing in the laboratory of Medical Genetics (UTUlab). In the CHARGE syndrome group, mutation was identified in 40.5% of the patients and in the Sotos syndrome group, in 34%, reflecting the use of the tests in routine diagnostics in differential diagnostics. In CHARGE syndrome, the low prevalence of structural aberrations was also confirmed. In 3p deletion syndrome, it was shown that small terminal deletions are not causative for the syndrome, and that testing with arraybased analysis provides a reliable estimate of the deletion size but benign copy number variants complicate result interpretation. During the development of the tests, it was discovered that finding an optimal molecular diagnostic strategy for a given syndrome is always a compromise between the sensitivity, specificity and feasibility of applying a new method. In addition, the clinical utility of the test should be considered prior to test development: sometimes a test performing well in a laboratory has limited utility for the patient, whereas a test performing poorly in the laboratory may have a great impact on the patient and their family. At present, the development of next generation sequencing methods is changing the concept of molecular diagnostics of rare diseases from single tests towards whole-genome analysis.

Genetic characterization of Brazilian bovine viral diarrhea virus isolates by partial nucleotide sequencing of the 5'-UTR region

Nineteen isolates of bovine viral diarrhea virus (BVDV) from Brazil were genetically characterized through partial nucleotide sequencing and analysis of the 5'UTR region. The isolates were grouped as BVDV-1 (11/19), BVDV-2 (6/19) or "atypical" pestivirus (2/19). Among the BVDV-1, eight isolates were classified as subgenotype BVDV-1a, whereas most (4 out of 6) BVDV-2 belonged to subgenotype 2b. Two isolates from aborted fetuses were not classified into any genetic group, being considered atypical BVDVs. Genetic diversity among Brazilian BVDV isolates may be responsible for vaccination and diag-nostic failure and therefore may influence the control strategies for BVDV infection in the country.

Sequencing and expression analysis of hepcidin mRNA in donkey (Equus asinus) liver

The hypoferremia that is observed during systemic inflammatory processes is mediated by hepcidin, which is a peptide that is mainly synthesized in the livers of several mammalian species. Hepcidin plays a key role in iron metabolism and in the innate immune system. It's up-regulation is particularly useful during acute inflammation, and it restricts the iron availability that is necessary for the growth of pathogenic microorganisms. In this study, the hepcidin mRNA of Equus asinus has been characterized, and the expression of donkey hepcidin in the liver has been determined. The donkey hepcidin sequence has an open reading frame (ORF) of 261 nucleotides, and the deduced corresponding protein sequence has 86 amino acids. The amino acid sequence of donkey hepcidin was most homologous to Equus caballus (98%). The mature donkey hepcidin sequence (25 amino acids) was 100% homologous to the equine mature hepcidin and has eight conserved cysteine residues that are found in all of the investigated hepcidin sequences. The expression profile of donkey hepcidin in the liver was high and was similar to the reference gene expression. The donkey hepcidin sequence was deposited in GenBankTM (HQ902884) and may be useful for additional studies on iron metabolism and the inflammatory process in this species.