Molecular and functional evolution of the fungal diterpene synthase genes.

BACKGROUND: Terpenes represent one of the largest and most diversified families of natural compounds and are used in numerous industrial applications. Terpene synthase (TPS) genes originated in bacteria as diterpene synthase (di-TPS) genes. They are also found in plant and fungal genomes. The recent availability of a large number of fungal genomes represents an opportunity to investigate how genes involved in diterpene synthesis were acquired by fungi, and to assess the consequences of this process on the fungal metabolism. RESULTS: In order to investigate the origin of fungal di-TPS, we implemented a search for potential fungal di-TPS genes and identified their presence in several unrelated Ascomycota and Basidiomycota species. The fungal di-TPS phylogenetic tree is function-related but is not associated with the phylogeny based on housekeeping genes. The lack of agreement between fungal and di-TPS-based phylogenies suggests the presence of Horizontal Gene Transfer (HGTs) events. Further evidence for HGT was provided by conservation of synteny of di-TPS and neighbouring genes in distantly related fungi. CONCLUSIONS: The results obtained here suggest that fungal di-TPSs originated from an ancient HGT event of a single di-TPS gene from a plant to a fungus in Ascomycota. In fungi, these di-TPSs allowed for the formation of clusters consisting in di-TPS, GGPPS and P450 genes to create functional clusters that were transferred between fungal species, producing diterpenes acting as hormones or toxins, thus affecting fungal development and pathogenicity.

An evaluation of prior influence on the predictive ability of Bayesian model averaging.

Model averaging is gaining popularity among ecologists for making inference and predictions. Methods for combining models include Bayesian model averaging (BMA) and Akaike's Information Criterion (AIC) model averaging. BMA can be implemented with different prior model weights, including the Kullback-Leibler prior associated with AIC model averaging, but it is unclear how the prior model weight affects model results in a predictive context. Here, we implemented BMA using the Bayesian Information Criterion (BIC) approximation to Bayes factors for building predictive models of bird abundance and occurrence in the Chihuahuan Desert of New Mexico. We examined how model predictive ability differed across four prior model weights, and how averaged coefficient estimates, standard errors and coefficients' posterior probabilities varied for 16 bird species. We also compared the predictive ability of BMA models to a best single-model approach. Overall, Occam's prior of parsimony provided the best predictive models. In general, the Kullback-Leibler prior, however, favored complex models of lower predictive ability. BMA performed better than a best single-model approach independently of the prior model weight for 6 out of 16 species. For 6 other species, the choice of the prior model weight affected whether BMA was better than the best single-model approach. Our results demonstrate that parsimonious priors may be favorable over priors that favor complexity for making predictions. The approach we present has direct applications in ecology for better predicting patterns of species' abundance and occurrence.

EVALUATING THE THRESHOLD DENSITY HYPOTHESIS FOR MOOSE (ALCES ALCES), WHITE-TAILED DEER (ODOCOILEUS VIRGINIANUS), AND PARELAPHOSTRONGYLUS TENUIS.

Despite the importance of the Parelaphostrongylus tenuis infection for moose (Alces alces) and white-tailed deer (Odocoileus virginianus) management, only one peer-reviewed study has evaluated the relationship between deer and moose densities and the potential for parasite-mediated competition between the species. Using polynomial-regression modeling, that study identified a deer-density threshold above which moose populations declined; however, the nature of the data and apparent outliers suggests the approach used to develop that threshold may not have been appropriate. We used the data from the original study to test whether alternative models, including linear models and negative binomial models would be less sensitive to outliers and could better explain that relationship. We found no evidence that moose density decreases as deer density increases. We concluded that, although the proposed moose-deer-P. tenuis relationship could be partially density dependent, additional factors, such as frequency dependence of disease transmission, gastropod abundance, and shared use of resources by moose and deer should also be considered.

Genome-wide computational prediction of tandem gene arrays: application in yeasts.

BACKGROUND: This paper describes an efficient in silico method for detecting tandem gene arrays (TGAs) in fully sequenced and compact genomes such as those of prokaryotes or unicellular eukaryotes. The originality of this method lies in the search of protein sequence similarities in the vicinity of each coding sequence, which allows the prediction of tandem duplicated gene copies independently of their functionality. RESULTS: Applied to nine hemiascomycete yeast genomes, this method predicts that 2% of the genes are involved in TGAs and gene relics are present in 11% of TGAs. The frequency of TGAs with degenerated gene copies means that a significant fraction of tandem duplicated genes follows the birth-and-death model of evolution. A comparison of sequence identity distributions between sets of homologous gene pairs shows that the different copies of tandem arrayed paralogs are less divergent than copies of dispersed paralogs in yeast genomes. It suggests that paralogs included in tandem structures are more recent or more subject to the gene conversion mechanism than other paralogs. CONCLUSION: The method reported here is a useful computational tool to provide a database of TGAs composed of functional or nonfunctional gene copies. Such a database has obvious applications in the fields of structural and comparative genomics. Notably, a detailed study of the TGA catalog will make it possible to tackle the fundamental questions of the origin and evolution of tandem gene clusters.

Complete mitochondrial genome sequence of the yeast Pichia farinosa and comparative analysis of closely related species.

Yeasts of the Pichia genus have been isolated from different natural environments. Phylogenies based on multigene sequence analysis have shown that the genus is polyphyletic. Some species of this genus are member of the CTG group. In order to have a better insight into the relationship among species assigned to the yeast genera Pichia into the CTG group, we first sequenced the mitochondrial genome of the osmotolerant yeast Pichia farinosa. We then compared this genome with mitochondrial genomes of yeasts of the CTG group. The P. farinosa mitochondrial DNA is a circular-mapping genome of 32,065 bp, which contains 43 genes transcribed from both strands. It contains a complete set of tRNAs, the small and the large rRNAs, as well as 14 protein-coding genes. Yeasts of the CTG group contain the same core of mitochondrial genes. Phylogenetic analysis based on mitochondrial sequences clearly shows that the CTG group is divided into two distinct clades: the first one contains diploid Candida species, whereas the second mainly contains haploid Pichia species. Moreover, this analysis provides clear evidence that Pichia farinosa and Pichia sorbitophila, which were known to be unique species, are two distinct species.

Population reduction by hunting helps control human-wildlife conflicts for a species that is a conservation success story.

Among the world's large Carnivores, American black bears (Ursus americanus) are the foremost conservation success story. Populations have been expanding across North America because the species is adaptable and tolerant of living near people, and because management agencies in the U.S. and Canada controlled hunting and other human-sources of mortality. As a result, human-black bear conflicts (damage to property, general nuisance, threat to human safety) have dramatically increased in some areas, making it urgently important to develop and deploy a variety of mitigation tools. Previous studies claimed that legal hunting did not directly reduce conflicts, but they did not evaluate whether hunting controlled conflicts via management of population size. Here, we compared temporal patterns of phoned-in complaints about black bears (total ~63,500) in Minnesota, USA, over 4 decades to corresponding bear population estimates: both doubled during the first decade. We also quantified natural bear foods, and found that large year-to-year fluctuations affected numbers of complaints; however, since this variation is due largely to weather, this factor cannot be managed. Complaints fell sharply when the management agency (1) shifted more responsibility for preventing and mitigating conflicts to the public; and (2) increased hunting pressure to reduce the bear population. This population reduction was more extreme than intended, however, and after hunting pressure was curtailed, population regrowth was slower than anticipated; consequently both population size and complaints remained at relatively low levels statewide for 2 decades (although with local hotspots). These long-term data indicated that conflicts can be kept in tolerable bounds by managing population size through hunting; but due to the bluntness of this instrument and deficiencies and uncertainties in monitoring and manipulating populations, it is wiser to maintain a population at a level where conflicts are socially-acceptable than try to reduce it once it is well beyond that point.

An evolutionary scenario for one of the largest yeast gene families.

The DUP gene family of Saccharomyces cerevisiae comprises 23 members that can be divided into two subfamilies--DUP240 and DUP380. The location of the DUP loci suggests that at least three mechanisms were responsible for their genomic dispersion: nonreciprocal translocation at chromosomal ends, tandem duplication and Ty-associated duplication. The data we present here suggest that these nonessential genes encode proteins that facilitate membrane trafficking processes. Dup240 proteins have three conserved domains (C1, C2 and C3) and two predicted transmembrane segments (H1 and H2). A direct repetition of the C1-H1-H2-C2 module is observed in Dup380p sequences. In this article, we propose an evolutionary model to account for the emergence of the two gene subfamilies.

Modeling habitat suitability for Greater Rheas based on satellite image texture.

Many wild species are affected by human activities occurring at broad spatial scales. For instance, in South America, habitat loss threatens Greater Rhea (Rhea americana) populations, making it important to model and map their habitat to better target conservation efforts. Spatially explicit habitat modeling is a powerful approach to understand and predict species occurrence and abundance. One problem with this approach is that commonly used land cover classifications do not capture the variability within a given land cover class that might constitute important habitat attribute information. Texture measures derived from remote sensing images quantify the variability in habitat features among and within habitat types; hence they are potentially a powerful tool to assess species-habitat relationships. Our goal was to explore the utility of texture measures for habitat modeling and to develop a habitat suitability map for Greater Rheas at the home range level in grasslands of Argentina. Greater Rhea group size obtained from aerial surveys was regressed against distance to roads, houses, and water, and land cover class abundance (dicotyledons, crops, grassland, forest, and bare soil), normalized difference vegetation index (NDVI), and selected first- and second-order texture measures derived from Landsat Thematic Mapper (TM) imagery. Among univariate models, Rhea group size was most strongly positively correlated with texture variables derived from near infrared reflectance measurement (TM band 4). The best multiple regression models explained 78% of the variability in Greater Rhea group size. Our results suggest that texture variables captured habitat heterogeneity that the conventional land cover classification did not detect. We used Greater Rhea group size as an indicator of habitat suitability; we categorized model output into different habitat quality classes. Only 16% of the study area represented high-quality habitat for Greater Rheas (group size > or =15). Our results stress the potential of image texture to capture within-habitat variability in habitat assessments, and the necessity to preserve the remaining natural habitat for Greater Rheas.

Spatio-temporal variation in prevalence and intensity of trematodes responsible for waterfowl die-offs in faucet snail-infested waterbodies of Minnesota, USA.

Several non-native trematodes hosted by the invasive Eurasian faucet snail, Bithynia tentaculata, have been causing die-offs of waterfowl in the Midwestern United States and Canada for several decades. Because of the potential implications of these die-offs on waterfowl in non-native settings, it is necessary to better understand the trematodes that cause the die-offs. Here, we studied the spatio-temporal dynamics of two trematodes, Cyathocotyle bushiensis and Sphaeridiotrema spp., known to infect waterfowl in northern Minnesota, USA, via their intermediate host, the faucet snail (Bithynia tentaculata). We studied prevalence (% of snails infected within a sample) and intensity (mean number of parasites per infected snail within a sample) of faucet snail infection with these two trematodes in small lakes, large lakes, ponds, and rivers in northern Minnesota in the spring, summer, and fall of 2011-2013. We tested whether parasite prevalence and infection intensity could be explained spatially (as a function of the abundance of faucet snails, average snail size, water depth, and proximity to known waterfowl groups) and temporally (across years and seasons) using generalized estimating equation models. The spatial and temporal patterns we observed varied within and among waterbodies. For both parasite species, parasite prevalence and intensity of infection were consistently higher in samples with larger snails and in deeper portions of the waterbodies. In Lake Winnibigoshish, prevalence was lower farther from the large waterfowl groups we observed, but the abundance of snails in a sample had no effect on prevalence or intensity of infection. Our findings help improve understanding of this multi-species system, but also illustrate the complexity of modeling the spatial and temporal dynamics of infections in waterbodies that are so variable in size, shape, waterfowl use, and function.

Expansion and contraction of the DUP240 multigene family in Saccharomyces cerevisiae populations.

The influence of duplicated sequences on chromosomal stability is poorly understood. To characterize chromosomal rearrangements involving duplicated sequences, we compared the organization of tandem repeats of the DUP240 gene family in 15 Saccharomyces cerevisiae strains of various origins. The DUP240 gene family consists of 10 members of unknown function in the reference strain S288C. Five DUP240 paralogs on chromosome I and two on chromosome VII are arranged as tandem repeats that are highly polymorphic in copy number and sequence. We characterized DNA sequences that are likely involved in homologous or nonhomologous recombination events and are responsible for intra- and interchromosomal rearrangements that cause the creation and disappearance of DUP240 paralogs. The tandemly repeated DUP240 genes seem to be privileged sites of gene birth and death.

Modelling avian biodiversity using raw, unclassified satellite imagery.

Applications of remote sensing for biodiversity conservation typically rely on image classifications that do not capture variability within coarse land cover classes. Here, we compare two measures derived from unclassified remotely sensed data, a measure of habitat heterogeneity and a measure of habitat composition, for explaining bird species richness and the spatial distribution of 10 species in a semi-arid landscape of New Mexico. We surveyed bird abundance from 1996 to 1998 at 42 plots located in the McGregor Range of Fort Bliss Army Reserve. Normalized Difference Vegetation Index values of two May 1997 Landsat scenes were the basis for among-pixel habitat heterogeneity (image texture), and we used the raw imagery to decompose each pixel into different habitat components (spectral mixture analysis). We used model averaging to relate measures of avian biodiversity to measures of image texture and spectral mixture analysis fractions. Measures of habitat heterogeneity, particularly angular second moment and standard deviation, provide higher explanatory power for bird species richness and the abundance of most species than measures of habitat composition. Using image texture, alone or in combination with other classified imagery-based approaches, for monitoring statuses and trends in biological diversity can greatly improve conservation efforts and habitat management.

The complete genome of Blastobotrys (Arxula) adeninivorans LS3 - a yeast of biotechnological interest.

BACKGROUND: The industrially important yeast Blastobotrys (Arxula) adeninivorans is an asexual hemiascomycete phylogenetically very distant from Saccharomyces cerevisiae. Its unusual metabolic flexibility allows it to use a wide range of carbon and nitrogen sources, while being thermotolerant, xerotolerant and osmotolerant. RESULTS: The sequencing of strain LS3 revealed that the nuclear genome of A. adeninivorans is 11.8 Mb long and consists of four chromosomes with regional centromeres. Its closest sequenced relative is Yarrowia lipolytica, although mean conservation of orthologs is low. With 914 introns within 6116 genes, A. adeninivorans is one of the most intron-rich hemiascomycetes sequenced to date. Several large species-specific families appear to result from multiple rounds of segmental duplications of tandem gene arrays, a novel mechanism not yet described in yeasts. An analysis of the genome and its transcriptome revealed enzymes with biotechnological potential, such as two extracellular tannases (Atan1p and Atan2p) of the tannic-acid catabolic route, and a new pathway for the assimilation of n-butanol via butyric aldehyde and butyric acid. CONCLUSIONS: The high-quality genome of this species that diverged early in Saccharomycotina will allow further fundamental studies on comparative genomics, evolution and phylogenetics. Protein components of different pathways for carbon and nitrogen source utilization were identified, which so far has remained unexplored in yeast, offering clues for further biotechnological developments. In the course of identifying alternative microorganisms for biotechnological interest, A. adeninivorans has already proved its strengthened competitiveness as a promising cell factory for many more applications.

Complete DNA sequence of Kuraishia capsulata illustrates novel genomic features among budding yeasts (Saccharomycotina).

The numerous yeast genome sequences presently available provide a rich source of information for functional as well as evolutionary genomics but unequally cover the large phylogenetic diversity of extant yeasts. We present here the complete sequence of the nuclear genome of the haploid-type strain of Kuraishia capsulata (CBS1993(T)), a nitrate-assimilating Saccharomycetales of uncertain taxonomy, isolated from tunnels of insect larvae underneath coniferous barks and characterized by its copious production of extracellular polysaccharides. The sequence is composed of seven scaffolds, one per chromosome, totaling 11.4 Mb and containing 6,029 protein-coding genes, ~13.5% of which being interrupted by introns. This GC-rich yeast genome (45.7%) appears phylogenetically related with the few other nitrate-assimilating yeasts sequenced so far, Ogataea polymorpha, O. parapolymorpha, and Dekkera bruxellensis, with which it shares a very reduced number of tRNA genes, a novel tRNA sparing strategy, and a common nitrate assimilation cluster, three specific features to this group of yeasts. Centromeres were recognized in GC-poor troughs of each scaffold. The strain bears MAT alpha genes at a single MAT locus and presents a significant degree of conservation with Saccharomyces cerevisiae genes, suggesting that it can perform sexual cycles in nature, although genes involved in meiosis were not all recognized. The complete absence of conservation of synteny between K. capsulata and any other yeast genome described so far, including the three other nitrate-assimilating species, validates the interest of this species for long-range evolutionary genomic studies among Saccharomycotina yeasts.

Insights into the life cycle of yeasts from the CTG clade revealed by the analysis of the Millerozyma (Pichia) farinosa species complex.

Among ascomycetous yeasts, the CTG clade is so-called because its constituent species translate CTG as serine instead of leucine. Though the biology of certain pathogenic species such as Candida albicans has been much studied, little is known about the life cycles of non-pathogen species of the CTG clade. Taking advantage of the recently obtained sequence of the biotechnological Millerozyma (Pichiasorbitophila) farinosa strain CBS 7064, we used MLST to better define phylogenic relationships between most of the Millerozyma farinosa strains available in public collections. This led to the constitution of four phylogenetic clades diverging from 8% to 15% at the DNA level and possibly constituting a species complex (M. farinosa) and to the proposal of two new species:Millerozyma miso sp. nov. CBS 2004(T) ( = CLIB 1230(T)) and Candida pseudofarinosa sp. nov.NCYC 386(T)( = CLIB 1231(T)). Further analysis showed that M. farinosa isolates exist as haploid and inter-clade hybrids. Despite the sequence divergence between the clades, secondary contacts after reproductive isolation were evidenced, as revealed by both introgression and mitochondria transfer between clades. We also showed that the inter-clade hybrids do sporulate to generate mainly viable vegetative diploid spores that are not the result of meiosis, and very rarely aneuploid spores possibly through the loss of heterozygosity during sporulation. Taken together, these results show that in this part of the CTG clade, non-Mendelian genetic exchanges occur at high rates through hybridization between divergent strains from distinct clades and subsequent massive loss of heterozygosity. This combination of mechanisms could constitute an alternative sexuality leading to an unsuspected biodiversity.

Pichia sorbitophila, an Interspecies Yeast Hybrid, Reveals Early Steps of Genome Resolution After Polyploidization.

Polyploidization is an important process in the evolution of eukaryotic genomes, but ensuing molecular mechanisms remain to be clarified. Autopolyploidization or whole-genome duplication events frequently are resolved in resulting lineages by the loss of single genes from most duplicated pairs, causing transient gene dosage imbalance and accelerating speciation through meiotic infertility. Allopolyploidization or formation of interspecies hybrids raises the problem of genetic incompatibility (Bateson-Dobzhansky-Muller effect) and may be resolved by the accumulation of mutational changes in resulting lineages. In this article, we show that an osmotolerant yeast species, Pichia sorbitophila, recently isolated in a concentrated sorbitol solution in industry, illustrates this last situation. Its genome is a mosaic of homologous and homeologous chromosomes, or parts thereof, that corresponds to a recently formed hybrid in the process of evolution. The respective parental contributions to this genome were characterized using existing variations in GC content. The genomic changes that occurred during the short period since hybrid formation were identified (e.g., loss of heterozygosity, unilateral loss of rDNA, reciprocal exchange) and distinguished from those undergone by the two parental genomes after separation from their common ancestor (i.e., NUMT (NUclear sequences of MiTochondrial origin) insertions, gene acquisitions, gene location movements, reciprocal translocation). We found that the physiological characteristics of this new yeast species are determined by specific but unequal contributions of its two parents, one of which could be identified as very closely related to an extant Pichia farinosa strain.

Tandem gene arrays, plastic chromosomal organizations.

This short article presents an overview of tandem gene arrays (TGAs) in hemiascomycete yeasts. In silico and in vivo analyses are combined to address structural, functional and evolutionary aspects of these particular chromosomal structures. Genomic instability of TGAs is discussed. We conclude that TGAs are generally dynamic regions of the genome in that they are the seats of chromosomal rearrangement events. In addition, they are often breeding grounds of new genes for a rapid adaptation of cells to demands of the environment.

Paleogenomics or the search for remnant duplicated copies of the yeast DUP240 gene family in intergenic areas.

Duplication, resulting in gene redundancy, is well known to be a driving force of evolutionary change. Gene families are therefore useful targets for approaching genome evolution. To address the gene death process, we examined the fate of the 10-member-large S288C DUP240 family in 15 Saccharomyces cerevisiae strains. Using an original three-step method of analysis reported here, both slightly and highly degenerate DUP240 copies, called pseudo-open reading frames (ORFs) and relics, respectively, were detected in strain S288C. It was concluded that two previously annotated ORFs correspond, in fact, to pseudo-ORFs and three additional relics were identified in intergenic areas. Comparative intraspecies analysis of these degenerate DUP240 loci revealed that the two pseudo-ORFs are present in a nondegenerate state in some other strains. This suggests that within a given gene family different loci are the target of the gene erasure process, which is therefore strain dependent. Besides, the variable positions observed indicate that the relic sequence may diverge faster than the flanking regions. All in all, this study shows that short conserved protein motifs provide a useful tool for detecting and accurately mapping degenerate gene remnants. The present results also highlight the strong contribution of comparative genomics for gene relic detection because the possibility of finding short conserved protein motifs in intergenic regions (IRs) largely depends on the choice of the most closely related paralog or ortholog. By mapping new genetic components in previously annotated IRs, our study constitutes a further refinement step in the crucial stage of genome annotation and provides a strategy for retracing ancient chromosomal reshaping events and, hence, for deciphering genome history.