Reproductive Capacity Evolves in Response to Ecology through Common Changes in Cell Number in Hawaiian Drosophila

Gene Protein Sequence Evolution Can Predict the Rapid Divergence of Ovariole Numbers in the Drosophila melanogaster Subgroup

Ovaries play key roles in fitness and evolution: they are essential female reproductive structures that develop and house the eggs in sexually reproducing animals. In Drosophila, the mature ovary contains multiple tubular egg-producing structures known as ovarioles. Ovarioles arise from somatic cellular structures in the larval ovary called terminal filaments (TFs), formed by TF cells and subsequently enclosed by sheath (SH) cells. As in many other insects, ovariole number per female varies extensively in Drosophila. At present, however, there is a striking gap of information on genetic mechanisms and evolutionary forces that shape the well-documented rapid interspecies divergence of ovariole numbers. To address this gap, here we studied genes associated with Drosophila melanogaster ovariole number or functions based on recent experimental and transcriptional datasets from larval ovaries, including TFs and SH cells, and assessed their rates and patterns of molecular evolution in five closely related species of the melanogaster subgroup that exhibit species-specific differences in ovariole numbers. From comprehensive analyses of protein sequence evolution (dN/dS), branch-site positive selection, expression specificity (tau), and phylogenetic regressions (phylogenetic generalized least squares), we report evidence of 42 genes that showed signs of playing roles in the genetic basis of interspecies evolutionary change of Drosophila ovariole number. These included the signaling genes upd2 and Ilp5 and extracellular matrix genes vkg and Col4a1, whose dN/dS predicted ovariole numbers among species. Together, we propose a model whereby a set of ovariole-involved gene proteins have an enhanced evolvability, including adaptive evolution, facilitating rapid shifts in ovariole number among Drosophila species.

Protective Geometry and Reproductive Anatomy as Candidate Determinants of Clutch Size Variation in Pentatomid Bugs

AbstractMany animals lay their eggs in clusters. Eggs on the periphery of clusters can be at higher risk of mortality. We asked whether the most commonly occurring clutch sizes in pentatomid bugs could result from geometrical arrangements that maximize the proportion of eggs in the cluster's interior. Although the most common clutch sizes do not correspond with geometric optimality, stink bugs do tend to lay clusters of eggs in shapes that protect increasing proportions of their offspring as clutch sizes increase. We also considered whether ovariole number, an aspect of reproductive anatomy that may be a fixed trait across many pentatomids, could explain observed distributions of clutch sizes. The most common clutch sizes across many species correspond with multiples of ovariole number. However, there are species with the same number of ovarioles that lay clutches of widely varying size, among which multiples of ovariole number are not overrepresented. In pentatomid bugs, reproductive anatomy appears to be more important than egg mass geometry in determining clutch size uniformity. In addition, our analysis demonstrates that groups of animals with little variation in ovariole number may nonetheless lay a broad range of clutch shapes and sizes.

Environmental microbes promote phenotypic plasticity in reproduction and sleep behaviour

The microbiome has been hypothesized as a driving force of phenotypic variation in host organisms that is capable of extending metabolic processes, altering development and in some cases, conferring novel functions that are critical for survival. Only a few studies have directly shown a causal role for the environmental microbiome in altering host phenotypic features. To assess the extent to which environmental microbes induce variation in host life-history traits and behaviour, we inoculated axenic Drosophila melanogaster with microbes isolated from drosophilid populations collected from two different field sites and generated two populations with distinct bacterial and fungal profiles. We show that microbes isolated from environmental sites with modest abiotic differences induce large variation in host reproduction, fatty acid levels, stress tolerance and sleep behaviour. Importantly, clearing microbes from each experimental population removed the phenotypic differences. The results support the causal role of environmental microbes as drivers of host phenotypic plasticity and potentially, rapid adaptation and evolution.

Retinoblastoma protein activity revealed by CRISPRi study of divergent Rbf1 and Rbf2 paralogs

Retinoblastoma tumor suppressor proteins regulate the key transition from G1 to S phase of the cell cycle. The mammalian Rb family comprises Rb, p107, and p130, with overlapping and unique roles in gene regulation. Drosophila experienced an independent gene duplication event, leading to the Rbf1 and Rbf2 paralogs. To uncover the significance of paralogy in the Rb family, we used CRISPRi. We engineered dCas9 fusions to Rbf1 and Rbf2, and deployed them to gene promoters in developing Drosophila tissue to study their relative impacts on gene expression. On some genes, both Rbf1 and Rbf2 mediate potent repression, in a highly distance-dependent manner. In other cases, the two proteins have different effects on phenotype and gene expression, indicating different functional potential. In a direct comparison of Rb activity on endogenous genes and transiently transfected reporters, we found that only qualitative, but not key quantitative aspects of repression were conserved, indicating that the native chromatin environment generates context-specific effects of Rb activity. Our study uncovers the complexity of Rb-mediated transcriptional regulation in a living organism, which is clearly impacted by the different promoter landscapes and the evolution of the Rb proteins themselves.

The evolution of ovary-biased gene expression in Hawaiian Drosophila

With detailed data on gene expression accessible from an increasingly broad array of species, we can test the extent to which our developmental genetic knowledge from model organisms predicts expression patterns and variation across species. But to know when differences in gene expression across species are significant, we first need to know how much evolutionary variation in gene expression we expect to observe. Here we provide an answer by analyzing RNAseq data across twelve species of Hawaiian Drosophilidae flies, focusing on gene expression differences between the ovary and other tissues. We show that over evolutionary time, there exists a cohort of ovary specific genes that is stable and that largely corresponds to described expression patterns from laboratory model Drosophila species. Our results also provide a demonstration of the prediction that, as phylogenetic distance increases, variation between species overwhelms variation between tissue types. Using ancestral state reconstruction of expression, we describe the distribution of evolutionary changes in tissue-biased expression, and use this to identify gains and losses of ovary-biased expression across these twelve species. We then use this distribution to calculate the evolutionary correlation in expression changes between genes, and demonstrate that genes with known interactions in D. melanogaster are significantly more correlated in their evolution than genes with no or unknown interactions. Finally, we use this correlation matrix to infer new networks of genes that share evolutionary trajectories, and we present these results as a dataset of new testable hypotheses about genetic roles and interactions in the function and evolution of the Drosophila ovary.

The early embryonic transcriptome of a Hawaiian Drosophila picture-wing fly shows evidence of altered gene expression and novel gene evolution

A massive adaptive radiation on the Hawaiian archipelago has produced approximately one-quarter of the fly species in the family Drosophilidae. The Hawaiian Drosophila clade has long been recognized as a model system for the study of both the ecology of island endemics and the evolution of developmental mechanisms, but relatively few genomic and transcriptomic datasets are available for this group. We present here a differential expression analysis of the transcriptional profiles of two highly conserved embryonic stages in the Hawaiian picture-wing fly Drosophila grimshawi. When we compared our results to previously published datasets across the family Drosophilidae, we identified cases of both gains and losses of gene representation in D. grimshawi, including an apparent delay in Hox gene activation. We also found a high expression of unannotated genes. Most transcripts of unannotated genes with open reading frames do not have identified homologs in non-Hawaiian Drosophila species, although the vast majority have sequence matches in genomes of other Hawaiian picture-wing flies. Some of these unannotated genes may have arisen from noncoding sequence in the ancestor of Hawaiian flies or during the evolution of the clade. Our results suggest that both the modified use of ancestral genes and the evolution of new ones may occur in rapid radiations.

Transcriptome Sequencing and Bioinformatics Analysis of Ovarian Tissues from Pomacea canaliculata in Guangdong and Hunan

In this study, the fecundity of Pomacea canaliculata was studied by collecting egg masses from Guangdong and Hunan using field egg collection and indoor propagation. Through high-throughput RNA sequencing (RNA-seq), we analyzed the ovarian tissue of the snails in Guangdong (G_O) and those in Hunan (H_O) using comparative analysis of transcription. Moreover, we used bioinformatics methods to screen the key pathways and genes that affect the fecundity of snails from the two locations. Results. The results showed that the absolute fecundity and weight-relative fecundity of Pomacea canaliculata in Guangdong were significantly higher than those in Hunan. We found 1,546 differential genes through differential gene screening (528 genes upregulated in snails from Guangdong and 1018 in snails from Hunan). The ribosomal signaling pathway and rpl23a, uba52 are critical pathways and essential genes that affect the fecundity of snails. Conclusions. The 27 differential genes in the ribosome signaling pathway, collected from H_O, were all downregulated. As a result, ovarian tissue protein synthesis is impaired, which is an important mechanism that affects snails' ability to reproduce.

Target of Rapamycin Drives Unequal Responses to Essential Amino Acid Depletion for Egg Laying in Drosophila Melanogaster

Nutrition shapes a broad range of life-history traits, ultimately impacting animal fitness. A key fitness-related trait, female fecundity is well known to change as a function of diet. In particular, the availability of dietary protein is one of the main drivers of egg production, and in the absence of essential amino acids egg laying declines. However, it is unclear whether all essential amino acids have the same impact on phenotypes like fecundity. Using a holidic diet, we fed adult female Drosophila melanogaster diets that contained all necessary nutrients except one of the 10 essential amino acids and assessed the effects on egg production. For most essential amino acids, depleting a single amino acid induced as rapid a decline in egg production as when there were no amino acids in the diet. However, when either methionine or histidine were excluded from the diet, egg production declined more slowly. Next, we tested whether GCN2 and TOR mediated this difference in response across amino acids. While mutations in GCN2 did not eliminate the differences in the rates of decline in egg laying among amino acid drop-out diets, we found that inhibiting TOR signalling caused egg laying to decline rapidly for all drop-out diets. TOR signalling does this by regulating the yolk-forming stages of egg chamber development. Our results suggest that amino acids differ in their ability to induce signalling via the TOR pathway. This is important because if phenotypes differ in sensitivity to individual amino acids, this generates the potential for mismatches between the output of a pathway and the animal’s true nutritional status.

Phylotranscriptomics reveals discordance in the phylogeny of Hawaiian Drosophila and Scaptomyza (Diptera: Drosophilidae)

Island radiations present natural laboratories for studying the evolutionary process. The Hawaiian Drosophilidae are one such radiation, with nearly 600 described species and substantial morphological and ecological diversification. These species are largely divided into a few major clades, but the relationship between clades remains uncertain. Here, we present new assembled transcriptomes from 12 species across these clades, and use these transcriptomes to resolve the base of the evolutionary radiation. We recover a new hypothesis for the relationship between clades, and demonstrate its support over previously published hypotheses. We then use the evolutionary radiation to explore dynamics of concordance in phylogenetic support, by analyzing the gene and site concordance factors for every possible topological combination of major groups. We show that high bootstrap values mask low evolutionary concordance, and we demonstrate that the most likely topology is distinct from the topology with the highest support across gene trees and from the topology with highest support across sites. We then combine all previously published genetic data for the group to estimate a time-calibrated tree for over 300 species of drosophilids. Finally, we digitize dozens of published Hawaiian Drosophilidae descriptions, and use this to pinpoint probable evolutionary shifts in reproductive ecology as well as body, wing, and egg size. We show that by examining the entire landscape of tree and trait space, we can gain a more complete understanding of how evolutionary dynamics play out across an island radiation.

Larger males facilitate population expansion in Ophraella communa

One of the most intriguing concepts in animal ecology is the reproductive advantages offered by larger body size, and the females prefer to mate with larger males to gain reproductive advantage. Currently, it is not clear how females recognize signs of male 'quality' and what mechanisms are involved in producing offspring with direct or indirect benefits. Our study aims to assess the preferences of females for males in Ophraella communa, determine the reproductive benefits and reveal the underlying mechanism behind this advantage. We demonstrate that male body size is an important determinant in the evolutionary process of O. communa, affecting female mate choice. Moreover, our study establishes that females prefer males with a larger body size, and this could further improve the developmental and reproductive fitness of their offspring. Finally, we focus on the seminal fluid proteins (SFPs) in O. communa, determine differentially expressed genes (i.e. OcACE, OcCBP and OcSFP) by analysing their proteomes and transcriptomes, and define the role of these SFPs-related genes through RNAi. Our study proved that the reproductive benefit of large males may be regulated by biased expression of crucial SFPs genes. The present study advances our understanding of the biological significance of preferential mating.

Repeated loss of variation in insect ovary morphology highlights the role of development in life-history evolution

The number of offspring an organism can produce is a key component of its evolutionary fitness and life history. Here we perform a test of the hypothesized trade-off between the number and size of offspring using thousands of descriptions of the number of egg-producing compartments in the insect ovary (ovarioles), a common proxy for potential offspring number in insects. We find evidence of a negative relationship between egg size and ovariole number when accounting for adult body size. However, in contrast to prior claims, we note that this relationship is not generalizable across all insect clades, and we highlight several factors that may have contributed to this size-number trade-off being stated as a general rule in previous studies. We reconstruct the evolution of the arrangement of cells that contribute nutrients and patterning information during oogenesis (nurse cells), and show that the diversification of ovariole number and egg size have both been largely independent of their presence or position within the ovariole. Instead, we show that ovariole number evolution has been shaped by a series of transitions between variable and invariant states, with multiple independent lineages evolving to have almost no variation in ovariole number. We highlight the implications of these invariant lineages on our understanding of the specification of ovariole number during development, as well as the importance of considering developmental processes in theories of life-history evolution.

Topology-driven protein-protein interaction network analysis detects genetic sub-networks regulating reproductive capacity

Understanding the genetic regulation of organ structure is a fundamental problem in developmental biology. Here, we use egg-producing structures of insect ovaries, called ovarioles, to deduce systems-level gene regulatory relationships from quantitative functional genetic analysis. We previously showed that Hippo signalling, a conserved regulator of animal organ size, regulates ovariole number in Drosophila melanogaster. To comprehensively determine how Hippo signalling interacts with other pathways in this regulation, we screened all known signalling pathway genes, and identified Hpo-dependent and Hpo-independent signalling requirements. Network analysis of known protein-protein interactions among screen results identified independent gene regulatory sub-networks regulating one or both of ovariole number and egg laying. These sub-networks predict involvement of previously uncharacterised genes with higher accuracy than the original candidate screen. This shows that network analysis combining functional genetic and large-scale interaction data can predict function of novel genes regulating development.