Promoters maintain their relative activity levels under different growth conditions

Most genes change expression levels across conditions, but it is unclear which of these changes represents specific regulation and what determines their quantitative degree. Here, we accurately measured activities of ～900 S. cerevisiae and ～1800 E. coli promoters using fluorescent reporters. We show that in both organisms 60–90% of promoters change their expression between conditions by a constant global scaling factor that depends only on the conditions and not on the promoter's identity. Quantifying such global effects allows precise characterization of specific regulation—promoters deviating from the global scale line. These are organized into few functionally related groups that also adhere to scale lines and preserve their relative activities across conditions. Thus, only several scaling factors suffice to accurately describe genome‐wide expression profiles across conditions. We present a parameter‐free passive resource allocation model that quantitatively accounts for the global scaling factors. It suggests that many changes in expression across conditions result from global effects and not specific regulation, and provides means for quantitative interpretation of expression profiles.     

Comparative genomics of two ‘Candidatus Accumulibacter' clades performing biological phosphorus removal

Members of the genus Candidatus Accumulibacter are important in many wastewater treatment systems performing enhanced biological phosphorus removal (EBPR). The Accumulibacter lineage can be subdivided phylogenetically into multiple clades, and previous work showed that these clades are ecologically distinct. The complete genome of Candidatus Accumulibacter phosphatis strain UW-1, a member of Clade IIA, was previously sequenced. Here, we report a draft genome sequence of Candidatus Accumulibacter spp. strain UW-2, a member of Clade IA, assembled following shotgun metagenomic sequencing of laboratory-scale bioreactor sludge. We estimate the genome to be 80–90% complete. Although the two clades share 16S rRNA sequence identity of >98.0%, we observed a remarkable lack of synteny between the two genomes. We identified 2317 genes shared between the two genomes, with an average nucleotide identity (ANI) of 78.3%, and accounting for 49% of genes in the UW-1 genome. Unlike UW-1, the UW-2 genome seemed to lack genes for nitrogen fixation and carbon fixation. Despite these differences, metabolic genes essential for denitrification and EBPR, including carbon storage polymer and polyphosphate metabolism, were conserved in both genomes. The ANI from genes associated with EBPR was statistically higher than that from genes not associated with EBPR, indicating a high selective pressure in EBPR systems. Further, we identified genomic islands of foreign origins including a near-complete lysogenic phage in the Clade IA genome. Interestingly, Clade IA appeared to be more phage susceptible based on it containing only a single Clustered Regularly Interspaced Short Palindromic Repeats locus as compared with the two found in Clade IIA. Overall, the comparative analysis provided a genetic basis to understand physiological differences and ecological niches of Accumulibacter populations, and highlights the importance of diversity in maintaining system functional resilience.     

Population specific fitness response of Drosophila subobscura to lead pollution

Abstract Differences in heavy metal tolerance among separate populations of the same species have often been interpreted as local adaptation. Persistence of differences after removing the stressor indicates that mechanisms responsible for the increased tolerance were genetically determined. Drosophila subobscura Collin (Diptera: Drosophilidae) populations were sampled from two localities with different history of heavy metal pollution, and reared for eight generations in the laboratory on a standard medium and on media with different concentrations of lead (Pb). To determine whether flies from different natural populations exposed to the Pb‐contaminated media in the laboratory show population specific variability in fitness components over generations, experimental groups with different concentrations of lead were assayed in three generations (F₂, F₅, and F₈) for fecundity, developmental time, and egg‐to‐adult viability. On the contaminated medium, fecundity was reduced in later generations and viability was increased, irrespective of the environmental origin of populations. For both populations, developmental time showed a tendency of slowing down on media with lead. Faster development was observed in later generations. Preadaptation to contamination, meaning higher fecundity, higher viability, and faster egg to adult development in all studied generations, was found in D. subobscura originating from the locality with a higher level of heavy metal pollution.     

Genome‐wide analysis of Nilaparvata lugens nymphal responses to high‐density and low‐quality rice hosts

The brown planthopper (BPH) Nilaparvata lugens is an economically important pest on rice plants. In this study, the higher population density and yellow‐ripe stage of rice plants were used to construct adverse survival conditions (ASC) against BPH nymphs. Simultaneously, the low population density and tillering stage of rice plants were used to establish a suitable survival condition (SSC) as a control. Solexa/Illumina sequencing was used to identify genes of BPH nymphs responding to ASC. Significantly longer duration development of BPH nymphs and significantly lower brachypterous ratio of BPH adults were observed by ASC compared with SSC. A total of 2 544 differentially expressed genes (DEGs) were obtained and analyzed by BLASTx, Gene Ontology and KEGG Orthology. Gene ontology analysis revealed that the DEGs were mainly involved in categories of cell, cell part, cellular process, binding, catalytic, organelle and metabolic processes. 1 138 DEGs having enzyme commission numbers were assigned to different metabolic pathways. The largest clusters were neurodegenerative diseases (137, 12.0%), followed by carbohydrate metabolism (113, 9.9%), amino acid metabolism (94, 8.3%), nucleotide metabolism (76, 6.7%), energy metabolism (64, 5.6%), translation (60, 5.3%), lipid metabolism (58, 5.1%), and folding, sorting and degradation (52, 4.6%). Expressing profile of 11 DEGs during eight nymphal developmental stages of BPH were analyzed by quantitative real‐time polymerase chain reaction. The 11 genes exhibited differential expression between ASC and SSC during at least one developmental stage. The DEGs identified in this study provide molecular proof of how BPH reconfigures its gene expression profile to adapt to overcrowding and low‐quality hosts.     

Confidence interval of intrinsic optimum temperature estimated using thermodynamic SSI model

The intrinsic optimum temperature for the development of ectotherms is one of the most important factors not only for their physiological processes but also for ecological and evolutional processes. The Sharpe–Schoolfield–Ikemoto (SSI) model succeeded in defining the temperature that can thermodynamically meet the condition that at a particular temperature the probability of an active enzyme reaching its maximum activity is realized. Previously, an algorithm was developed by Ikemoto (Tropical malaria does not mean hot environments. Journal of Medical Entomology, 45, 963–969) to estimate model parameters, but that program was computationally very time consuming. Now, investigators can use the SSI model more easily because a full automatic computer program was designed by Shi et al. (A modified program for estimating the parameters of the SSI model. Environmental Entomology, 40, 462–469). However, the statistical significance of the point estimate of the intrinsic optimum temperature for each ectotherm has not yet been determined. Here, we provided a new method for calculating the confidence interval of the estimated intrinsic optimum temperature by modifying the approximate bootstrap confidence intervals method. For this purpose, it was necessary to develop a new program for a faster estimation of the parameters in the SSI model, which we have also done.     

Protein array of Coxiella burnetii probed with Q fever sera

Coxiella burnetii is the etiological agent of Q fever.To identify its major seroreactive proteins,a subgenomic protein array was developed.A total of 101 assumed virulence-associated recombinant proteins of C.burnetii were probed with sera from mice experimentally infected with C.burnetii and sera from Q fever patients.Sixteen proteins were recognized as major seroreactive antigens by the mouse sera.Seven of these 16 proteins reacted positively with at least 45% of Q fever patient sera.Notably,HspB had the highest fluorescence intensity value and positive frequency of all the proteins on the array when probed with both Q fever patient sera and mouse sera.These results suggest that these seven major seroreactive proteins,particularly HspB,are potential serodiagnostic and subunit vaccine antigens of Q fever.     

SNP discovery using Paired‐End RAD‐tag sequencing on pooled genomic DNA of Sisymbrium austriacum (Brassicaceae)

Single nucleotide polymorphisms SNPs are rapidly replacing anonymous markers in population genomic studies, but their use in non model organisms is hampered by the scarcity of cost‐effective approaches to uncover genome‐wide variation in a comprehensive subset of individuals. The screening of one or only a few individuals induces ascertainment bias. To discover SNPs for a population genomic study of the Pyrenean rocket (Sisymbrium austriacum subsp. chrysanthum), we undertook a pooled RAD‐PE (Restriction site Associated DNA Paired‐End sequencing) approach. RAD tags were generated from the PstI‐digested pooled genomic DNA of 12 individuals sampled across the species distribution range and paired‐end sequenced using Illumina technology to produce ~24.5 Mb of sequences, covering ~7% of the specie's genome. Sequences were assembled into ~76 000 contigs with a mean length of 323 bp (N₅₀ = 357 bp, sequencing depth = 24x). In all, >15 000 SNPs were called, of which 47% were annotated in putative genic regions based on homology with the Arabidopsis thaliana genome. Gene ontology (GO) slim categorization demonstrated that the identified SNPs covered extant genic variation well. The validation of 300 SNPs on a larger set of individuals using a KASPar assay underpinned the utility of pooled RAD‐PE as an inexpensive genome‐wide SNP discovery technique (success rate: 87%). In addition to SNPs, we discovered >600 putative SSR markers.