Genesis and expansion of metazoan transcription factor gene classes

We know little about the genomic events that led to the advent of a multicellular grade of organization in animals, one of the most dramatic transitions in evolution. Metazoan multicellularity is correlated with the evolution of embryogenesis, which presumably was underpinned by a gene regulatory network reliant on the differential activation of signaling pathways and transcription factors. Many transcription factor genes that play critical roles in bilaterian development largely appear to have evolved before the divergence of cnidarian and bilaterian lineages. In contrast, sponges seem to have a more limited suite of transcription factors, suggesting that the developmental regulatory gene repertoire changed markedly during early metazoan evolution. Using whole- genome information from the sponge Amphimedon queenslandica, a range of eumetazoans, and the choanoflagellate Monosiga brevicollis, we investigate the genesis and expansion of homeobox, Sox, T- box, and Fox transcription factor genes. Comparative analyses reveal that novel transcription factor domains ( such as Paired, POU, and T- box) arose very early in metazoan evolution, prior to the separation of extant metazoan phyla but after the divergence of choanoflagellate and metazoan lineages. Phylogenetic analyses indicate that transcription factor classes then gradually expanded at the base of Metazoa before the bilaterian radiation, with each class following a different evolutionary trajectory. Based on the limited number of transcription factors in the Amphimedon genome, we infer that the genome of the metazoan last common ancestor included fewer gene members in each class than are present in extant eumetazoans. Transcription factor orthologues present in sponge, cnidarian, and bilaterian genomes may represent part of the core metazoan regulatory network underlying the origin of animal development and multicellularity.

Evolution of microbial virulence: the benefits of stress

Although genome sequencing of microbial pathogens has shed light on the evolution of virulence, the drivers of the gain and loss of genes and of pathogenicity islands (gene clusters), which contribute to the emergence of new disease outbreaks, are unclear. Recent experiments with the bean pathogen Pseudomonas syringae pv. phaseolicola illustrate how exposure to resistance mechanisms acts as the driving force for genome reorganization. Here we argue that the antimicrobial conditions generated by host defences can accelerate the generation of genome rearrangements that provide selective advantages to the invading microbe. Similar exposure to environmental stress outside the host could also drive the horizontal gene transfer that has led to the evolution of pathogenicity towards both animals and plants.

N-nitroso compound exposure-associated transcriptomic profiles are indicative of an increased risk for colorectal cancer

Endogenous formation of N-nitroso compounds (NOCs), which are known animal carcinogens, could contribute to human carcinogenesis but definitive evidence is still lacking. To investigate the relevance of NOCs in human colorectal cancer (CRC) development, we analyzed whole genome gene expression modifications in human colon biopsies in relation to fecal NOC exposure. We had a particular interest in patients suffering from intestinal inflammation as this may stimulate endogenous NOC formation, and consequently predispose to CRC risk. Inflammatory bowel disease (IBD) patients diagnosed with ulcerative colitis and irritable bowel syndrome patients without inflammation, serving as controls, were therefore recruited. Fecal NOC were demonstrated in the majority of subjects. By associating gene expression levels of all subjects to fecal NOC levels, we identified a NOC exposure-associated transcriptomic response that suggests that physiological NOC concentrations may potentially induce genotoxic responses and chromatin modifications in human colon tissue, both of which are linked to carcinogenicity. In a network analysis, chromatin modifications were linked to 11 significantly modulated histone genes, pointing towards a possible epigenetic mechanism that may be relevant in comprehending NOC-induced carcinogenesis. In addition, pro-inflammatory transcriptomic modifications were identified in visually non-inflamed regions of the IBD colon. However, fecal NOC levels were slightly but not significantly increased in IBD patients, suggesting that inflammation did not strongly stimulate NOC formation. We conclude that NOC exposure is associated with gene expression modifications in the human colon that may suggest a potential role of these compounds in CRC development.

A comparison of the Thlaspi caerulescens and Thlaspi arvense shoot transcriptomes

Whole-genome transcriptome profiling is revealing how biological systems are regulated at the transcriptional level. This study reports the development of a robust method to profile and compare the transcriptomes of two nonmodel plant species, Thlaspi caerulescens, a zinc (Zn) hyperaccumulator, and Thlaspi arvense, a nonhyperaccumulator, using Affymetrix Arabidopsis thaliana ATH1-121501 GeneChip (R) arrays (Affymetrix, Santa Clara, CA, USA). Transcript abundance was quantified in the shoots of agar- and compost-grown plants of both species. Analyses were optimized using a genomic DNA (gDNA)-based probe-selection strategy based on the hybridization efficiency of Thlaspi gDNA with corresponding A. thaliana probes. In silico alignments of GeneChip (R) probes with Thlaspi gene sequences, and quantitative real-time PCR, confirmed the validity of this approach. Approximately 5000 genes were differentially expressed in the shoots of T. caerulescens compared with T. arvense, including genes involved in Zn transport and compartmentalization. Future functional analyses of genes identified as differentially expressed in the shoots of these closely related species will improve our understanding of the molecular mechanisms of Zn hyperaccumulation.

Non-replication of the association between 5HTTLPR and response to psychological therapy for child anxiety disorders

Background We previously reported an association between 5HTTLPR genotype and outcome following cognitive–behavioural therapy (CBT) in child anxiety (Cohort 1). Children homozygous for the low-expression short-allele showed more positive outcomes. Other similar studies have produced mixed results, with most reporting no association between genotype and CBT outcome. Aims To replicate the association between 5HTTLPR and CBT outcome in child anxiety from the Genes for Treatment study (GxT Cohort 2, n = 829). Method Logistic and linear mixed effects models were used to examine the relationship between 5HTTLPR and CBT outcomes. Mega-analyses using both cohorts were performed. Results There was no significant effect of 5HTTLPR on CBT outcomes in Cohort 2. Mega-analyses identified a significant association between 5HTTLPR and remission from all anxiety disorders at follow-up (odds ratio 0.45, P = 0.014), but not primary anxiety disorder outcomes. Conclusions The association between 5HTTLPR genotype and CBT outcome did not replicate. Short-allele homozygotes showed more positive treatment outcomes, but with small, non-significant effects. Future studies would benefit from utilising whole genome approaches and large, homogenous samples.

Isolation of novel multipotent neural crest-derived stem cells from adult human inferior turbinate

Adult human neural crest-derived stem cells (NCSCs) are of extraordinary high plasticity and promising candidates for the use in regenerative medicine. Here we describe for the first time a novel neural crest-derived stem cell population within the respiratory epithelium of human adult inferior turbinate. In contrast to superior and middle turbinates, high amounts of source material could be isolated from human inferior turbinates. Using minimally-invasive surgery methods isolation is efficient even in older patients. Within their endogenous niche, inferior turbinate stem cells (ITSCs) expressed high levels of nestin, p75(NTR), and S100. Immunoelectron microscopy using anti-p75 antibodies displayed that ITSCs are of glial origin and closely related to nonmyelinating Schwann cells. Cultivated ITSCs were positive for nestin and S100 and the neural crest markers Slug and SOX10. Whole genome microarray analysis showed pronounced differences to human ES cells in respect to pluripotency markers OCT4, SOX2, LIN28, and NANOG, whereas expression of WDR5, KLF4, and c-MYC was nearly similar. ITSCs were able to differentiate into cells with neuro-ectodermal and mesodermal phenotype. Additionally ITSCs are able to survive and perform neural crest typical chain migration in vivo when transplanted into chicken embryos. However ITSCs do not form teratomas in severe combined immunodeficient mice. Finally, we developed a separation strategy based on magnetic cell sorting of p75(NTR) positive ITSCs that formed larger neurospheres and proliferated faster than p75(NTR) negative ITSCs. Taken together our study describes a novel, readily accessible source of multipotent human NCSCs for potential cell-replacement therapy.

Transcriptome analysis of nicotine-exposed cells from the brainstem of neonate spontaneously hypertensive and Wistar Kyoto rats

In this study, the effects of nicotine on global gene expression of cultured cells from the brainstem of spontaneously hypertensive rat (SHR) and normotensive Wistar Kyoto (WKY) rats were evaluated using whole-genome oligoarrays. We found that nicotine may act differentially on the gene expression profiles of SHR and WKY. The influence of strain was present in 321 genes that were differentially expressed in SHR as compared with WKY brainstem cells independently of the nicotine treatment. A total of 146 genes had their expression altered in both strains after nicotine exposure. Interaction between nicotine treatment and the strain was observed to affect the expression of 229 genes that participate in cellular pathways related to neurotransmitter secretion, intracellular trafficking and cell communication, and are possibly involved in the phenotypic differentiation between SHR and WKY rats, including hypertension. Further characterization of their function in hypertension development is warranted. The Pharmacogenomics Journal (2010) 10, 134-160; doi:10.1038/tpj.2009.42; published online 15 September 2009

Genomic imbalances associated with mullerian aplasia

Background: Aplasia of the mullerian ducts leads to absence of the uterine corpus, uterine cervix, and upper (superior) vagina. Patients with mullerian aplasia (MA) often exhibit additional clinical features such as renal, vertebral and cardiac defects. A number of different syndromes have been associated with MA, and in most cases its aetiology remains poorly understood. Objective and methods: 14 syndromic patients with MA and 46, XX G-banded karyotype were screened for DNA copy number changes by similar to 1 Mb whole genome bacterial artificial chromosome (BAC) array based comparative genomic hybridisation (CGH). The detected alterations were validated by an independent method and further mapped by high resolution oligo-arrays. Results: Submicroscopic genomic imbalances affecting the 1q21.1, 17q12, 22q11.21, and Xq21.31 chromosome regions were detected in four probands. Presence of the alterations in the normal mother of one patient suggests incomplete penetrance and/or variable expressivity. Conclusion: 4 of the 14 patients (29%) were found to have cryptic genomic alterations. The imbalances on 22q11.21 support recent findings by us and others that alterations in this chromosome region may result in impairment of mullerian duct development. The remaining imbalances indicate involvement of previously unknown chromosome regions in MA, and point specifically to LHX1 and KLHL4 as candidate genes.

New Trypanosoma (Duttonella) vivax genotypes from tsetse flies in East Africa

Salivarian trypanosomes pose a substantial threat to livestock, but their full diversity is not known. To survey trypanosomes carried by tsetse in Tanzania, DNA samples from infected proboscides of Glossina pallidipes and G. swynnertoni were identified using fluorescent fragment length barcoding (FFLB), which discriminates species by size polymorphisms in multiple regions of the ribosomal RNA locus. FELLB identified the trypanosomes in 65 of 105 (61.9%) infected proboscides, revealing 9 mixed infections. Of 7 different FFLB profiles, 2 were similar but not identical to reference West African Trypanosoma vivax; 5 other profiles belonged to known species also identified in fly midguts. Phylogenetic analysis of the glycosomal glyceraldehyde phosphate dehydrogenase gene revealed that the Tanzanian T. vivax samples fell into 2 distinct groups, both outside the main chide of African and South American T. vivax. These new T. vivax genotypes were common and widespread in tsetse in Tanzania. The T. brucei-like trypanosome previously described from tsetse midguts was also found in 2 proboscides, demonstrating a salivarian transmission route. Investigation of mammalian host range and pathogenicity will reveal the importance of these new trypanosomes for the epidemiology and control of animal trypanosomiasis in East Africa.

Novel Statistical Methods in Quantitative Genetics : Modeling Genetic Variance for Quantitative Trait Loci Mapping and Genomic Evaluation

This thesis develops and evaluates statistical methods for different types of genetic analyses, including quantitative trait loci (QTL) analysis, genome-wide association study (GWAS), and genomic evaluation. The main contribution of the thesis is to provide novel insights in modeling genetic variance, especially via random effects models. In variance component QTL analysis, a full likelihood model accounting for uncertainty in the identity-by-descent (IBD) matrix was developed. It was found to be able to correctly adjust the bias in genetic variance component estimation and gain power in QTL mapping in terms of precision.  Double hierarchical generalized linear models, and a non-iterative simplified version, were implemented and applied to fit data of an entire genome. These whole genome models were shown to have good performance in both QTL mapping and genomic prediction. A re-analysis of a publicly available GWAS data set identified significant loci in Arabidopsis that control phenotypic variance instead of mean, which validated the idea of variance-controlling genes.  The works in the thesis are accompanied by R packages available online, including a general statistical tool for fitting random effects models (hglm), an efficient generalized ridge regression for high-dimensional data (bigRR), a double-layer mixed model for genomic data analysis (iQTL), a stochastic IBD matrix calculator (MCIBD), a computational interface for QTL mapping (qtl.outbred), and a GWAS analysis tool for mapping variance-controlling loci (vGWAS).

Xylella fastidiosa: An in vivo system to study possible survival strategies within citrus xylem vessels based on global gene expression analysis

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Xylella fastidiosa gene expression analysis by DNA microarrays

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Genetic linkage maps of chicken chromosomes 6, 7, 8, 11 and 13 from a Brazilian resource population

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Genomic content and new insights on the origin of the B chromosome of the cichlid fish Astatotilapia latifasciata

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Comparative RH Maps of the River Buffalo and Bovine Y Chromosomes

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)