Sex-Specific Expression of Alternative Transcripts in \u3ci\u3eDrosophila\u3c/i\u3e

Background: Many genes produce multiple transcripts due to alternative splicing or utilization of alternative transcription initiation/termination sites. This \u27transcriptome expansion\u27 is thought to increase phenotypic complexity by allowing a single locus to produce several functionally distinct proteins. However, sex, genetic and developmental variation in the representation of alternative transcripts has never been examined systematically. Here, we describe a genome-wide analysis of sex-specific expression of alternative transcripts in Drosophila melanogaster. Results: We compared transcript profiles in males and females from eight Drosophila lines (OregonR and 2b, and 6 RIL) using a newly designed 60-mer oligonucleotide microarray that allows us to distinguish a large proportion of alternative transcripts. The new microarray incorporates 7,207 oligonucleotides, satisfying stringent binding and specificity criteria that target both the common and the unique regions of 2,768 multi-transcript genes, as well as 12,912 oligonucleotides that target genes with a single known transcript. We estimate that up to 22% of genes that produce multiple transcripts show a sex-specific bias in the representation of alternative transcripts. Sexual dimorphism in overall transcript abundance was evident for 53% of genes. The X chromosome contains a significantly higher proportion of genes with female-biased transcription than the autosomes. However, genes on the X chromosome are no more likely to have a sexual bias in alternative transcript representation than autosomal genes. Conclusion: Widespread sex-specific expression of alternative transcripts in Drosophila suggests that a new level of sexual dimorphism at the molecular level exists

Sex-specific expression of alternative transcripts in Drosophila

BACKGROUND: Many genes produce multiple transcripts due to alternative splicing or utilization of alternative transcription initiation/termination sites. This 'transcriptome expansion' is thought to increase phenotypic complexity by allowing a single locus to produce several functionally distinct proteins. However, sex, genetic and developmental variation in the representation of alternative transcripts has never been examined systematically. Here, we describe a genome-wide analysis of sex-specific expression of alternative transcripts in Drosophila melanogaster. RESULTS: We compared transcript profiles in males and females from eight Drosophila lines (OregonR and 2b, and 6 RIL) using a newly designed 60-mer oligonucleotide microarray that allows us to distinguish a large proportion of alternative transcripts. The new microarray incorporates 7,207 oligonucleotides, satisfying stringent binding and specificity criteria that target both the common and the unique regions of 2,768 multi-transcript genes, as well as 12,912 oligonucleotides that target genes with a single known transcript. We estimate that up to 22% of genes that produce multiple transcripts show a sex-specific bias in the representation of alternative transcripts. Sexual dimorphism in overall transcript abundance was evident for 53% of genes. The X chromosome contains a significantly higher proportion of genes with female-biased transcription than the autosomes. However, genes on the X chromosome are no more likely to have a sexual bias in alternative transcript representation than autosomal genes. CONCLUSION: Widespread sex-specific expression of alternative transcripts in Drosophila suggests that a new level of sexual dimorphism at the molecular level exists

A Search for Parent-of-Origin Effects on Honey Bee Gene Expression

Parent-specific gene expression (PSGE) is little known outside of mammals and plants. PSGE occurs when the expression level of a gene depends on whether an allele was inherited from the mother or the father. Kin selection theory predicts that there should be extensive PSGE in social insects because social insect parents can gain inclusive fitness benefits by silencing parental alleles in female offspring. We searched for evidence of PSGE in honey bees using transcriptomes from reciprocal crosses between European and Africanized strains. We found 46 transcripts with significant parent-of-origin effects on gene expression, many of which overexpressed the maternal allele. Interestingly, we also found a large proportion of genes showing a bias toward maternal alleles in only one of the reciprocal crosses. These results indicate that PSGE may occur in social insects. The nonreciprocal effects could be largely driven by hybrid incompatibility between these strains. Future work will help to determine if these are indeed parent-of-origin effects that can modulate inclusive fitness benefits

Integration of hybridization-based markers (overgos) into physical maps for comparative and evolutionary explorations in the genus Oryza and in Sorghum

BACKGROUND: With the completion of the genome sequence for rice (Oryza sativa L.), the focus of rice genomics research has shifted to the comparison of the rice genome with genomes of other species for gene cloning, breeding, and evolutionary studies. The genus Oryza includes 23 species that shared a common ancestor 8–10 million years ago making this an ideal model for investigations into the processes underlying domestication, as many of the Oryza species are still undergoing domestication. This study integrates high-throughput, hybridization-based markers with BAC end sequence and fingerprint data to construct physical maps of rice chromosome 1 orthologues in two wild Oryza species. Similar studies were undertaken in Sorghum bicolor, a species which diverged from cultivated rice 40–50 million years ago. RESULTS: Overgo markers, in conjunction with fingerprint and BAC end sequence data, were used to build sequence-ready BAC contigs for two wild Oryza species. The markers drove contig merges to construct physical maps syntenic to rice chromosome 1 in the wild species and provided evidence for at least one rearrangement on chromosome 1 of the O. sativa versus Oryza officinalis comparative map. When rice overgos were aligned to available S. bicolor sequence, 29% of the overgos aligned with three or fewer mismatches; of these, 41% gave positive hybridization signals. Overgo hybridization patterns supported colinearity of loci in regions of sorghum chromosome 3 and rice chromosome 1 and suggested that a possible genomic inversion occurred in this syntenic region in one of the two genomes after the divergence of S. bicolor and O. sativa. CONCLUSION: The results of this study emphasize the importance of identifying conserved sequences in the reference sequence when designing overgo probes in order for those probes to hybridize successfully in distantly related species. As interspecific markers, overgos can be used successfully to construct physical maps in species which diverged less than 8 million years ago, and can be used in a more limited fashion to examine colinearity among species which diverged as much as 40 million years ago. Additionally, overgos are able to provide evidence of genomic rearrangements in comparative physical mapping studies

The inference of gray whale (Eschrichtius robustus) historical population attributes from whole-genome sequences

Commercial whaling caused extensive demographic declines in many great whale species, including gray whales that were extirpated from the Atlantic Ocean and dramatically reduced in the Pacific Ocean. The Eastern Pacific gray whale has recovered since the 1982 ban on commercial whaling, but the Western Pacific gray whale-once considered possibly extinct-consists of only about 200 individuals and is considered critically endangered by some international authorities. Herein, we use whole-genome sequencing to investigate the demographic history of gray whales from the Pacific and use environmental niche modelling to make predictions about future gene flow.Our sequencing efforts and habitat niche modelling indicate that: i) western gray whale effective population sizes have declined since the last glacial maximum; ii) contemporary gray whale genomes, both eastern and western, harbor less autosomal nucleotide diversity than most other marine mammals and megafauna; iii) the extent of inbreeding, as measured by autozygosity, is greater in the Western Pacific than in the Eastern Pacific populations; and iv) future climate change is expected to open new migratory routes for gray whales.Our results indicate that gray whale genomes contain low nucleotide diversity and have been subject to both historical and recent inbreeding. Population sizes over the last million years likely peaked about 25,000 years before present and have declined since then. Our niche modelling suggests that novel migratory routes may develop within the next century and if so this could help retain overall genetic diversity, which is essential for adaption and successful recovery in light of global environmental change and past exploitation

Data from: New insights into the phylogenetics and population structure of the prairie falcon (Falco mexicanus)

Background: Management requires a robust understanding of between- and within-species genetic variability, however such data are still lacking in many species. For example, although multiple population genetics studies of the peregrine falcon (Falco peregrinus) have been conducted, no similar studies have been done of the closely-related prairie falcon (F. mexicanus) and it is unclear how much genetic variation and population structure exists across the species’ range. Furthermore, the phylogenetic relationship of F. mexicanus relative to other falcon species is contested. We utilized a genomics approach (i.e., genome sequencing and assembly followed by single nucleotide polymorphism genotyping) to rapidly address these gaps in knowledge. Results: We sequenced the genome of a single female prairie falcon and generated a 1.17 Gb (gigabases) draft genome assembly. We generated maximum likelihood phylogenetic trees using complete mitochondrial genomes as well as nuclear protein-coding genes. This process provided evidence that F. mexicanus is an outgroup to the clade that includes the peregrine falcon and members of the subgenus Hierofalco. We annotated > 16,000 genes and almost 600,000 high-quality single nucleotide polymorphisms (SNPs) in the nuclear genome, providing the raw material for a SNP assay design featuring > 140 gene-associated markers and a molecular-sexing marker. We subsequently genotyped ~ 100 individuals from California (including the San Francisco East Bay Area, Pinnacles National Park and the Mojave Desert) and Idaho (Snake River Birds of Prey National Conservation Area). We tested for population structure and found evidence that individuals sampled in California and Idaho represent a single panmictic population. Conclusions: Our study illustrates how genomic resources can rapidly shed light on genetic variability in understudied species and resolve phylogenetic relationships. Furthermore, we found evidence of a single, randomly mating population of prairie falcons across our sampling locations. Prairie falcons are highly mobile and relatively rare long-distance dispersal events may promote gene flow throughout the range. As such, California’s prairie falcons might be managed as a single population, indicating that management actions undertaken to benefit the species at the local level have the potential to influence the species as a whole

Falco mexicanus genome assembly scaffolds

The Falco mexicanus draft genome was generated from paired-end, mate-paired and long reads using ABySS vs. 1.5.2 and SPAdes 3.1.1. This file includes the scaffolds greater than 10 kb in length

Additional file 1: of The inference of gray whale (Eschrichtius robustus) historical population attributes from whole-genome sequences

Table S1. Information on raw reads filtering statistics. Paired-end libraries were sequenced on an Illumina HiSeq 2500. Table S2. Environmental variables used in AQUAMAPS to generate maps of suitable habitat for gray whales during the Holocene. Table S3. The D–test statistic evaluates the number (n) of ABBA and BABA sites (D = (nABBA - nBABA) / (nABBA + nBABA)) and D < 0 means that P1 is more closely related to P3 than to P2, whereas D > 0 indicates that P2 is more closely related to P3 than P1. The significance of the D test was evaluated with a Z-score, where |Z-scores| > 3 was used as the critical value for a significant test. Figure S1. Inferred effective population sizes (Ne) over time. Estimates are averages based on 11 autosomal scaffolds larger than 30 Mb. A substitution rate of a) 10 × 10− 10 bp− 1 year− 1 and b) 1.5 × 10− 10 bp− 1 year− 1 were used. (DOCX 446 kb

Falco mexicanus mitochondrial genome assembly

The F. mexicanus mitochondrial genome was assembled using MITObim vs. 1.6. The genome assembly is 17,117 bp in length and characterized by 13 protein-coding genes, two ribosomal subunit genes (rRNA), 22 transfer RNA genes and a control region

Falco mexicanus genome annotation - protein sequences

16,320 genes were annotated in the Falco mexicanus genome, this is the .fasta file with predicted protein sequences