A complex genomic architecture underlies reproductive isolation in a North American oriole hybrid zone.

Natural hybrid zones provide powerful opportunities for identifying the mechanisms that facilitate and inhibit speciation. Documenting the extent of genomic admixture allows us to discern the architecture of reproductive isolation through the identification of isolating barriers. This approach is particularly powerful for characterizing the accumulation of isolating barriers in systems exhibiting varying levels of genomic divergence. Here, we use a hybrid zone between two species-the Baltimore (Icterus galbula) and Bullock's (I. bullockii) orioles-to investigate this architecture of reproductive isolation. We combine whole genome re-sequencing with data from an additional 313 individuals amplityped at ancestry-informative markers to characterize fine-scale patterns of admixture, and to quantify links between genes and the plumage traits. On a genome-wide scale, we document several putative barriers to reproduction, including elevated peaks of divergence above a generally high genomic baseline, a large putative inversion on the Z chromosome, and complex interactions between melanogenesis-pathway candidate genes. Concordant and coincident clines for these different genomic regions further suggest the coupling of pre- and post-mating barriers. Our findings of complex and coupled interactions between pre- and post-mating barriers suggest a relatively rapid accumulation of barriers between these species, and they demonstrate the complexities of the speciation process.

Interspecific forced copulations generate most hybrids in broadly sympatric ducks.

Although rare, hybrids are more common in broadly sympatric waterfowl than in any other avian family; yet, the behavioral ecology explaining their generation has remained controversial. Leading hypotheses are forced interspecific copulations, mis-imprinting caused by mixed broods, and scarcity of conspecific mates. Using a large sample of hybrid ducks solicited from North American hunters we evaluated these hypotheses by genetically determining the mother and father species of F1 hybrids. Based on abundances in areas where their breeding ranges overlap, the frequency of hybrids varied greatly from expectations, with hybrids between species within recently derived clades being much more frequent than those between more divergent clades. Forced copulations, as measured by large phallus-length asymmetries between parentals, strongly predicted the father species of most F1 hybrids. Thus, most Anas acuta x A. platyrhynchos (Northern Pintail x Mallard) F1s were sired by A. acuta, and most A. platyrhynchos x Mareca strepera (Mallard x Gadwall) F1s were sired by A. platyrhynchos. Siring asymmetries were consistent with phallus length asymmetries in five additional parental combinations, but none had samples large enough to be individually statistically significant. The exception to this trend was our sample of nine A. platyrhynchos x Mareca americana (Mallard x Gadwall) F1s, for which a large phallus asymmetry failed to predict the father species. Hybrids were rare in brood parasitic species, suggesting mis-imprinting to be an unlikely cause of most hybrids; however, our samples of hybrids from regular brood parasites were inadequate to strongly address this hypothesis. We could test the scarcity of mates hypothesis for only a single hybrid combination and it contradicted our prediction: most F1 M. Penelope x M. americana (Eurasian x American Wigeon) were sired by M. penelope, strongly contradicting our prediction that female M. penelope wintering in enormous flocks of M. americana (American Wigeon) on the west coast of North America would have difficulty finding conspecific mates. In general, our results support interspecific forced copulations as the predominant behavioral mechanism generating hybrids in North temperate waterfowl.

Genomics of rapid ecological divergence and parallel adaptation in four tidal marsh sparrows.

Theory suggests that different taxa having colonized a similar, challenging environment will show parallel or lineage-specific adaptations to shared selection pressures, but empirical examples of parallel evolution in independent taxa are exceedingly rare. We employed comparative genomics to identify parallel and lineage-specific responses to selection within and among four species of North American sparrows that represent four independent, post-Pleistocene colonization events by an ancestral, upland subspecies and a derived salt marsh specialist. We identified multiple cases of parallel adaptation in these independent comparisons following salt marsh colonization, including selection of 12 candidate genes linked to osmoregulation. In addition to detecting shared genetic targets of selection across multiple comparisons, we found many novel, species-specific signatures of selection, including evidence of selection of loci associated with both physiological and behavioral mechanisms of osmoregulation. Demographic reconstructions of all four species highlighted their recent divergence and small effective population sizes, as expected given their rapid radiation into saline environments. Our results highlight the interplay of both shared and lineage-specific selection pressures in the colonization of a biotically and abiotically challenging habitat and confirm theoretical expectations that steep environmental clines can drive repeated and rapid evolutionary diversification in birds.

Double-digest RAD sequencing outperforms microsatellite loci at assigning paternity and estimating relatedness: A proof of concept in a highly promiscuous bird.

Information on genetic relationships among individuals is essential to many studies of the behaviour and ecology of wild organisms. Parentage and relatedness assays based on large numbers of single nucleotide polymorphism (SNP) loci hold substantial advantages over the microsatellite markers traditionally used for these purposes. We present a double-digest restriction site-associated DNA sequencing (ddRAD-seq) analysis pipeline that, as such, simultaneously achieves the SNP discovery and genotyping steps and which is optimized to return a statistically powerful set of SNP markers (typically 150-600 after stringent filtering) from large numbers of individuals (up to 240 per run). We explore the trade-offs inherent in this approach through a set of experiments in a species with a complex social system, the variegated fairy-wren (Malurus lamberti) and further validate it in a phylogenetically broad set of other bird species. Through direct comparisons with a parallel data set from a robust panel of highly variable microsatellite markers, we show that this ddRAD-seq approach results in substantially improved power to discriminate among potential relatives and considerably more precise estimates of relatedness coefficients. The pipeline is designed to be universally applicable to all bird species (and with minor modifications to many other taxa), to be cost- and time-efficient, and to be replicable across independent runs such that genotype data from different study periods can be combined and analysed as field samples are accumulated.

Modulation of extracytoplasmic function (ECF) sigma factor promoter selectivity by spacer region sequence.

The ability of bacteria to adapt to stress depends on the conditional expression of specific sets of genes. Bacillus subtilis encodes seven extracytoplasmic function (ECF) sigma (σ) factors that regulate functions important for survival under conditions eliciting cell envelope stress. Of these, four have been studied in detail: σM, σW, σX and σV. These four σ factors recognize overlapping sets of promoters, although the sequences that determine this overlapping recognition are incompletely understood. A major role in promoter selectivity has been ascribed to the core -10 and -35 promoter elements. Here, we demonstrate that a homopolymeric T-tract motif, proximal to the -35 element, functions in combination with the core promoter sequences to determine selectivity for ECF sigma factors. This motif is most critical for promoter activation by σV, and contributes variably to activation by σM, σX and σW. We propose that this motif, which is a feature of the deduced promoter consensus for a subset of ECF σ factors from many species, imparts intrinsic DNA curvature to influence promoter activity. The differential effect of this region among ECF σ factors thereby provides a mechanism to modulate the nature and extent of regulon overlap.

The ECF sigma factor, PSPTO_1043, in Pseudomonas syringae pv. tomato DC3000 is induced by oxidative stress and regulates genes involved in oxidative stress response.

The bacterial plant pathogen Pseudomonas syringae adapts to changes in the environment by modifying its gene expression profile. In many cases, the response is mediated by the activation of extracytoplasmic function (ECF) sigma factors that direct RNA polymerase to transcribe specific sets of genes. In this study we focus on PSPTO_1043, one of ten ECF sigma factors in P. syringae pv. tomato DC3000 (DC3000). PSPTO_1043, together with PSPTO_1042, encode an RpoERsp/ChrR-like sigma/anti-sigma factor pair. Although this gene pair is unique to the P. syringae group among the pseudomonads, homologous genes can be found in photosynthetic genera such as Rhodospirillum, Thalassospira, Phaeospirillum and Parvibaculum. Using ChIP-Seq, we detected 137 putative PSPTO_1043 binding sites and identified a likely promoter motif. We characterized 13 promoter candidates, six of which regulate genes that appear to be found only in P. syringae. PSPTO_1043 responds to the presence of singlet oxygen (1O2) and tert-butyl hydroperoxide (tBOOH) and several of the genes regulated by PSPTO_1043 appear to be involved in response to oxidative stress.

Virulence of Pseudomonas syringae pv. tomato DC3000 Is Influenced by the Catabolite Repression Control Protein Crc.

Pseudomonas syringae infects diverse plant species and is widely used as a model system in the study of effector function and the molecular basis of plant diseases. Although the relationship between bacterial metabolism, nutrient acquisition, and virulence has attracted increasing attention in bacterial pathology, it is largely unexplored in P. syringae. The Crc (catabolite repression control) protein is a putative RNA-binding protein that regulates carbon metabolism as well as a number of other factors in the pseudomonads. Here, we show that deletion of crc increased bacterial swarming motility and biofilm formation. The crc mutant showed reduced growth and symptoms in Arabidopsis and tomato when compared with the wild-type strain. We have evidence that the crc mutant shows delayed hypersensitive response (HR) when infiltrated into Nicotiana benthamiana and tobacco. Interestingly, the crc mutant was more susceptible to hydrogen peroxide, suggesting that, in planta, the mutant may be sensitive to reactive oxygen species generated during pathogen-associated molecular pattern-triggered immunity (PTI). Indeed, HR was further delayed when PTI-induced tissues were challenged with the crc mutant. The crc mutant did not elicit an altered PTI response in plants compared with the wild-type strain. We conclude that Crc plays an important role in growth and survival during infection.

Disruption of the carA gene in Pseudomonas syringae results in reduced fitness and alters motility.

BACKGROUND: Pseudomonas syringae infects diverse plant species and is widely used in the study of effector function and the molecular basis of disease. Although the relationship between bacterial metabolism, nutrient acquisition and virulence has attracted increasing attention in bacterial pathology, there is limited knowledge regarding these studies in Pseudomonas syringae. The aim of this study was to investigate the function of the carA gene and the small RNA P32, and characterize the regulation of these transcripts. RESULTS: Disruption of the carA gene (ΔcarA) which encodes the predicted small chain of carbamoylphosphate synthetase, resulted in arginine and pyrimidine auxotrophy in Pseudomonas syringae pv. tomato DC3000. Complementation with the wild type carA gene was able to restore growth to wild-type levels in minimal medium. Deletion of the small RNA P32, which resides immediately upstream of carA, did not result in arginine or pyrimidine auxotrophy. The expression of carA was influenced by the concentrations of both arginine and uracil in the medium. When tested for pathogenicity, ΔcarA showed reduced fitness in tomato as well as Arabidopsis when compared to the wild-type strain. In contrast, mutation of the region encoding P32 had minimal effect in planta. ΔcarA also exhibited reduced motility and increased biofilm formation, whereas disruption of P32 had no impact on motility or biofilm formation. CONCLUSIONS: Our data show that carA plays an important role in providing arginine and uracil for growth of the bacteria and also influences other factors that are potentially important for growth and survival during infection. Although we find that the small RNA P32 and carA are co-transcribed, P32 does not play a role in the phenotypes that carA is required for, such as motility, cell attachment, and virulence. Additionally, our data suggests that pyrimidines may be limited in the apoplastic space of the plant host tomato.

Plumage Genes and Little Else Distinguish the Genomes of Hybridizing Warblers.

When related taxa hybridize extensively, their genomes may become increasingly homogenized over time. This mixing via hybridization creates conservation challenges when it reduces genetic or phenotypic diversity and when it endangers previously distinct species via genetic swamping [1]. However, hybridization also facilitates admixture mapping of traits that distinguish each species and the associated genes that maintain distinctiveness despite ongoing gene flow [2]. We address these dual aspects of hybridization in the golden-winged/blue-winged warbler complex, two phenotypically divergent warblers that are indistinguishable using traditional molecular markers and that draw substantial conservation attention [3-5]. Whole-genome comparisons show that differentiation is extremely low: only six small genomic regions exhibit strong differences. Four of these divergence peaks occur in proximity to genes known to be involved in feather development or pigmentation: agouti signaling protein (ASIP), follistatin (FST), ecodysplasin (EDA), wingless-related integration site (Wnt), and beta-carotene oxygenase 2 (BCO2). Throat coloration-the most striking plumage difference between these warblers-is perfectly associated with the promoter region of agouti, and genotypes at this locus obey simple Mendelian recessive inheritance of the black-throated phenotype characteristic of golden-winged warblers. The more general pattern of genomic similarity between these warblers likely results from a protracted period of hybridization, contradicting the broadly accepted hypothesis that admixture results from solely anthropogenic habitat change in the past two centuries [4]. Considered in concert, these results are relevant to both the genetic architecture of avian feather pigmentation and the evolutionary history and conservation challenges associated with these declining songbirds.

Analysis of the small RNA spf in the plant pathogen Pseudomonas syringae pv. tomato strain DC3000.

Bacteria contain small non-coding RNAs (ncRNAs) that are typically responsible for altering transcription, translation or mRNA stability. ncRNAs are important because they often regulate virulence factors and susceptibility to various stresses. Here, the regulation of a recently described ncRNA of Pseudomonas syringae DC3000, spot 42 (now referred to as spf), was investigated. A putative RpoE binding site was identified upstream of spf in strain DC3000. RpoE is shown to regulate the expression of spf. Also, deletion of spf results in increased sensitivity to hydrogen peroxide compared with the wild-type strain, suggesting that spf plays a role in susceptibility to oxidative stress. Furthermore, expression of alg8 is shown to be influenced by spf, suggesting that this ncRNA plays a role in alginate biosynthesis. Structural and comparative genomic analyses show this ncRNA is well conserved among the pseudomonads. The findings provide new information on the regulation and role of this ncRNA in P. syringae.

CrcZ and CrcX regulate carbon source utilization in Pseudomonas syringae pathovar tomato strain DC3000.

Small non-coding RNAs (ncRNAs) are important components of many regulatory pathways in bacteria and play key roles in regulating factors important for virulence. Carbon catabolite repression control is modulated by small RNAs (crcZ or crcZ and crcY) in Pseudomonas aeruginosa and Pseudomonas putida. In this study, we demonstrate that expression of crcZ and crcX (formerly designated psr1 and psr2, respectively) is dependent upon RpoN together with the two-component system CbrAB, and is influenced by the carbon source present in the medium in the model plant pathogen Pseudomonas syringae pv tomato DC3000. The distribution of the members of the Crc ncRNA family was also determined by screening available genomic sequences of the Pseudomonads. Interestingly, variable numbers of the Crc family members exist in Pseudomonas genomes. The ncRNAs are comprised of three main subfamilies, named CrcZ, CrcX and CrcY. Most importantly the CrcX subfamily appears to be unique to all P. syringae strains sequenced to date.

Analysis of the small RNA P16/RgsA in the plant pathogen Pseudomonas syringae pv. tomato strain DC3000.

Bacteria contain small non-coding RNAs (ncRNAs) that are responsible for altering transcription, translation or mRNA stability. ncRNAs are important because they regulate virulence factors and susceptibility to various stresses. Here, the regulation of a recently described ncRNA of Pseudomonas syringae pv. tomato DC3000, P16, was investigated. We determined that RpoS regulates the expression of P16. We found that deletion of P16 results in increased sensitivity to hydrogen peroxide compared to the wild-type strain, suggesting that P16 plays a role in the bacteria's susceptibility to oxidative stress. Additionally the P16 mutant displayed enhanced resistance to heat stress. Our findings provide new information on the regulation and role of this ncRNA in P. syringae.

Regulons of three Pseudomonas syringae pv. tomato DC3000 iron starvation sigma factors.

Pseudomonas syringae pv. tomato DC3000 contains genes for 15 sigma factors. The majority are members of the extracytoplasmic function class of sigma factors, including five that belong to the iron starvation subgroup. In this study, we identified the genes controlled by three iron starvation sigma factors. Their regulons are composed of a small number of genes likely to be involved in iron uptake.

An extracytoplasmic function sigma factor-mediated cell surface signaling system in Pseudomonas syringae pv. tomato DC3000 regulates gene expression in response to heterologous siderophores.

The diversity of regulatory systems encoded by bacteria provides an indication of the variety of stresses and interactions that these organisms encounter in nature. We have been investigating how the plant pathogen Pseudomonas syringae pv. tomato DC3000 responds to iron limitation and have focused on the iron starvation (IS) sigma factors to identify regulon members and to explore the mechanistic details of genetic control for this class of regulators. In the study described in this report, we used chromatin immunoprecipitation paired with high-throughput sequencing (ChIP-Seq) to screen the genome for locations associated with binding of the P. syringae IS sigma factor PSPTO_1203. We used multiple methods to demonstrate differential regulation of two genes identified in the ChIP-Seq screen and characterize the promoter elements that facilitate PSPTO_1203-dependent regulation. The genes regulated by PSPTO_1203 encode a TonB-dependent transducer (PSPTO_1206) and a cytoplasmic membrane protein (PSPTO_2145), which is located in the P. syringae pyoverdine cluster. Additionally, we identified siderophores that induce the activity of PSPTO_1203 and used this information to investigate the functional components of the signal transduction cascade.

Characterization of the Fur regulon in Pseudomonas syringae pv. tomato DC3000.

The plant pathogen Pseudomonas syringae pv. tomato DC3000 (DC3000) is found in a wide variety of environments and must monitor and respond to various environmental signals such as the availability of iron, an essential element for bacterial growth. An important regulator of iron homeostasis is Fur (ferric uptake regulator), and here we present the first study of the Fur regulon in DC3000. Using chromatin immunoprecipitation followed by massively parallel sequencing (ChIP-seq), 312 chromosomal regions were highly enriched by coimmunoprecipitation with a C-terminally tagged Fur protein. Integration of these data with previous microarray and global transcriptome analyses allowed us to expand the putative DC3000 Fur regulon to include genes both repressed and activated in the presence of bioavailable iron. Using nonradioactive DNase I footprinting, we confirmed Fur binding in 41 regions, including upstream of 11 iron-repressed genes and the iron-activated genes encoding two bacterioferritins (PSPTO_0653 and PSPTO_4160), a ParA protein (PSPTO_0855), and a two-component system (TCS) (PSPTO_3382 to PSPTO_3380).

Genome-wide identification of transcriptional start sites in the plant pathogen Pseudomonas syringae pv. tomato str. DC3000.

RNA-Seq has provided valuable insights into global gene expression in a wide variety of organisms. Using a modified RNA-Seq approach and Illumina's high-throughput sequencing technology, we globally identified 5'-ends of transcripts for the plant pathogen Pseudomonas syringae pv. tomato str. DC3000. A substantial fraction of 5'-ends obtained by this method were consistent with results obtained using global RNA-Seq and 5'RACE. As expected, many 5'-ends were positioned a short distance upstream of annotated genes. We also captured 5'-ends within intergenic regions, providing evidence for the expression of un-annotated genes and non-coding RNAs, and detected numerous examples of antisense transcription, suggesting additional levels of complexity in gene regulation in DC3000. Importantly, targeted searches for sequence patterns in the vicinity of 5'-ends revealed over 1200 putative promoters and other regulatory motifs, establishing a broad foundation for future investigations of regulation at the genomic and single gene levels.

Transcriptome analysis of Pseudomonas syringae identifies new genes, noncoding RNAs, and antisense activity.

To fully understand how bacteria respond to their environment, it is essential to assess genome-wide transcriptional activity. New high-throughput sequencing technologies make it possible to query the transcriptome of an organism in an efficient unbiased manner. We applied a strand-specific method to sequence bacterial transcripts using Illumina's high-throughput sequencing technology. The resulting sequences were used to construct genome-wide transcriptional profiles. Novel bioinformatics analyses were developed and used in combination with proteomics data for the qualitative classification of transcriptional activity in defined regions. As expected, most transcriptional activity was consistent with predictions from the genome annotation. Importantly, we identified and confirmed transcriptional activity in areas of the genome inconsistent with the annotation and in unannotated regions. Further analyses revealed potential RpoN-dependent promoter sequences upstream of several noncoding RNAs (ncRNAs), suggesting a role for these ncRNAs in RpoN-dependent phenotypes. We were also able to validate a number of transcriptional start sites, many of which were consistent with predicted promoter motifs. Overall, our approach provides an efficient way to survey global transcriptional activity in bacteria and enables rapid discovery of specific areas in the genome that merit further investigation.

In vitro mutagenesis of Bacillus subtilis by using a modified Tn7 transposon with an outward-facing inducible promoter.

A Tn7 donor plasmid, pTn7SX, was constructed for use with the model gram-positive bacterium Bacillus subtilis. This new mini-Tn7, mTn7SX, contains a spectinomycin resistance cassette and an outward-facing, xylose-inducible promoter, thereby allowing for the regulated expression of genes downstream of the transposon. We demonstrate that mTn7SX inserts are obtained at a high frequency and occur randomly throughout the B. subtilis genome. The utility of this system was demonstrated by the selection of mutants with increased resistance to the antibiotic fosfomycin or duramycin.

The yydFGHIJ operon of Bacillus subtilis encodes a peptide that induces the LiaRS two-component system.

The Bacillus subtilis LiaRS two-component system (TCS) responds to perturbations of the cell envelope induced by lipid II-interacting antibiotics, such as vancomycin, ramoplanin, nisin, and bacitracin. Here, we characterize Tn7-generated mutations that induce the liaRS TCS. In addition to insertions in liaF, a known negative regulator of the LiaRS TCS, we identified two disruptions in the last two genes of the yydFGHIJ operon. This operon is predicted to encode a 49-amino-acid peptide (YydF), a modification enzyme (YydG), a membrane-embedded protease (YydH), and an ATP-binding cassette (ABC) transporter (YydIJ). Genome sequence comparisons suggest that the yydFGHIJ operon may have been acquired by horizontal transfer. Inactivation of the YydIJ transporter resulted in increased expression from the LiaR-dependent P(liaI) promoter only in the presence of the yydFGH genes. Cells harboring the complete yydFGHIJ operon induced LiaR activity in cocultured cells lacking either this transporter or the complete operon. These results suggest that this operon is involved in the synthesis and export of a modified peptide (YydF*) that elicits cell envelope stress sensed by the LiaRS TCS.

LiaRS-dependent gene expression is embedded in transition state regulation in Bacillus subtilis.

Maintaining envelope integrity is crucial for the survival of any bacterial cell, especially those living in a complex and ever-changing habitat such as the soil ecosystem. The LiaRS two-component system is part of the regulatory network orchestrating the cell-envelope stress response in Bacillus subtilis. It responds to perturbations of the cell envelope, especially the presence of antibiotics that interfere with the lipid II cycle, such as bacitracin or vancomycin. LiaRS-dependent regulation is strictly repressed by the membrane protein LiaF in the absence of inducing conditions. Here, it is shown that the LiaR-dependent liaI promoter is induced at the onset of stationary phase without addition of exogenous stresses. Its activity is embedded in the complex regulatory cascade governing adaptation at the onset of stationary phase. The liaI promoter is directly repressed by the transition state regulator AbrB and responds indirectly to the activity of Spo0A, the master regulator of sporulation. The activity of the liaI promoter is therefore tightly regulated by at least five regulators to ensure an appropriate level of liaIH expression.