Search

According to archaeological records, chickpea (Cicer arietinum) was first domesticated in the Fertile Crescent about 10,000 years BP. Its subsequent diversification in Middle East, South Asia, Ethiopia, and the Western Mediterranean, however, remains obscure and cannot be resolved using only archeological and historical evidence. Moreover, chickpea has two market types: "desi" and "kabuli," for which the geographic origin is a matter of debate. To decipher chickpea history, we took the genetic data from 421 chickpea landraces unaffected by the green revolution and tested complex historical hypotheses of chickpea migration and admixture on two hierarchical spatial levels: within and between major regions of cultivation. For chickpea migration within regions, we developed popdisp, a Bayesian model of population dispersal from a regional representative center toward the sampling sites that considers geographical proximities between sites. This method confirmed that chickpea spreads within each geographical region along optimal geographical routes rather than by simple diffusion and estimated representative allele frequencies for each region. For chickpea migration between regions, we developed another model, migadmi, that takes allele frequencies of populations and evaluates multiple and nested admixture events. Applying this model to desi populations, we found both Indian and Middle Eastern traces in Ethiopian chickpea, suggesting the presence of a seaway from South Asia to Ethiopia. As for the origin of kabuli chickpeas, we found significant evidence for its origin from Turkey rather than Central Asia.

Enteric methane emissions are the single largest source of direct greenhouse gas emissions (GHG) in beef and dairy value chains and a substantial contributor to anthropogenic methane emissions globally. In late 2019, the World Wildlife Fund (WWF), the Advanced Research Projects Agency-Energy (ARPA-E) and the Foundation for Food and Agriculture Research (FFAR) convened approximately 50 stakeholders representing research and production of seaweeds, animal feeds, dairy cattle, and beef and dairy foods to discuss challenges and opportunities associated with the use of seaweed-based ingredients to reduce enteric methane emissions. This Perspective article describes the considerations identified by the workshop participants and suggests next steps for the further development and evaluation of seaweed-based feed ingredients as enteric methane mitigants. Although numerous compounds derived from sources other than seaweed have been identified as having enteric methane mitigation potential, these mitigants are outside the scope of this article.

Enteric methane emissions are the single largest source of direct greenhouse gas emissions (GHG) in beef and dairy value chains and a substantial contributor to anthropogenic methane emissions globally. In late 2019, the World Wildlife Fund (WWF), the Advanced Research Projects Agency-Energy (ARPA-E) and the Foundation for Food and Agriculture Research (FFAR) convened approximately 50 stakeholders representing research and production of seaweeds, animal feeds, dairy cattle, and beef and dairy foods to discuss challenges and opportunities associated with the use of seaweed-based ingredients to reduce enteric methane emissions. This Perspective article describes the considerations identified by the workshop participants and suggests next steps for the further development and evaluation of seaweed-based feed ingredients as enteric methane mitigants. Although numerous compounds derived from sources other than seaweed have been identified as having enteric methane mitigation potential, these mitigants are outside the scope of this article.

A defining challenge of the 21st century is meeting the nutritional demands of the growing human population, under a scenario of limited land and water resources and under the specter of climate change. The Vavilov seed bank contains numerous landraces collected nearly a hundred years ago, and thus may contain 'genetic gems' with the potential to enhance modern breeding efforts. Here, we analyze 407 landraces, sampled from major historic centers of chickpea cultivation and secondary diversification. Genome-Wide Association Studies (GWAS) conducted on both phenotypic traits and bioclimatic variables at landraces sampling sites as extended phenotypes resulted in 84 GWAS hits associated to various regions. The novel haploblock-based test identified haploblocks enriched for single nucleotide polymorphisms (SNPs) associated with phenotypes and bioclimatic variables. Subsequent bi-clustering of traits sharing enriched haploblocks underscored both non-random distribution of SNPs among several haploblocks and their association with multiple traits. We hypothesize that these clusters of pleiotropic SNPs represent co-adapted genetic complexes to a range of environmental conditions that chickpea experienced during domestication and subsequent geographic radiation. Linking genetic variation to phenotypic data and a wealth of historic information preserved in historic seed banks are the keys for genome-based and environment-informed breeding intensification.

Gene regulatory networks control the complex programs that drive development. Deciphering the connections between transcription factors (TFs) and target genes is challenging, in part because TFs bind to thousands of places in the genome but control expression through a subset of these binding events. We hypothesize that we can combine natural variation of expression levels and predictions of TF binding sites to identify TF targets. We gather RNA-seq data from 71 genetically distinct F1 Drosophila melanogaster embryos and calculate the correlations between TF and potential target genes' expression levels, which we call "regulatory strength." To separate direct and indirect TF targets, we hypothesize that direct TF targets will have a preponderance of binding sites in their upstream regions. Using 14 TFs active during embryogenesis, we find that 12 TFs showed a significant correlation between their binding strength and regulatory strength on downstream targets, and 10 TFs showed a significant correlation between the number of binding sites and the regulatory effect on target genes. The general roles, e.g. bicoid's role as an activator, and the particular interactions we observed between our TFs, e.g. twist's role as a repressor of sloppy paired and odd paired, generally coincide with the literature.

Noncoding DNA sequences, which play various roles in gene expression and regulation, are under evolutionary pressure. Gene regulation requires specific protein-DNA binding events, and our previous studies showed that both DNA sequence and shape readout are employed by transcription factors (TFs) to achieve DNA binding specificity. By investigating the shape-disrupting properties of single nucleotide polymorphisms (SNPs) in human regulatory regions, we established a link between disruptive local DNA shape changes and loss of specific TF binding. Furthermore, we described cases where disease-associated SNPs may alter TF binding through DNA shape changes. This link led us to hypothesize that local DNA shape within and around TF binding sites is under selection pressure. To verify this hypothesis, we analyzed SNP data derived from 216 natural strains of Drosophila melanogaster. Comparing SNPs located in functional and nonfunctional regions within experimentally validated cis-regulatory modules (CRMs) from D. melanogaster that are active in the blastoderm stage of development, we found that SNPs within functional regions tended to cause smaller DNA shape variations. Furthermore, SNPs with higher minor allele frequency were more likely to result in smaller DNA shape variations. The same analysis based on a large number of SNPs in putative CRMs of the D. melanogaster genome derived from DNase I accessibility data confirmed these observations. Taken together, our results indicate that common SNPs in functional regions tend to maintain DNA shape, whereas shape-disrupting SNPs are more likely to be eliminated through purifying selection.

Domesticated species are impacted in unintended ways during domestication and breeding. Changes in the nature and intensity of selection impart genetic drift, reduce diversity, and increase the frequency of deleterious alleles. Such outcomes constrain our ability to expand the cultivation of crops into environments that differ from those under which domestication occurred. We address this need in chickpea, an important pulse legume, by harnessing the diversity of wild crop relatives. We document an extreme domestication-related genetic bottleneck and decipher the genetic history of wild populations. We provide evidence of ancestral adaptations for seed coat color crypsis, estimate the impact of environment on genetic structure and trait values, and demonstrate variation between wild and cultivated accessions for agronomic properties. A resource of genotyped, association mapping progeny functionally links the wild and cultivated gene pools and is an essential resource chickpea for improvement, while our methods inform collection of other wild crop progenitor species.

Gene regulatory networks control the complex programs that drive development. Deciphering the connections between transcription factors (TFs) and target genes is challenging, in part because TFs bind to thousands of places in the genome but control expression through a subset of these binding events. We hypothesize that we can combine natural variation of expression levels and predictions of TF binding sites to identify TF targets. We gather RNA-seq data from 71 genetically distinct F1 Drosophila melanogaster embryos and calculate the correlations between TF and potential target genes' expression levels, which we call "regulatory strength." To separate direct and indirect TF targets, we hypothesize that direct TF targets will have a preponderance of binding sites in their upstream regions. Using 14 TFs active during embryogenesis, we find that 12 TFs showed a significant correlation between their binding strength and regulatory strength on downstream targets, and 10 TFs showed a significant correlation between the number of binding sites and the regulatory effect on target genes. The general roles, e.g. bicoid's role as an activator, and the particular interactions we observed between our TFs, e.g. twist's role as a repressor of sloppy paired and odd paired, generally coincide with the literature.

Noncoding DNA sequences, which play various roles in gene expression and regulation, are under evolutionary pressure. Gene regulation requires specific protein-DNA binding events, and our previous studies showed that both DNA sequence and shape readout are employed by transcription factors (TFs) to achieve DNA binding specificity. By investigating the shape-disrupting properties of single nucleotide polymorphisms (SNPs) in human regulatory regions, we established a link between disruptive local DNA shape changes and loss of specific TF binding. Furthermore, we described cases where disease-associated SNPs may alter TF binding through DNA shape changes. This link led us to hypothesize that local DNA shape within and around TF binding sites is under selection pressure. To verify this hypothesis, we analyzed SNP data derived from 216 natural strains of Drosophila melanogaster. Comparing SNPs located in functional and nonfunctional regions within experimentally validated cis-regulatory modules (CRMs) from D. melanogaster that are active in the blastoderm stage of development, we found that SNPs within functional regions tended to cause smaller DNA shape variations. Furthermore, SNPs with higher minor allele frequency were more likely to result in smaller DNA shape variations. The same analysis based on a large number of SNPs in putative CRMs of the D. melanogaster genome derived from DNase I accessibility data confirmed these observations. Taken together, our results indicate that common SNPs in functional regions tend to maintain DNA shape, whereas shape-disrupting SNPs are more likely to be eliminated through purifying selection.

Domesticated species are impacted in unintended ways during domestication and breeding. Changes in the nature and intensity of selection impart genetic drift, reduce diversity, and increase the frequency of deleterious alleles. Such outcomes constrain our ability to expand the cultivation of crops into environments that differ from those under which domestication occurred. We address this need in chickpea, an important pulse legume, by harnessing the diversity of wild crop relatives. We document an extreme domestication-related genetic bottleneck and decipher the genetic history of wild populations. We provide evidence of ancestral adaptations for seed coat color crypsis, estimate the impact of environment on genetic structure and trait values, and demonstrate variation between wild and cultivated accessions for agronomic properties. A resource of genotyped, association mapping progeny functionally links the wild and cultivated gene pools and is an essential resource chickpea for improvement, while our methods inform collection of other wild crop progenitor species.

The Vavilov Institute of Plant Genetic Resources (VIR), in St. Petersburg, Russia, houses a unique genebank, with historical collections of landraces. When they were collected, the geographical distribution and genetic diversity of most crops closely reflected their historical patterns of cultivation established over the preceding millennia. We employed a combination of genomics, computational biology and phenotyping to characterize VIR's 147 chickpea accessions from Turkey and Ethiopia, representing chickpea's center of origin and a major location of secondary diversity. Genotyping by sequencing identified 14,059 segregating polymorphisms and genome-wide association studies revealed 28 GWAS hits in potential candidate genes likely to affect traits of agricultural importance. The proportion of polymorphisms shared among accessions is a strong predictor of phenotypic resemblance, and of environmental similarity between historical sampling sites. We found that 20 out of 28 polymorphisms, associated with multiple traits, including days to maturity, plant phenology, and yield-related traits such as pod number, localized to chromosome 4. We hypothesize that selection and introgression via inadvertent hybridization between more and less advanced morphotypes might have resulted in agricultural improvement genes being aggregated to genomic 'agro islands', and in genotype-to-phenotype relationships resembling widespread pleiotropy.

The Vavilov Institute of Plant Genetic Resources (VIR), in St. Petersburg, Russia, houses a unique genebank, with historical collections of landraces. When they were collected, the geographical distribution and genetic diversity of most crops closely reflected their historical patterns of cultivation established over the preceding millennia. We employed a combination of genomics, computational biology and phenotyping to characterize VIR's 147 chickpea accessions from Turkey and Ethiopia, representing chickpea's center of origin and a major location of secondary diversity. Genotyping by sequencing identified 14,059 segregating polymorphisms and genome-wide association studies revealed 28 GWAS hits in potential candidate genes likely to affect traits of agricultural importance. The proportion of polymorphisms shared among accessions is a strong predictor of phenotypic resemblance, and of environmental similarity between historical sampling sites. We found that 20 out of 28 polymorphisms, associated with multiple traits, including days to maturity, plant phenology, and yield-related traits such as pod number, localized to chromosome 4. We hypothesize that selection and introgression via inadvertent hybridization between more and less advanced morphotypes might have resulted in agricultural improvement genes being aggregated to genomic 'agro islands', and in genotype-to-phenotype relationships resembling widespread pleiotropy.

High soil salinity negatively influences plant growth and yield. Some taxa have evolved mechanisms for avoiding or tolerating elevated soil salinity, which can be modulated by the environment experienced by parents or offspring. We tested the contribution of the parental and offspring environments on salinity adaptation and their potential underlying mechanisms. In a two-generation greenhouse experiment, we factorially manipulated salinity concentrations for genotypes of Medicago truncatula that were originally collected from natural populations that differed in soil salinity. To compare population level adaptation to soil salinity and to test the potential mechanisms involved we measured two aspects of plant performance, reproduction and vegetative biomass, and phenological and physiological traits associated with salinity avoidance and tolerance. Saline-origin populations had greater biomass and reproduction under saline conditions than non-saline populations, consistent with local adaptation to saline soils. Additionally, parental environmental exposure to salt increased this difference in performance. In terms of environmental effects on mechanisms of salinity adaptation, parental exposure to salt spurred phenological differences that facilitated salt avoidance, while offspring exposure to salt resulted in traits associated with greater salt tolerance. Non-saline origin populations expressed traits associated with greater growth in the absence of salt while, for saline adapted populations, the ability to maintain greater performance in saline environments was also associated with lower growth potential in the absence of salt. Plastic responses induced by parental and offspring environments in phenology, leaf traits, and gas exchange contribute to salinity adaptation in M. truncatula. The ability of plants to tolerate environmental stress, such as high soil salinity, is likely modulated by a combination of parental effects and within-generation phenotypic plastic

High soil salinity negatively influences plant growth and yield. Some taxa have evolved mechanisms for avoiding or tolerating elevated soil salinity, which can be modulated by the environment experienced by parents or offspring. We tested the contribution of the parental and offspring environments on salinity adaptation and their potential underlying mechanisms. In a two-generation greenhouse experiment, we factorially manipulated salinity concentrations for genotypes of Medicago truncatula that were originally collected from natural populations that differed in soil salinity. To compare population level adaptation to soil salinity and to test the potential mechanisms involved we measured two aspects of plant performance, reproduction and vegetative biomass, and phenological and physiological traits associated with salinity avoidance and tolerance. Saline-origin populations had greater biomass and reproduction under saline conditions than non-saline populations, consistent with local adaptation to saline soils. Additionally, parental environmental exposure to salt increased this difference in performance. In terms of environmental effects on mechanisms of salinity adaptation, parental exposure to salt spurred phenological differences that facilitated salt avoidance, while offspring exposure to salt resulted in traits associated with greater salt tolerance. Non-saline origin populations expressed traits associated with greater growth in the absence of salt while, for saline adapted populations, the ability to maintain greater performance in saline environments was also associated with lower growth potential in the absence of salt. Plastic responses induced by parental and offspring environments in phenology, leaf traits, and gas exchange contribute to salinity adaptation in M. truncatula. The ability of plants to tolerate environmental stress, such as high soil salinity, is likely modulated by a combination of parental effects and within-generation phenotypic plastic

BackgroundAs our world becomes warmer, agriculture is increasingly impacted by rising soil salinity and understanding plant adaptation to salt stress can help enable effective crop breeding. Salt tolerance is a complex plant phenotype and we know little about the pathways utilized by naturally tolerant plants. Legumes are important species in agricultural and natural ecosystems, since they engage in symbiotic nitrogen-fixation, but are especially vulnerable to salinity stress.ResultsOur studies of the model legume Medicago truncatula in field and greenhouse settings demonstrate that Tunisian populations are locally adapted to saline soils at the metapopulation level and that saline origin genotypes are less impacted by salt than non-saline origin genotypes; these populations thus likely contain adaptively diverged alleles. Whole genome resequencing of 39 wild accessions reveals ongoing migration and candidate genomic regions that assort non-randomly with soil salinity. Consistent with natural selection acting at these sites, saline alleles are typically rare in the range-wide species' gene pool and are also typically derived relative to the sister species M. littoralis. Candidate regions for adaptation contain genes that regulate physiological acclimation to salt stress, such as abscisic acid and jasmonic acid signaling, including a novel salt-tolerance candidate orthologous to the uncharacterized gene AtCIPK21. Unexpectedly, these regions also contain biotic stress genes and flowering time pathway genes. We show that flowering time is differentiated between saline and non-saline populations and may allow salt stress escape.ConclusionsThis work nominates multiple potential pathways of adaptation to naturally stressful environments in a model legume. These candidates point to the importance of both tolerance and avoidance in natural legume populations. We have uncovered several promising targets that could be used to breed for enhanced salt tolerance in crop legumes

BackgroundAs our world becomes warmer, agriculture is increasingly impacted by rising soil salinity and understanding plant adaptation to salt stress can help enable effective crop breeding. Salt tolerance is a complex plant phenotype and we know little about the pathways utilized by naturally tolerant plants. Legumes are important species in agricultural and natural ecosystems, since they engage in symbiotic nitrogen-fixation, but are especially vulnerable to salinity stress.ResultsOur studies of the model legume Medicago truncatula in field and greenhouse settings demonstrate that Tunisian populations are locally adapted to saline soils at the metapopulation level and that saline origin genotypes are less impacted by salt than non-saline origin genotypes; these populations thus likely contain adaptively diverged alleles. Whole genome resequencing of 39 wild accessions reveals ongoing migration and candidate genomic regions that assort non-randomly with soil salinity. Consistent with natural selection acting at these sites, saline alleles are typically rare in the range-wide species' gene pool and are also typically derived relative to the sister species M. littoralis. Candidate regions for adaptation contain genes that regulate physiological acclimation to salt stress, such as abscisic acid and jasmonic acid signaling, including a novel salt-tolerance candidate orthologous to the uncharacterized gene AtCIPK21. Unexpectedly, these regions also contain biotic stress genes and flowering time pathway genes. We show that flowering time is differentiated between saline and non-saline populations and may allow salt stress escape.ConclusionsThis work nominates multiple potential pathways of adaptation to naturally stressful environments in a model legume. These candidates point to the importance of both tolerance and avoidance in natural legume populations. We have uncovered several promising targets that could be used to breed for enhanced salt tolerance in crop legumes

Background Allotetraploids carry pairs of diverged homoeologs for most genes. With the genome doubled in size, the number of putative interactions is enormous. This poses challenges on how to coordinate the two disparate genomes, and creates opportunities by enhancing the phenotypic variation. New combinations of alleles co-adapt and respond to new environmental pressures. Three stages of the allopolyploidization process - parental species divergence, hybridization, and genome duplication - have been well analyzed. The last stage of evolutionary adjustments remains mysterious. Results Homoeolog-specific retention and use were analyzed in Arabidopsis suecica (As), a species derived from A. thaliana (At) and A. arenosa (Aa) in a single event 12,000 to 300,000 years ago. We used 405,466 diagnostic features on tiling microarrays to recognize At and Aa contributions to the As genome and transcriptome: 324 genes lacked Aa contributions and 614 genes lacked At contributions within As. In leaf tissues, 3,458 genes preferentially expressed At homoeologs while 4,150 favored Aa homoeologs. These patterns were validated with resequencing. Genes with preferential use of Aa homoeologs were enriched for expression functions, consistent with the dominance of Aa transcription. Heterologous networks - mixed from At and Aa transcripts - were underrepresented. Conclusions Thousands of deleted and silenced homoeologs in the genome of As were identified. Since heterologous networks may be compromised by interspecies incompatibilities, these networks evolve co-biases, expressing either only Aa or only At homoeologs. This progressive change towards predominantly pure parental networks might contribute to phenotypic variability and plasticity, and enable the species to exploit a larger range of environments.

Background Allotetraploids carry pairs of diverged homoeologs for most genes. With the genome doubled in size, the number of putative interactions is enormous. This poses challenges on how to coordinate the two disparate genomes, and creates opportunities by enhancing the phenotypic variation. New combinations of alleles co-adapt and respond to new environmental pressures. Three stages of the allopolyploidization process - parental species divergence, hybridization, and genome duplication - have been well analyzed. The last stage of evolutionary adjustments remains mysterious. Results Homoeolog-specific retention and use were analyzed in Arabidopsis suecica (As), a species derived from A. thaliana (At) and A. arenosa (Aa) in a single event 12,000 to 300,000 years ago. We used 405,466 diagnostic features on tiling microarrays to recognize At and Aa contributions to the As genome and transcriptome: 324 genes lacked Aa contributions and 614 genes lacked At contributions within As. In leaf tissues, 3,458 genes preferentially expressed At homoeologs while 4,150 favored Aa homoeologs. These patterns were validated with resequencing. Genes with preferential use of Aa homoeologs were enriched for expression functions, consistent with the dominance of Aa transcription. Heterologous networks - mixed from At and Aa transcripts - were underrepresented. Conclusions Thousands of deleted and silenced homoeologs in the genome of As were identified. Since heterologous networks may be compromised by interspecies incompatibilities, these networks evolve co-biases, expressing either only Aa or only At homoeologs. This progressive change towards predominantly pure parental networks might contribute to phenotypic variability and plasticity, and enable the species to exploit a larger range of environments.

Genes with overlapping expression and function may gradually diverge despite retaining some common functions. To test whether such genes show distinct patterns of molecular evolution within species, we examined sequence variation at the bric à brac (bab) locus of Drosophila melanogaster. This locus is composed of two anciently duplicated paralogs, bab1 and bab2, which are involved in patterning the adult abdomen, legs, and ovaries. We have sequenced the 148 kb genomic region spanning the bab1 and bab2 genes from 94 inbred lines of D. melanogaster sampled from a single location. Two non-coding regions, one in each paralog, appear to be under selection. The strongest evidence of directional selection is found in a region of bab2 that has no known functional role. The other region is located in the bab1 paralog and is known to contain a cis-regulatory element that controls sex-specific abdominal pigmentation. The coding region of bab1 appears to be under stronger functional constraint than the bab2 coding sequences. Thus, the two paralogs are evolving under different selective regimes in the same natural population, illuminating the different evolutionary trajectories of partially redundant duplicate genes.

Genes with overlapping expression and function may gradually diverge despite retaining some common functions. To test whether such genes show distinct patterns of molecular evolution within species, we examined sequence variation at the bric à brac (bab) locus of Drosophila melanogaster. This locus is composed of two anciently duplicated paralogs, bab1 and bab2, which are involved in patterning the adult abdomen, legs, and ovaries. We have sequenced the 148 kb genomic region spanning the bab1 and bab2 genes from 94 inbred lines of D. melanogaster sampled from a single location. Two non-coding regions, one in each paralog, appear to be under selection. The strongest evidence of directional selection is found in a region of bab2 that has no known functional role. The other region is located in the bab1 paralog and is known to contain a cis-regulatory element that controls sex-specific abdominal pigmentation. The coding region of bab1 appears to be under stronger functional constraint than the bab2 coding sequences. Thus, the two paralogs are evolving under different selective regimes in the same natural population, illuminating the different evolutionary trajectories of partially redundant duplicate genes.

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 and2b, 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.

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.

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 and2b, 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.