Methyl CpG Binding Domain Ultra‐Sequencing: a novel method for identifying inter‐individual and cell‐type‐specific variation in DNA methylation

Experience‐dependent changes in DNA methylation can exert profound effects on neuronal function and behaviour. A single learning event can induce a variety of DNA modifications within the neuronal genome, some of which may be common to all individuals experiencing the event, whereas others may occur in a subset of individuals. Variations in experience‐induced DNA methylation may subsequently confer increased vulnerability or resilience to the development of neuropsychiatric disorders. However, the detection of experience‐dependent changes in DNA methylation in the brain has been hindered by the interrogation of heterogeneous cell populations, regional differences in epigenetic states and the use of pooled tissue obtained from multiple individuals. Methyl CpG Binding Domain Ultra‐Sequencing (MBD Ultra‐Seq) overcomes current limitations on genome‐wide epigenetic profiling by incorporating fluorescence‐activated cell sorting and sample‐specific barcoding to examine cell‐type‐specific CpG methylation in discrete brain regions of individuals. We demonstrate the value of this method by characterizing differences in 5‐methylcytosine (5mC) in neurons and non‐neurons of the ventromedial prefrontal cortex of individual adult C57BL/6 mice, using as little as 50 ng of genomic DNA per sample. We find that the neuronal methylome is characterized by greater CpG methylation as well as the enrichment of 5mC within intergenic loci. In conclusion, MBD Ultra‐Seq is a robust method for detecting DNA methylation in neurons derived from discrete brain regions of individual animals. This protocol will facilitate the detection of experience‐dependent changes in DNA methylation in a variety of behavioural paradigms and help identify aberrant experience‐induced DNA methylation that may underlie risk and resiliency to neuropsychiatric disease.     

NGS‐based targeted resequencing identified rare subtypes of albinism: Providing accurate molecular diagnosis for Japanese patients with albinism

Albinism, which is commonly inherited as an autosomal recessive trait, is characterized by a reduction or absence of melanin in the eyes, skin, and hair. To date, more than 20 causal genes for albinism have been identified; thus, the accurate diagnosis of albinism requires next‐generation sequencing (NGS). In this study, we analyzed 46 patients who tested negative for oculocutaneous albinism (OCA)1–4 and Hermansky‐Pudlak syndrome (HPS)1 based on conventional analysis, in addition to 28 new Japanese patients, using NGS‐based targeted resequencing. We identified a genetic background for albinism in 18 of the 46 patients (39%), who were previously tested negative according to the conventional analysis. In addition, we unveiled a genetic predisposition toward albinism in 23 of the 28 new patients (82%). We identified six patients with rare subtypes of albinism, including HPS3, HPS4, and HPS6, and found 12 novel pathological mutations in albinism‐related genes. Furthermore, most patients who were not diagnosed with albinism by the NGS analysis showed mild manifestations of albinism without apparent eye symptoms and harbored only one heterozygous mutation, occasionally in combination with skin‐color associated gene variants.     

Targeted multiplex next‐generation sequencing: advances in techniques of mitochondrial and nuclear DNA sequencing for population genomics

Next‐generation sequencing (NGS) is emerging as an efficient and cost‐effective tool in population genomic analyses of nonmodel organisms, allowing simultaneous resequencing of many regions of multi‐genomic DNA from multiplexed samples. Here, we detail our synthesis of protocols for targeted resequencing of mitochondrial and nuclear loci by generating indexed genomic libraries for multiplexing up to 100 individuals in a single sequencing pool, and then enriching the pooled library using custom DNA capture arrays. Our use of DNA sequence from one species to capture and enrich the sequencing libraries of another species (i.e. cross‐species DNA capture) indicates that efficient enrichment occurs when sequences are up to about 12% divergent, allowing us to take advantage of genomic information in one species to sequence orthologous regions in related species. In addition to a complete mitochondrial genome on each array, we have included between 43 and 118 nuclear loci for low‐coverage sequencing of between 18 kb and 87 kb of DNA sequence per individual for single nucleotide polymorphisms discovery from 50 to 100 individuals in a single sequencing lane. Using this method, we have generated a total of over 500 whole mitochondrial genomes from seven cetacean species and green sea turtles. The greater variation detected in mitogenomes relative to short mtDNA sequences is helping to resolve genetic structure ranging from geographic to species‐level differences. These NGS and analysis techniques have allowed for simultaneous population genomic studies of mtDNA and nDNA with greater genomic coverage and phylogeographic resolution than has previously been possible in marine mammals and turtles.     

Multi‐Omics Reveals Impact of Cysteine Feed Concentration and Resulting Redox Imbalance on Cellular Energy Metabolism and Specific Productivity in CHO Cell Bioprocessing

Chinese hamster ovary (CHO) cells are currently the primary host cell lines used in biotherapeutic manufacturing of monoclonal antibodies (mAbs) and other biopharmaceuticals. Cellular energy metabolism and endoplasmic reticulum (ER) stress are known to greatly impact cell growth, viability, and specific productivity of a biotherapeutic; but the molecular mechanisms are not fully understood. The authors previously employed multi‐omics profiling to investigate the impact of a reduction in cysteine (Cys) feed concentration in a fed‐batch process and found that disruption of the redox balance led to a substantial decline in cell viability and titer. Here, the multi‐omics findings are expanded, and the impact redox imbalance has on ER stress, mitochondrial homeostasis, and lipid metabolism is explored. The reduced Cys feed activates the amino acid response (AAR), increases mitochondrial stress, and initiates gluconeogenesis. Multi‐omics analysis reveals that together, ER stress and AAR signaling shift the cellular energy metabolism to rely primarily on anaplerotic reactions, consuming amino acids and producing lactate, to maintain energy generation. Furthermore, the pathways are demonstrated in which this shift in metabolism leads to a substantial decline in specific productivity and altered mAb glycosylation. Through this work, meaningful bioprocess markers and targets for genetic engineering are identified.     

Development and characterization of polymorphic microsatellite markers for Neolitsea sericea using Illumina paired‐end draft sequencing data

Simple sequence repeat (SSR) markers were developed and characterized for Neolitsea sericea (Bl.) Koidz. (Lauraceae). Out of 196 designed primer pairs, a total of 144 pairs showed amplification, of which 44 had clear and stable chromatograms. Polymorphism of these 44 loci was tested using 32 individuals sampled from a single population of N. sericea. The number of alleles and the polymorphism information content varied from 3 to 12 and 0.271 to 0.853, respectively. A significant departure from the Hardy‐Weinberg equilibrium was observed in one of the 44 loci. These SSR markers are useful for population genetic studies and parentage analysis in N. sericea, which is one of the most common evergreen species in coastal Pinus thunbergii forests in central‐western Japan.     

Forest understory plant and soil microbial response to an experimentally induced drought and heat‐pulse event: the importance of maintaining the continuum

Drought duration and intensity are expected to increase with global climate change. How changes in water availability and temperature affect the combined plant–soil–microorganism response remains uncertain. We excavated soil monoliths from a beech (Fagus sylvatica L.) forest, thus keeping the understory plant–microbe communities intact, imposed an extreme climate event, consisting of drought and/or a single heat‐pulse event, and followed microbial community dynamics over a time period of 28 days. During the treatment, we labeled the canopy with ¹³CO₂ with the goal of (i) determining the strength of plant–microbe carbon linkages under control, drought, heat and heat–drought treatments and (ii) characterizing microbial groups that are tightly linked to the plant–soil carbon continuum based on ¹³C‐labeled PLFAs. Additionally, we used 16S rRNA sequencing of bacteria from the Ah horizon to determine the short‐term changes in the active microbial community. The treatments did not sever within‐plant transport over the experiment, and carbon sinks belowground were still active. Based on the relative distribution of labeled carbon to roots and microbial PLFAs, we determined that soil microbes appear to have a stronger carbon sink strength during environmental stress. High‐throughput sequencing of the 16S rRNA revealed multiple trajectories in microbial community shifts within the different treatments. Heat in combination with drought had a clear negative effect on microbial diversity and resulted in a distinct shift in the microbial community structure that also corresponded to the lowest level of label found in the PLFAs. Hence, the strongest changes in microbial abundances occurred in the heat–drought treatment where plants were most severely affected. Our study suggests that many of the shifts in the microbial communities that we might expect from extreme environmental stress will result from the plant–soil–microbial dynamics rather than from direct effects of drought and heat on soil microbes alone.     

The easy road to genome‐wide medium density SNP screening in a non‐model species: development and application of a 10 K SNP‐chip for the house sparrow (Passer domesticus)

With the advent of next generation sequencing, new avenues have opened to study genomics in wild populations of non‐model species. Here, we describe a successful approach to a genome‐wide medium density Single Nucleotide Polymorphism (SNP) panel in a non‐model species, the house sparrow (Passer domesticus), through the development of a 10 K Illumina iSelect HD BeadChip. Genomic DNA and cDNA derived from six individuals were sequenced on a 454 GS FLX system and generated a total of 1.2 million sequences, in which SNPs were detected. As no reference genome exists for the house sparrow, we used the zebra finch (Taeniopygia guttata) reference genome to determine the most likely position of each SNP. The 10 000 SNPs on the SNP‐chip were selected to be distributed evenly across 31 chromosomes, giving on average one SNP per 100 000 bp. The SNP‐chip was screened across 1968 individual house sparrows from four island populations. Of the original 10 000 SNPs, 7413 were found to be variable, and 99% of these SNPs were successfully called in at least 93% of all individuals. We used the SNP‐chip to demonstrate the ability of such genome‐wide marker data to detect population sub‐division, and compared these results to similar analyses using microsatellites. The SNP‐chip will be used to map Quantitative Trait Loci (QTL) for fitness‐related phenotypic traits in natural populations.     

QDD version 3.1: a user‐friendly computer program for microsatellite selection and primer design revisited: experimental validation of variables determining genotyping success rate

Microsatellite marker development has been greatly simplified by the use of high‐throughput sequencing followed by in silico microsatellite detection and primer design. However, the selection of markers designed by the existing pipelines depends either on arbitrary criteria, or older studies on PCR success. Based on wet laboratory experiments, we have identified the following factors that are most likely to influence genotyping success rate: alignment score between the primers and the amplicon; the distance between primers and microsatellites; the length of the PCR product; target region complexity and the number of reads underlying the sequence. The QDD pipeline has been modified to include these most pertinent factors in the output to help the selection of markers. Furthermore, new features are also included in the present version: (i) not only raw sequencing reads are accepted as input, but also contigs, allowing the analysis of assembled high‐coverage data; (ii) input data can be both in fasta and fastq format to facilitate the use of Illumina and IonTorrent reads; (iii) A comparison to known transposable elements allows their detection; (iv) A contamination check can be carried out by BLASTing potential markers against the nucleotide (nt) database of NCBI; (v) QDD3 is now also available imbedded into a virtual machine making installation easier and operating system independent. It can be used both on command‐line version as well as integrated into a Galaxy server, providing a user‐friendly interface, as well as the possibility to utilize a large variety of NGS tools.     

Host‐associated divergence in sympatric host races of the leaf beetle Lochmaea capreae: implications for local adaptation and reproductive isolation

Ecological specialization is widely recognized as a major determinant of the emergence and maintenance of biodiversity. We studied two critical facets of specialization – local adaptation and habitat choice – in the host races of the leaf beetle Lochmaea capreae on willow and birch. Our results revealed that there is asymmetric disruptive selection for host use traits, and host races achieved different adaptive sets of life history traits through association with their host plant. Beetles from each host race exhibited food and oviposition preference for their own host plant. Reciprocal transplant displayed significant variation in host acceptance and performance: all families from the willow race rejected the alternative host plant before initiation of feeding and all died on this host plant. By contrast, all families from the birch race accepted willow for feeding, but they consumed less and performed less well. Intriguingly, families that performed well on birch also performed well on willow, suggesting positive genetic correlation rather than genetic trade‐offs. Our results suggest that the major proximal determinant of host specialization in the willow race is the behavioural acceptance of a plant rather than the toxicity of the food resource. However, in the birch race a combination of behavioural host acceptance and performance may play a role in specialization. Our study sheds light on the mechanisms by which divergent host adaptation might influence the evolution of reproductive isolation between herbivorous populations.