Multi-tissue RNA-Seq Analysis and Long-read-based Genome Assembly Reveal Complex Sex-specific Gene Regulation and Molecular Evolution in the Manila Clam.

The molecular factors and gene regulation involved in sex determination and gonad differentiation in bivalve molluscs are unknown. It has been suggested that doubly uniparental inheritance (DUI) of mitochondria may be involved in these processes in species such as the ubiquitous and commercially relevant Manila clam, Ruditapes philippinarum. We present the first long-read-based de novo genome assembly of a Manila clam, and a RNA-Seq multi-tissue analysis of 15 females and 15 males. The highly contiguous genome assembly was used as reference to investigate gene expression, alternative splicing, sequence evolution, tissue-specific co-expression networks, and sexual contrasting SNPs. Differential expression (DE) and differential splicing (DS) analyses revealed sex-specific transcriptional regulation in gonads, but not in somatic tissues. Co-expression networks revealed complex gene regulation in gonads, and genes in gonad-associated modules showed high tissue specificity. However, male gonad-associated modules showed contrasting patterns of sequence evolution and tissue specificity. One gene set was related to the structural organization of male gametes and presented slow sequence evolution but high pleiotropy, whereas another gene set was enriched in reproduction-related processes and characterized by fast sequence evolution and tissue specificity. Sexual contrasting SNPs were found in genes overrepresented in mitochondrial-related functions, providing new candidates for investigating the relationship between mitochondria and sex in DUI species. Together, these results increase our understanding of the role of DE, DS, and sequence evolution of sex-specific genes in an understudied taxon. We also provide resourceful genomic data for studies regarding sex diagnosis and breeding in bivalves.

Revealing the Draft Genome Sequence of Bradyrhizobium sp. Strain USDA 3458, an Effective Symbiotic Diazotroph Isolated from Cowpea (Vigna unguiculata) Genotype IT82E-16.

Pairing plants with plant growth-promoting bacteria is critical to the future of agriculture. Bradyrhizobium sp. strain USDA 3458 isolated from Vigna unguiculata (cowpea) paired with cowpea genotype IT82E-16 represents a novel combination in arid regions. Here, we report the draft genome sequence of strain USDA 3458.

Elucidation of the Genome of Bradyrhizobium sp. Strain USDA 3456, a Historic Agricultural Diazotroph from Cowpea (Vigna unguiculata).

Bradyrhizobium sp. strain USDA 3456 is a historic strain from the United States Department of Agriculture (USDA) Agricultural Research Service (ARS) National Rhizobium Germplasm Collection isolated from Vigna unguiculata (cowpea) in 1966. Strain USDA 3456 has been utilized in global agricultural applications, including improving soil nitrogen fertility. The draft genome sequence here provides a genetic reference of a novel diazotroph.

Association mapping reveals novel serpentine adaptation gene clusters in a population of symbiotic Mesorhizobium.

The genetic variants that underlie microbial environmental adaptation are key components of models of microbial diversification. Characterizing adaptive variants and the pangenomic context in which they evolve remains a frontier in understanding how microbial diversity is generated. The genomics of rhizobium adaptation to contrasting soil environments is ecologically and agriculturally important because these bacteria are responsible for half of all current biologically fixed nitrogen, yet they live the majority of their lives in soil. Our study uses whole-genome sequencing to describe the pan-genome of a focal clade of wild mesorhizobia that show contrasting levels of nickel adaptation despite high relatedness (99.8% identity at 16S). We observe ecotypic specialization within an otherwise genomically cohesive population, rather than finding distinct specialized bacterial lineages in contrasting soil types. This finding supports recent reports that heterogeneous environments impose selection that maintains differentiation only at a small fraction of the genome. Our work further uses a genome-wide association study to propose candidate genes for nickel adaptation. Several candidates show homology to genetic systems involved in nickel tolerance and one cluster of candidates correlates perfectly with soil origin, which validates our approach of ascribing genomic variation to adaptive divergence.

A cost-effective method for high-throughput construction of illumina sequencing libraries.

Despite the plummeting cost of next-generation sequencing, the preparation of sequencing libraries using commercially available kits still remains expensive. The cost can be prohibitive for large-scale comparative or experimental studies, where hundreds to thousands of samples need to be analyzed. The increasing use of multiplexing dozens to hundreds of samples underscores the urgent need to develop a cost-effective and time-efficient high-throughput method for library preparation. By optimizing and scaling down the steps in library construction and using commonly available reagents, the protocol described here allows for the preparation of DNA libraries in a 96-well format using no specialized equipment and at a substantial savings in both reagent cost and personnel hours. Utilizing this optimized high-throughput format results in a 10-fold cost reduction, compared to commercially available kits, making per library or pooled sample costs ∼$12.60-14.90 for individually prepared libraries and ∼$8.60-10.60 for pooled libraries with individual barcodes; both techniques allow for up to 144 samples to be pooled on a single lane with the barcodes tested herein.

Widespread selection across coding and noncoding DNA in the pea aphid genome.

Genome-wide patterns of diversity and selection are critical measures for understanding how evolution has shaped the genome. Yet, these population genomic estimates are available for only a limited number of model organisms. Here we focus on the population genomics of the pea aphid (Acyrthosiphon pisum). The pea aphid is an emerging model system that exhibits a range of intriguing biological traits not present in classic model systems. We performed low-coverage genome resequencing of 21 clonal pea aphid lines collected from alfalfa host plants in North America to characterize genome-wide patterns of diversity and selection. We observed an excess of low-frequency polymorphisms throughout coding and noncoding DNA, which we suggest is the result of a founding event and subsequent population expansion in North America. Most gene regions showed lower levels of Tajima's D than synonymous sites, suggesting that the majority of the genome is not evolving neutrally but rather exhibits significant constraint. Furthermore, we used the pea aphid's unique manner of X-chromosome inheritance to assign genomic scaffolds to either autosomes or the X chromosome. Comparing autosomal vs. X-linked sequence variation, we discovered that autosomal genes show an excess of low frequency variants indicating that purifying selection acts more efficiently on the X chromosome. Overall, our results provide a critical first step in characterizing the genetic diversity and evolutionary pressures on an aphid genome.

Somatic sex-specific transcriptome differences in Drosophila revealed by whole transcriptome sequencing.

BACKGROUND: Understanding animal development and physiology at a molecular-biological level has been advanced by the ability to determine at high resolution the repertoire of mRNA molecules by whole transcriptome resequencing. This includes the ability to detect and quantify rare abundance transcripts and isoform-specific mRNA variants produced from a gene.The sex hierarchy consists of a pre-mRNA splicing cascade that directs the production of sex-specific transcription factors that specify nearly all sexual dimorphism. We have used deep RNA sequencing to gain insight into how the Drosophila sex hierarchy generates somatic sex differences, by examining gene and transcript isoform expression differences between the sexes in adult head tissues. RESULTS: Here we find 1,381 genes that differ in overall expression levels and 1,370 isoform-specific transcripts that differ between males and females. Additionally, we find 512 genes not regulated downstream of transformer that are significantly more highly expressed in males than females. These 512 genes are enriched on the × chromosome and reside adjacent to dosage compensation complex entry sites, which taken together suggests that their residence on the × chromosome might be sufficient to confer male-biased expression. There are no transcription unit structural features, from a set of features, that are robustly significantly different in the genes with significant sex differences in the ratio of isoform-specific transcripts, as compared to random isoform-specific transcripts, suggesting that there is no single molecular mechanism that generates isoform-specific transcript differences between the sexes, even though the sex hierarchy is known to include three pre-mRNA splicing factors. CONCLUSIONS: We identify thousands of genes that show sex-specific differences in overall gene expression levels, and identify hundreds of additional genes that have differences in the abundance of isoform-specific transcripts. No transcription unit structural feature was robustly enriched in the sex-differentially expressed transcript isoforms. Additionally, we found that many genes with male-biased expression were enriched on the × chromosome and reside adjacent to dosage compensation entry sites, suggesting that differences in sex chromosome composition contributes to dimorphism in gene expression. Taken together, this study provides new insight into the molecular underpinnings of sexual differentiation.

Genome-wide analysis of a long-term evolution experiment with Drosophila.

Experimental evolution systems allow the genomic study of adaptation, and so far this has been done primarily in asexual systems with small genomes, such as bacteria and yeast. Here we present whole-genome resequencing data from Drosophila melanogaster populations that have experienced over 600 generations of laboratory selection for accelerated development. Flies in these selected populations develop from egg to adult ∼20% faster than flies of ancestral control populations, and have evolved a number of other correlated phenotypes. On the basis of 688,520 intermediate-frequency, high-quality single nucleotide polymorphisms, we identify several dozen genomic regions that show strong allele frequency differentiation between a pooled sample of five replicate populations selected for accelerated development and pooled controls. On the basis of resequencing data from a single replicate population with accelerated development, as well as single nucleotide polymorphism data from individual flies from each replicate population, we infer little allele frequency differentiation between replicate populations within a selection treatment. Signatures of selection are qualitatively different than what has been observed in asexual species; in our sexual populations, adaptation is not associated with 'classic' sweeps whereby newly arising, unconditionally advantageous mutations become fixed. More parsimonious explanations include 'incomplete' sweep models, in which mutations have not had enough time to fix, and 'soft' sweep models, in which selection acts on pre-existing, common genetic variants. We conclude that, at least for life history characters such as development time, unconditionally advantageous alleles rarely arise, are associated with small net fitness gains or cannot fix because selection coefficients change over time.

High nucleotide divergence in developmental regulatory genes contrasts with the structural elements of olfactory pathways in caenorhabditis.

Almost all organismal function is controlled by pathways composed of interacting genetic components. The relationship between pathway structure and the evolution of individual pathway components is not completely understood. For the nematode Caenorhabditis elegans, chemosensory pathways regulate critical aspects of an individual's life history and development. To help understand how olfaction evolves in Caenorhabditis and to examine patterns of gene evolution within transduction pathways in general, we analyzed nucleotide variation within and between species across two well-characterized olfactory pathways, including regulatory genes controlling the fate of the cells in which the pathways are expressed. In agreement with previous studies, we found much higher levels of polymorphism within C. remanei than within the related species C. elegans and C. briggsae. There are significant differences in the rates of nucleotide evolution for genes across the two pathways but no particular association between evolutionary rate and gene position, suggesting that the evolution of functional pathways must be considered within the context of broader gene network structure. However, developmental regulatory genes show both higher levels of divergence and polymorphism than the structural genes of the pathway. These results show that, contrary to the emerging paradigm in the evolution of development, important structural changes can accumulate in transcription factors.

Rapid and cost-effective polymorphism identification and genotyping using restriction site associated DNA (RAD) markers.

Restriction site associated DNA (RAD) tags are a genome-wide representation of every site of a particular restriction enzyme by short DNA tags. Most organisms segregate large numbers of DNA sequence polymorphisms that disrupt restriction sites, which allows RAD tags to serve as genetic markers spread at a high density throughout the genome. Here, we demonstrate the applicability of RAD markers for both individual and bulk-segregant genotyping. First, we show that these markers can be identified and typed on pre-existing microarray formats. Second, we present a method that uses RAD marker DNA to rapidly produce a low-cost microarray genotyping resource that can be used to efficiently identify and type thousands of RAD markers. We demonstrate the utility of the former approach by using a tiling path array for the fruit fly to map a recombination breakpoint, and the latter approach by creating and using an enriched RAD marker array for the threespine stickleback. The high number of RAD markers enabled localization of a previously identified region, as well as a second region also associated with the lateral plate phenotype. Taken together, our results demonstrate that RAD markers, and the method to develop a RAD marker microarray resource, allow high-throughput, high-resolution genotyping in both model and nonmodel systems.