Search

Abstract Background The comparison of complete genomes has revealed surprisingly large numbers of conserved non-protein-coding (CNC) DNA regions. However, the biological function of CNC remains elusive. CNC differ in two aspects from conserved protein-coding regions. They are not conserved across phylum boundaries, and they do not contain readily detectable sub-domains. Here we characterize the persistence length and time of CNC and conserved protein-coding regions in the vertebrate and insect lineages. Results The persistence length is the length of a genome region over which a certain level of sequence identity is consistently maintained. The persistence time is the evolutionary period during which a conserved region evolves under the same selective constraints. Our main findings are: (i) Insect genomes contain 1.60 times less conserved information than vertebrates; (ii) Vertebrate CNC have a higher persistence length than conserved coding regions or insect CNC; (iii) CNC have shorter persistence times as compared to conserved coding regions in both lineages. Conclusion Higher persistence length of vertebrate CNC indicates that the conserved information in vertebrates and insects is organized in functional elements of different lengths. These findings might be related to the higher morphological complexity of vertebrates and give clues about the structure of active CNC elements. Shorter persistence time might explain the previously puzzling observations of highly conserved CNC within each phylum, and of a lack of conservation between phyla. It suggests that CNC divergence might be a key factor in vertebrate evolution. Further evolutionary studies will help to relate individual CNC to specific developmental processes.

Abstract Background High throughput gene expression profiling (GEP) is becoming a routine technique in life science laboratories. With experimental designs that repeatedly span thousands of genes and hundreds of samples, relying on a dedicated database infrastructure is no longer an option. GEP technology is a fast moving target, with new approaches constantly broadening the field diversity. This technology heterogeneity, compounded by the informatics complexity of GEP databases, means that software developments have so far focused on mainstream techniques, leaving less typical yet established techniques such as Nylon microarrays at best partially supported. Results MAF (MicroArray Facility) is the laboratory database system we have developed for managing the design, production and hybridization of spotted microarrays. Although it can support the widely used glass microarrays and oligo-chips, MAF was designed with the specific idiosyncrasies of Nylon based microarrays in mind. Notably single channel radioactive probes, microarray stripping and reuse, vector control hybridizations and spike-in controls are all natively supported by the software suite. MicroArray Facility is MIAME supportive and dynamically provides feedback on missing annotations to help users estimate effective MIAME compliance. Genomic data such as clone identifiers and gene symbols are also directly annotated by MAF software using standard public resources. The MAGE-ML data format is implemented for full data export. Journalized database operations (audit tracking), data anonymization, material traceability and user/project level confidentiality policies are also managed by MAF. Conclusion MicroArray Facility is a complete data management system for microarray producers and end-users. Particular care has been devoted to adequately model Nylon based microarrays. The MAF system, developed and implemented in both private and academic environments, has proved a robust solution for shared facilities and industry service providers alike.

The climate fluctuations of the Quaternary shaped the movement of species in and out of glacial refugia. In Europe, the majority of species followed one of the described traditional postglacial recolonization routes from the southern peninsulas towards the north. Like most organisms, barn owls are assumed to have colonized the British Isles by crossing over Doggerland, a land bridge that connected Britain to northern Europe. However, while they are dark rufous in northern Europe, barn owls in the British Isles are conspicuously white, a contrast that could suggest selective forces are at play on the islands. Yet, our analysis of known candidate genes involved in coloration found no signature of selection. Instead, using whole genome sequences and species distribution modelling, we found that owls colonised the British Isles soon after the last glaciation, directly from a white coloured refugium in the Iberian Peninsula, before colonising northern Europe. They would have followed a hitherto unknown post-glacial colonization route to the Isles over a westwards path of suitable habitat in now submerged land in the Bay of Biscay, thus not crossing Doggerland. As such, they inherited the white colour of their Iberian founders and maintained it through low gene flow with the mainland that prevents the import of rufous alleles. Thus, we contend that neutral processes probably explain this contrasting white colour compared to continental owls. With the barn owl being a top predator, we expect future research will show this unanticipated route was used by other species from its paleo community., (© 2021 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Dioecious plants vary in whether their sex chromosomes are heteromorphic or homomorphic, but even homomorphic sex chromosomes may show divergence between homologues in the non-recombining, sex-determining region (SDR). Very little is known about the SDR of these species, which might represent particularly early stages of sex-chromosome evolution. Here, we assess the size and content of the SDR of the diploid dioecious herb Mercurialis annua , a species with homomorphic sex chromosomes and mild Y-chromosome degeneration. We used RNA sequencing (RNAseq) to identify new Y-linked markers for M. annua. Twelve of 24 transcripts showing male-specific expression in a previous experiment could be amplified by polymerase chain reaction (PCR) only from males, and are thus likely to be Y-linked. Analysis of genome-capture data from multiple populations of M. annua pointed to an additional six male-limited (and thus Y-linked) sequences. We used these markers to identify and sequence 17 sex-linked bacterial artificial chromosomes (BACs), which form 11 groups of non-overlapping sequences, covering a total sequence length of about 1.5 Mb. Content analysis of this region suggests that it is enriched for repeats, has low gene density, and contains few candidate sex-determining genes. The BACs map to a subset of the sex-linked region of the genetic map, which we estimate to be at least 14.5 Mb. This is substantially larger than estimates for other dioecious plants with homomorphic sex chromosomes, both in absolute terms and relative to their genome sizes. Our data provide a rare, high-resolution view of the homomorphic Y chromosome of a dioecious plant.

Candida glabrata is an important fungal pathogen which develops rapid antifungal resistance in treated patients. It is known that azole treatments lead to antifungal resistance in this fungal species and that multidrug efflux transporters are involved in this process. Specific mutations in the transcriptional regulator PDR1 result in upregulation of the transporters. In addition, we showed that the PDR1 mutations can contribute to enhance virulence in animal models. In this study, we were interested to compare genomes of two specific C. glabrata -related isolates, one of which was azole susceptible (DSY562) while the other was azole resistant (DSY565). DSY565 contained a PDR1 mutation (L280F) and was isolated after a time-lapse of 50 d of azole therapy. We expected that genome comparisons between both isolates could reveal additional mutations reflecting host adaptation or even additional resistance mechanisms. The PacBio technology used here yielded 14 major contigs (sizes 0.18-1.6 Mb) and mitochondrial genomes from both DSY562 and DSY565 isolates that were highly similar to each other. Comparisons of the clinical genomes with the published CBS138 genome indicated important genome rearrangements, but not between the clinical strains. Among the unique features, several retrotransposons were identified in the genomes of the investigated clinical isolates. DSY562 and DSY565 each contained a large set of adhesin-like genes (101 and 107, respectively), which exceed by far the number of reported adhesins (63) in the CBS138 genome. Comparison between DSY562 and DSY565 yielded 17 nonsynonymous SNPs (among which the was the expected PDR1 mutation) as well as small size indels in coding regions (11) but mainly in adhesin-like genes. The genomes contained a DNA mismatch repair allele of MSH2 known to be involved in the so-called hyper-mutator phenotype of this yeast species and the number of accumulated mutations between both clinical isolates is consistent with the presence of a MSH2 defect. In conclusion, this study is the first to compare genomes of C. glabrata sequential clinical isolates using the PacBio technology as an approach. The genomes of these isolates taken in the same patient at two different time points exhibited limited variations, even if submitted to the host pressure., (Copyright © 2017 Vale-Silva et al.)

Mycoplasma mycoidessubsp.mycoidesis the etiologic agent of contagious bovine pleuropneumonia. We report here the complete genome sequence of the strain T1/44, which is widely used as a live vaccine in Africa., (Copyright © 2016 Gourgues et al.)

The surrounding capsule of Streptococcus pneumoniae has been identified as a major virulence factor and is targeted by pneumococcal conjugate vaccines (PCV). However, nonencapsulated S. pneumoniae (non-Ec-Sp) have also been isolated globally, mainly in carriage studies. It is unknown if non-Ec-Sp evolve sporadically, if they have high antibiotic nonsusceptiblity rates and a unique, specific gene content. Here, whole-genome sequencing of 131 non-Ec-Sp isolates sourced from 17 different locations around the world was performed. Results revealed a deep-branching classic lineage that is distinct from multiple sporadic lineages. The sporadic lineages clustered with a previously sequenced, global collection of encapsulated S. pneumoniae (Ec-Sp) isolates while the classic lineage is comprised mainly of the frequently identified multilocus sequences types (STs) ST344 (n = 39) and ST448 (n = 40). All ST344 and nine ST448 isolates had high nonsusceptiblity rates to β-lactams and other antimicrobials. Analysis of the accessory genome reveals that the classic non-Ec-Sp contained an increased number of mobile elements, than Ec-Sp and sporadic non-Ec-Sp. Performing adherence assays to human epithelial cells for selected classic and sporadic non-Ec-Sp revealed that the presence of a integrative conjugative element (ICE) results in increased adherence to human epithelial cells (P = 0.005). In contrast, sporadic non-Ec-Sp lacking the ICE had greater growth in vitro possibly resulting in improved fitness. In conclusion, non-Ec-Sp isolates from the classic lineage have evolved separately. They have spread globally, are well adapted to nasopharyngeal carriage and are able to coexist with Ec-Sp. Due to continued use of PCV, non-Ec-Sp may become more prevalent., (© The Author(s) 2014. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

We present and validate BlastR, a method for efficiently and accurately searching non-coding RNAs. Our approach relies on the comparison of di-nucleotides using BlosumR, a new log-odd substitution matrix. In order to use BlosumR for comparison, we recoded RNA sequences into protein-like sequences. We then showed that BlosumR can be used along with the BlastP algorithm in order to search non-coding RNA sequences. Using Rfam as a gold standard, we benchmarked this approach and show BlastR to be more sensitive than BlastN. We also show that BlastR is both faster and more sensitive than BlastP used with a single nucleotide log-odd substitution matrix. BlastR, when used in combination with WU-BlastP, is about 5% more accurate than WU-BlastN and about 50 times slower. The approach shown here is equally effective when combined with the NCBI-Blast package. The software is an open source freeware available from www.tcoffee.org/blastr.html.

Ants have evolved very complex societies and are key ecosystem members. Some ants, such as the fire ant Solenopsis invicta, are also major pests. Here, we present a draft genome of S. invicta, assembled from Roche 454 and Illumina sequencing reads obtained from a focal haploid male and his brothers. We used comparative genomic methods to obtain insight into the unique features of the S. invicta genome. For example, we found that this genome harbors four adjacent copies of vitellogenin. A phylogenetic analysis revealed that an ancestral vitellogenin gene first underwent a duplication that was followed by possibly independent duplications of each of the daughter vitellogenins. The vitellogenin genes have undergone subfunctionalization with queen- and worker-specific expression, possibly reflecting differential selection acting on the queen and worker castes. Additionally, we identified more than 400 putative olfactory receptors of which at least 297 are intact. This represents the largest repertoire reported so far in insects. S. invicta also harbors an expansion of a specific family of lipid-processing genes, two putative orthologs to the transformer/feminizer sex differentiation gene, a functional DNA methylation system, and a single putative telomerase ortholog. EST data indicate that this S. invicta telomerase ortholog has at least four spliceforms that differ in their use of two sets of mutually exclusive exons. Some of these and other unique aspects of the fire ant genome are likely linked to the complex social behavior of this species.

Profiling microRNA (miRNA) expression is of widespread interest given the critical role of miRNAs in many cellular functions. Profiling can be achieved via hybridization-based (microarrays), sequencing-based, or amplification-based (quantitative reverse transcription-PCR, qPCR) technologies. Among these, microarrays face the significant challenge of accurately distinguishing between mature and immature miRNA forms, and different vendors have developed different methods to meet this challenge. Here we measure differential miRNA expression using the Affymetrix, Agilent, and Illumina microarray platforms, as well as qPCR (Applied Biosystems) and ultra high-throughput sequencing (Illumina). We show that the differential expression measurements are more divergent when the three types of microarrays are compared than when the Agilent microarray, qPCR, and sequencing technology measurements are compared, which exhibit a good overall concordance.

Background: The comparison of complete genomes has revealed surprisingly large numbers of conserved non-protein-coding (CNC) DNA regions. However, the biological function of CNC remains elusive. CNC differ in two aspects from conserved protein-coding regions. They are not conserved across phylum boundaries, and they do not contain readily detectable sub-domains. Here we characterize the persistence length and time of CNC and conserved protein-coding regions in the vertebrate and insect lineages., Results: The persistence length is the length of a genome region over which a certain level of sequence identity is consistently maintained. The persistence time is the evolutionary period during which a conserved region evolves under the same selective constraints. Our main findings are: (i) Insect genomes contain 1.60 times less conserved information than vertebrates; (ii) Vertebrate CNC have a higher persistence length than conserved coding regions or insect CNC; (iii) CNC have shorter persistence times as compared to conserved coding regions in both lineages., Conclusion: Higher persistence length of vertebrate CNC indicates that the conserved information in vertebrates and insects is organized in functional elements of different lengths. These findings might be related to the higher morphological complexity of vertebrates and give clues about the structure of active CNC elements. Shorter persistence time might explain the previously puzzling observations of highly conserved CNC within each phylum, and of a lack of conservation between phyla. It suggests that CNC divergence might be a key factor in vertebrate evolution. Further evolutionary studies will help to relate individual CNC to specific developmental processes.

Background: High throughput gene expression profiling (GEP) is becoming a routine technique in life science laboratories. With experimental designs that repeatedly span thousands of genes and hundreds of samples, relying on a dedicated database infrastructure is no longer an option.GEP technology is a fast moving target, with new approaches constantly broadening the field diversity. This technology heterogeneity, compounded by the informatics complexity of GEP databases, means that software developments have so far focused on mainstream techniques, leaving less typical yet established techniques such as Nylon microarrays at best partially supported., Results: MAF (MicroArray Facility) is the laboratory database system we have developed for managing the design, production and hybridization of spotted microarrays. Although it can support the widely used glass microarrays and oligo-chips, MAF was designed with the specific idiosyncrasies of Nylon based microarrays in mind. Notably single channel radioactive probes, microarray stripping and reuse, vector control hybridizations and spike-in controls are all natively supported by the software suite. MicroArray Facility is MIAME supportive and dynamically provides feedback on missing annotations to help users estimate effective MIAME compliance. Genomic data such as clone identifiers and gene symbols are also directly annotated by MAF software using standard public resources. The MAGE-ML data format is implemented for full data export. Journalized database operations (audit tracking), data anonymization, material traceability and user/project level confidentiality policies are also managed by MAF., Conclusion: MicroArray Facility is a complete data management system for microarray producers and end-users. Particular care has been devoted to adequately model Nylon based microarrays. The MAF system, developed and implemented in both private and academic environments, has proved a robust solution for shared facilities and industry service providers alike.

Breast cancer is a complex genetic disease characterized by the accumulation of multiple molecular alterations. The resulting clinical heterogeneity makes current therapeutic strategies-based on clinicopathlogical factors-less than perfectly adapted to each patient. Today, DNA arrays, by allowing the simultaneous and quantitative analysis of the mRNA expression levels of thousands of genes in a single assay, provide novel tools to tackle this complexity. Potential applications are multiple in the cancer field and the first research results are promising. Using home-made DNA arrays in an approach easily compatible with academic research-nylon support and radioactive detection-we identified a predictor set of 23 genes whose expression patterns differentiated two groups of breast cancer patients with different survival after adjuvant chemotherapy. We then validated and further extended these results in a larger, independent and homogeneous series of poor prognosis primary breast cancers treated with adjuvant anthracyclin-based chemotherapy. We confirmed the prognostic classification provided by the 23-gene set predictor. We then improved the predictor set and refined the classification by sorting the tumors into three classes with significantly different long-term survival. These results show the potential of the technology with an accessible approach for academic research teams. We also showed that nylon DNA arrays with radioactive detection are associated with excellent sensitivity, an advantage in clinical situations where the amount of available material is limited.

The diversity of the T cell repertoire of mature T splenocytes is generated, in the thymus, by pairing of alpha and beta variable domains of the alpha beta TCR and by the rearrangements of various gene segments encoding these domains. In the periphery, it results from competition between various T cell subpopulations including recent thymic migrants and long-lived T cells. Quantitative data on the actual size of the T cell repertoire are lacking. Using PCR methods and extensive sequencing, we have measured for the first time the size of the TCR-alpha beta repertoire of naive mouse T splenocytes. There are 5-8 x 105 different nucleotide sequences of BV chains in the whole spleen of young adult mice. We have also determined the size of the BV repertoire in a subpopulation of AV2+ T splenocytes, which allows us to provide a minimum estimate of the alpha beta repertoire. We find that the mouse spleen harbors about 2 x 106 clones of about 10 cells each. This figure, although orders of magnitude smaller than the maximum theoretical diversity (estimated up to 1015), is still large enough to maintain a high functional diversity.

The positive selection of CD4+ T cells requires the expression of major histocompatibility complex (MHC) class II molecules in the thymus, but the role of self-peptides complexed to class II molecules is still a matter of debate. Recently, it was observed that transgenic mice expressing a single peptide-MHC class II complex positively select significant numbers of diverse CD4+ T cells in the thymus. However, the number of selected T cell specificities has not been evaluated so far. Here, we have sequenced 700 junctional complementarity determining regions 3 (CDR3) from T cell receptors (TCRs) carrying Vbeta11-Jbeta1.1 or Vbeta12-Jbeta1.1 rearrangements. We found that a single peptide-MHC class II complex positively selects at least 10(5) different Vbeta rearrangements. Our data yield a first evaluation of the size of the T cell repertoire. In addition, they provide evidence that the single Ealpha52-68-I-Ab complex skews the amino acid frequency in the TCR CDR3 loop of positively selected T cells. A detailed analysis of CDR3 sequences indicates that a fraction of the beta chain repertoire bears the imprint of the selecting self-peptide.