Non-model Organisms Research Articles

Tick hemocytes have a pleiotropic role in microbial infection and arthropod fitness

Uncovering the complexity of systems in non-model organisms is critical for understanding arthropod immunology. Prior efforts have mostly focused on Dipteran insects, which only account for a subset of existing arthropod species in nature. Here we use and develop advanced techniques to describe immune cells (hemocytes) from the clinically relevant tick Ixodes scapularis at a single-cell resolution. We observe molecular alterations in hemocytes upon feeding and infection with either the Lyme disease spirochete Borrelia burgdorferi or the rickettsial agent Anaplasma phagocytophilum. We reveal hemocyte clusters exhibiting defined signatures related to immunity, metabolism, and proliferation. Depletion of phagocytic hemocytes affects hemocytin and astakine levels, two I. scapularis hemocyte markers, impacting blood-feeding, molting behavior, and bacterial acquisition. Mechanistically, astakine alters hemocyte proliferation, whereas hemocytin affects the c-Jun N-terminal kinase (JNK) signaling pathway in I. scapularis. Altogether, we discover a role for tick hemocytes in immunophysiology and provide a valuable resource for comparative biology in arthropods.

Constructing the metabolic network of wheat kernels based on structure-guided chemical modification and multi-omics data

Metabolic network construction plays a pivotal role in unraveling the regulatory mechanism of biological activities, although it often proves to be challenging and labor-intensive, particularly with non-model organisms. In this study, we develop a computational approach that employs reaction models based on the structure-guided chemical modification and related compounds to construct a metabolic network in wheat. This construction results in a comprehensive structure-guided network, including 625 identified metabolites and additional 333 putative reactions compared with the Kyoto Encyclopedia of Genes and Genomes database. Using a combination of gene annotation, reaction classification, structure similarity, and correlations from transcriptome and metabolome analysis, a total of 229 potential genes related to these reactions are identified within this network. To validate the network, the functionality of a hydroxycinnamoyltransferase (TraesCS3D01G314900) for the synthesis of polyphenols and a rhamnosyltransferase (TraesCS2D01G078700) for the modification of flavonoids are verified through in vitro enzymatic studies and wheat mutant tests, respectively. Our research thus supports the utility of structure-guided chemical modification as an effective tool in identifying causal candidate genes for constructing metabolic networks and further in metabolomic genetic studies.

BSXplorer: analytical framework for exploratory analysis of BS-seq data

BackgroundBisulfite sequencing detects and quantifies DNA methylation patterns, contributing to our understanding of gene expression regulation, genome stability maintenance, conservation of epigenetic mechanisms across divergent taxa, epigenetic inheritance and, eventually, phenotypic variation. Graphical representation of methylation data is crucial in exploring epigenetic regulation on a genome-wide scale in both plants and animals. This is especially relevant for non-model organisms with poorly annotated genomes and/or organisms where genome sequences are not yet assembled on chromosome level. Despite being a technology of choice to profile DNA methylation for many years now there are surprisingly few lightweight and robust standalone tools available for efficient graphical analysis of data in non-model systems. This significantly limits evolutionary studies and agrigenomics research. BSXplorer is a tool specifically developed to fill this gap and assist researchers in explorative data analysis and in visualising and interpreting bisulfite sequencing data more easily.ResultsBSXplorer provides in-depth graphical analysis of sequencing data encompassing (a) profiling of methylation levels in metagenes or in user-defined regions using line plots and heatmaps, generation of summary statistics charts, (b) enabling comparative analyses of methylation patterns across experimental samples, methylation contexts and species, and (c) identification of modules sharing similar methylation signatures at functional genomic elements. The tool processes methylation data quickly and offers API and CLI capabilities, along with the ability to create high-quality figures suitable for publication.ConclusionsBSXplorer facilitates efficient methylation data mining, contrasting and visualization, making it an easy-to-use package that is highly useful for epigenetic research.

CRISPR/Cas9-mediated knockout of scarlet gene produces eye color mutants in the soybean looper, Chrysodeixis includens.

The CRISPR/Cas9 technology has greatly progressed research on non-model organisms, demonstrating successful applications in genome editing for various insects. However, its utilization in the case of the soybean looper, Chrysodeixis includens, a notable pest affecting soybean crops, has not been explored due to constraints such as limited genomic information and the embryonic microinjection technique. This study presents successful outcomes in generating heritable knockout mutants for a pigment transporter gene, scarlet, in C. includens through CRISPR/Cas9-mediated mutagenesis. The scarlet locus identified in the genome assembly of C. includens consists of 14 exons, with a coding sequence extending for 1,986 bp. Two single guide RNAs (sgRNAs) were designed to target the first exon of scarlet. Microinjection of these two sgRNAs along with the Cas9 protein into fresh embryos resulted in the successful production of variable phenotypes, particularly mutant eyes. The observed mutation rate accounted for about 16%. Genotype analysis revealed diverse indel mutations at the target site, presumably originating from double-strand breaks followed by the nonhomologous end joining repair, leading to a premature stop codon due to frame shift. Single-pair mating of the mutant moths produced G1 offspring, and the establishment of a homozygous mutant strain occurred in G2. The mutant moths exhibited lightly greenish or yellowish compound eyes in both sexes, confirming the involvement of scarlet in pigmentation in C. includens. Notably, the CRISPR/Cas9-mediated genome editing technique serves as a visible phenotypic marker, demonstrating its proof-of-concept applicability in C. includens, as other pigment transporter genes have been utilized as visible markers to establish genetic control for various insects. These results provide the first successful case that the CRISPR/Cas9 method effectively induces mutations in C. includes, an economically important soybean insect pest.

Annotated genome and transcriptome of the endangered Caribbean mountainous star coral (Orbicella faveolata) using PacBio long-read sequencing

Long-read sequencing is revolutionizing de-novo genome assemblies, with continued advancements making it more readily available for previously understudied, non-model organisms. Stony corals are one such example, with long-read de-novo genome assemblies now starting to be publicly available, opening the door for a wide array of ‘omics-based research. Here we present a new de-novo genome assembly for the endangered Caribbean star coral, Orbicella faveolata, using PacBio circular consensus reads. Our genome assembly improved the contiguity (51 versus 1,933 contigs) and complete and single copy BUSCO orthologs (93.6% versus 85.3%, database metazoa_odb10), compared to the currently available reference genome generated using short-read methodologies. Our new de-novo assembled genome also showed comparable quality metrics to other coral long-read genomes. Telomeric repeat analysis identified putative chromosomes in our scaffolded assembly, with these repeats at either one, or both ends, of scaffolded contigs. We identified 32,172 protein coding genes in our assembly through use of long-read RNA sequencing (ISO-seq) of additional O. faveolata fragments exposed to a range of abiotic and biotic treatments, and publicly available short-read RNA-seq data. With anthropogenic influences heavily affecting O. faveolata, as well as its increasing incorporation into reef restoration activities, this updated genome resource can be used for population genomics and other ‘omics analyses to aid in the conservation of this species.

Preliminary analysis of pathways and their implications during salinity stress in abalone

Transcriptome sequencing has offered immense opportunities to study non-model organisms. Abalone is an important marine mollusk that encounters harsh environmental conditions in its natural habitat and under aquaculture conditions; hence, research that increases molecular information to understand abalone physiology and stress response is noteworthy. Accordingly, the study used transcriptome sequencing of the gill tissues of abalone exposed to low salinity stress. The aim is to explore some enriched pathways during salinity stress and the crosstalk and functions of the genes involved in the candidate biological processes for future further analysis of their expression patterns.The data suggest that abalone genes such as YAP/TAZ, Myc, Nkd, and Axin (involved in the Hippo signaling pathway) and PI3K/Akt, SHC, and RTK (involved in the Ras signaling pathways) might mediate growth and development. Thus, deregulation of the Hippo and Ras pathways by salinity stress could be a possible mechanism by which unfavorable salinities influence growth in abalone. Furthermore, PEPCK, GYS, and PLC genes (mediating the Glucagon signaling pathway) might be necessary for glucose homeostasis, reproduction, and abalone meat sensory qualities; hence, a need to investigate how they might be influenced by environmental stress. Genes such as MYD88, IRAK1/4, JNK, AP-1, and TRAF6 (mediating the MAPK signaling pathway) could be useful in understanding abalone's innate immune response to environmental stresses. Finally, the aminoacyl-tRNA biosynthesis pathway hints at the mechanism by which new raw materials for protein biosynthesis are mobilized for physiological processes and how abalone might respond to this process during salinity stress. Low salinity clearly regulated genes in these pathways in a time-dependent manner, as hinted by the heat maps. In the future, qRT-PCR verification and in-depth study of the various genes and proteins discussed would provide enormous molecular information resources for the abalone biology.

XPORT ENTRAP: A droplet microfluidic platform for enhanced DNA transfer between microbial species

A significant hurdle for the widespread implementation and use of synthetic biology is the challenge of highly efficient introduction of DNA into microorganisms. This is especially a barrier for the utilization of non-model organisms and/or novel chassis species for a variety of applications, ranging from molecular biology to biotechnology and biomanufacturing applications. Common approaches to episomal and chromosomal gene editing, which employ techniques such as chemical competence and electroporation, are typically only amenable to a small subset of microbial species while leaving the vast majority of microorganisms in nature genetically inaccessible. To address this challenge, we have employed the previously described B. subtilis broad-host conjugation strain, XPORT, which was modularly designed for loading DNA cargo and conjugating such DNA into recalcitrant microbes. In this current work, we have leveraged and adapted the XPORT strain for use in a droplet microfluidic platform to enable increased efficiency of conjugation-based DNA transfer. The system named DNA ENTRAP (DNA ENhanced TRAnsfer Platform) utilizes cell-encapsulated water-in-oil emulsion droplets as pico-liter-volume bioreactors that allows controlled contacts between the donor and receiver cells within the emulsion bioreactor. This allowed enhanced XPORT-mediated genetic transfer over the current benchtop XPORT process, demonstrated using two different Bacillus subtilis strains (donor and receiver), as well as increased throughput (e.g., number of successfully conjugated cells) due to the automated assay steps inherent to microfluidic lab-on-a-chip systems. DNA ENTRAP paves the way for a streamlined automation of culturing and XPORT-mediated genetic transfer processes as well as future high-throughput cell engineering and screening applications.

Development of the Static and Dynamic Gene Expression Regulation Toolkit in Pseudomonas chlororaphis.

The advancement of metabolic engineering and synthetic biology has promoted in-depth research on the nonmodel microbial metabolism, and the potential of nonmodel organisms in industrial biotechnology is becoming increasingly evident. The nonmodel organism Pseudomonas chlororaphis is a safe plant growth promoting bacterium for the production of phenazine compounds; however, its application is seriously hindered due to the lack of an effective gene expression precise regulation toolkit. In this study, we constructed a library of 108 promoter-5'-UTR (PUTR) and characterized them through fluorescent protein detection. Then, 6 PUTRs with stable low, intermediate, and high intensities were further characterized by report genes lacZ encoding β-galactosidase from Escherichia coli K12 and phzO encoding PCA monooxygenase from P. chlororaphis GP72 and thus developed as a static gene expression regulation system. Furthermore, the stable and high-intensity expressed PMOK_RS0128085UTR was fused with the LacO operator to construct an IPTG-induced plasmid, and a self-induced plasmid was constructed employing the high-intensity PMOK_RS0116635UTR regulated by cell density, resulting in a dynamic gene expression regulation system. In summary, this study established two sets of static and dynamic regulatory systems for P. chlororaphis, providing an effective toolkit for fine-tuning gene expression and reprograming the metabolism flux.

MINNO: An Open Source Software for Refining Metabolic Networks and Investigating Complex Network Activity Using Empirical Metabolomics Data.

Metabolomics is a powerful tool for uncovering biochemical diversity in a wide range of organisms. Metabolic network modeling is commonly used to frame metabolomics data in the context of a broader biological system. However, network modeling of poorly characterized nonmodel organisms remains challenging due to gene homology mismatches which lead to network architecture errors. To address this, we developed the Metabolic Interactive Nodular Network for Omics (MINNO), a web-based mapping tool that uses empirical metabolomics data to refine metabolic networks. MINNO allows users to create, modify, and interact with metabolic pathway visualizations for thousands of organisms, in both individual and multispecies contexts. Herein, we illustrate the use of MINNO in elucidating the metabolic networks of understudied species, such as those of the Borrelia genus, which cause Lyme and relapsing fever diseases. Using a hybrid genomics-metabolomics modeling approach, we constructed species-specific metabolic networks for threeBorrelia species. Using these empirically refined networks, we were able to metabolically differentiate these species via their nucleotide metabolism, which cannot be predicted from genomic networks. Additionally, using MINNO, we identified 18 missing reactions from the KEGG database, of which nine were supported by the primary literature. These examples illustrate the use of metabolomics for the empirical refining of genetically constructed networks and show how MINNO can be used to study nonmodel organisms.

CryptoCEN: A Co-Expression Network for Cryptococcus neoformans reveals novel proteins involved in DNA damage repair.

Elucidating gene function is a major goal in biology, especially among non-model organisms. However, doing so is complicated by the fact that molecular conservation does not always mirror functional conservation, and that complex relationships among genes are responsible for encoding pathways and higher-order biological processes. Co-expression, a promising approach for predicting gene function, relies on the general principal that genes with similar expression patterns across multiple conditions will likely be involved in the same biological process. For Cryptococcus neoformans, a prevalent human fungal pathogen greatly diverged from model yeasts, approximately 60% of the predicted genes in the genome lack functional annotations. Here, we leveraged a large amount of publicly available transcriptomic data to generate a C. neoformans Co-Expression Network (CryptoCEN), successfully recapitulating known protein networks, predicting gene function, and enabling insights into the principles influencing co-expression. With 100% predictive accuracy, we used CryptoCEN to identify 13 new DNA damage response genes, underscoring the utility of guilt-by-association for determining gene function. Overall, co-expression is a powerful tool for uncovering gene function, and decreases the experimental tests needed to identify functions for currently under-annotated genes.

In vitro proliferation of Mytilus edulis male germ cell progenitors.

Our understanding of basic cellular processes has mostly been provided by mammalian cell culture, and by some non-mammalian vertebrate and few invertebrate cell culture models. Developing reliable culture conditions for non-model organisms is essential to allow investigation of more unusual cellular processes. Here, we investigate how cells isolated from different tissues of the marine mussel Mytilus edulis thrive and survive in vitro in the hope of establishing a suitable laboratory model for the investigation of cellular mechanisms specific to these bivalve mollusks. We found that cells dissociated from mantle tissue attached to the culture vessels and proliferated well in vitro, whereas cells isolated from gills, although remaining viable, did not maintain divisions over three to four weeks in culture. We used antibodies against the germ-line marker DEAD-box helicase 4 (DDX4), also known as VASA, and the epithelial cell marker cytokeratin to distinguish different cell types in culture. DDX4-positive cells were predominant in 25-day-old cultures from male mantles. Cells from other tissues remained in low numbers and did not seem to change in composition over time. Overall, the culture conditions described here allow an efficient selection of male germ cells that could be used to study specific cellular mechanisms in vitro.

Mining NCBI Sequence Read Archive Database: An Untapped Source of Organelle Genomes for Taxonomic and Comparative Genomics Research

The NCBI SRA database is constantly expanding due to the large amount of genomic and transcriptomic data from various organisms generated by next-generation sequencing, and re-searchers worldwide regularly deposit new data into the database. This high-coverage genomic and transcriptomic information can be re-evaluated regardless of the original research subject. The database-deposited NGS data can offer valuable insights into the genomes of organelles, particularly for non-model organisms. Here, we developed an automated bioinformatics workflow called “OrgaMiner”, designed to unveil high-quality mitochondrial and chloroplast genomes by data mining the NCBI SRA database. OrgaMiner, a Python-based pipeline, automatically orchestrates various tools to extract, assemble, and annotate organelle genomes for non-model organisms without available organelle genome sequences but with data in the NCBI SRA. To test the usability and feasibility of the pipeline, “mollusca” was selected as a keyword, and 76 new mitochondrial genomes were de novo assembled and annotated automatically without writing one single code. The applicability of the pipeline can be expanded to identify organelles in diverse invertebrate, vertebrate, and plant species by simply specifying the taxonomic name. OrgaMiner provides an easy-to-use, end-to-end solution for biologists mainly working with taxonomy and population genetics.

Food Waste from Campus Dining Hall as a Potential Feedstock for 2,3-Butanediol Production via Non-Sterilized Fermentation.

Food waste is a major issue that is increasingly affecting our environment. More than one-third of food is wasted, resulting in over $400 billion in losses to the U.S. economy. While composting and other small recycling practices are encouraged from person-to-person, it is not enough to balance the net loss of 80 million tons per year. Currently, one of the most promising routes for reducing food waste is through microbial fermentation, which can convert the waste into valuable bioproducts. Among the compounds produced from fermentation, 2,3-butanediol (2,3-BDO) has gained interest recently due to its molecular structure as a building block for many other derivatives used in perfumes, synthetic rubber, fumigants, antifreeze agents, fuel additives, and pharmaceuticals. Waste feedstocks, such as food waste, are a potential source of renewable energy due to their lack of cost and availability. Food waste also possesses microbial requirements for growth such as carbohydrates, proteins, fats, and more. However, food waste is highly inconsistent and the variability in composition may hinder its ability to be a stable source for bioproducts such as 2,3-BDO. This current study focuses specifically on post-consumer food waste and how 2,3-BDO can be produced through a non-model organism, Bacillus licheniformis YNP5-TSU during non-sterile fermentation. From the dining hall at Tennessee State University, 13 food waste samples were collected over a 6-month period and the compositional analysis was performed. On average, these samples consisted of fat (19.7%), protein (18.7%), ash (4.8%), fiber (3.4%), starch (27.1%), and soluble sugars (20.9%) on a dry basis with an average moisture content of 34.7%. Food waste samples were also assessed for their potential production of 2,3-BDO during non-sterile thermophilic fermentation, resulting in a max titer of 12.12 g/L and a 33% g/g yield of 2,3-BDO/carbohydrates. These findings are promising and can lead to the better understanding of food waste as a defined feedstock for 2,3-BDO and other fermentation end-products.

Genetic basis and expression of ventral colour in polymorphic common lizards.

Colour is an important visual cue that can correlate with sex, behaviour, life history or ecological strategies, and has evolved divergently and convergently across animal lineages. Its genetic basis in non-model organisms is rarely known, but such information is vital for determining the drivers and mechanisms of colour evolution. Leveraging genetic admixture in a rare contact zone between oviparous and viviparous common lizards (Zootoca vivipara), we show that females (N = 558) of the two otherwise morphologically indistinguishable reproductive modes differ in their ventral colouration (from pale to vibrant yellow) and intensity of melanic patterning. We find no association between female colouration and reproductive investment, and no evidence for selection on colour. Using a combination of genetic mapping and transcriptomic evidence, we identified two candidate genes associated with ventral colour differentiation, DGAT2 and PMEL. These are genes known to be involved in carotenoid metabolism and melanin synthesis respectively. Ventral melanic spots were associated with two genomic regions, including a SNP close to protein tyrosine phosphatase (PTP) genes. Using genome re-sequencing data, our results show that fixed coding mutations in the candidate genes cannot account for differences in colouration. Taken together, our findings show that the evolution of ventral colouration and its associations across common lizard lineages is variable. A potential genetic mechanism explaining the flexibility of ventral colouration may be that colouration in common lizards, but also across squamates, is predominantly driven by regulatory genetic variation.

High-fidelity (repeat) consensus sequences from short reads using combined read clustering and assembly

BackgroundDespite the many cheap and fast ways to generate genomic data, good and exact genome assembly is still a problem, with especially the repeats being vastly underrepresented and often misassembled. As short reads in low coverage are already sufficient to represent the repeat landscape of any given genome, many read cluster algorithms were brought forward that provide repeat identification and classification. But how can trustworthy, reliable and representative repeat consensuses be derived from unassembled genomes?ResultsHere, we combine methods from repeat identification and genome assembly to derive these robust consensuses. We test several use cases, such as (1) consensus building from clustered short reads of non-model genomes, (2) from genome-wide amplification setups, and (3) specific repeat-centred questions, such as the linked vs. unlinked arrangement of ribosomal genes. In all our use cases, the derived consensuses are robust and representative. To evaluate overall performance, we compare our high-fidelity repeat consensuses to RepeatExplorer2-derived contigs and check, if they represent real transposable elements as found in long reads. Our results demonstrate that it is possible to generate useful, reliable and trustworthy consensuses from short reads by a combination from read cluster and genome assembly methods in an automatable way.ConclusionWe anticipate that our workflow opens the way towards more efficient and less manual repeat characterization and annotation, benefitting all genome studies, but especially those of non-model organisms.

One-step knock-in of two antimicrobial peptide transgenes at multiple loci of catfish by CRISPR/Cas9-mediated multiplex genome engineering

CRISPR/Cas9-mediated multiplex genome editing (MGE) conventionally uses multiple single-guide RNAs (sgRNAs) for gene-targeted mutagenesis via the non-homologous end joining (NHEJ) pathway. MGE has been proven to be highly efficient for functional gene disruption/knockout (KO) at multiple loci in mammalian cells or organisms. However, in the absence of a DNA donor, this approach is limited to small indels without transgene integration. Here, we establish the linear double-stranded DNA (dsDNA) and double-cut plasmid (dcPlasmid) combination-assisted MGE in channel catfish (Ictalurus punctatus), allowing combinational deletion mutagenesis and transgene knock-in (KI) at multiple sites through NHEJ/homology-directed repair (HDR) pathway in parallel. In this study, we used single-sgRNA-based genome editing (ssGE) and multi-sgRNA-based MGE (msMGE) to replace the luteinizing hormone (lh) and melanocortin-4 receptor (mc4r) genes with the cathelicidin (As-Cath) transgene and the myostatin (two target sites: mstn1, mstn2) gene with the cecropin (Cec) transgene, respectively. A total of 9000 embryos were microinjected from three families, and 1004 live fingerlings were generated and analyzed. There was no significant difference in hatchability (all P > 0.05) and fry survival (all P > 0.05) between ssGE and msMGE. Compared to ssGE, CRISPR/Cas9-mediated msMGE assisted by the mixture of dsDNA and dcPlasmid donors yielded a higher knock-in (KI) efficiency of As-Cath (19.93 %, [59/296] vs. 12.96 %, [45/347]; P = 0.018) and Cec (22.97 %, [68/296] vs. 10.80 %, [39/361]; P = 0.003) transgenes, respectively. The msMGE strategy can be used to generate transgenic fish carrying two transgenes at multiple loci. In addition, double and quadruple mutant individuals can be produced with high efficiency (36.3 % ∼ 71.1 %) in one-step microinjection. In conclusion, we demonstrated that the CRISPR/Cas9-mediated msMGE allows the one-step generation of simultaneous insertion of the As-Cath and Cec transgenes at four sites, and the simultaneous disruption of the lh, mc4r, mstn1 and mstn2 alleles. This msMGE system, aided by the mixture donors, promises to pioneer a new dimension in the drive and selection of multiple designated traits in other non-model organisms.

GeneMiner: A tool for extracting phylogenetic markers from next-generation sequencing data.

The advancement of next-generation sequencing (NGS) technologies has been revolutionary for the field of evolutionary biology. This technology has led to an abundance of available genomes and transcriptomes for researchers to mine. Specifically, researchers can mine for various types of molecular markers that are vital for phylogenetic, evolutionary and ecological studies. Numerous tools have been developed to extract these molecular markers from NGS data. However, due to an insufficient number of well-annotated reference genomes for non-model organisms, it remains challenging to obtain these markers accurately and efficiently. Here, we present GeneMiner, an improved and expanded version of our previous tool, Easy353. GeneMiner combines the reference-guided de Bruijn graph assembly with seed self-discovery and greedy extension. Additionally, it includes a verification step using a parameter-bootstrap method to reduce the pitfalls associated with using a relatively distant reference. Our results, using both experimental and simulation data, showed GeneMiner can accurately acquire phylogenetic molecular markers for plants using transcriptomic, genomic and other NGS data. GeneMiner is designed to be user-friendly, fast and memory-efficient. Further, it is compatible with Linux, Windows and macOS. All source codes are publicly available on GitHub (https://github.com/sculab/GeneMiner) and Gitee (https://gitee.com/sculab/GeneMiner) for easy accessibility and transparency.

Screening H3 Histone Acetylation in a Wild Bird, the House Sparrow (Passer Domesticus).

Epigenetic mechanisms are increasingly understood to have major impacts across ecology. However, one molecular epigenetic mechanism, DNA methylation, currently dominates the literature. A second mechanism, histone modification, is likely important to ecologically relevant phenotypes and thus warrants investigation, especially because molecular interplay between methylation and histone acetylation can strongly affect gene expression. There are a limited number of histone acetylation studies on non-model organisms, yet those that exist show that it can impact gene expression and phenotypic plasticity. Wild birds provide an excellent system to investigate histone acetylation, as free-living individuals must rapidly adjust to environmental change. Here, we screen histone acetylation in the house sparrow (Passer domesticus); we studied this species because DNA methylation was important in the spread of this bird globally. This species has one of the broadest geographic distributions in the world, and part of this success is related to the way that it uses methylation to regulate its gene expression. Here, we verify that a commercially available assay that was developed for mammals can be used in house sparrows. We detected high variance in histone acetylation among individuals in both liver and spleen tissue. Further, house sparrows with higher epigenetic potential in the Toll Like Receptor-4 (TLR-4) promoter (i.e., CpG content) had higher histone acetylation in liver. Also, there was a negative correlation between histone acetylation in spleen and TLR-4 expression. In addition to validating a method for measuring histone acetylation in wild songbirds, this study also shows that histone acetylation is related to epigenetic potential and gene expression, adding a new study option for ecological epigenetics.

Accurate Bayesian inference of sex chromosome karyotypes and sex-linked scaffolds from low-depth sequencing data.

The identification of sex-linked scaffolds and the genetic sex of individuals, i.e. their sex karyotype, is a fundamental step in population genomic studies. If sex-linked scaffolds are known, single individuals may be sexed based on read counts of next-generation sequencing data. If both sex-linked scaffolds as well as sex karyotypes are unknown, as is often the case for non-model organisms, they have to be jointly inferred. For both cases, current methods rely on arbitrary thresholds, which limits their power for low-depth data. In addition, most current methods are limited to euploid sex karyotypes (XX and XY). Here we develop BeXY, a fully Bayesian method to jointly infer the posterior probabilities for each scaffold to be autosomal, X- or Y-linked and for each individual to be any of the sex karyotypes XX, XY, X0, XXX, XXY, XYY and XXYY. If the sex-linked scaffolds are known, it also identifies autosomal trisomies and estimates the sex karyotype posterior probabilities for single individuals. As we show with downsampling experiments, BeXY has higher power than all existing methods. It accurately infers the sex karyotype of ancient human samples with as few as 20,000 reads and accurately infers sex-linked scaffolds from data sets of just a handful of samples or with highly imbalanced sex ratios, also in the case of low-quality reference assemblies. We illustrate the power of BeXY by applying it to both whole-genome shotgun and target enrichment sequencing data of ancient and modern humans, as well as several non-model organisms.

Measuring Mitotic Spindle and Microtubule Dynamics in Marine Embryos and Non-model Organisms.

During eukaryotic cell division a microtubule-based structure, the mitotic spindle, aligns and segregates chromosomes between daughter cells. Understanding how this cellular structure is assembled and coordinated in space and in time requires measuring microtubule dynamics and visualizing spindle assembly with high temporal and spatial resolution. Visualization is often achieved by the introduction and the detection of molecular probes and fluorescence microscopy. Microtubules and mitotic spindles are highly conserved across eukaryotes; however, several technical limitations have restricted these investigations to only a few species. The ability to monitor microtubule and chromosome choreography in a wide range of species is fundamental to reveal conserved mechanisms or unravel unconventional strategies that certain forms of life have developed to ensure faithful partitioning of chromosomes during cell division. Here, we describe a technique based on injection of purified proteins that enables the visualization of microtubules and chromosomes with a high contrast in several divergent marine embryos. We also provide analysis methods and tools to extract microtubule dynamics and monitor spindle assembly. These techniques can be adapted to a wide variety of species in order to measure microtubule dynamics and spindle assembly kinetics when genetic tools are not available or in parallel to the development of such techniques in non-model organisms.