Antagonistic interactions subdue inter‐species green‐beard cooperation in bacteria

Cooperation can be favoured through the green‐beard mechanism, where a set of linked genes encodes both a cooperative trait and a phenotypic marker (green beard), which allows carriers of the trait to selectively direct cooperative acts to other carriers. In theory, the green‐beard mechanism should favour cooperation even when interacting partners are totally unrelated at the genome level. Here, we explore such an extreme green‐beard scenario between two unrelated bacterial species—Pseudomonas aeruginosa and Burkholderia cenocepacia, which share a cooperative locus encoding the public good pyochelin (an iron‐scavenging siderophore) and its cognate receptor (green beard) required for iron–pyochelin uptake. We show that pyochelin, when provided in cell‐free supernatants, can be mutually exchanged between species and provide fitness benefits under iron limitation. However, in co‐culture we observed that these cooperative benefits vanished and communities were dominated by P. aeruginosa, regardless of strain background and species starting frequencies. Our results further suggest that P. aeruginosa engages in interference competition to suppress B. cenocepacia, indicating that inter‐species conflict arising from dissimilarities at the genome level overrule the aligned cooperative interests at the pyochelin locus. Thus, green‐beard cooperation is subdued by competition, indicating that interspecific siderophore cooperation is difficult to evolve and to be maintained.     

Additional Psychometric Evidence and Construct Validity for a Revised Preferences Scale of Morningness

The Preferences Scale of morningness was originally conceived as a two‐factor, 11‐item model. More recently, a two‐factor, six‐item model has been developed and supported in two independent samples. These competing models were examined using structural equation modeling in a mixed student and working sample (n=120). The results supported the two factor, six‐item model (χ²(8, N=120)=10.84, p>0.05) as best fitting the sample data. The two factors explained 59% of the variance and Cronbach's alpha was 0.71. Significant differences (p<0.01) in self‐reported driving ability between morning and evening types were obtained by time‐of‐day, providing some evidence of construct validity. Morning types rated their performance as better in the morning hours and evening types rated their performance as better at night. We then examined whether self‐rated performance is subject to some degree of bias. For both morning and evening types, there was a tendency for those scoring high on self‐deception to rate their driving ability as better, but these differences were not significant. Overall, these findings suggest that self‐ratings are suitable for use in determining construct validity. Recommendations for future studies are made.     

Fates of Evolutionarily Distinct,Plastid‐type Glyceraldehyde 3‐phosphate Dehydrogenase Genes in Kareniacean Dinoflagellates

The ancestral kareniacean dinoflagellate has undergone tertiary endosymbiosis, in which the original plastid is replaced by a haptophyte endosymbiont. During this plastid replacement, the endosymbiont genes were most likely flowed into the host dinoflagellate genome (endosymbiotic gene transfer or EGT). Such EGT may have generated the redundancy of functionally homologous genes in the host genome—one has resided in the host genome prior to the haptophyte endosymbiosis, while the other transferred from the endosymbiont genome. However, it remains to be well understood how evolutionarily distinct but functionally homologous genes were dealt in the dinoflagellate genomes bearing haptophyte‐derived plastids. To model the gene evolution after EGT in plastid replacement, we here compared the characteristics of the two evolutionally distinct genes encoding plastid‐type glyceraldehyde 3‐phosphate dehydrogenase (GAPDH) in Karenia brevis and K. mikimotoi bearing haptophyte‐derived tertiary plastids: “gapC1h” acquired from the haptophyte endosymbiont and “gapC1p” inherited from the ancestral dinoflagellate. Our experiments consistently and clearly demonstrated that, in the two species examined, the principal plastid‐type GAPDH is encoded by gapC1h rather than gapC1p. We here propose an evolutionary scheme resolving the EGT‐derived redundancy of genes involved in plastid function and maintenance in the nuclear genomes of dinoflagellates that have undergone plastid replacements. Although K. brevis and K. mikimotoi are closely related to each other, the statuses of the two evolutionarily distinct gapC1 genes in the two Karenia species correspond to different steps in the proposed scheme.     

De novo genome assembly of the stress tolerant forest species Casuarina equisetifolia provides insight into secondary growth

Casuarina equisetifolia (C. equisetifolia), a conifer‐like angiosperm with resistance to typhoon and stress tolerance, is mainly cultivated in the coastal areas of Australasia. C. equisetifolia, making it a valuable model to study secondary growth associated genes and stress‐tolerance traits. However, the genome sequence is unavailable and therefore wood‐associated growth rate and stress resistance at the molecular level is largely unexplored. We therefore constructed a high‐quality draft genome sequence of C. equisetifolia by a combination of Illumina second‐generation sequencing reads and Pacific Biosciences single‐molecule real‐time (SMRT) long reads to advance the investigation of this species. Here, we report the genome assembly, which contains approximately 300 megabases (Mb) and scaffold size of N50 is 1.06 Mb. Additionally, gene annotation, assisted by a combination of prediction and RNA‐seq data, generated 29 827 annotated protein‐coding genes and 1983 non‐coding genes, respectively. Furthermore, we found that the total number of repetitive sequences account for one‐third of the genome assembly. Here we also construct the genome‐wide map of DNA modification, such as two novel forms N⁶‐adenine (6mA) and N4‐methylcytosine (4mC) at the level of single‐nucleotide resolution using single‐molecule real‐time (SMRT) sequencing. Interestingly, we found that 17% of 6mA modification genes and 15% of 4mC modification genes also included alternative splicing events. Finally, we investigated cellulose, hemicellulose, and lignin‐related genes, which were associated with secondary growth and contained different DNA modifications. The high‐quality genome sequence and annotation of C. equisetifolia in this study provide a valuable resource to strengthen our understanding of the diverse traits of trees.     

Genome plasticity in response to stress in Tetrahymena thermophila: selective and reversible chromosome amplification and paralogous expansion of metallothionein genes

Extreme stress situations can induce genetic variations including genome reorganization. In ciliates like Tetrahymena thermophila, the approximately 45‐fold ploidy of the somatic macronucleus may enable adaptive responses that depend on genome plasticity. To identify potential genome‐level adaptations related to metal toxicity, we isolated three Tetrahymena thermophila strains after an extended adaptation period to extreme metal concentrations (Cd²⁺, Cu²⁺ or Pb²⁺). In the Cd‐adapted strain, we found a approximately five‐fold copy number increase of three genes located in the same macronuclear chromosome, including two metallothionein genes, MTT1 and MTT3. The apparent amplification of this macronuclear chromosome was reversible and reproducible, depending on the presence of environmental metal. We also analysed three knockout (KO) and/or knockdown (KD) strains for MTT1 and/or MTT5. In the MTT5KD strain, we found at least two new genes arising from paralogous expansion of MTT1, which encode truncated variants of MTT1. The expansion can be explained by a model based on somatic recombination between MTT1 genes on pairs of macronuclear chromosomes. At least two of the new paralogs are transcribed and upregulated in response to Cd²⁺. Altogether, we have thus identified two distinct mechanisms, both involving genomic plasticity in the polyploid macronucleus that may represent adaptive responses to metal‐related stress.     

Facilitation among plants can accelerate density‐dependent mortality and steepen self‐thinning lines in stressful environments

The speed and slope of plant self‐thinning are all affected by plant–plant interactions across environmental gradients. Possible mechanisms driving the self‐thinning dynamics include the relative strength of root versus shoot competition, and the interplay between competition and facilitation. Although these mechanisms often act in concert, their relative importance has not yet been fully explored. We used both a one‐layer and a two‐layer zone‐of‐influence (ZOI) model to examine how competition and facilitation drive self‐thinning across stress gradients. As a development of the traditional ZOI model, the two‐layer version explicitly models shoot and root growth and neighbor interactions, and thus the overall size‐symmetry of competition is regulated by the relative strength of root versus shoot competition. One‐layer model simulations revealed that increasingly asymmetric competition accelerated thinning, and steepened (slope ranged from about –1 to –4/3) and lowered self‐thinning lines. Stress slowed down density‐dependent mortality considerably when competition was not completely symmetric. Stress significantly decreased the self‐thinning intercept, while facilitation simply counteracted stress effects. Both stress and facilitation showed little effect on the slope. In the two‐layer model, both stress and facilitation affected mortality in the same way as in the one‐layer version when competition was not completely symmetric. Different from the one‐layer model, the two‐layer version showed that the effects of stress and facilitation on the self‐thinning slope were mediated by the asymmetry of competition. As stress increased, the overall asymmetry of competition shifted from asymmetric to symmetric due to increased relative strength of root competition. High stress thus dramatically flattened self‐thinning lines, whereas the inclusion of facilitation counteracted stress and led to steeper self‐thinning lines. Our two‐layer model is based on the current knowledge of plant–plant interactions, and better represents ecological realities. It can help elaborate experiments for testing the role of competition and facilitation in driving plant population dynamics.     

MADS-box out of the black box

The compelling elegance of using genome‐wide scans to detect the signature of selection is difficult to resist, but is countered by the low demonstrated efficacy of pinpointing the actual genes and traits that are the targets of selection in nonmodel species. While the difficulty of going from a suggestive signature to a functional nucleotide polymorphism should not prevent researchers from using genome scans, it does lessen their long‐term utility within and across study systems. In a new study published in this issue of Molecular Ecology ( Mariac et al. 2011 ), researchers have gone a long way towards increasing the relevance of genome‐wide scans for selection via two approaches: (i) they tailored the markers used in the scan to target a family of developmental genes that were good candidates for controlling a trait of interest and (ii) they used an independent mapping population to confirm the association of the gene with polymorphism in the trait of interest. All of this was completed in the nonmodel system of pearl millet (Pennisetum glaucum) and may provide a road map for other researchers hoping to pin down solid candidate genes for selected traits in natural or cultivated systems. Outside of these broad methodological innovations, the paper specifically focuses on a trait (flowering time) that varies across an environmental gradient (rainfall). This environmental gradient potentially serves as a model for environmental change over time, and allele frequencies at the gene can therefore be used to track how populations of pearl millet will adapt to future climate shifts at the genetic level.     

Flip‐flop organization in the chloroplast genome of Capsosiphon fulvescens (Ulvophyceae,Chlorophyta)

To better understand organelle genome evolution of the ulvophycean green alga Capsosiphon fulvescens, we sequenced and characterized its complete chloroplast genome. The circular chloroplast genome was 111,561 bp in length with 31.3% GC content that contained 108 genes including 77 protein‐coding genes, two copies of rRNA operons, and 27 tRNAs. In this analysis, we found the two types of isoform, called heteroplasmy, were likely caused by a flip‐flop organization. The flip‐flop mechanism may have caused structural variation and gene conversion in the chloroplast genome of C. fulvescens. In a phylogenetic analysis based on all available ulvophycean chloroplast genome data, including a new C. fulvescens genome, we found three major conflicting signals for C. fulvescens and its sister taxon Pseudoneochloris marina within 70 individual genes: (i) monophyly with Ulotrichales, (ii) monophyly with Ulvales, and (iii) monophyly with the clade of Ulotrichales and Ulvales. Although the 70‐gene concatenated phylogeny supported monophyly with Ulvales for both species, these complex phylogenetic signals of individual genes need further investigations using a data‐rich approach (i.e., organelle genome data) from broader taxon sampling.     

Recent selection for self‐compatibility in a population of Leavenworthia alabamica

The evolution of self‐compatibility (SC) is the first step in the evolutionary transition in plants from outcrossing enforced by self‐incompatibility (SI) to self‐fertilization. In the Brassicaceae, SI is controlled by alleles of two tightly linked genes at the S‐locus. Despite permitting inbreeding, mutations at the S‐locus leading to SC may be selected if they provide reproductive assurance and/or gain a transmission advantage in a population when SC plants self‐ and outcross. Positive selection can leave a genomic signature in the regions physically linked to the focus of selection when selection has occurred recently. From an SC population of Leavenworthia alabamica with a known nonfunctional mutation at the S‐locus, we collected sequence data from a ～690 Kb region surrounding the S‐locus, as well as from regions not linked to the S‐locus. To test for recent positive selection acting at the S‐locus, we examined polymorphism and the site‐frequency spectra. Using forward simulations, we demonstrate that recent selection of the strength expected for SC at a locus formerly under balancing selection can generate patterns similar to those seen in our empirical data.     

Nuclear–mitochondrial epistasis: a gene's eye view of genomic conflict

We use population genetic models to investigate the cooperative and conflicting synergistic fitness effects between genes from the nucleus and the mitochondrion. By varying fitness parameters, we examine the scope for conflict relative to cooperation among genomes and the utility of the “gene's eye view” analytical approach, which is based on the marginal average fitness of specific alleles. Because sexual conflict can maintain polymorphism of mitochondrial haplotypes, we can explore two types of evolutionary conflict (genomic and sexual) with one epistatic model. We find that the nuclear genetic architecture (autosomal, X‐linked, or Z‐linked) and the mating system change the regions of parameter space corresponding to the evolution by sexual and genomic conflict. For all models, regardless of conflict or cooperation, we find that population mean fitness increases monotonically as evolution proceeds. Moreover, we find that the process of gene frequency change with positive, synergistic fitnesses is self‐accelerating, as the success of an allele in one genome or in one sex increases the frequency of the interacting allele upon which its success depends. This results in runaway evolutionary dynamics caused by the positive intergenomic associations generated by selection. An inbreeding mating system tends to further accelerate these runaway dynamics because it maintains favorable host–symbiont or male–female gene combinations. In contrast, where conflict predominates, the success of an allele in one genome or in one sex diminishes the frequency of the corresponding allele in the other, resulting in considerably slower evolutionary dynamics. The rate of change of mean fitness is also much faster with positive, synergistic fitnesses and much slower where conflict is predominant. Consequently, selection rapidly fixes cooperative gene combinations, while leaving behind a slowing evolving residue of conflicting gene combinations at mutation–selection balance. We discuss how an emphasis on marginal fitness averages may obscure the interdependence of allelic fitness across genomes, making the evolutionary trajectories appear independent of one another when they are not.     

GENOME‐WIDE INVESTIGATION OF REPRODUCTIVE ISOLATION IN EXPERIMENTAL LINEAGES AND NATURAL SPECIES OF NEUROSPORA: IDENTIFYING CANDIDATE REGIONS BY MICROARRAY‐BASED GENOTYPING AND MAPPING

Inherent incompatibilities between genetic components from genomes of different species may cause intrinsic reproductive isolation. In evolution experiments designed to instigate speciation in laboratory populations of the filamentous fungus Neurospora, we previously discovered a pair of incompatibility loci (dfe and dma) that interact negatively to cause severe defects in sexual reproduction. Here we show that the dfe‐dma incompatibility also is a significant cause of genetic isolation between two naturally occurring species of Neurospora (N. crassa and N. intermedia). The strong incompatibility interaction has a simple genetic basis (two biallelic loci) and antagonistic epistasis occurs between heterospecific alleles only, consistent with the Dobzhansky–Muller model of genic incompatibility. We developed microarray‐based, restriction‐site associated DNA (RAD) markers that identified ～1500 polymorphisms between the genomes of the two species, and constructed the first interspecific physical map of Neurospora. With this new mapping resource, the approximate genomic locations of the incompatibility loci were determined using three different approaches: genome scanning, bulk‐segregant analyses, and introgression. These population, quantitative, and classical genetics methods concordantly identified two candidate regions, narrowing the search for each incompatibility locus to only ～2% of the nuclear genome. This study demonstrates how advances in high‐throughput, genome‐wide genotyping can be applied to mapping reproductive isolation genes and speciation research.     

A negative genetic interaction map in isogenic cancer cell lines reveals cancer cell vulnerabilities

Improved efforts are necessary to define the functional product of cancer mutations currently being revealed through large‐scale sequencing efforts. Using genome‐scale pooled shRNA screening technology, we mapped negative genetic interactions across a set of isogenic cancer cell lines and confirmed hundreds of these interactions in orthogonal co‐culture competition assays to generate a high‐confidence genetic interaction network of differentially essential or differential essentiality (DiE) genes. The network uncovered examples of conserved genetic interactions, densely connected functional modules derived from comparative genomics with model systems data, functions for uncharacterized genes in the human genome and targetable vulnerabilities. Finally, we demonstrate a general applicability of DiE gene signatures in determining genetic dependencies of other non‐isogenic cancer cell lines. For example, the PTEN^?/? DiE genes reveal a signature that can preferentially classify PTEN‐dependent genotypes across a series of non‐isogenic cell lines derived from the breast, pancreas and ovarian cancers. Our reference network suggests that many cancer vulnerabilities remain to be discovered through systematic derivation of a network of differentially essential genes in an isogenic cancer cell model.     

The limited role of differential fractionation in genome content variation and function in maize (Zea mays L.) inbred lines

Maize is a diverse paleotetraploid species with considerable presence/absence variation and copy number variation. One mechanism through which presence/absence variation can arise is differential fractionation. Fractionation refers to the loss of duplicate gene pairs from one of the maize subgenomes during diploidization. Differential fractionation refers to non‐shared gene loss events between individuals following a whole‐genome duplication event. We investigated the prevalence of presence/absence variation resulting from differential fractionation in the syntenic portion of the genome using two whole‐genome de novo assemblies of the inbred lines B73 and PH207. Between these two genomes, syntenic genes were highly conserved with less than 1% of syntenic genes being subject to differential fractionation. The few variably fractionated syntenic genes that were identified are unlikely to contribute to functional phenotypic variation, as there is a significant depletion of these genes in annotated gene sets. In further comparisons of 60 diverse inbred lines, non‐syntenic genes were six times more likely to be variable than syntenic genes, suggesting that comparisons among additional genome assemblies are not likely to result in the discovery of large‐scale presence/absence variation among syntenic genes.     

An integrated analysis to predict micro‐RNAs targeting both stemness and metastasis in breast cancer stem cells

Several evidences support the idea that a small population of tumour cells representing self‐renewal potential are involved in initiation, maintenance, metastasis, and outcomes of cancer therapy. Elucidation of microRNAs/genes regulatory networks activated in cancer stem cells (CSCs) is necessary for the identification of new targets for cancer therapy. The aim of the present study was to predict the miRNAs pattern, which can target both metastasis and self‐renewal pathways using integration of literature and data mining. For this purpose, mammospheres derived from MCF‐7, MDA‐MB231, and MDA‐MB468 were used as breast CSCs model. They had higher migration, invasion, and colony formation potential, with increasing in stemness‐ and EMT‐related genes expression. Our results determined that miR‐204, ‐200c, ‐34a, and ‐10b contemporarily could target both self‐renewal and EMT pathways. This core regulatory of miRNAs could increase the survival rate of breast invasive carcinoma via up‐regulation of OCT4, SOX2, KLF4, c‐MYC, NOTCH1, SNAI1, ZEB1, and CDH2 and down‐regulation of CDH1. The majority of those target genes were involved in the regulation of pluripotency, MAPK, WNT, Hedgehog, p53, and transforming growth factor β pathways. Hence, this study provides novel insights for targeting core regulatory of miRNAs in breast CSCs to target both self‐renewal and metastasis potential and eradication of breast cancer.     

Organoid technologies meet genome engineering

Three‐dimensional (3D) stem cell differentiation cultures recently emerged as a novel model system for investigating human embryonic development and disease progression in vitro, complementing existing animal and two‐dimensional (2D) cell culture models. Organoids, the 3D self‐organizing structures derived from pluripotent or somatic stem cells, can recapitulate many aspects of structural organization and functionality of their in vivo organ counterparts, thus holding great promise for biomedical research and translational applications. Importantly, faithful recapitulation of disease and development processes relies on the ability to modify the genomic contents in organoid cells. The revolutionary genome engineering technologies, CRISPR/Cas9 in particular, enable investigators to generate various reporter cell lines for prompt validation of specific cell lineages as well as to introduce disease‐associated mutations for disease modeling. In this review, we provide historical overviews, and discuss technical considerations, and potential future applications of genome engineering in 3D organoid models.     

Genome streamlining via complete loss of introns has occurred multiple times in lichenized fungal mitochondria

Reductions in genome size and complexity are a hallmark of obligate symbioses. The mitochondrial genome displays clear examples of these reductions, with the ancestral alpha‐proteobacterial genome size and gene number having been reduced by orders of magnitude in most descendent modern mitochondrial genomes. Here, we examine patterns of mitochondrial evolution specifically looking at intron size, number, and position across 58 species from 21 genera of lichenized Ascomycete fungi, representing a broad range of fungal diversity and niches. Our results show that the cox1gene always contained the highest number of introns out of all the mitochondrial protein‐coding genes, that high intron sequence similarity (>90%) can be maintained between different genera, and that lichens have undergone at least two instances of complete, genome‐wide intron loss consistent with evidence for genome streamlining via loss of parasitic, noncoding DNA, in Phlyctis boliviensisand Graphis lineola. Notably, however, lichenized fungi have not only undergone intron loss but in some instances have expanded considerably in size due to intron proliferation (e.g., Alectoria fallacina and Parmotrema neotropicum), even between closely related sister species (e.g., Cladonia). These results shed light on the highly dynamic mitochondrial evolution that is occurring in lichens and suggest that these obligate symbiotic organisms are in some cases undergoing recent, broad‐scale genome streamlining via loss of protein‐coding genes as well as noncoding, parasitic DNA elements.     

A genome‐wide association study suggests new candidate genes for milk production traits in Chinese Holstein cattle

A genome‐wide association study (GWAS) was conducted on 15 milk production traits in Chinese Holstein. The experimental population consisted of 445 cattle, each genotyped by the GGP (GeneSeek genomic profiling)‐BovineLD V3 SNP chip, which had 26 151 public SNPs in its manifest file. After data cleaning, 20 326 SNPs were retained for the GWAS. The phenotypes were estimated breeding values of traits, provided by a public dairy herd improvement program center that had been collected once a month for 3 years. Two statistical models, a fixed‐effect linear regression model and a mixed‐effect linear model, were used to estimate the association effects of SNPs on each of the phenotypes. Genome‐wide significant and suggestive thresholds were set at 2.46E‐06 and 4.95E‐05 respectively. The two statistical models concurrently identified two genome‐wide significant (P < 0.05) SNPs on milk production traits in this Chinese Holstein population. The positional candidate genes, which were the ones closest to these two identified SNPs, were EEF2K (eukaryotic elongation factor 2 kinase) and KLHL1 (kelch like family member 1). These two genes could serve as new candidate genes for milk yield and lactation persistence, yet their roles need to be verified in further function studies.