Aspergillus fumigatus is a human fungal pathogen that causes aspergillosis in immunocompromised hosts. A. fumigatus is believed to be exposed to diverse environmental stresses in the host cells. The ...adaptation mechanisms are critical for infections in human bodies. Transcriptional networks in response to diverse environmental challenges remain to be elucidated. To gain insights into the adaptation to environmental stresses in A. fumigatus mycelia, we conducted time series transcriptome analyses.
With the aid of RNA-seq, we explored the global gene expression profiles of mycelia in A. fumigatus upon exposure to diverse environmental changes, including heat, superoxide, and osmotic stresses. From the perspective of global transcriptomes, transient responses to superoxide and osmotic stresses were observed while responses to heat stresses were gradual. We identified the stress-responsive genes for particular stresses, and the 266 genes whose expression levels drastically fluctuated upon exposure to all tested stresses. Among these, the 77 environmental stress response genes are conserved in S. cerevisiae, suggesting that these genes might be more general prerequisites for adaptation to environmental stresses. Finally, we revealed the strong correlations among expression profiles of genes related to 'rRNA processing'.
The time series transcriptome analysis revealed the stress-responsive genes underlying the adaptation mechanisms in A. fumigatus mycelia. These results will shed light on the regulatory networks underpinning the adaptation of the filamentous fungi.
Infectious disease is the single greatest threat to taxa such as amphibians (chytrid fungus), bats (white nose syndrome), Tasmanian devils (devil facial tumor disease), and black-footed ferrets ...(canine distemper virus, plague). Although understanding the genetic basis to disease susceptibility is important for the long-term persistence of these groups, most research has been limited to major-histocompatibility and Toll-like receptor genes. To better understand the genetic basis of infectious disease susceptibility in a species of conservation concern, we sequenced all known/predicted immune response genes (i.e., the immunomes) in 16 Florida gopher tortoises, Gopherus polyphemus. All tortoises produced antibodies against Mycoplasma agassizii (an etiologic agent of infectious upper respiratory tract disease; URTD) and, at the time of sampling, either had (n = 10) or lacked (n = 6) clinical signs.
We found several variants associated with URTD clinical status in complement and lectin genes, which may play a role in Mycoplasma immunity. Thirty-five genes deviated from neutrality according to Tajima's D. These genes were enriched in functions relating to macromolecule and protein modifications, which are vital to immune system functioning.
These results are suggestive of genetic differences that might contribute to disease severity, a finding that is consistent with other mycoplasmal diseases. This has implications for management because tortoises across their range may possess genetic variation associated with a more severe response to URTD. More generally: 1) this approach demonstrates that a broader consideration of immune genes is better able to identify important variants, and; 2) this data pipeline can be adopted to identify alleles associated with disease susceptibility or resistance in other taxa, and therefore provide information on a population's risk of succumbing to disease, inform translocations to increase genetic variation for disease resistance, and help to identify potential treatments.
In Nile tilapia sex determination is governed by a male heterogametic system XX/XY either on LG1 or LG23. The latter carries a Y-specific duplicate of the amh gene, which is a testis-determining ...factor. Allelic variants in the amh gene demonstrated to be major triggers for autosomal and temperature-dependent sex reversal. Further, QTL on LG23 and LG20 show a temperature-responsiveness with influence on the phenotypic sex relative to the sex chromosomes. Here we present a ddRADseq based approach to identify genomic regions that show unusual large differentiation in terms of fixation index (F
) between temperature-treated pseudomales and non-masculinized females using a comparative genome-scan. Genome-wide associations were identified for the temperature-dependent sex using a genetically all-female population devoid of amh-ΔY.
Twenty-two thousand three hundred ninety-two SNPs were interrogated for the comparison of temperature-treated pseudomales and females, which revealed the largest differentiation on LG23. Outlier F
-values (0.35-0.44) were determined for six SNPs in the genomic interval (9,190,077-11,065,693) harbouring the amh gene (9,602,693-9,605,808), exceeding the genome-wide low F
of 0.013. Association analysis with a set of 9104 selected SNPs confirmed that the same genomic region on LG23 exerts a significant effect on the temperature-dependent sex.
This study highlights the role of LG23 in sex determination, harbouring major determinants for temperature-dependent sex reversal in Nile tilapia. Furthermore F
outlier detection proves a powerful tool for detection of sex-determining regions in fish genomes.
Polyketide synthase (PKS) catalyzes the biosynthesis of polyketides, which are structurally and functionally diverse natural products in microorganisms and plants. Here, we have analyzed available ...full genome sequences of microscopic and macroscopic algae for the presence of type I PKS genes.
Type I PKS genes are present in 15 of 32 analyzed algal species. In chlorophytes, large proteins in the MDa range are predicted in most sequenced species, and PKSs with free-standing acyltransferase domains (trans-AT PKSs) predominate. In a phylogenetic tree, PKS sequences from different algal phyla form clades that are distinct from PKSs from other organisms such as non-photosynthetic protists or cyanobacteria. However, intermixing is observed in some cases, for example polyunsaturated fatty acid (PUFA) and glycolipid synthases of various origins. Close relationships between type I PKS modules from different species or between modules within the same multimodular enzyme were identified, suggesting module duplications during evolution of algal PKSs. In contrast to type I PKSs, nonribosomal peptide synthetases (NRPSs) are relatively rare in algae (occurrence in 7 of 32 species).
Our phylogenetic analysis of type I PKSs in algae supports an evolutionary scenario whereby integrated AT domains were displaced to yield trans-AT PKSs. Together with module duplications, the displacement of AT domains may constitute a major mechanism of PKS evolution in algae. This study advances our understanding of the diversity of eukaryotic PKSs and their evolutionary trajectories.
Sequence alignments are the starting point for most evolutionary and comparative analyses. Full genome sequences can be compared to study patterns of within and between species variation. Genome ...sequence alignments are complex structures containing information such as coordinates, quality scores and synteny structure, which are stored in Multiple Alignment Format (MAF) files. Processing these alignments therefore involves parsing and manipulating typically large MAF files in an efficient way.
MafFilter is a command-line driven program written in C++ that enables the processing of genome alignments stored in the Multiple Alignment Format in an efficient and extensible manner. It provides an extensive set of tools which can be parametrized and combined by the user via option files. We demonstrate the software's functionality and performance on several biological examples covering Primate genomics and fungal population genomics. Example analyses involve window-based alignment filtering, feature extractions and various statistics, phylogenetics and population genomics calculations.
MafFilter is a highly efficient and flexible tool to analyse multiple genome alignments. By allowing the user to combine a large set of available methods, as well as designing his/her own, it enables the design of custom data filtering and analysis pipelines for genomic studies. MafFilter is an open source software available at http://bioweb.me/maffilter.
The advances and decreasing economical cost of whole genome sequencing (WGS), will soon make this technology available for routine infectious disease epidemiology. In epidemiological studies, ...outbreak isolates have very little diversity and require extensive genomic analysis to differentiate and classify isolates. One of the successfully and broadly used methods is analysis of single nucletide polymorphisms (SNPs). Currently, there are different tools and methods to identify SNPs including various options and cut-off values. Furthermore, all current methods require bioinformatic skills. Thus, we lack a standard and simple automatic tool to determine SNPs and construct phylogenetic tree from WGS data.
Here we introduce snpTree, a server for online-automatic SNPs analysis. This tool is composed of different SNPs analysis suites, perl and python scripts. snpTree can identify SNPs and construct phylogenetic trees from WGS as well as from assembled genomes or contigs. WGS data in fastq format are aligned to reference genomes by BWA while contigs in fasta format are processed by Nucmer. SNPs are concatenated based on position on reference genome and a tree is constructed from concatenated SNPs using FastTree and a perl script. The online server was implemented by HTML, Java and python script.The server was evaluated using four published bacterial WGS data sets (V. cholerae, S. aureus CC398, S. Typhimurium and M. tuberculosis). The evaluation results for the first three cases was consistent and concordant for both raw reads and assembled genomes. In the latter case the original publication involved extensive filtering of SNPs, which could not be repeated using snpTree.
The snpTree server is an easy to use option for rapid standardised and automatic SNP analysis in epidemiological studies also for users with limited bioinformatic experience. The web server is freely accessible at http://www.cbs.dtu.dk/services/snpTree-1.0/.
Gene variants within regulatory regions are thought to be major contributors of the variation of complex traits/diseases. Genome wide association studies (GWAS), have identified scores of genetic ...variants that appear to contribute to human disease risk. However, most of these variants do not appear to be functional. Thus, the significance of the association may be brought up by still unknown mechanisms or by linkage disequilibrium (LD) with functional polymorphisms. In the present study, focused on functional variants related with the binding of microRNAs (miR), we utilized SNP data, including newly released 1000 Genomes Project data to perform a genome-wide scan of SNPs that abrogate or create miR recognition element (MRE) seed sites (MRESS).
We identified 2723 SNPs disrupting, and 22295 SNPs creating MRESSs. We estimated the percent of SNPs falling within both validated (5%) and predicted conserved MRESSs (3%). We determined 87 of these MRESS SNPs were listed in GWAS association studies, or in strong LD with a GWAS SNP, and may represent the functional variants of identified GWAS SNPs. Furthermore, 39 of these have evidence of co-expression of target mRNA and the predicted miR. We also gathered previously published eQTL data supporting a functional role for four of these SNPs shown to associate with disease phenotypes. Comparison of FST statistics (a measure of population subdivision) for predicted MRESS SNPs against non MRESS SNPs revealed a significantly higher (P = 0.0004) degree of subdivision among MRESS SNPs, suggesting a role for these SNPs in environmentally driven selection.
We have demonstrated the potential of publicly available resources to identify high priority candidate SNPs for functional studies and for disease risk prediction.
To identify the genes involved in the development of low temperature (LT) tolerance in hexaploid wheat, we examined the global changes in expression in response to cold of the 55,052 potentially ...unique genes represented in the Affymetrix Wheat Genome microarray. We compared the expression of genes in winter-habit (winter Norstar and winter Manitou) and spring-habit (spring Manitou and spring Norstar)) cultivars, wherein the locus for the vernalization gene Vrn-A1 was swapped between the parental winter Norstar and spring Manitou in the derived near-isogenic lines winter Manitou and spring Norstar. Global expression of genes in the crowns of 3-leaf stage plants cold-acclimated at 6°C for 0, 2, 14, 21, 38, 42, 56 and 70 days was examined.
Analysis of variance of gene expression separated the samples by genetic background and by the developmental stage before or after vernalization saturation was reached. Using gene-specific ANOVA we identified 12,901 genes (at p < 0.001) that change in expression with respect to both genotype and the duration of cold-treatment. We examined in more detail a subset of these genes (2,771) where expression was highly influenced by the interaction between these two main factors. Functional assignments using GO annotations showed that genes involved in transport, oxidation-reduction, and stress response were highly represented. Clustering based on the pattern of transcript accumulation identified genes that were up or down-regulated by cold-treatment. Our data indicate that the cold-sensitive lines can up-regulate known cold-responsive genes comparable to that of cold-hardy lines. The levels of expression of these genes were highly influenced by the initial rate and the duration of the gene's response to cold. We show that the Vrn-A1 locus controls the duration of gene expression but not its initial rate of response to cold treatment. Furthermore, we provide evidence that Ta.Vrn-A1 and Ta.Vrt1 originally hypothesized to encode for the same gene showed different patterns of expression and therefore are distinct.
This study provides novel insight into the underlying mechanisms that regulate the expression of cold-responsive genes in wheat. The results support the developmental model of LT tolerance gene regulation and demonstrate the complex genotype by environment interactions that determine LT adaptation in winter annual cereals.
The conventional variant calling of pathogenic alleles in exome and genome sequencing requires the presence of the non-pathogenic alleles as genome references. This hinders the correct identification ...of variants with minor and/or pathogenic reference alleles warranting additional approaches for variant calling.
More than 26,000 Exome Aggregation Consortium (ExAC) variants have a minor reference allele including variants with known ClinVar disease alleles. For instance, in a number of variants related to clotting disorders, the phenotype-associated allele is a human genome reference allele (rs6025, rs6003, rs1799983, and rs2227564 using the assembly hg19). We highlighted how the current variant calling standards miss homozygous reference disease variants in these sites and provided a bioinformatic panel that can be used to screen these variants using commonly available variant callers. We present exome sequencing results from an individual with venous thrombosis to emphasize how pathogenic alleles in clinically relevant variants escape variant calling while non-pathogenic alleles are detected.
This article highlights the importance of specialized variant calling strategies in clinical variants with minor reference alleles especially in the context of personal genomes and exomes. We provide here a simple strategy to screen potential disease-causing variants when present in homozygous reference state.
PDF