Search

Study Question: Which transcriptomic alterations in mid-luteal endometrial scratch biopsies, taken prior to the assisted reproductive treatment (ART) treatment cycle are associated with unsuccessful pregnancy?, Summary Answer: Dysregulated interleukin-17 (IL-17) pathway components are demonstrated in women who fail to become pregnant after ART., What Is Known Already: Implantation failure is now recognised as a critical factor in unexplained infertility and may be an important component of failed ART., Study Design, Size, Duration: Using a prospective longitudinal study design, 29 nulliparous women with unexplained infertility undergoing ART were recruited between October 2016 and February 2018. Mid-luteal stage endometrium and matched serum samples were collected, and patients underwent a single embryo transfer in the subsequent cycle. RNA-seq analysis of endometrial biopsies was performed on the discovery cohort (n = 20)., Participants/materials, Setting, Methods: Gene set enrichment analysis of the differentially expressed genes (DEGs) was performed. Endometrium and serum were then prepared for IL-17A analysis by ELISA., Main Results and the Role of Chance: There were 204 differentially expressed protein-coding genes identified in tissue from women who became pregnant (n = 9) compared with tissue from women who failed to become pregnant (n = 11) (false discovery rate; P < 0.05). Of the 204 DEGs, 166 were decreased while 38 were increased in the pregnant compared to the non-pregnant groups. Gene set enrichment analysis of the DEGs identified an over-representation of IL-17 and Pl3K-Akt signalling pathways. All the DEGs within the IL-17 signalling pathway (MMP3, MMP1, IL1β, LCN2, S100A9 and FOSL1) demonstrated decreased expression in the pregnant group. Serum IL-17 protein levels were increased in the non-pregnant discovery cohort (n = 11) and these findings were confirmed a validation cohort (n = 9)., Limitations, Reasons for Caution: Limitations of our study include the cohort size and the lack of aneuploidy data for the embryos; however, all embryos transferred were single good or top-quality blastocysts., Wider Implications of the Findings: These findings demonstrate dysregulated IL-17 pathway components in women who fail to become pregnant after ART. Elevated serum levels of the pro-inflammatory cytokine IL-17 may predict failure of ART in women with unexplained infertility. Future trials of anti-IL-17 therapies in this cohort warrant further investigation., Study Funding/competing Interest(s): Funding from the UCD Wellcome Institutional Strategic Support Fund, which was financed jointly by University College Dublin and the SFI-HRB-Wellcome Biomedical Research Partnership (ref 204844/Z/16/Z), is acknowledged. The authors have no competing interests., Trial Registration Number: NA., (© The Author(s) 2020. Published by Oxford University Press on behalf of European Society of Human Reproduction and Embryology. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Background: Although many of the genic features in Mycobacterium abscessus have been fully validated, a comprehensive understanding of the regulatory elements remains lacking. Moreover, there is little understanding of how the organism regulates its transcriptomic profile, enabling cells to survive in hostile environments. Here, to computationally infer the gene regulatory network for Mycobacterium abscessus we propose a novel statistical computational modelling approach: BayesIan gene regulatory Networks inferreD via gene coExpression and compaRative genomics (BINDER). In tandem with derived experimental coexpression data, the property of genomic conservation is exploited to probabilistically infer a gene regulatory network in Mycobacterium abscessus.Inference on regulatory interactions is conducted by combining 'primary' and 'auxiliary' data strata. The data forming the primary and auxiliary strata are derived from RNA-seq experiments and sequence information in the primary organism Mycobacterium abscessus as well as ChIP-seq data extracted from a related proxy organism Mycobacterium tuberculosis. The primary and auxiliary data are combined in a hierarchical Bayesian framework, informing the apposite bivariate likelihood function and prior distributions respectively. The inferred relationships provide insight to regulon groupings in Mycobacterium abscessus., Results: We implement BINDER on data relating to a collection of 167,280 regulator-target pairs resulting in the identification of 54 regulator-target pairs, across 5 transcription factors, for which there is strong probability of regulatory interaction., Conclusions: The inferred regulatory interactions provide insight to, and a valuable resource for further studies of, transcriptional control in Mycobacterium abscessus, and in the family of Mycobacteriaceae more generally. Further, the developed BINDER framework has broad applicability, useable in settings where computational inference of a gene regulatory network requires integration of data sources derived from both the primary organism of interest and from related proxy organisms.

Background: Mycobacterium abscessus subsp. abscessus (MAB) is a highly drug resistant mycobacterium and the most common respiratory pathogen among the rapidly growing non-tuberculous mycobacteria. MAB is also one of the most deadly of the emerging cystic fibrosis (CF) pathogens requiring prolonged treatment with multiple antibiotics. In addition to its "mycobacterial" virulence genes, the genome of MAB harbours a large accessory genome, presumably acquired via lateral gene transfer including homologs shared with the CF pathogens Pseudomonas aeruginosa and Burkholderia cepacia. While multiple genome sequences are available there is little functional genomics data available for this important pathogen., Results: We report here the first multi-omics approach to characterize the primary transcriptome, coding potential and potential regulatory regions of the MAB genome utilizing differential RNA sequencing (dRNA-seq), RNA-seq, Ribosome profiling and LC-MS proteomics. In addition we attempt to address the genomes contribution to the molecular systems that underlie MAB's adaptation and persistence in the human host through an examination of MABs transcriptional response to a number of clinically relevant conditions. These include hypoxia, exposure to sub-inhibitory concentrations of antibiotics and growth in an artificial sputum designed to mimic the conditions within the cystic fibrosis lung., Conclusions: Our integrated map provides the first comprehensive view of the primary transcriptome of MAB and evidence for the translation of over one hundred new short open reading frames (sORFs). Our map will act as a resource for ongoing functional genomics characterization of MAB and our transcriptome data from clinically relevant stresses informs our understanding of MAB's adaptation to life in the CF lung. MAB's adaptation to growth in artificial CF sputum highlights shared metabolic strategies with other CF pathogens including P. aeruginosa and mirrors the transcriptional responses that lead to persistence in mycobacteria. These strategies include an increased reliance on amino acid metabolism, and fatty acid catabolism and highlights the relevance of the glyoxylate shunt to growth in the CF lung. Our data suggests that, similar to what is seen in chronically persisting P. aeruginosa, progression towards a biofilm mode of growth would play a more prominent role in a longer-term MAB infection. Finally, MAB's transcriptional response to antibiotics highlights the role of antibiotic modifications enzymes, active transport and the evolutionarily conserved WhiB7 driven antibiotic resistance regulon.

To compare gene expression among bovine tissues, large bovine RNA-seq datasets were used, comprising 280 samples from 10 different bovine tissues (uterine endometrium, granulosa cells, theca cells, cervix, embryos, leucocytes, liver, hypothalamus, pituitary, muscle) and generating 260 Gbases of data. Twin approaches were used: an information-theoretic analysis of the existing annotated transcriptome to identify the most tissue-specific genes and a de-novo transcriptome annotation to evaluate general features of the transcription landscape. Expression was detected for 97% of the Ensembl transcriptome with at least one read in one sample and between 28% and 66% at a level of 10 tags per million (TPM) or greater in individual tissues. Over 95% of genes exhibited some level of tissue-specific gene expression. This was mostly due to different levels of expression in different tissues rather than exclusive expression in a single tissue. Less than 1% of annotated genes exhibited a highly restricted tissue-specific expression profile and approximately 2% exhibited classic housekeeping profiles. In conclusion, it is the combined effects of the variable expression of large numbers of genes (73%-93% of the genome) and the specific expression of a small number of genes (<1% of the transcriptome) that contribute to determining the outcome of the function of individual tissues.

Background: Mycobacterium abscessus (MAB) is an emerging pathogen causing pulmonary infections in those with inflammatory lung disorders, such as Cystic Fibrosis (CF), and is associated with the highest fatality rate among rapidly growing mycobacteria (RGM). Phenotypically, MAB manifests as either a Smooth (MAB-S) or a Rough (MAB-R) morphotype, which differ in their levels of cell wall glycopeptidolipids (GPLs) and in their pathogenicity in vivo. As one of the primary immune cells encountered by MAB, we sought to examine the early transcriptional events within macrophages, following infection with both MAB-S or MAB-R., Results: We sampled the transcriptomes (mRNA and miRNA) of THP-1-derived macrophages infected with both MAB-R and MAB-S at 1, 4 and 24 h post-infection (hpi) using RNA-seq. A core set of 606 genes showed consistent expression profiles in response to both morphotypes, corresponding to the early transcriptional response to MAB. The core response is type I Interferon (IFN)-driven, involving the NF-κB and MAPK signaling pathways with concomitant pro-inflammatory cytokine production, and network analysis identified STAT1, EGR1, and SRC as key hub and bottleneck genes. MAB-S elicited a more robust transcriptional response at both the mRNA and miRNA levels, which was reflected in higher cytokine levels in culture supernatants. The transcriptional profiles of macrophages infected with both morphotypes were highly correlated, however, and a direct comparison identified few genes to distinguish them. Most of the induced miRNAs have previously been associated with mycobacterial infection and overall miRNA expression patterns were similarly highly correlated between the morphotypes., Conclusions: The report here details the first whole transcriptome analysis of the early macrophage response to MAB infection. The overall picture at the early stages of macrophage infection is similar to that of other mycobacteria, reflected in a core type I IFN and pro-inflammatory cytokine response. Large-scale transcriptional differences in the host response to the different MAB morphotypes are not evident in the early stages of infection, however the subset of genes with distinct expression profiles suggest potentially interesting differences in internal trafficking of MAB within macrophages.

Background: Domestication of the now-extinct wild aurochs, Bos primigenius, gave rise to the two major domestic extant cattle taxa, B. taurus and B. indicus. While previous genetic studies have shed some light on the evolutionary relationships between European aurochs and modern cattle, important questions remain unanswered, including the phylogenetic status of aurochs, whether gene flow from aurochs into early domestic populations occurred, and which genomic regions were subject to selection processes during and after domestication. Here, we address these questions using whole-genome sequencing data generated from an approximately 6,750-year-old British aurochs bone and genome sequence data from 81 additional cattle plus genome-wide single nucleotide polymorphism data from a diverse panel of 1,225 modern animals., Results: Phylogenomic analyses place the aurochs as a distinct outgroup to the domestic B. taurus lineage, supporting the predominant Near Eastern origin of European cattle. Conversely, traditional British and Irish breeds share more genetic variants with this aurochs specimen than other European populations, supporting localized gene flow from aurochs into the ancestors of modern British and Irish cattle, perhaps through purposeful restocking by early herders in Britain. Finally, the functions of genes showing evidence for positive selection in B. taurus are enriched for neurobiology, growth, metabolism and immunobiology, suggesting that these biological processes have been important in the domestication of cattle., Conclusions: This work provides important new information regarding the origins and functional evolution of modern cattle, revealing that the interface between early European domestic populations and wild aurochs was significantly more complex than previously thought.

Related species are often used to understand the molecular underpinning of virulence through examination of a shared set of biological features attributable to a core genome of orthologous genes. An important but insufficiently studied issue, however, is the extent to which the regulatory architectures are similarly conserved. A small number of studies have compared the primary transcriptomes of different bacterial species, but few have compared closely related species with clearly divergent evolutionary histories. We addressed the impact of differing modes of evolution within the genus Mycobacterium through comparison of the primary transcriptome of M. marinum with that of a closely related lineage, M. bovis. Both are thought to have evolved from an ancestral generalist species, with M. bovis and other members of the M. tuberculosis complex having subsequently undergone downsizing of their genomes during the transition to obligate pathogenicity. M. marinum, in contrast, has retained a large genome, appropriate for an environmental organism, and is a broad-host-range pathogen. We also examined changes over a shorter evolutionary time period through comparison of the primary transcriptome of M. bovis with that of another member of the M. tuberculosis complex (M. tuberculosis) which possesses an almost identical genome but maintains a distinct host preference. Importance: Our comparison of the transcriptional start site (TSS) maps of M. marinum and M. bovis uncovers a pillar of conserved promoters, noncoding RNA (NCRNA), and a genome-wide signal in the -35 promoter regions of both species. We identify evolutionarily conserved transcriptional attenuation and highlight its potential contribution to multidrug resistance mediated through the transcriptional regulator whiB7. We show that a species population history is reflected in its transcriptome and posit relaxed selection as the main driver of an abundance of canonical -10 promoter sites in M. bovis relative to M. marinum. It appears that transcriptome composition in mycobacteria is driven primarily by the availability of such sites and that their frequencies diverge significantly across the mycobacterial clade. Finally, through comparison of M. bovis and M. tuberculosis, we illustrate that single nucleotide polymorphism (SNP)-driven promoter differences likely underpin many of the transcriptional differences between M. tuberculosis complex lineages., (Copyright © 2014 Dinan et al.)

The rise and spread of antibiotic resistance is among the most severe challenges facing modern medicine. Despite this fact, attempts to develop novel classes of antibiotic have been largely unsuccessful. The traditional mechanisms by which antibiotics work are subject to relatively rapid bacterial resistance via mutation, and hence have a limited period of efficacy. One promising strategy to ameliorate this problem is to shift from the use of chemical compounds targeting protein structures and processes to a new era of RNA-based therapeutics. RNA-mediated regulation (riboregulation) has evolved naturally in bacteria and is therefore a highly efficient means by which gene expression can be manipulated. Here, we describe recent advances toward the development of effective anti-bacterial therapies, which operate through various strategies centered on RNA.

Background: Mycobacterium bovis, the causative agent of bovine tuberculosis, is an intracellular pathogen that can persist inside host macrophages during infection via a diverse range of mechanisms that subvert the host immune response. In the current study, we have analysed and compared the transcriptomes of M. bovis-infected monocyte-derived macrophages (MDM) purified from six Holstein-Friesian females with the transcriptomes of non-infected control MDM from the same animals over a 24 h period using strand-specific RNA sequencing (RNA-seq). In addition, we compare gene expression profiles generated using RNA-seq with those previously generated by us using the high-density Affymetrix® GeneChip® Bovine Genome Array platform from the same MDM-extracted RNA., Results: A mean of 7.2 million reads from each MDM sample mapped uniquely and unambiguously to single Bos taurus reference genome locations. Analysis of these mapped reads showed 2,584 genes (1,392 upregulated; 1,192 downregulated) and 757 putative natural antisense transcripts (558 upregulated; 119 downregulated) that were differentially expressed based on sense and antisense strand data, respectively (adjusted P-value ≤ 0.05). Of the differentially expressed genes, 694 were common to both the sense and antisense data sets, with the direction of expression (i.e. up- or downregulation) positively correlated for 693 genes and negatively correlated for the remaining gene. Gene ontology analysis of the differentially expressed genes revealed an enrichment of immune, apoptotic and cell signalling genes. Notably, the number of differentially expressed genes identified from RNA-seq sense strand analysis was greater than the number of differentially expressed genes detected from microarray analysis (2,584 genes versus 2,015 genes). Furthermore, our data reveal a greater dynamic range in the detection and quantification of gene transcripts for RNA-seq compared to microarray technology., Conclusions: This study highlights the value of RNA-seq in identifying novel immunomodulatory mechanisms that underlie host-mycobacterial pathogen interactions during infection, including possible complex post-transcriptional regulation of host gene expression involving antisense RNA.

Background: The Amoebozoa constitute one of the primary divisions of eukaryotes, encompassing taxa of both biomedical and evolutionary importance, yet its genomic diversity remains largely unsampled. Here we present an analysis of a whole genome assembly of Acanthamoeba castellanii (Ac) the first representative from a solitary free-living amoebozoan., Results: Ac encodes 15,455 compact intron-rich genes, a significant number of which are predicted to have arisen through inter-kingdom lateral gene transfer (LGT). A majority of the LGT candidates have undergone a substantial degree of intronization and Ac appears to have incorporated them into established transcriptional programs. Ac manifests a complex signaling and cell communication repertoire, including a complete tyrosine kinase signaling toolkit and a comparable diversity of predicted extracellular receptors to that found in the facultatively multicellular dictyostelids. An important environmental host of a diverse range of bacteria and viruses, Ac utilizes a diverse repertoire of predicted pattern recognition receptors, many with predicted orthologous functions in the innate immune systems of higher organisms., Conclusions: Our analysis highlights the important role of LGT in the biology of Ac and in the diversification of microbial eukaryotes. The early evolution of a key signaling facility implicated in the evolution of metazoan multicellularity strongly argues for its emergence early in the Unikont lineage. Overall, the availability of an Ac genome should aid in deciphering the biology of the Amoebozoa and facilitate functional genomic studies in this important model organism and environmental host.

Mucus within the cervical canal represents a hormonally regulated barrier that reconciles the need to exclude the vaginal microflora from the uterus during progesterone dominance, while permitting sperm transport at estrus. Its characteristics change during the estrous cycle to facilitate these competing functional requirements. Hydrated mucin glycoproteins synthesized by the endocervical epithelium form the molecular scaffold of this mucus. This study uses the bovine cervix as a model to examine functional groups of genes related to mucin biosynthesis and mucus production over the periestrous period when functional changes in cervical barrier function are most prominent. Cervical tissue samples were collected from 30 estrus synchronized beef heifers. Animals were slaughtered in groups starting 12 h after the withdrawal of intravaginal progesterone releasing devices (controlled internal drug releases) until 7 days postonset of estrus (luteal phase). Subsequent groupings represented proestrus, early estrus, late estrus, metestrus, and finally the early luteal phase. Tissues were submitted to next generation RNA-seq transcriptome analysis. We identified 114 genes associated with biosynthesis and intracellular transport of mucins, and postsecretory modifications of cervical; 53 of these genes showed at least a twofold change in one or more experimental group in relation to onset of estrus, and the differences between groups were significant (P < 0.05). The majority of these genes showed the greatest alteration in their expression in the 48 h postestrus and luteal phase groups.

Polycomb group proteins are repressive chromatin modifiers with essential roles in metazoan development, cellular differentiation and cell fate maintenance. How Polycomb proteins access active chromatin to confer transcriptional silencing during lineage transitions remains unclear. Here we show that the Polycomb repressive complex 2 (PRC2) component PHF19 binds trimethylated histone H3 Lys36 (H3K36me3), a mark of active chromatin, via its Tudor domain. PHF19 associates with the H3K36me3 demethylase NO66, and it is required to recruit the PRC2 complex and NO66 to stem cell genes during differentiation, leading to PRC2-mediated trimethylation of histone H3 Lys27 (H3K27), loss of H3K36me3 and transcriptional silencing. We propose a model whereby PHF19 functions during mouse embryonic stem cell differentiation to transiently bind the H3K36me3 mark via its Tudor domain, forming essential contact points that allow recruitment of PRC2 and H3K36me3 demethylase activity to active gene loci during their transition to a Polycomb-repressed state.

The aims of this study were to 1) identify the earliest transcriptional response of the bovine endometrium to the presence of the conceptus (using RNAseq), 2) investigate if these genes are regulated by interferon tau (IFNT) in vivo, and 3) determine if they are predictive of the pregnancy status of postpartum dairy cows. RNAseq identified 459 differentially expressed genes (DEGs) between pregnant and cyclic endometria on day 16. Quantitative real-time PCR analysis of selected genes revealed PARP12, ZNFX1, HERC6, IFI16, RNF213, and DDX58 expression increased in pregnant compared with cyclic endometria on day 16 and were directly upregulated by intrauterine infusion of IFNT in vivo for 2 h (P < 0.05). On day 13 following estrous endometrial expression of nine genes increased [ARHGAP1, MGC127874, LIMS2, TBC1D1, FBXL7, C25H16orf71, LOC507810, ZSWIM4, and one novel gene (ENSBTAT00000050193)] and seven genes decreased (SERBP1, SRGAP2, AL7A1, TBK1, F2RL2, MGC128929, and WBSCR17; P < 0.05) in pregnant compared with cyclic heifers. Of these DEGs, significant differences in expression between pregnant and cyclic endometria were maintained on day 16 for F2RL2, LIMS2, LOC507810, MGC127874, TBC1D1, WBSCR17, and ZSWIM4 (P < 0.05) both their expression was not directly regulated by IFNT in vivo. Analysis of the expression of selected interferon-stimulated genes in blood samples from postpartum dairy cows revealed a significant increase (P < 0.05) in expression of ZXFX1, PARP12, SAMD9, and HERC6 on day 18 following artificial insemination in cows subsequently confirmed pregnant compared with cyclic controls. In conclusion, RNAseq identified a number of novel pregnancy-associated genes in the endometrium of cattle during early pregnancy that are not regulated by IFNT in vivo. In addition, a number of genes that are directly regulated by short term exposure to IFNT in vivo are differentially expressed on day 18 following estrus detection in the blood of postpartum dairy cows depending on their pregnancy status.

Clinical isolates of Staphylococcus aureus can express biofilm phenotypes promoted by the major cell wall autolysin and the fibronectin-binding proteins or the icaADBC-encoded polysaccharide intercellular adhesin/poly-N-acetylglucosamine (PIA/PNAG). Biofilm production in methicillin-susceptible S. aureus (MSSA) strains is typically dependent on PIA/PNAG whereas methicillin-resistant isolates express an Atl/FnBP-mediated biofilm phenotype suggesting a relationship between susceptibility to β-lactam antibiotics and biofilm. By introducing the methicillin resistance gene mecA into the PNAG-producing laboratory strain 8325-4 we generated a heterogeneously resistant (HeR) strain, from which a homogeneous, high-level resistant (HoR) derivative was isolated following exposure to oxacillin. The HoR phenotype was associated with a R₆₀₂H substitution in the DHHA1 domain of GdpP, a recently identified c-di-AMP phosphodiesterase with roles in resistance/tolerance to β-lactam antibiotics and cell envelope stress. Transcription of icaADBC and PNAG production were impaired in the 8325-4 HoR derivative, which instead produced a proteinaceous biofilm that was significantly inhibited by antibodies against the mecA-encoded penicillin binding protein 2a (PBP2a). Conversely excision of the SCCmec element in the MRSA strain BH1CC resulted in oxacillin susceptibility and reduced biofilm production, both of which were complemented by mecA alone. Transcriptional activity of the accessory gene regulator locus was also repressed in the 8325-4 HoR strain, which in turn was accompanied by reduced protease production and significantly reduced virulence in a mouse model of device infection. Thus, homogeneous methicillin resistance has the potential to affect agr- and icaADBC-mediated phenotypes, including altered biofilm expression and virulence, which together are consistent with the adaptation of healthcare-associated MRSA strains to the antibiotic-rich hospital environment in which they are frequently responsible for device-related infections in immuno-compromised patients.

Second generation sequencing has prompted a number of groups to re-interrogate the transcriptomes of several bacterial and archaeal species. One of the central findings has been the identification of complex networks of small non-coding RNAs that play central roles in transcriptional regulation in all growth conditions and for the pathogen's interaction with and survival within host cells. Legionella pneumophila is a gram-negative facultative intracellular human pathogen with a distinct biphasic lifestyle. One of its primary environmental hosts in the free-living amoeba Acanthamoeba castellanii and its infection by L. pneumophila mimics that seen in human macrophages. Here we present analysis of strand specific sequencing of the transcriptional response of L. pneumophila during exponential and post-exponential broth growth and during the replicative and transmissive phase of infection inside A. castellanii. We extend previous microarray based studies as well as uncovering evidence of a complex regulatory architecture underpinned by numerous non-coding RNAs. Over seventy new non-coding RNAs could be identified; many of them appear to be strain specific and in configurations not previously reported. We discover a family of non-coding RNAs preferentially expressed during infection conditions and identify a second copy of 6S RNA in L. pneumophila. We show that the newly discovered putative 6S RNA as well as a number of other non-coding RNAs show evidence for antisense transcription. The nature and extent of the non-coding RNAs and their expression patterns suggests that these may well play central roles in the regulation of Legionella spp. specific traits and offer clues as to how L. pneumophila adapts to its intracellular niche. The expression profiles outlined in the study have been deposited into Genbank's Gene Expression Omnibus (GEO) database under the series accession GSE27232.

Background: Digital gene expression profiling was used to characterize the assembly of genes expressed in equine skeletal muscle and to identify the subset of genes that were differentially expressed following a ten-month period of exercise training. The study cohort comprised seven Thoroughbred racehorses from a single training yard. Skeletal muscle biopsies were collected at rest from the gluteus medius at two time points: T(1) - untrained, (9 +/- 0.5 months old) and T(2) - trained (20 +/- 0.7 months old)., Results: The most abundant mRNA transcripts in the muscle transcriptome were those involved in muscle contraction, aerobic respiration and mitochondrial function. A previously unreported over-representation of genes related to RNA processing, the stress response and proteolysis was observed. Following training 92 tags were differentially expressed of which 74 were annotated. Sixteen genes showed increased expression, including the mitochondrial genes ACADVL, MRPS21 and SLC25A29 encoded by the nuclear genome. Among the 58 genes with decreased expression, MSTN, a negative regulator of muscle growth, had the greatest decrease.Functional analysis of all expressed genes using FatiScan revealed an asymmetric distribution of 482 Gene Ontology (GO) groups and 18 KEGG pathways. Functional groups displaying highly significant (P < 0.0001) increased expression included mitochondrion, oxidative phosphorylation and fatty acid metabolism while functional groups with decreased expression were mainly associated with structural genes and included the sarcoplasm, laminin complex and cytoskeleton., Conclusion: Exercise training in Thoroughbred racehorses results in coordinate changes in the gene expression of functional groups of genes related to metabolism, oxidative phosphorylation and muscle structure.

Background: The derivation of domestic cattle from the extinct wild aurochs (Bos primigenius) has been well-documented by archaeological and genetic studies. Genetic studies point towards the Neolithic Near East as the centre of origin for Bos taurus, with some lines of evidence suggesting possible, albeit rare, genetic contributions from locally domesticated wild aurochsen across Eurasia. Inferences from these investigations have been based largely on the analysis of partial mitochondrial DNA sequences generated from modern animals, with limited sequence data from ancient aurochsen samples. Recent developments in DNA sequencing technologies, however, are affording new opportunities for the examination of genetic material retrieved from extinct species, providing new insight into their evolutionary history. Here we present DNA sequence analysis of the first complete mitochondrial genome (16,338 base pairs) from an archaeologically-verified and exceptionally-well preserved aurochs bone sample., Methodology: DNA extracts were generated from an aurochs humerus bone sample recovered from a cave site located in Derbyshire, England and radiocarbon-dated to 6,738+/-68 calibrated years before present. These extracts were prepared for both Sanger and next generation DNA sequencing technologies (Illumina Genome Analyzer). In total, 289.9 megabases (22.48%) of the post-filtered DNA sequences generated using the Illumina Genome Analyzer from this sample mapped with confidence to the bovine genome. A consensus B. primigenius mitochondrial genome sequence was constructed and was analysed alongside all available complete bovine mitochondrial genome sequences., Conclusions: For all nucleotide positions where both Sanger and Illumina Genome Analyzer sequencing methods gave high-confidence calls, no discrepancies were observed. Sequence analysis reveals evidence of heteroplasmy in this sample and places this mitochondrial genome sequence securely within a previously identified aurochsen haplogroup (haplogroup P), thus providing novel insights into pre-domestic patterns of variation. The high proportion of authentic, endogenous aurochs DNA preserved in this sample bodes well for future efforts to determine the complete genome sequence of a wild ancestor of domestic cattle.

In this report, we describe the glutathione transferase (GST) gene family in the dengue vector Aedes aegypti and suggest a novel role for a new class of mosquito GSTs. Twenty-six GST genes are present in Ae. aegypti, two of which are alternatively spliced to give a total of 29 transcripts for cytosolic GSTs. The six classes identified in other insect species are all represented and, as in Anopheles gambiae, the majority of the mosquito GSTs belong to the insect-specific Delta and Epsilon classes with eight members each. Sixteen secure 1:1 orthologs were identified between GSTs in Ae. aegypti and An. gambiae, but only four of these have recognisable orthologs in Drosophila melanogaster. Three mosquito-specific GSTs were identified which did not belong to any previously recognised GST classes. One of these, GSTx2, has been previously implicated in conferring 1,1,1-trichloro-2,2-bis-(p-chlorophenyl)ethane (DDT) resistance in Ae. aegypti from South America. However, we found no evidence for increased levels of this GST protein in DDT/pyrethroid-resistant populations from Thailand. Furthermore, we show that the recombinant GSTX2-2 protein is unable to metabolise DDT. Interestingly, GSTX2-2 showed an affinity for hematin, and this, together with the restricted distribution of this class to haematophagous insects, may indicate a role for these enzymes in protecting mosquitoes against heme toxicity during blood feeding.

Background: Aedes aegypti is the principal vector of yellow fever and dengue viruses throughout the tropical world. To provide a set of manually curated and annotated sequences from the Ae. aegypti genome, 14 mapped bacterial artificial chromosome (BAC) clones encompassing 1.57 Mb were sequenced, assembled and manually annotated using a combination of computational gene-finding, expressed sequence tag (EST) matches and comparative protein homology. PCR and sequencing were used to experimentally confirm expression and sequence of a subset of these transcripts., Results: Of the 51 manual annotations, 50 and 43 demonstrated a high level of similarity to Anopheles gambiae and Drosophila melanogaster genes, respectively. Ten of the 12 BAC sequences with more than one annotated gene exhibited synteny with the A. gambiae genome. Putative transcripts from eight BAC clones were found in multiple copies (two copies in most cases) in the Aedes genome assembly, which point to the probable presence of haplotype polymorphisms and/or misassemblies., Conclusion: This study not only provides a benchmark set of manually annotated transcripts for this genome that can be used to assess the quality of the auto-annotation pipeline and the assembly, but it also looks at the effect of a high repeat content on the genome assembly and annotation pipeline.

The genome sequence of the protistan parasite Entamoeba histolytica HM-1:IMSS has been completed recently. Among the findings has been a unique organisation for the tRNA genes in this organism. Forty-two of the tRNA isoacceptor types are encoded in tandem arrays that vary in unit length from 490 to 1775 basepairs and contain from 1 to 5 tRNA genes. In three cases a 5S RNA gene is also present in the unit. An estimated 10% of the genome is made up of these arrays. Interspersed between RNA-encoding sequences are short tandem repeats that are polymorphic between isolates and, in some cases, within isolates. The number and organisation of tRNA genes in E. histolytica is unprecedented. In addition to encoding the tRNAs of the organism we propose that the arrays may fulfil a structural role in the genome.

Acanthamoeba castellanii is a free-living amoeba found in soil, freshwater, and marine environments and an important predator of bacteria. Acanthamoeba castellanii is also an opportunistic pathogen of clinical interest, responsible for several distinct diseases in humans. In order to provide a genomic platform for the study of this ubiquitous and important protist, we generated a sequence survey of approximately 0.5 x coverage of the genome. The data predict that A. castellanii exhibits a greater biosynthetic capacity than the free-living Dictyostelium discoideum and the parasite Entamoeba histolytica, providing an explanation for the ability of A. castellanii to inhabit a diversity of environments. Alginate lyase may provide access to bacteria within biofilms by breaking down the biofilm matrix, and polyhydroxybutyrate depolymerase may facilitate utilization of the bacterial storage compound polyhydroxybutyrate as a food source. Enzymes for the synthesis and breakdown of cellulose were identified, and they likely participate in encystation and excystation as in D. discoideum. Trehalose-6-phosphate synthase is present, suggesting that trehalose plays a role in stress adaptation. Detection and response to a number of stress conditions is likely accomplished with a large set of signal transduction histidine kinases and a set of putative receptor serine/threonine kinases similar to those found in E. histolytica. Serine, cysteine and metalloproteases were identified, some of which are likely involved in pathogenicity.

The sequencing of the genome of Entamoeba histolytica has allowed a reconstruction of its metabolic pathways, many of which are unusual for a eukaryote. Based on the genome sequence, it appears that amino acids may play a larger role than previously thought in energy metabolism, with roles in both ATP synthesis and NAD regeneration. Arginine decarboxylase may be involved in survival of E. histolytica during its passage through the stomach. The usual pyrimidine synthesis pathway is absent, but a partial pyrimidine degradation pathway could be part of a novel pyrimidine synthesis pathway. Ribonucleotide reductase was not found in the E. histolytica genome, but it was found in the close relatives Entamoeba invadens and Entamoeba moshkovskii, suggesting a recent loss from E. histolytica. The usual eukaryotic glucose transporters are not present, but members of a prokaryotic monosaccharide transporter family are present.

Rab proteins are ubiquitous small GTP-binding proteins that form a highly conserved family and regulate vesicular trafficking. Recent completion of the genome of the enteric protozoan parasite Entamoeba histolytica enabled us to identify an extremely large number (>90) of putative Rab genes. Multiple alignment and phylogenic analysis of amebic, human, and yeast Rab showed that only 22 amebic Rab proteins including EhRab1, EhRab2, EhRab5, EhRab7, EhRab8, EhRab11, and EhRab21 showed significant similarity to Rab from other organisms. The 69 remaining amebic Rab proteins showed only moderate similarity (<40% identity) to Rab proteins from other organisms. Approximately one-third of Rab proteins including Rab7, Rab11, and RabC form 15 subfamilies, which contain up to nine isoforms. Approximately 70% of amebic Rab genes contain single or multiple introns, and this proportion is significantly higher than that of common genes in this organism. Twenty-five Rabs possess an atypical carboxyl terminus such as CXXX, XCXX, XXCX, XXXC, and no cysteine. We propose annotation of amebic Rab genes and discuss biological significance of this extraordinary diversity of EhRab proteins in this organism.

Cryptococcus neoformans is a basidiomycetous yeast ubiquitous in the environment, a model for fungal pathogenesis, and an opportunistic human pathogen of global importance. We have sequenced its approximately 20-megabase genome, which contains approximately 6500 intron-rich gene structures and encodes a transcriptome abundant in alternatively spliced and antisense messages. The genome is rich in transposons, many of which cluster at candidate centromeric regions. The presence of these transposons may drive karyotype instability and phenotypic variation. C. neoformans encodes unique genes that may contribute to its unusual virulence properties, and comparison of two phenotypically distinct strains reveals variation in gene content in addition to sequence polymorphisms between the genomes.

The mosquito, Aedes aegypti, is the primary, worldwide arthropod vector for the yellow fever and dengue viruses. As it is also one of the most tractable mosquito species for laboratory studies, it has been and remains one of the most intensively studied arthropod species. This has resulted in the development of detailed genetic and physical maps for Ae. aegypti and considerable insight into its genome organization. The research community is well-advanced in developing important molecular tools that will facilitate a whole genome sequencing effort. This includes generation of BAC clone end sequences, physical mapping of selected BAC clones and generation of EST sequences. Whole genome sequence information for Ae. aegypti will provide important insight into mosquito chromosome evolution and allow for the identification of genes and gene function. These functions may be common to all mosquitoes or perhaps unique to individual species, possibly specific to host-seeking and blood-feeding behaviors, as well as the innate immune response to pathogens encountered during blood-feeding. This information will be invaluable to the global effort to develop novel strategies for preventing arthropod-borne disease transmission.

Cysteine proteases are known to be important pathogenicity factors of the protozoan parasite Entamoeba histolytica. So far, a total of eight genes coding for cysteine proteases have been identified in E. histolytica, two of which are absent in the closely related nonpathogenic species E. dispar. However, present knowledge is restricted to enzymes expressed during in vitro cultivation of the parasite, which might represent only a subset of the entire repertoire. Taking advantage of the current E. histolytica genome-sequencing efforts, we analyzed databases containing more than 99% of all ameba gene sequences for the presence of cysteine protease genes. A total of 20 full-length genes was identified (including all eight genes previously reported), which show 10 to 86% sequence identity. The various genes obviously originated from two separate ancestors since they form two distinct clades. Despite cathepsin B-like substrate specificities, all of the ameba polypeptides are structurally related to cathepsin L-like enzymes. None of the previously described enzymes but 7 of the 12 newly identified proteins are unique compared to cathepsins of higher eukaryotes in that they are predicted to have transmembrane or glycosylphosphatidylinositol anchor attachment domains. Southern blot analysis revealed that orthologous sequences for all of the newly identified proteases are present in E. dispar. Interestingly, the majority of the various cysteine protease genes are not expressed in E. histolytica or E. dispar trophozoites during in vitro cultivation. Therefore, it is likely that at least some of these enzymes are required for infection of the human host and/or for completion of the parasite life cycle.

A Bacterial Artificial Chromosome (BAC) genomic DNA library of Anopheles gambiae, the major human malaria vector in sub-Saharan Africa, was constructed and characterized. This library (ND-TAM) is composed of 30,720 BAC clones in eighty 384-well plates. The estimated average insert size of the library is 133 kb, with an overall genome coverage of approximately 14-fold. The ends of approximately two-thirds of the clones in the library were sequenced, yielding 32,340 pair-mate ends. A statistical analysis (G-test) of the results of PCR screening of the library indicated a random distribution of BACs in the genome, although one gap encompassing the white locus on the X-chromosome was identified. Furthermore, combined with another previously constructed BAC library (ND-1), ~2,000 BACs have been physically mapped by polytene chromosomal in situ hybridization. These BAC end pair mates and physically mapped BACs have been useful for both the assembly of a fully sequenced A. gambiae genome and for linking the assembled sequence to the three polytene chromosomes. This ND-TAM library is now publicly available at both http://www.malaria.mr4.org/mr4pages/index.html/ and http://hbz.tamu.edu/, providing a valuable resource to the mosquito research community.

Entamoeba histolytica and Spironucleus barkhanus have genes that encode short iron-dependent hydrogenases (Fe-hydrogenases), even though these protists lack hydrogenosomes. To understand better the biochemistry of the protist Fe-hydrogenases, we prepared a recombinant E. histolytica short Fe-hydrogenase and measured its activity in vitro. A Giardia lamblia gene encoding a short Fe-hydrogenase was identified from shotgun genomic sequences, and RT-PCR showed that cultured entamoebas and giardias transcribe short Fe-hydrogenase mRNAs. A second E. histolytica gene, which encoded a long Fe-hydrogenase, was identified from shotgun genomic sequences. Phylogenetic analyses suggested that the short Fe-hydrogenase genes of entamoeba and diplomonads share a common ancestor, while the long Fe-hydrogenase gene of entamoeba appears to have been laterally transferred from a bacterium. These results are discussed in the context of competing ideas for the origins of genes encoding fermentation enzymes of these protists.

Giardia lamblia and Entamoeba histolytica are amitochondriate, microaerophilic protists which use fermentation enzymes like those of bacteria to survive anaerobic conditions within the intestinal lumen. Genes encoding fermentation enzymes and related electron transport peptides (e.g., ferredoxins) in giardia organisms and amebae are hypothesized to be derived from either an ancient anaerobic eukaryote (amitochondriate fossil hypothesis), a mitochondrial endosymbiont (hydrogen hypothesis), or anaerobic bacteria (lateral transfer hypothesis). The goals here were to complete the molecular characterization of giardial and amebic fermentation enzymes and to determine the origins of the genes encoding them, when possible. A putative giardia [2Fe-2S]ferredoxin which had a hypothetical organelle-targeting sequence at its N terminus showed similarity to mitochondrial ferredoxins and the hydrogenosomal ferredoxin of Trichomonas vaginalis (another luminal protist). However, phylogenetic trees were star shaped, with weak bootstrap support, so we were unable to confirm or rule out the endosymbiotic origin of the giardia [2Fe-2S]ferredoxin gene. Putative giardial and amebic 6-kDa ferredoxins, ferredoxin-nitroreductase fusion proteins, and oxygen-insensitive nitroreductases each tentatively supported the lateral transfer hypothesis. Although there were not enough sequences to perform meaningful phylogenetic analyses, the unique common occurrence of these peptides and enzymes in giardia organisms, amebae, and the few anaerobic prokaryotes suggests the possibility of lateral transfer. In contrast, there was more robust phylogenetic evidence for the lateral transfer of G. lamblia genes encoding an NADH oxidase from a gram-positive coccus and a microbial group 3 alcohol dehydrogenase from thermoanaerobic prokaryotes. In further support of lateral transfer, the G. lamblia NADH oxidase and adh3 genes appeared to have an evolutionary history distinct from those of E. histolytica.

Short introns occur in numerous protist lineages, but there are no reports of intervening sequences in the protists Giardia lamblia and Trichomonas vaginalis, which may represent the deepest known branches in the eukaryotic line of descent. We have discovered a 35-bp spliceosomal intron in a gene encoding a putative [2Fe-2S] ferredoxin of G. lamblia. The Giardia intron contains a canonical splice site at its 3' end (AG), a noncanonical splice site at its 5' end (CT), and a branch point sequence that fits the yeast consensus sequence of TACTAAC except for the first nucleotide (AACTAAC). We have also identified several G. lamblia genes with spliceosomal peptides, including homologues of eukaryote-specific spliceosomal peptides (Prp8 and Prp11), several DExH-box RNA-helicases that have homologues in eubacteria, but serve essential functions in the splicing of introns in eukaryotes, and 11 predicted archaebacteria-like Sm and like-Sm core peptides, which coat small nuclear RNAs. Phylogenetic analyses show the Giardia Sm core peptides are the products of multiple, ancestral gene duplications followed by divergence, but they retain strong similarity to Sm and like-Sm peptides of other eukaryotes. Although we have documented only a single intron in Giardia, it likely has other introns and fully functional, spliceosomal machinery. If introns were added during eukaryotic evolution (the introns-late hypothesis), then these results push back the date of this event before the branching of G. lamblia.

The 2,160,837-base pair genome sequence of an isolate of Streptococcus pneumoniae, a Gram-positive pathogen that causes pneumonia, bacteremia, meningitis, and otitis media, contains 2236 predicted coding regions; of these, 1440 (64%) were assigned a biological role. Approximately 5% of the genome is composed of insertion sequences that may contribute to genome rearrangements through uptake of foreign DNA. Extracellular enzyme systems for the metabolism of polysaccharides and hexosamines provide a substantial source of carbon and nitrogen for S. pneumoniae and also damage host tissues and facilitate colonization. A motif identified within the signal peptide of proteins is potentially involved in targeting these proteins to the cell surface of low-guanine/cytosine (GC) Gram-positive species. Several surface-exposed proteins that may serve as potential vaccine candidates were identified. Comparative genome hybridization with DNA arrays revealed strain differences in S. pneumoniae that could contribute to differences in virulence and antigenicity.

Paget's disease of bone (PDB) is a common disorder characterized by focal areas of increased and disorganized osteoclastic bone resorption, leading to bone pain, deformity, pathological fracture, and an increased risk of osteosarcoma. Genetic factors play an important role in the pathogenesis of Paget's disease. In some families, the disease has been found to be linked to a susceptibility locus on chromosome 18q21-22, which also contains the gene responsible for familial expansile osteolysis (FEO)--a rare bone dysplasia with many similarities to Paget's disease. Insertion mutations of the TNFRSF11A gene encoding Receptor Activator of NF kappa B (RANK) have recently been found to be responsible for FEO and rare cases of early onset familial Paget's disease. Loss of heterozygosity (LOH) affecting the PDB/FEO critical region has also been described in osteosarcomas suggesting that TNFRSF11A might also be involved in the development of osteosarcoma. In order to investigate the possible role of TNFRSF11A in the pathogenesis of Paget's disease and osteosarcoma, we conducted mutation screening of the TNFRSF11A gene in patients with familial and sporadic Paget's disease as well as DNA extracted from Pagetic bone lesions, an osteosarcoma arising in Pagetic bone and six osteosarcoma cell lines. No specific abnormalities of the TNFRSF11A gene were identified in a Pagetic osteosarcoma, the osteosarcoma cell lines, DNA extracted from Pagetic bone lesions, or DNA extracted from peripheral blood in patients with familial or sporadic Paget's disease including several individuals with early onset Paget's disease. These data indicate that TNFRSF11A mutations contribute neither to the vast majority of cases of sporadic or familial PDB, nor to the development of osteosarcoma.

TIGRFAMs is a collection of protein families featuring curated multiple sequence alignments, hidden Markov models and associated information designed to support the automated functional identification of proteins by sequence homology. We introduce the term 'equivalog' to describe members of a set of homologous proteins that are conserved with respect to function since their last common ancestor. Related proteins are grouped into equivalog families where possible, and otherwise into protein families with other hierarchically defined homology types. TIGRFAMs currently contains over 800 protein families, available for searching or downloading at www.tigr.org/TIGRFAMs. Classification by equivalog family, where achievable, complements classification by orthology, superfamily, domain or motif. It provides the information best suited for automatic assignment of specific functions to proteins from large-scale genome sequencing projects.

We have determined the detailed molecular structure and evolution of an alpha satellite junction from human chromosome 16p11. The analysis reveals that the alpha satellite sequence bordering the transition lacks higher-order structure and that the non-alpha satellite portion consists of a mosaic of duplicated segments of complex evolutionary origin. The 16p11 junction was formed recently (5-10 million years ago) by the duplication and transposition of genomic segments from Xq28 and 4q24. Once this mosaic structure was formed, a larger complex was spread among multiple pericentromeric regions. This resulted in the formation of large (>62 kb) paralogous segments that share a high degree ( approximately 97%) of sequence similarity. Both phylogenetic and comparative analyses indicate that these pericentromeric-directed duplications occurred around the time of the divergence of the human, gorilla and chimpanzee lineages, resulting in the subtle restructuring of the primate genome among these species. The available data suggest that such chimeric structures are a general property of several different human chromosomes near their alpha satellite junctions.

Several publicly funded large-scale sequencing efforts have been initiated with the goal of completing the first reference human genome sequence by the year 2005. Here we present the results of analysis of 11.8 Mb of genomic sequence from chromosome 16. The apparent gene density varies throughout the region, but the number of genes predicted (84) suggests that this is a gene-poor region. This result may also suggest that the total number of human genes is likely to be at the lower end of published estimates. One of the most interesting aspects of this region of the genome is the presence of highly homologous, recently duplicated tracts of sequence distributed throughout the p-arm. Such duplications have implications for mapping and gene analysis as well as the predisposition to recurrent chromosomal structural rearrangements associated with genetic disease., (Copyright 1999 Academic Press.)