FOAM (Functional Ontology Assignments for Metagenomes): a Hidden Markov Model (HMM) database with environmental focus

A new functional gene database, FOAM (Functional Ontology Assignments for Metagenomes), was developed to screen environmental metagenomic sequence datasets. FOAM provides a new functional ontology dedicated to classify gene functions relevant to environmental microorganisms based on Hidden Markov Models (HMMs). Sets of aligned protein sequences (i.e. ‘profiles’) were tailored to a large group of target KEGG Orthologs (KOs) from which HMMs were trained. The alignments were checked and curated to make them specific to the targeted KO. Within this process, sequence profiles were enriched with the most abundant sequences available to maximize the yield of accurate classifier models. An associated functional ontology was built to describe the functional groups and hierarchy. FOAM allows the user to select the target search space before HMM-based comparison steps and to easily organize the results into different functional categories and subcategories. FOAM is publicly available at http://portal.nersc.gov/project/m1317/FOAM/.     

The Three Receptor Tyrosine Kinases c-KIT,VEGFR2 and PDGFRα, Closely Spaced at 4q12, Show Increased Protein Expression in Triple-Negative Breast Cancer

Background

Triple-negative breast cancer (TNBC) is a heterogeneous subgroup of breast cancer with poor prognosis and no targeted therapy available. Receptor tyrosine kinases (RTKs) are emerging targets in anticancer therapy and many RTK-inhibiting drugs are currently being developed. The aim of this study was to elucidate if there is a correlation between the protein expression of three RTKs c-KIT, VEGFR2 and PDGFRα, their gene copy number, and prognosis in TNBC compared to non-TNBC.

Methods

Tumor tissue samples from patients diagnosed with primary breast cancer were stained with immunohistochemistry (IHC) for protein assessment, and with fluorescence in situ hybridization (FISH) for gene copy number determination. Breast cancer mortality (BCM), measured from the date of surgery to death, was used as endpoint.

Results

The cohort included 464 patients, out of which 34 (7.3%) had a TNBC. High expression of the three RTKs was more common in TNBC compared to non-TNBC: c-KIT 49% vs. 10% (P<0.001), PDGFRα 32% vs. 19% (P = 0.07) and VEGFR2 32% vs. 6% (P<0.001). The odds ratio (OR) of c-KIT, VEGFR2 and PDGFRα positivity, adjusted for tumor characteristics, was 6.8, 3.6 and 1.3 times higher for TNBC than for non-TNBC. 73.5% of the TNBC had high expression of at least one of the three investigated receptors, compared to 30.0% of the non-TNBC (P<0.001). Survival analysis showed no significant difference in BCM for TNBC patients with high vs. low c-KIT, PDGFRα or VEGFR2 protein expression. 193 (42%) tumors were evaluated with FISH. No correlation was seen between increased gene copy number and TNBC, or between increased gene copy number and high protein expression of the RTK.

Conclusion

c-KIT, VEGFR2 and PDGFRα show higher protein expression in TNBC compared to non-TNBC. Further investigation clarifying the importance of these RTKs in TNBC is encouraged, as they are possible targets for anticancer therapy.     

Reduced Motivation in the BACHD Rat Model of Huntington Disease Is Dependent on the Choice of Food Deprivation Strategy

Huntington disease (HD) is an inherited neurodegenerative disease characterized by motor, cognitive, psychiatric and metabolic symptoms. Animal models of HD show phenotypes that can be divided into similar categories, with the metabolic phenotype of certain models being characterized by obesity. Although interesting in terms of modeling metabolic symptoms of HD, the obesity phenotype can be problematic as it might confound the results of certain behavioral tests. This concerns the assessment of cognitive function in particular, as tests for such phenotypes are often based on food depriving the animals and having them perform tasks for food rewards. The BACHD rat is a recently established animal model of HD, and in order to ensure that behavioral characterization of these rats is done in a reliable way, a basic understanding of their physiology is needed. Here, we show that BACHD rats are obese and suffer from discrete developmental deficits. When assessing the motivation to lever push for a food reward, BACHD rats were found to be less motivated than wild type rats, although this phenotype was dependent on the food deprivation strategy. Specifically, the phenotype was present when rats of both genotypes were deprived to 85% of their respective free-feeding body weight, but not when deprivation levels were adjusted in order to match the rats'' apparent hunger levels. The study emphasizes the importance of considering metabolic abnormalities as a confounding factor when performing behavioral characterization of HD animal models.     

Genetic Inhibition of Phosphorylation of the Translation Initiation Factor eIF2α Does Not Block Aβ-Dependent Elevation of BACE1 and APP Levels or Reduce Amyloid Pathology in a Mouse Model of Alzheimer’s Disease

β-site amyloid precursor protein (APP) cleaving enzyme 1 (BACE1) initiates the production of β-amyloid (Aβ), the major constituent of amyloid plaques in Alzheimer’s disease (AD). BACE1 is elevated ∼2–3 fold in AD brain and is concentrated in dystrophic neurites near plaques, suggesting BACE1 elevation is Aβ−dependent. Previously, we showed that phosphorylation of the translation initiation factor eIF2α de-represses translation of BACE1 mRNA following stress such as energy deprivation. We hypothesized that stress induced by Aβ might increase BACE1 levels by the same translational mechanism involving eIF2α phosphorylation. To test this hypothesis, we used three different genetic strategies to determine the effects of reducing eIF2α phosphorylation on Aβ-dependent BACE1 elevation in vitro and in vivo: 1) a two-vector adeno-associated virus (AAV) system to express constitutively active GADD34, the regulatory subunit of PP1c eIF2α phosphatase; 2) a non-phosphorylatable eIF2α S51A knockin mutation; 3) a BACE1-YFP transgene lacking the BACE1 mRNA 5′ untranslated region (UTR) required for eIF2α translational regulation. The first two strategies were used in primary neurons and 5XFAD transgenic mice, while the third strategy was employed only in 5XFAD mice. Despite very effective reduction of eIF2α phosphorylation in both primary neurons and 5XFAD brains, or elimination of eIF2α-mediated regulation of BACE1-YFP mRNA translation in 5XFAD brains, Aβ-dependent BACE1 elevation was not decreased. Additionally, robust inhibition of eIF2α phosphorylation did not block Aβ-dependent APP elevation in primary neurons, nor did it reduce amyloid pathology in 5XFAD mice. We conclude that amyloid-associated BACE1 elevation is not caused by translational de-repression via eIF2α phosphorylation, but instead appears to involve a post-translational mechanism. These definitive genetic results exclude a role for eIF2α phosphorylation in Aβ-dependent BACE1 and APP elevation. We suggest a vicious pathogenic cycle wherein Aβ42 toxicity induces peri-plaque BACE1 and APP accumulation in dystrophic neurites leading to exacerbated Aβ production and plaque progression.     

The Light-Harvesting Chlorophyll a/b Binding Proteins Lhcb1 and Lhcb2 Play Complementary Roles during State Transitions in Arabidopsis

Photosynthetic light harvesting in plants is regulated by phosphorylation-driven state transitions: functional redistributions of the major trimeric light-harvesting complex II (LHCII) to balance the relative excitation of photosystem I and photosystem II. State transitions are driven by reversible LHCII phosphorylation by the STN7 kinase and PPH1/TAP38 phosphatase. LHCII trimers are composed of Lhcb1, Lhcb2, and Lhcb3 proteins in various trimeric configurations. Here, we show that despite their nearly identical amino acid composition, the functional roles of Lhcb1 and Lhcb2 are different but complementary. Arabidopsis thaliana plants lacking only Lhcb2 contain thylakoid protein complexes similar to wild-type plants, where Lhcb2 has been replaced by Lhcb1. However, these do not perform state transitions, so phosphorylation of Lhcb2 seems to be a critical step. In contrast, plants lacking Lhcb1 had a more profound antenna remodeling due to a decrease in the amount of LHCII trimers influencing thylakoid membrane structure and, more indirectly, state transitions. Although state transitions are also found in green algae, the detailed architecture of the extant seed plant light-harvesting antenna can now be dated back to a time after the divergence of the bryophyte and spermatophyte lineages, but before the split of the angiosperm and gymnosperm lineages more than 300 million years ago.     

Balancing selection and heterozygote advantage in major histocompatibility complex loci of the bottlenecked Finnish wolf population

Maintaining effective immune response is an essential factor in the survival of small populations. One of the most important immune gene regions is the highly polymorphic major histocompatibility complex (MHC). We investigated how a population bottleneck and recovery have influenced the diversity and selection in three MHC class II loci, DLA‐DRB1, DLA‐DQA1 and DLA‐DQB1, in the Finnish wolf population. We studied the larger Russian Karelian wolf population for comparison and used 17 microsatellite markers as reference loci. The Finnish and Karelian wolf populations did not differ substantially in their MHC diversities ( = 0.047, P = 0.377), but differed in neutral microsatellite diversities ( = 0.148, P = 0.008). MHC allele frequency distributions in the Finnish population were more even than expected under neutrality, implying balancing selection. In addition, an excess of nonsynonymous compared to synonymous polymorphisms indicated historical balancing selection. We also studied association between helminth (Trichinella spp. and Echinococcus canadensis) prevalence and MHC diversity at allele and SNP level. MHC‐heterozygous wolves were less often infected by Trichinella spp. and carriers of specific MHC alleles, SNP haplotypes and SNP alleles had less helminth infections. The associated SNP haplotypes and alleles were shared by different MHC alleles, which emphasizes the necessity of single‐nucleotide‐level association studies also in MHC. Here, we show that strong balancing selection has had similar effect on MHC diversities in the Finnish and Russian Karelian wolf populations despite significant genetic differentiation at neutral markers and small population size in the Finnish population.     

Insights into Conifer Giga-Genomes

Insights from sequenced genomes of major land plant lineages have advanced research in almost every aspect of plant biology. Until recently, however, assembled genome sequences of gymnosperms have been missing from this picture. Conifers of the pine family (Pinaceae) are a group of gymnosperms that dominate large parts of the world’s forests. Despite their ecological and economic importance, conifers seemed long out of reach for complete genome sequencing, due in part to their enormous genome size (20–30 Gb) and the highly repetitive nature of their genomes. Technological advances in genome sequencing and assembly enabled the recent publication of three conifer genomes: white spruce (Picea glauca), Norway spruce (Picea abies), and loblolly pine (Pinus taeda). These genome sequences revealed distinctive features compared with other plant genomes and may represent a window into the past of seed plant genomes. This Update highlights recent advances, remaining challenges, and opportunities in light of the publication of the first conifer and gymnosperm genomes.Conifers are the most widely distributed group of gymnosperms, with 600 to 630 species in 69 genera, including 220 to 250 species of the Pinaceae family (Wang and Ran, 2014). Coniferous forests cover an estimated 39% of the world’s forests (Armenise et al., 2012). Conifers dominate many natural and planted forests in the northern hemisphere and are also planted as exotics for commercial forestry in the southern hemisphere. The importance of conifers for global ecosystem services, their value for forestry-dependent economies, and their contrasting biology with angiosperms are major drivers behind efforts to understand the complex structure, functions, and evolution of their genomes. However, owing to their nonmodel system attributes (i.e. slow-growing and long-lived life history traits), extremely large genome size (Fig. 1), and repeat-rich genome sequence with repeats mostly in the form of transposable elements, no reports of a conifer genome assembly, or any gymnosperm genome for that matter (Soltis and Soltis, 2013), were published until recently. Following early releases of the white spruce (Picea glauca) and loblolly pine (Pinus taeda) genome sequences in public databases (e.g. National Center for Biotechnology Information and http://dendrome.ucdavis.edu/treegenes/), a series of articles described the first conifer genome assemblies for Norway spruce (Picea abies; Nystedt et al., 2013) and interior white spruce, a genetic admix of white spruce (Birol et al., 2013) and loblolly pine (Neale et al., 2014; Zimin et al., 2014). Norway spruce is a prominent forest tree in northern Europe. White spruce is a dominant tree species across the large Canadian forest landscape. Loblolly pine dominates commercial forestry in the southeastern United States. White spruce, Norway spruce, and loblolly pine represent some of the most economically important conifers worldwide, and they are the subjects of important tree improvement/breeding programs (Mullin et al., 2011). This Update highlights significant insights obtained from these genomes as well as some ongoing challenges and recent developments in conifer genomics.Open in a separate windowFigure 1.Size and assembly of conifer genomes compared with other plant genomes. Genome size is plotted against the number of scaffolds divided by the haploid chromosome number for a range of plant species. As such, an assembly that reconstructs a genome with perfect contiguity will have a value of 1, and values greater than 1 represent increasing genome fragmentation. Genome assemblies that utilized Sanger sequencing either in full or in part are represented as white circles. Assemblies constructed using only next generation sequencing technologies are represented as black circles. Both axes are plotted on a log₁₀ scale. With the exception of Populus tremula, Hordeum vulgare, and the three conifer genomes, all genomes were obtained from the Phytozome resource (version 10; http://phytozome.jgi.doe.gov/). The early release draft assembly of P. tremula was obtained from the PopGenIE.org FTP resource (ftp://popgenie.org/popgenie/UPSC_genomes/UPSC_Draft_Assemblies/Current/Genome/) and H. vulgare ‘Morex’ from the Munich Information Center for Protein Sequences barley genome database FTP resource (ftp://ftpmips.helmholtz-muenchen.de/plants/barley/public_data/sequences/). The conifer genomes are detailed by Birol et al. (2013), Nystedt et al. (2013), and Zimin et al. (2014).     

Occurrence of mixotrophic flagellates in relation to bacterioplankton production, light regime and availability of inorganic nutrients in unproductive lakes with differing humic contents

1. Field data from five unproductive Swedish lakes were used to investigate the occurrence of mixotrophic flagellates in relation to bacterioplankton, autotrophic phytoplankton, heterotrophic flagellates and abiotic environmental factors. Three different sources of data were used: (i) a 3‐year study (1995–97) of the humic Lake Örträsket, (ii) seasonal measurements from five lakes with widely varying dissolved organic carbon (DOC) concentrations, and (iii) whole lake enrichment experiments with inorganic nutrients and organic carbon. 2. Mixotrophic flagellates usually dominated over autotrophic phytoplankton in Lake Örträsket in early summer, when both bacterial production and light levels were high. Comparative data from the five lakes demonstrated that the ratio between the biomasses of mixotrophic flagellates and autotrophic phytoplankton (the M/A‐ratio) was positively correlated to bacterioplankton production, but not to the light regime. Whole lake carbon addition (white sugar) increased bacterial biomass, and production, reduced the biomass of autotrophs by a factor of 16, and increased the M/A‐ratio from 0.03 to 3.4. Collectively, the results indicate that the dominance of mixotrophs among phytoplankton was positively related to bacterioplankton production. 3. Whole lake fertilisation with nitrogen (N) and phosphorus (P) demonstrated that the obligate autotrophic phytoplankton was limited by N. N‐addition increased the biomass of the autotrophic phytoplankton but had no effect on mixotrophic flagellates or bacteria, and the M/A‐ratio decreased from 1.2 to 0.6 after N‐enrichment. Therefore, we suggest that bacteria under natural conditions, by utilising allochthonous DOC as an energy and carbon source, are able to outcompete autotrophs for available inorganic nutrients. Consequently, mixotrophic flagellates can become the dominant phytoplankters when phagotrophy permits them to use nutrients stored in bacterial biomass. 4. In Lake Örträsket, the biomass of mixotrophs was usually higher than the biomass of heterotrophs during the summer. This dominance could not be explained by higher grazing rates among the mixotrophs. Instead, ratios between mixotrophic and heterotrophic biomass (the M/H‐ratio) were positively related to light availability. Therefore, we suggest that photosynthesis can enable mixotrophic flagellates to outcompete heterotrophic flagellates.     

The effect of chain length on the melting temperature and size of dialkyldimethylammonium bromide vesicles

Differential scanning calorimetry (DSC) and dynamic light scattering (DLS) were used to obtain the gel to liquid-crystalline phase transition temperature (Tm) and the apparent hydrodynamic radius (Rh) of spontaneously formed cationic vesicles of dialkyldimethylammonium bromide salts (CnH2n+1)2(CH3)2N+.Br-, with varying chain lengths. The preparation of cationic vesicles from aqueous solution of these surfactants, for n=12, 14, 16 and 18 (DDAB, DTDAB, DHDAB and DODAB, respectively), requires the knowledge of the surfactant gel to liquid-crystalline phase transition temperature, or melting temperature (Tm) since below this temperature these surfactants are poorly or not soluble in water. That series of cationic surfactants has been widely investigated as vesicle-forming surfactants, although C12 and C18, DDAB and DODAB are by far the most investigated from this series. The dependence of Tm of these surfactants on the number n of carbons in the surfactant tails is reported. The Tm obtained by DSC increases non-linearly with n, and the vesicle apparent radius Rh is about the same for DHDAB and DODAB, but much smaller for DDAB.     

Fatty acid binding protein 4 in human skeletal muscle

The mechanisms that regulate intramyocellular triglycerol (IMTG) storage and mobilization are largely unknown. However, during the last decades several intracellular fatty acid binding proteins (FABPs) have been identified. FABP3 is the dominating FABP in skeletal muscle. Expression of additional FABPs is suggested from findings in FABP3-null mutated mice. In the present study, our aims were to investigate if FABP4 is expressed within skeletal muscle fibers and if FABP3 and FABP4 are more abundant in skeletal muscle fibers in endurance-trained than in control subjects. We show that FABP4 protein is expressed within the skeletal muscle fibers and that FABP4 mRNA and protein are more abundant in the endurance trained subjects. Still, FABP4 is markedly less expressed than FABP3, which is the generally accepted dominating FABP in skeletal muscle tissue.