Cis-effects on gene expression in the human prenatal brain associated with genetic risk for neuropsychiatric disorders.

The majority of common risk alleles identified for neuropsychiatric disorders reside in noncoding regions of the genome and are therefore likely to impact gene regulation. However, the genes that are primarily affected and the nature and developmental timing of these effects remain unclear. Given the hypothesized role for early neurodevelopmental processes in these conditions, we here define genetic predictors of gene expression in the human fetal brain with which we perform transcriptome-wide association studies (TWASs) of attention deficit hyperactivity disorder (ADHD), autism spectrum disorder, bipolar disorder, major depressive disorder, and schizophrenia. We identify prenatal cis-regulatory effects on 63 genes and 166 individual transcripts associated with genetic risk for these conditions. We observe pleiotropic effects of expression predictors for a number of genes and transcripts, including those of decreased DDHD2 expression in association with risk for schizophrenia and bipolar disorder, increased expression of a ST3GAL3 transcript with risk for schizophrenia and ADHD, and increased expression of an XPNPEP3 transcript with risk for schizophrenia, bipolar disorder, and major depression. For the protocadherin alpha cluster genes PCDHA7 and PCDHA8, we find that predictors of low expression are associated with risk for major depressive disorder while those of higher expression are associated with risk for schizophrenia. Our findings support a role for altered gene regulation in the prenatal brain in susceptibility to various neuropsychiatric disorders and prioritize potential risk genes for further neurobiological investigation.

Assessing population structure and host specialization in lichenized cyanobacteria.

Coevolutionary theory predicts that the distribution of obligately symbiotic organisms will be determined by the dispersal ability and ecological range of both partners. We examined this prediction for lichen-forming fungi that form obligate symbioses with cyanobacteria. We compared genotypes of both partners of 250 lichens collected at multiple spatial scales in British Columbia, Canada. Multilocus sequence data collected from a subset of 128 of the specimens were used to determine the degree of recombination within the cyanobacterial populations. We found that six distinct clusters of cyanobacterial genotypes are distributed throughout the known global phylogeny of the genus Nostoc, and that each appears to be evolving clonally. Fungal specialization is high, with each species associating with either one or two of the cyanobacterial clusters, while cyanobacterial specialization varies, with clusters associating with between one and 12 different fungal species. Specialization also varies geographically, with some combinations restricted to a single site despite the availability of both partners elsewhere. Photobiont association patterns are determined by a combination of genetically based specificity, spatial population structure, and ecological factors and cannot be easily predicted by photobiont dispersal syndromes.

Extensive remodeling of the Pseudomonas syringae pv. avellanae type III secretome associated with two independent host shifts onto hazelnut.

BACKGROUND: Hazelnut (Corylus avellana) decline disease in Greece and Italy is caused by the convergent evolution of two distantly related lineages of Pseudomonas syringae pv. avellanae (Pav). We sequenced the genomes of three Pav isolates to determine if their convergent virulence phenotype had a common genetic basis due to either genetic exchange between lineages or parallel evolution. RESULTS: We found little evidence for horizontal transfer (recombination) of genes between Pav lineages, but two large genomic islands (GIs) have been recently acquired by one of the lineages. Evolutionary analyses of the genes encoding type III secreted effectors (T3SEs) that are translocated into host cells and are important for both suppressing and eliciting defense responses show that the two Pav lineages have dramatically different T3SE profiles, with only two shared putatively functional T3SEs. One Pav lineage has undergone unprecedented secretome remodeling, including the acquisition of eleven new T3SEs and the loss or pseudogenization of 15, including five of the six core T3SE families that are present in the other Pav lineage. Molecular dating indicates that divergence within both of the Pav lineages predates their observation in the field. This suggest that both Pav lineages have been cryptically infecting hazelnut trees or wild relatives for many years, and that the emergence of hazelnut decline in the 1970s may have been due to changes in agricultural practice. CONCLUSIONS: These data show that divergent lineages of P. syringae can converge on identical disease etiology on the same host plant using different virulence mechanisms and that dramatic shifts in the arsenal of T3SEs can accompany disease emergence.

HopAS1 recognition significantly contributes to Arabidopsis nonhost resistance to Pseudomonas syringae pathogens.

• Plant immunity is activated by sensing either conserved microbial signatures, called pathogen/microbe-associated molecular patterns (P/MAMPs), or specific effectors secreted by pathogens. However, it is not known why most microbes are nonpathogenic in most plant species. • Nonhost resistance (NHR) consists of multiple layers of innate immunity and protects plants from the vast majority of potentially pathogenic microbes. Effector-triggered immunity (ETI) has been implicated in race-specific disease resistance. However, the role of ETI in NHR is unclear. • Pseudomonas syringae pv. tomato (Pto) T1 is pathogenic in tomato (Solanum lycopersicum) yet nonpathogenic in Arabidopsis. Here, we show that, in addition to the type III secretion system (T3SS)-dependent effector (T3SE) avrRpt2, a second T3SE of Pto T1, hopAS1, triggers ETI in nonhost Arabidopsis. • hopAS1 is broadly present in P. syringae strains, contributes to virulence in tomato, and is quantitatively required for Arabidopsis NHR to Pto T1. Strikingly, all tested P. syringae strains that are pathogenic in Arabidopsis carry truncated hopAS1 variants of forms, demonstrating that HopAS1-triggered immunity plays an important role in Arabidopsis NHR to a broad-range of P. syringae strains.

Identification and Phenotypic Characterization of Hsp90 Phosphorylation Sites That Modulate Virulence Traits in the Major Human Fungal Pathogen Candida albicans.

The highly conserved, ubiquitous molecular chaperone Hsp90 is a key regulator of cellular proteostasis and environmental stress responses. In human pathogenic fungi, which kill more than 1.6 million patients each year worldwide, Hsp90 governs cellular morphogenesis, drug resistance, and virulence. Yet, our understanding of the regulatory mechanisms governing fungal Hsp90 function remains sparse. Post-translational modifications are powerful components of nature's toolbox to regulate protein abundance and function. Phosphorylation in particular is critical in many cellular signaling pathways and errant phosphorylation can have dire consequences for the cell. In the case of Hsp90, phosphorylation affects its stability and governs its interactions with co-chaperones and clients. Thereby modulating the cell's ability to cope with environmental stress. Candida albicans, one of the leading human fungal pathogens, causes ~750,000 life-threatening invasive infections worldwide with unacceptably high mortality rates. Yet, it remains unknown if and how Hsp90 phosphorylation affects C. albicans virulence traits. Here, we show that phosphorylation of Hsp90 is critical for expression of virulence traits. We combined proteomics, molecular evolution analyses and structural modeling with molecular biology to characterize the role of Hsp90 phosphorylation in this non-model pathogen. We demonstrated that phosphorylation negatively affects key virulence traits, such as the thermal stress response, morphogenesis, and drug susceptibility. Our results provide the first record of a specific Hsp90 phosphorylation site acting as modulator of fungal virulence. Post-translational modifications of Hsp90 could prove valuable in future exploitations as antifungal drug targets.

Tobacco Hornworm (Manduca sexta) caterpillars as a novel host model for the study of fungal virulence and drug efficacy.

The two leading yeast pathogens of humans, Candida albicans and Cryptococcus neoformans, cause systemic infections in >1.4 million patients worldwide with mortality rates approaching 75%. It is thus imperative to study fungal virulence mechanisms, efficacy of antifungal drugs, and host response pathways. While this is commonly done in mammalian models, which are afflicted by ethical and practical concerns, invertebrate models, such as wax moth larvae and nematodes have been introduced over the last two decades. To complement existing invertebrate host models, we developed fifth instar caterpillars of the Tobacco Hornworm moth Manduca sexta as a novel host model. These caterpillars can be maintained at 37°C, are suitable for injections with defined amounts of yeast cells, and are susceptible to the most threatening yeast pathogens, including C. albicans, C. neoformans, C. auris, and C. glabrata. Importantly, fungal burden can be assessed daily throughout the course of infection in a single caterpillar's feces and hemolymph. Infected caterpillars can be rescued by treatment with antifungal drugs. Notably, these animals are large enough for weight to provide a reliable and reproducible measure of fungal disease and to facilitate host tissue-specific expression analyses. M. sexta caterpillars combine a suite of parameters that make them suitable for the study of fungal virulence.

Novel Insight Into the Etiology of Autism Spectrum Disorder Gained by Integrating Expression Data With Genome-wide Association Statistics.

BACKGROUND: A recent genome-wide association study (GWAS) of autism spectrum disorder (ASD) (ncases = 18,381, ncontrols = 27,969) has provided novel opportunities for investigating the etiology of ASD. Here, we integrate the ASD GWAS summary statistics with summary-level gene expression data to infer differential gene expression in ASD, an approach called transcriptome-wide association study (TWAS). METHODS: Using FUSION software, ASD GWAS summary statistics were integrated with predictors of gene expression from 16 human datasets, including adult and fetal brains. A novel adaptation of established statistical methods was then used to test for enrichment within candidate pathways and specific tissues and at different stages of brain development. The proportion of ASD heritability explained by predicted expression of genes in the TWAS was estimated using stratified linkage disequilibrium score regression. RESULTS: This study identified 14 genes as significantly differentially expressed in ASD, 13 of which were outside of known genome-wide significant loci (±500 kb). XRN2, a gene proximal to an ASD GWAS locus, was inferred to be significantly upregulated in ASD, providing insight into the functional consequence of this associated locus. One novel transcriptome-wide significant association from this study is the downregulation of PDIA6, which showed minimal evidence of association in the GWAS, and in gene-based analysis using MAGMA. Predicted gene expression in this study accounted for 13.0% of the total ASD single nucleotide polymorphism heritability. CONCLUSIONS: This study has implicated several genes as significantly up/downregulated in ASD, providing novel and useful information for subsequent functional studies. This study also explores the utility of TWAS-based enrichment analysis and compares TWAS results with a functionally agnostic approach.

Light-induced dynamic structural color by intracellular 3D photonic crystals in brown algae.

Natural photonic crystals are responsible for strong reflectance at selective wavelengths in different natural systems. We demonstrate that intracellular opal-like photonic crystals formed from lipids within photosynthetic cells produce vivid structural color in the alga Cystoseira tamariscifolia. The reflectance of the opaline vesicles is dynamically responsive to environmental illumination. The structural color is present in low light-adapted samples, whereas higher light levels produce a slow disappearance of the structural color such that it eventually vanishes completely. Once returned to low-light conditions, the color re-emerges. Our results suggest that these complex intracellular natural photonic crystals are responsive to environmental conditions, changing their packing structure reversibly, and have the potential to manipulate light for roles beyond visual signaling.

Expression quantitative trait loci in the developing human brain and their enrichment in neuropsychiatric disorders.

BACKGROUND: Genetic influences on gene expression in the human fetal brain plausibly impact upon a variety of postnatal brain-related traits, including susceptibility to neuropsychiatric disorders. However, to date, there have been no studies that have mapped genome-wide expression quantitative trait loci (eQTL) specifically in the human prenatal brain. RESULTS: We performed deep RNA sequencing and genome-wide genotyping on a unique collection of 120 human brains from the second trimester of gestation to provide the first eQTL dataset derived exclusively from the human fetal brain. We identify high confidence cis-acting eQTL at the individual transcript as well as whole gene level, including many mapping to a common inversion polymorphism on chromosome 17q21. Fetal brain eQTL are enriched among risk variants for postnatal conditions including attention deficit hyperactivity disorder, schizophrenia, and bipolar disorder. We further identify changes in gene expression within the prenatal brain that potentially mediate risk for neuropsychiatric traits, including increased expression of C4A in association with genetic risk for schizophrenia, increased expression of LRRC57 in association with genetic risk for bipolar disorder, and altered expression of multiple genes within the chromosome 17q21 inversion in association with variants influencing the personality trait of neuroticism. CONCLUSIONS: We have mapped eQTL operating in the human fetal brain, providing evidence that these confer risk to certain neuropsychiatric disorders, and identifying gene expression changes that potentially mediate susceptibility to these conditions.

Repeated evolution of net venation and fleshy fruits among monocots in shaded habitats confirms a priori predictions: evidence from an ndhF phylogeny.

We present a well-resolved, highly inclusive phylogeny for monocots, based on ndhF sequence variation, and use it to test a priori hypotheses that net venation and vertebrate-dispersed fleshy fruits should undergo concerted convergence, representing independent but often concurrent adaptations to shaded conditions. Our data demonstrate that net venation arose at least 26 times and was lost eight times over the past 90 million years; fleshy fruits arose at least 21 times and disappeared 11 times. Both traits show a highly significant pattern of concerted convergence (p<10(-9)), arising 16 times and disappearing four times in tandem. This phenomenon appears driven by even stronger tendencies for both traits to evolve in shade and be lost in open habitats (p<10(-13)-10(-29)). These patterns are among the strongest ever demonstrated for evolutionary convergence in individual traits and the predictability of evolution, and the strongest evidence yet uncovered for concerted convergence. The rate of adaptive shifts per taxon has declined exponentially over the past 90 million years, as expected when large-scale radiations fill adaptive zones.

Genomic and Gene-Expression Comparisons among Phage-Resistant Type-IV Pilus Mutants of Pseudomonas syringae pathovar phaseolicola.

Pseudomonas syringae pv. phaseolicola (Pph) is a significant bacterial pathogen of agricultural crops, and phage Φ6 and other members of the dsRNA virus family Cystoviridae undergo lytic (virulent) infection of Pph, using the type IV pilus as the initial site of cellular attachment. Despite the popularity of Pph/phage Φ6 as a model system in evolutionary biology, Pph resistance to phage Φ6 remains poorly characterized. To investigate differences between phage Φ6 resistant Pph strains, we examined genomic and gene expression variation among three bacterial genotypes that differ in the number of type IV pili expressed per cell: ordinary (wild-type), non-piliated, and super-piliated. Genome sequencing of non-piliated and super-piliated Pph identified few mutations that separate these genotypes from wild type Pph--and none present in genes known to be directly involved in type IV pilus expression. Expression analysis revealed that 81.1% of gene ontology (GO) terms up-regulated in the non-piliated strain were down-regulated in the super-piliated strain. This differential expression is particularly prevalent in genes associated with respiration--specifically genes in the tricarboxylic acid cycle (TCA) cycle, aerobic respiration, and acetyl-CoA metabolism. The expression patterns of the TCA pathway appear to be generally up and down-regulated, in non-piliated and super-piliated Pph respectively. As pilus retraction is mediated by an ATP motor, loss of retraction ability might lead to a lower energy draw on the bacterial cell, leading to a different energy balance than wild type. The lower metabolic rate of the super-piliated strain is potentially a result of its loss of ability to retract.

Next-generation genomics of Pseudomonas syringae.

The first wave of Pseudomonas syringae next-generation genomic studies has revealed insights into host-specific virulence and immunity, genome dynamics and evolution, and genetic and metabolic specialization. These studies have further enhanced our understanding of type III effector diversity, identified an atypical type III secretion system (T3SS) in a new clade of nonpathogenic P. syringae, identified metabolic pathways common to pathogens of woody hosts and revealed extensive genomic diversity among strains that infect common hosts. In general, these discoveries have illustrated the utility of draft genome sequencing for quickly and economically identifying candidate loci for more refined genetic and functional analyses.

Evolution of plant pathogenesis in Pseudomonas syringae: a genomics perspective.

The phytopathogenic bacterium Pseudomonas syringae causes serious diseases in a wide range of important crop plants, with recent severe outbreaks on the New Zealand kiwifruit crop and among British horse chestnut trees. Next-generation genome sequencing of over 25 new strains has greatly broadened our understanding of how this species adapts to a diverse range of plant hosts. Not unexpectedly, the genomes were found to be highly dynamic, and extensive polymorphism was found in the distribution of type III secreted effectors (T3SEs) and other virulence-associated genes, even among strains within the same pathovar. An underexplored area brought to light by these data is the specific metabolic adaptations required for growth on woody hosts. These studies provide a tremendous wealth of candidates for more refined functional characterization, which is greatly enhancing our ability to disentangle the web of host-pathogen interactions that determine disease outcomes.

Use of low-coverage, large-insert, short-read data for rapid and accurate generation of enhanced-quality draft Pseudomonas genome sequences.

Next-generation genomic technology has both greatly accelerated the pace of genome research as well as increased our reliance on draft genome sequences. While groups such as the Genomics Standards Consortium have made strong efforts to promote genome standards there is a still a general lack of uniformity among published draft genomes, leading to challenges for downstream comparative analyses. This lack of uniformity is a particular problem when using standard draft genomes that frequently have large numbers of low-quality sequencing tracts. Here we present a proposal for an "enhanced-quality draft" genome that identifies at least 95% of the coding sequences, thereby effectively providing a full accounting of the genic component of the genome. Enhanced-quality draft genomes are easily attainable through a combination of small- and large-insert next-generation, paired-end sequencing. We illustrate the generation of an enhanced-quality draft genome by re-sequencing the plant pathogenic bacterium Pseudomonas syringae pv. phaseolicola 1448A (Pph 1448A), which has a published, closed genome sequence of 5.93 Mbp. We use a combination of Illumina paired-end and mate-pair sequencing, and surprisingly find that de novo assemblies with 100x paired-end coverage and mate-pair sequencing with as low as low as 2-5x coverage are substantially better than assemblies based on higher coverage. The rapid and low-cost generation of large numbers of enhanced-quality draft genome sequences will be of particular value for microbial diagnostics and biosecurity, which rely on precise discrimination of potentially dangerous clones from closely related benign strains.

Assessing reproductive isolation in highly diverse communities of the lichen-forming fungal genus peltigera.

The lichen-forming fungal genus Peltigera includes a number of species that are extremely widespread, both geographically and ecologically. However, morphological variability has lead to doubts about the distinctness of some species, and it has been suggested that hybridization is common in nature. We examined species boundaries by looking for evidence of hybridization and gene flow among seven described species collected at five sites in British Columbia, Canada. We found no evidence of gene flow or hybridization between described species, with fixed differences between species for two or more of the three loci examined. Reproductive isolation did not reflect a solely clonal mode of reproduction as there was evidence of ongoing gene flow within species. In addition, we found five undescribed species that were reproductively isolated, although there was evidence of ongoing or historical gene flow between two of the new species. These results indicate that the genus Peltigera is more diverse in western North America than originally perceived, and that morphological variability is due largely to the presence of undescribed species rather than hybridization or intraspecific variation.

Moths that vector a plant pathogen also transport endophytic fungi and mycoparasitic antagonists.

Claviceps paspali, a common fungal pathogen of Paspalum grasses, attracts moth vectors by producing sugary exudates in the grass florets it infects. These exudates also support mycoparasitic Fusarium species that may negatively influence C. paspali fitness. We examined the potential for moths on which C. paspali depends to also transmit mycoparasitic Fusarium and fungal endophytes, which inhabit asymptomatic plant tissue and may influence host susceptibility to pathogens. We quantified infections by C. paspali, Fusarium spp., and endophytic fungi associated with Paspalum spp. at focal sites in the southeastern USA and used data from the nuclear internal transcribed spacer (ITS rDNA) to compare communities of plant-associated and moth-borne fungi. ITS sequences of moth-borne fungi were identical to reference sequences of mycoparasitic Fusarium heterosporum and to three distinct endophytic fungi isolated from Paspalum species. Our results demonstrate an unexpected overlap of fungal communities between disparate locations and among plant species and plant tissues, and suggest an unexpected role of moths, which vector a plant pathogen, to transmit other guilds of fungi. In turn, the potential for insects to transmit plant pathogens as well as mycoparasites and endophytic fungi suggests complex interactions underlying a commonly observed grass-pathogen system.

Correlated evolution of self-incompatibility and clonal reproduction in Solanum (Solanaceae).

It has been suggested that clonality provides reproductive assurance in cross-fertilizing species subject to pollen limitation, relieving one of the main selective pressures favoring the evolution of self-fertilization. According to this hypothesis, cross-fertilizing species subject to pollen limitation should often be clonal. Here, we investigated the association between clonality and a genetic mechanism enforcing outcrossing, self-incompatibility, in Solanum (Solanaceae). We collected self-incompatibility and clonality information on 87 species, and looked for an association between these two traits. To account for the contribution of shared evolutionary history to this association, we incorporated phylogenetic information from chloroplast (NADH dehydrogenase subunit F) sequence data. We found that self-incompatibility is strongly associated with clonal reproduction: all self-incompatible species reproduce clonally, while the absence of clonality is widespread among self-compatible taxa. The observed correlation persists after taking into account shared phylogenetic history, assumptions about the evolutionary history of self-incompatibility, uncertainty associated with phylogeny estimation, and associations with life history (annual/perennial). Our results are consistent with the hypothesis that clonality provides reproductive assurance, and suggest that the consequences of clonal growth in the evolution of plant reproductive strategies may be more significant than previously thought.

Fungal community analysis by large-scale sequencing of environmental samples.

Fungi are an important and diverse component of soil communities, but these communities have proven difficult to study in conventional biotic surveys. We evaluated soil fungal diversity at two sites in a temperate forest using direct isolation of small-subunit and internal transcribed spacer (ITS) rRNA genes by PCR and high-throughput sequencing of cloned fragments. We identified 412 sequence types from 863 fungal ITS sequences, as well as 112 ITS sequences from other eukaryotic microorganisms. Equal proportions of Basidiomycota and Ascomycota sequences were present in both the ITS and small-subunit libraries, while members of other fungal phyla were recovered at much lower frequencies. Many sequences closely matched sequences from mycorrhizal, plant-pathogenic, and saprophytic fungi. Compositional differences were observed among samples from different soil depths, with mycorrhizal species predominating deeper in the soil profile and saprophytic species predominating in the litter layer. Richness was consistently lowest in the deepest soil horizon samples. Comparable levels of fungal richness have been observed following traditional specimen-based collecting and culturing surveys, but only after much more extensive sampling. The high rate at which new sequence types were recovered even after sampling 863 fungal ITS sequences and the dominance of fungi in our libraries relative to other eukaryotes suggest that the abundance and diversity of fungi in forest soils may be much higher than previously hypothesized. All sequences were deposited in GenBank, with accession numbers AY 969316 to AY 970290 for the ITS sequences and AY 969135 to AY 969315 for the SSU sequences.

Inference of higher-order relationships in the cycads from a large chloroplast data set.

We investigated higher-order relationships in the cycads, an ancient group of seed-bearing plants, by examining a large portion of the chloroplast genome from seven species chosen to exemplify our current understanding of taxonomic diversity in the order. The regions considered span approximately 13.5 kb of unaligned data per taxon, and comprise a diverse range of coding sequences, introns and intergenic spacers dispersed throughout the plastid genome. Our results provide substantial support for most of the inferred backbone of cycad phylogeny, and weak evidence that the sister-group of the cycads among living seed plants is Ginkgo biloba. Cycas (representing Cycadaceae) is the sister-group of the remaining cycads; Dioon is part of the next most basal split. Two of the three commonly recognized families of cycads (Zamiaceae and Stangeriaceae) are not monophyletic; Stangeria is embedded within Zamiaceae, close to Zamia and Ceratozamia, and not closely allied to the other genus of Stangeriaceae, Bowenia. In contrast to the other seed plants, cycad chloroplast genomes share two features with Ginkgo: a reduced rate of evolution and an elevated transition:transversion ratio. We demonstrate that the latter aspect of their molecular evolution is unlikely to have affected inference of cycad relationships in the context of seed-plant wide analyses.