Next-Generation Phylogeography: A Targeted Approach for Multilocus Sequencing of Non-Model Organisms

The field of phylogeography has long since realized the need and utility of incorporating nuclear DNA (nDNA) sequences into analyses. However, the use of nDNA sequence data, at the population level, has been hindered by technical laboratory difficulty, sequencing costs, and problematic analytical methods dealing with genotypic sequence data, especially in non-model organisms. Here, we present a method utilizing the 454 GS-FLX Titanium pyrosequencing platform with the capacity to simultaneously sequence two species of sea star (Meridiastra calcar and Parvulastra exigua) at five different nDNA loci across 16 different populations of 20 individuals each per species. We compare results from 3 populations with traditional Sanger sequencing based methods, and demonstrate that this next-generation sequencing platform is more time and cost effective and more sensitive to rare variants than Sanger based sequencing. A crucial advantage is that the high coverage of clonally amplified sequences simplifies haplotype determination, even in highly polymorphic species. This targeted next-generation approach can greatly increase the use of nDNA sequence loci in phylogeographic and population genetic studies by mitigating many of the time, cost, and analytical issues associated with highly polymorphic, diploid sequence markers

Optical Dissection of Neural Circuits Responsible for Drosophila Larval Locomotion with Halorhodopsin

Halorhodopsin (NpHR), a light-driven microbial chloride pump, enables silencing of neuronal function with superb temporal and spatial resolution. Here, we generated a transgenic line of Drosophila that drives expression of NpHR under control of the Gal4/UAS system. Then, we used it to dissect the functional properties of neural circuits that regulate larval peristalsis, a continuous wave of muscular contraction from posterior to anterior segments. We first demonstrate the effectiveness of NpHR by showing that global and continuous NpHR-mediated optical inhibition of motor neurons or sensory feedback neurons induce the same behavioral responses in crawling larvae to those elicited when the function of these neurons are inhibited by Shibirets, namely complete paralyses or slowed locomotion, respectively. We then applied transient and/or focused light stimuli to inhibit the activity of motor neurons in a more temporally and spatially restricted manner and studied the effects of the optical inhibition on peristalsis. When a brief light stimulus (1–10 sec) was applied to a crawling larva, the wave of muscular contraction stopped transiently but resumed from the halted position when the light was turned off. Similarly, when a focused light stimulus was applied to inhibit motor neurons in one or a few segments which were about to be activated in a dissected larva undergoing fictive locomotion, the propagation of muscular constriction paused during the light stimulus but resumed from the halted position when the inhibition (>5 sec) was removed. These results suggest that (1) Firing of motor neurons at the forefront of the wave is required for the wave to proceed to more anterior segments, and (2) The information about the phase of the wave, namely which segment is active at a given time, can be memorized in the neural circuits for several seconds

Effector diversification within compartments of the Leptosphaeria maculans genome affected by Repeat-Induced Point mutations

Fungi are of primary ecological, biotechnological and economic importance. Many fundamental biological processes that are shared by animals and fungi are studied in fungi due to their experimental tractability. Many fungi are pathogens or mutualists and are model systems to analyse effector genes and their mechanisms of diversification. In this study, we report the genome sequence of the phytopathogenic ascomycete Leptosphaeria maculans and characterize its repertoire of protein effectors. The L. maculans genome has an unusual bipartite structure with alternating distinct guanine and cytosine-equilibrated and adenine and thymine (AT)-rich blocks of homogenous nucleotide composition. The AT-rich blocks comprise one-third of the genome and contain effector genes and families of transposable elements, both of which are affected by repeat-induced point mutation, a fungal-specific genome defence mechanism. This genomic environment for effectors promotes rapid sequence diversification and underpins the evolutionary potential of the fungus to adapt rapidly to novel host-derived constraints

Evolution of Linked Avirulence Effectors in Leptosphaeria maculans Is Affected by Genomic Environment and Exposure to Resistance Genes in Host Plants

Brassica napus (canola) cultivars and isolates of the blackleg fungus, Leptosphaeria maculans interact in a ‘gene for gene’ manner whereby plant resistance (R) genes are complementary to pathogen avirulence (Avr) genes. Avirulence genes encode proteins that belong to a class of pathogen molecules known as effectors, which includes small secreted proteins that play a role in disease. In Australia in 2003 canola cultivars with the Rlm1 resistance gene suffered a breakdown of disease resistance, resulting in severe yield losses. This was associated with a large increase in the frequency of virulence alleles of the complementary avirulence gene, AvrLm1, in fungal populations. Surprisingly, the frequency of virulence alleles of AvrLm6 (complementary to Rlm6) also increased dramatically, even though the cultivars did not contain Rlm6. In the L. maculans genome, AvrLm1 and AvrLm6 are linked along with five other genes in a region interspersed with transposable elements that have been degenerated by Repeat-Induced Point (RIP) mutations. Analyses of 295 Australian isolates showed deletions, RIP mutations and/or non-RIP derived amino acid substitutions in the predicted proteins encoded by these seven genes. The degree of RIP mutations within single copy sequences in this region was proportional to their proximity to the degenerated transposable elements. The RIP alleles were monophyletic and were present only in isolates collected after resistance conferred by Rlm1 broke down, whereas deletion alleles belonged to several polyphyletic lineages and were present before and after the resistance breakdown. Thus, genomic environment and exposure to resistance genes in B. napus has affected the evolution of these linked avirulence genes in L. maculans

Epistatic Association Mapping in Homozygous Crop Cultivars

The genetic dissection of complex traits plays a crucial role in crop breeding. However, genetic analysis and crop breeding have heretofore been performed separately. In this study, we designed a new approach that integrates epistatic association analysis in crop cultivars with breeding by design. First, we proposed an epistatic association mapping (EAM) approach in homozygous crop cultivars. The phenotypic values of complex traits, along with molecular marker information, were used to perform EAM. In our EAM, all the main-effect quantitative trait loci (QTLs), environmental effects, QTL-by-environment interactions and QTL-by-QTL interactions were included in a full model and estimated by empirical Bayes approach. A series of Monte Carlo simulations was performed to confirm the reliability of the new method. Next, the information from all detected QTLs was used to mine novel alleles for each locus and to design elite cross combination. Finally, the new approach was adopted to dissect the genetic basis of seed length in 215 soybean cultivars obtained, by stratified random sampling, from 6 geographic ecotypes in China. As a result, 19 main-effect QTLs and 3 epistatic QTLs were identified, more than 10 novel alleles were mined and 3 elite parental combinations, such as Daqingdou and Zhengzhou790034, were predicted

Policing the legume-Rhizobium symbiosis: a critical test of partner choice

In legume-Rhizobium symbioses, specialised soil bacteria fix atmospheric nitrogen in return for carbon. However, ineffective strains can arise, making discrimination essential. Discrimination can occur via partner choice, where legumes prevent ineffective strains from entering, or via sanctioning, where plants provide fewer resources. Several studies have inferred that legumes exercise partner choice, but the rhizobia compared were not otherwise isogenic. To test when and how plants discriminate ineffective strains we developed sets of fixing and non-fixing strains that differed only in the expression of nifH - essential for nitrogen fixation - and could be visualised using marker genes. We show that the plant is unable to select against the non-fixing strain at the point of entry, but that non-fixing nodules are sanctioned. We also used the technique to characterise mixed nodules (containing both a fixing and a non-fixing strain), whose frequency could be predicted using a simple diffusion model. We discuss that sanctioning is likely to evolve in preference to partner choice in any symbiosis where partner quality cannot be adequately assessed until goods or services are actively exchanged

The Repeatability of Adaptive Radiation During Long-Term Experimental Evolution of Escherichia coli in a Multiple Nutrient Environment

Adaptive radiations occur when a species diversifies into different ecological specialists due to competition for resources and trade-offs associated with the specialization. The evolutionary outcome of an instance of adaptive radiation cannot generally be predicted because chance (stochastic events) and necessity (deterministic events) contribute to the evolution of diversity. With increasing contributions of chance, the degree of parallelism among different instances of adaptive radiations and the predictability of an outcome will decrease. To assess the relative contributions of chance and necessity during adaptive radiation, we performed a selection experiment by evolving twelve independent microcosms of Escherichia coli for 1000 generations in an environment that contained two distinct resources. Specialization to either of these resources involves strong trade-offs in the ability to use the other resource. After selection, we measured three phenotypic traits: 1) fitness, 2) mean colony size, and 3) colony size diversity. We used fitness relative to the ancestor as a measure of adaptation to the selective environment; changes in colony size as a measure of the evolution of new resource specialists because colony size has been shown to correlate with resource specialization; and colony size diversity as a measure of the evolved ecological diversity. Resource competition led to the rapid evolution of phenotypic diversity within microcosms. Measurements of fitness, colony size, and colony size diversity within and among microcosms showed that the repeatability of adaptive radiation was high, despite the evolution of genetic variation within microcosms. Consistent with the observation of parallel evolution, we show that the relative contributions of chance are far smaller and less important than effects due to adaptation for the traits investigated. The two-resource environment imposed similar selection pressures in independent populations and promoted parallel phenotypic adaptive radiations in all independently evolved microcosms

Comprehensive annotation of the Parastagonospora nodorum reference genome using next-generation genomics, transcriptomics and proteogenomics

Parastagonospora nodorum, the causal agent of Septoria nodorum blotch (SNB), is an economically important pathogen of wheat (Triticum spp.), and a model for the study of necrotrophic pathology and genome evolution. The reference P. nodorum strain SN15 was the first Dothideomycete with a published genome sequence, and has been used as the basis for comparison within and between species. Here we present an updated reference genome assembly with corrections of SNP and indel errors in the underlying genome assembly from deep resequencing data as well as extensive manual annotation of gene models using transcriptomic and proteomic sources of evidence (https://github.com/robsyme/Parastagonospora_nodorum_SN15). The updated assembly and annotation includes 8,366 genes with modified protein sequence and 866 new genes. This study shows the benefits of using a wide variety of experimental methods allied to expert curation to generate a reliable set of gene models

On the phylogeny of Mustelidae subfamilies: analysis of seventeen nuclear non-coding loci and mitochondrial complete genomes

<p>Abstract</p> <p>Background</p> <p>Mustelidae, as the largest and most-diverse family of order Carnivora, comprises eight subfamilies. Phylogenetic relationships among these Mustelidae subfamilies remain argumentative subjects in recent years. One of the main reasons is that the mustelids represent a typical example of rapid evolutionary radiation and recent speciation event. Prior investigation has been concentrated on the application of different mitochondrial (mt) sequence and nuclear protein-coding data, herein we employ 17 nuclear non-coding loci (>15 kb), in conjunction with mt complete genome data (>16 kb), to clarify these enigmatic problems.</p> <p>Results</p> <p>The combined nuclear intron and mt genome analyses both robustly support that Taxidiinae diverged first, followed by Melinae. Lutrinae and Mustelinae are grouped together in all analyses with strong supports. The position of Helictidinae, however, is enigmatic because the mt genome analysis places it to the clade uniting Lutrinae and Mustelinae, whereas the nuclear intron analysis favores a novel view supporting a closer relationship of Helictidinae to Martinae. This finding emphasizes a need to add more data and include more taxa to resolve this problem. In addition, the molecular dating provides insights into the time scale of the origin and diversification of the Mustelidae subfamilies. Finally, the phylogenetic performances and limits of nuclear introns and mt genes are discussed in the context of Mustelidae phylogeny.</p> <p>Conclusion</p> <p>Our study not only brings new perspectives on the previously obscured phylogenetic relationships among Mustelidae subfamilies, but also provides another example demonstrating the effectiveness of nuclear non-coding loci for reconstructing evolutionary histories in a group that has undergone rapid bursts of speciation.</p