Population genomic response to geographic gradients by widespread and endemic fishes of the Arabian Peninsula.

Genetic structure within marine species may be driven by local adaptation to their environment, or alternatively by historical processes, such as geographic isolation. The gulfs and seas bordering the Arabian Peninsula offer an ideal setting to examine connectivity patterns in coral reef fishes with respect to environmental gradients and vicariance. The Red Sea is characterized by a unique marine fauna, historical periods of desiccation and isolation, as well as environmental gradients in salinity, temperature, and primary productivity that vary both by latitude and by season. The adjacent Arabian Sea is characterized by a sharper environmental gradient, ranging from extensive coral cover and warm temperatures in the southwest, to sparse coral cover, cooler temperatures, and seasonal upwelling in the northeast. Reef fish, however, are not confined to these seas, with some Red Sea fishes extending varying distances into the northern Arabian Sea, while their pelagic larvae are presumably capable of much greater dispersal. These species must therefore cope with a diversity of conditions that invoke the possibility of steep clines in natural selection. Here, we test for genetic structure in two widespread reef fish species (a butterflyfish and surgeonfish) and eight range-restricted butterflyfishes across the Red Sea and Arabian Sea using genome-wide single nucleotide polymorphisms. We performed multiple matrix regression with randomization analyses on genetic distances for all species, as well as reconstructed scenarios for population subdivision in the species with signatures of isolation. We found that (a) widespread species displayed more genetic subdivision than regional endemics and (b) this genetic structure was not correlated with contemporary environmental parameters but instead may reflect historical events. We propose that the endemic species may be adapted to a diversity of local conditions, but the widespread species are instead subject to ecological filtering where different combinations of genotypes persist under divergent ecological regimes.

Exploration of the Biosynthetic Potential of the Populus Microbiome.

Natural products (NPs) isolated from bacteria have dramatically advanced human society, especially in medicine and agriculture. The rapidity and ease of genome sequencing have enabled bioinformatics-guided NP discovery and characterization. As a result, NP potential and diversity within a complex community, such as the microbiome of a plant, are rapidly expanding areas of scientific exploration. Here, we assess biosynthetic diversity in the Populus microbiome by analyzing both bacterial isolate genomes and metagenome samples. We utilize the fully sequenced genomes of isolates from the Populus root microbiome to characterize a subset of organisms for NP potential. The more than 3,400 individual gene clusters identified in 339 bacterial isolates, including 173 newly sequenced organisms, were diverse across NP types and distinct from known NP clusters. The ribosomally synthesized and posttranslationally modified peptides were both widespread and divergent from previously characterized molecules. Lactones and siderophores were prevalent in the genomes, suggesting a high level of communication and pressure to compete for resources. We then consider the overall bacterial diversity and NP variety of metagenome samples compared to the sequenced isolate collection and other plant microbiomes. The sequenced collection, curated to reflect the phylogenetic diversity of the Populus microbiome, also reflects the overall NP diversity trends seen in the metagenomic samples. In our study, only about 1% of all clusters from sequenced isolates were positively matched to a previously characterized gene cluster, suggesting a great opportunity for the discovery of novel NPs involved in communication and control in the Populus root microbiome. IMPORTANCE The plant root microbiome is one of the most diverse and abundant biological communities known. Plant-associated bacteria can have a profound effect on plant growth and development, and especially on protection from disease and environmental stress. These organisms are also known to be a rich source of antibiotic and antifungal drugs. In order to better understand the ways bacterial communities influence plant health, we evaluated the diversity and uniqueness of the natural product gene clusters in bacteria isolated from poplar trees. The complex molecule clusters are abundant, and the majority are unique, suggesting a great potential to discover new molecules that could not only affect plant health but also could have applications as antibiotic agents.

Genus-Wide Assessment of Lignocellulose Utilization in the Extremely Thermophilic Genus Caldicellulosiruptor by Genomic, Pangenomic, and Metagenomic Analyses.

Metagenomic data from Obsidian Pool (Yellowstone National Park, USA) and 13 genome sequences were used to reassess genus-wide biodiversity for the extremely thermophilic Caldicellulosiruptor The updated core genome contains 1,401 ortholog groups (average genome size for 13 species = 2,516 genes). The pangenome, which remains open with a revised total of 3,493 ortholog groups, encodes a variety of multidomain glycoside hydrolases (GHs). These include three cellulases with GH48 domains that are colocated in the glucan degradation locus (GDL) and are specific determinants for microcrystalline cellulose utilization. Three recently sequenced species, Caldicellulosiruptor sp. strain Rt8.B8 (renamed here Caldicellulosiruptor morganii), Thermoanaerobacter cellulolyticus strain NA10 (renamed here Caldicellulosiruptor naganoensis), and Caldicellulosiruptor sp. strain Wai35.B1 (renamed here Caldicellulosiruptor danielii), degraded Avicel and lignocellulose (switchgrass). C. morganii was more efficient than Caldicellulosiruptor bescii in this regard and differed from the other 12 species examined, both based on genome content and organization and in the specific domain features of conserved GHs. Metagenomic analysis of lignocellulose-enriched samples from Obsidian Pool revealed limited new information on genus biodiversity. Enrichments yielded genomic signatures closely related to that of Caldicellulosiruptor obsidiansis, but there was also evidence for other thermophilic fermentative anaerobes (Caldanaerobacter, Fervidobacterium, Caloramator, and Clostridium). One enrichment, containing 89.8% Caldicellulosiruptor and 9.7% Caloramator, had a capacity for switchgrass solubilization comparable to that of C. bescii These results refine the known biodiversity of Caldicellulosiruptor and indicate that microcrystalline cellulose degradation at temperatures above 70°C, based on current information, is limited to certain members of this genus that produce GH48 domain-containing enzymes.IMPORTANCE The genus Caldicellulosiruptor contains the most thermophilic bacteria capable of lignocellulose deconstruction, which are promising candidates for consolidated bioprocessing for the production of biofuels and bio-based chemicals. The focus here is on the extant capability of this genus for plant biomass degradation and the extent to which this can be inferred from the core and pangenomes, based on analysis of 13 species and metagenomic sequence information from environmental samples. Key to microcrystalline hydrolysis is the content of the glucan degradation locus (GDL), a set of genes encoding glycoside hydrolases (GHs), several of which have GH48 and family 3 carbohydrate binding module domains, that function as primary cellulases. Resolving the relationship between the GDL and lignocellulose degradation will inform efforts to identify more prolific members of the genus and to develop metabolic engineering strategies to improve this characteristic.

Russulaceae: a new genomic dataset to study ecosystem function and evolutionary diversification of ectomycorrhizal fungi with their tree associates.

The family Russulaceae is considered an iconic lineage of mostly mushroom-forming basidiomycetes due to their importance as edible mushrooms in many parts of the world, and their ubiquity as ectomycorrhizal symbionts in both temperate and tropical forested biomes. Although much research has been focused on this group, a comprehensive or cohesive synthesis by which to understand the functional diversity of the group has yet to develop. Interest in ectomycorrhizal fungi, of which Russulaceae is a key lineage, is prodigious due to the important roles they play as plant root mutualists in ecosystem functioning, global carbon sequestration, and a potential role in technology development toward environmental sustainability. As one of the most species-diverse ectomycorrhizal lineages, the Russulaceae has recently been the focus of a dense sampling and genome sequencing initiative with the Joint Genome Institute aimed at untangling their functional roles and testing whether functional niche specialization exists for independent lineages of ectomycorrhizal fungi. Here we present a review of important studies on this group to contextualize what we know about its members' evolutionary history and ecosystem functions, as well as to generate hypotheses establishing the Russulaceae as a valuable experimental system.

Draft genome of an iconic Red Sea reef fish, the blacktail butterflyfish (Chaetodon austriacus): current status and its characteristics.

Butterflyfish are among the most iconic of the coral reef fishes and represent a model system to study general questions of biogeography, evolution and population genetics. We assembled and annotated the genome sequence of the blacktail butterflyfish (Chaetodon austriacus), an Arabian region endemic species that is reliant on coral reefs for food and shelter. Using available bony fish (superclass Osteichthyes) genomes as a reference, a total of 28 926 high-quality protein-coding genes were predicted from 13 967 assembled scaffolds. The quality and completeness of the draft genome of C. austriacus suggest that it has the potential to serve as a resource for studies on the co-evolution of reef fish adaptations to the unique Red Sea environment, as well as a comparison of gene sequences between closely related congeneric species of butterflyfish distributed more broadly across the tropical Indo-Pacific.

Using a butterflyfish genome as a general tool for RAD-Seq studies in specialized reef fish.

Data from a large-scale restriction site-associated DNA sequencing (RAD-Seq) study of nine butterflyfish species in the Red Sea and Arabian Sea provided a means to test the utility of a recently published draft genome (Chaetodon austriacus) and assess apparent bias in this method of isolating nuclear loci. We here processed double-digest restriction site-associated DNA (ddRAD) sequencing data to identify single nucleotide polymorphism (SNP) markers and their associated function with and without our reference genome to see whether it improves the quality of RAD-Seq. Our analyses indicate (i) a modest gap between the number of nonannotated versus annotated SNPs across all species, (ii) an advantage of using genomic resources for closely related but not distantly related butterflyfish species based on the ability to assign putative gene function to SNPs and (iii) an enrichment of genes among sister butterflyfish taxa related to calcium transmembrane transport and binding. The latter result highlights the potential for this approach to reveal insights into adaptive mechanisms in populations inhabiting challenging coral reef environments such as the Red Sea, Arabian Sea and Arabian Gulf with further study.

Seascape genetics along environmental gradients in the Arabian Peninsula: insights from ddRAD sequencing of anemonefishes.

Understanding the processes that shape patterns of genetic structure across space is a central aim of landscape genetics. However, it remains unclear how geographical features and environmental variables shape gene flow, particularly for marine species in large complex seascapes. Here, we evaluated the genomic composition of the two-band anemonefish Amphiprion bicinctus across its entire geographical range in the Red Sea and Gulf of Aden, as well as its close relative, Amphiprion omanensis endemic to the southern coast of Oman. Both the Red Sea and the Arabian Sea are complex and environmentally heterogeneous marine systems that provide an ideal scenario to address these questions. Our findings confirm the presence of two genetic clusters previously reported for A. bicinctus in the Red Sea. Genetic structure analyses suggest a complex seascape configuration, with evidence of both isolation by distance (IBD) and isolation by environment (IBE). In addition to IBD and IBE, genetic structure among sites was best explained when two barriers to gene flow were also accounted for. One of these coincides with a strong oligotrophic-eutrophic gradient at around 16-20˚N in the Red Sea. The other agrees with a historical bathymetric barrier at the straight of Bab al Mandab. Finally, these data support the presence of interspecific hybrids at an intermediate suture zone at Socotra and indicate complex patterns of genomic admixture in the Gulf of Aden with evidence of introgression between species. Our findings highlight the power of recent genomic approaches to resolve subtle patterns of gene flow in marine seascapes.

Viral metagenomics: analysis of begomoviruses by illumina high-throughput sequencing.

Traditional DNA sequencing methods are inefficient, lack the ability to discern the least abundant viral sequences, and ineffective for determining the extent of variability in viral populations. Here, populations of single-stranded DNA plant begomoviral genomes and their associated beta- and alpha-satellite molecules (virus-satellite complexes) (genus, Begomovirus; family, Geminiviridae) were enriched from total nucleic acids isolated from symptomatic, field-infected plants, using rolling circle amplification (RCA). Enriched virus-satellite complexes were subjected to Illumina-Next Generation Sequencing (NGS). CASAVA and SeqMan NGen programs were implemented, respectively, for quality control and for de novo and reference-guided contig assembly of viral-satellite sequences. The authenticity of the begomoviral sequences, and the reproducibility of the Illumina-NGS approach for begomoviral deep sequencing projects, were validated by comparing NGS results with those obtained using traditional molecular cloning and Sanger sequencing of viral components and satellite DNAs, also enriched by RCA or amplified by polymerase chain reaction. As the use of NGS approaches, together with advances in software development, make possible deep sequence coverage at a lower cost; the approach described herein will streamline the exploration of begomovirus diversity and population structure from naturally infected plants, irrespective of viral abundance. This is the first report of the implementation of Illumina-NGS to explore the diversity and identify begomoviral-satellite SNPs directly from plants naturally-infected with begomoviruses under field conditions.

Activities and specificities of homodimeric TALENs in Saccharomyces cerevisiae.

The development of highly efficient genome engineering reagents is of paramount importance to launch the next wave of biotechnology. TAL effectors have been developed as an adaptable DNA binding scaffold that can be engineered to bind to any user-defined sequence. Thus, TAL-based DNA binding modules have been used to generate chimeric proteins for a variety of targeted genome modifications across eukaryotic species. For example, TAL effectors fused to the catalytic domain of FokI endonuclease (TALENs) were used to generate site-specific double strand breaks (DSBs), the repair of which can be harnessed to dictate user-desired, genome-editing outcomes. To cleave DNA, FokI endonuclease must dimerize which can be achieved using a pair of TALENs that bind to the DNA targeted in a tail-to-tail orientation with proper spacing allowing the dimer formation. Because TALENs binding to DNA are dependent on their repeat sequences and nucleotides binding specificities, homodimers and heterodimers binding can be formed. In the present study, we used several TALEN monomers with increased repeats binding degeneracy to allow homodimer formation at increased number of genomic loci. We assessed their binding specificities and genome modification activities. Our results indicate that homodimeric TALENs could be used to modify the yeast genome in a site-specific manner and their binding to the promoter regions might modulate the expression of target genes. Taken together, our data indicate that homodimeric TALENs could be used to achieve different engineering possibilities of biotechnological applications and that their transcriptional modulations need to be considered when analyzing their phenotypic effects.

Simplified method for predicting a functional class of proteins in transcription factor complexes.

BACKGROUND: Initiation of transcription is essential for most of the cellular responses to environmental conditions and for cell and tissue specificity. This process is regulated through numerous proteins, their ligands and mutual interactions, as well as interactions with DNA. The key such regulatory proteins are transcription factors (TFs) and transcription co-factors (TcoFs). TcoFs are important since they modulate the transcription initiation process through interaction with TFs. In eukaryotes, transcription requires that TFs form different protein complexes with various nuclear proteins. To better understand transcription regulation, it is important to know the functional class of proteins interacting with TFs during transcription initiation. Such information is not fully available, since not all proteins that act as TFs or TcoFs are yet annotated as such, due to generally partial functional annotation of proteins. In this study we have developed a method to predict, using only sequence composition of the interacting proteins, the functional class of human TF binding partners to be (i) TF, (ii) TcoF, or (iii) other nuclear protein. This allows for complementing the annotation of the currently known pool of nuclear proteins. Since only the knowledge of protein sequences is required in addition to protein interaction, the method should be easily applicable to many species. RESULTS: Based on experimentally validated interactions between human TFs with different TFs, TcoFs and other nuclear proteins, our two classification systems (implemented as a web-based application) achieve high accuracies in distinguishing TFs and TcoFs from other nuclear proteins, and TFs from TcoFs respectively. CONCLUSION: As demonstrated, given the fact that two proteins are capable of forming direct physical interactions and using only information about their sequence composition, we have developed a completely new method for predicting a functional class of TF interacting protein partners with high precision and accuracy.

Rapid and highly efficient construction of TALE-based transcriptional regulators and nucleases for genome modification.

Transcription activator-like effectors (TALEs) can be used as DNA-targeting modules by engineering their repeat domains to dictate user-selected sequence specificity. TALEs have been shown to function as site-specific transcriptional activators in a variety of cell types and organisms. TALE nucleases (TALENs), generated by fusing the FokI cleavage domain to TALE, have been used to create genomic double-strand breaks. The identity of the TALE repeat variable di-residues, their number, and their order dictate the DNA sequence specificity. Because TALE repeats are nearly identical, their assembly by cloning or even by synthesis is challenging and time consuming. Here, we report the development and use of a rapid and straightforward approach for the construction of designer TALE (dTALE) activators and nucleases with user-selected DNA target specificity. Using our plasmid set of 100 repeat modules, researchers can assemble repeat domains for any 14-nucleotide target sequence in one sequential restriction-ligation cloning step and in only 24 h. We generated several custom dTALEs and dTALENs with new target sequence specificities and validated their function by transient expression in tobacco leaves and in vitro DNA cleavage assays, respectively. Moreover, we developed a web tool, called idTALE, to facilitate the design of dTALENs and the identification of their genomic targets and potential off-targets in the genomes of several model species. Our dTALE repeat assembly approach along with the web tool idTALE will expedite genome-engineering applications in a variety of cell types and organisms including plants.