Chromosomal inversions and ecotypic differentiation in Anopheles gambiae: the perspective from whole‐genome sequencing

The molecular mechanisms and genetic architecture that facilitate adaptive radiation of lineages remain elusive. Polymorphic chromosomal inversions, due to their recombination‐reducing effect, are proposed instruments of ecotypic differentiation. Here, we study an ecologically diversifying lineage of Anopheles gambiae, known as the Bamako chromosomal form based on its unique complement of three chromosomal inversions, to explore the impact of these inversions on ecotypic differentiation. We used pooled and individual genome sequencing of Bamako, typical (non‐Bamako) An. gambiae and the sister species Anopheles coluzzii to investigate evolutionary relationships and genomewide patterns of nucleotide diversity and differentiation among lineages. Despite extensive shared polymorphism and limited differentiation from the other taxa, Bamako clusters apart from the other taxa, and forms a maximally supported clade in neighbour‐joining trees based on whole‐genome data (including inversions) or solely on collinear regions. Nevertheless, F_ST outlier analysis reveals that the majority of differentiated regions between Bamako and typical An. gambiae are located inside chromosomal inversions, consistent with their role in the ecological isolation of Bamako. Exceptionally differentiated genomic regions were enriched for genes implicated in nervous system development and signalling. Candidate genes associated with a selective sweep unique to Bamako contain substitutions not observed in sympatric samples of the other taxa, and several insecticide resistance gene alleles shared between Bamako and other taxa segregate at sharply different frequencies in these samples. Bamako represents a useful window into the initial stages of ecological and genomic differentiation from sympatric populations in this important group of malaria vectors.     

PSMC analysis of effective population sizes in molecular ecology and its application to black‐and‐white Ficedula flycatchers

Climatic fluctuations during the Quaternary period governed the demography of species and contributed to population differentiation and ultimately speciation. Studies of these past processes have previously been hindered by a lack of means and genetic data to model changes in effective population size (N_e) through time. However, based on diploid genome sequences of high quality, the recently developed pairwise sequentially Markovian coalescent (PSMC) can estimate trajectories of changes in N_e over considerable time periods. We applied this approach to resequencing data from nearly 200 genomes of four species and several populations of the Ficedula species complex of black‐and‐white flycatchers. N_e curves of Atlas, collared, pied and semicollared flycatcher converged 1–2 million years ago (Ma) at an N_e of ≈ 200 000, likely reflecting the time when all four species last shared a common ancestor. Subsequent separate N_e trajectories are consistent with lineage splitting and speciation. All species showed evidence of population growth up until 100–200 thousand years ago (kya), followed by decline and then start of a new phase of population expansion. However, timing and amplitude of changes in N_e differed among species, and for pied flycatcher, the temporal dynamics of N_e differed between Spanish birds and central/northern European populations. This cautions against extrapolation of demographic inference between lineages and calls for adequate sampling to provide representative pictures of the coalescence process in different species or populations. We also empirically evaluate criteria for proper inference of demographic histories using PSMC and arrive at recommendations of using sequencing data with a mean genome coverage of ≥18X, a per‐site filter of ≥10 reads and no more than 25% of missing data.     

Rapid identification of lettuce seed germination mutants by bulked segregant analysis and whole genome sequencing

Lettuce (Lactuca sativa) seeds exhibit thermoinhibition, or failure to complete germination when imbibed at warm temperatures. Chemical mutagenesis was employed to develop lettuce lines that exhibit germination thermotolerance. Two independent thermotolerant lettuce seed mutant lines, TG01 and TG10, were generated through ethyl methanesulfonate mutagenesis. Genetic and physiological analyses indicated that these two mutations were allelic and recessive. To identify the causal gene(s), we applied bulked segregant analysis by whole genome sequencing. For each mutant, bulked DNA samples of segregating thermotolerant (mutant) seeds were sequenced and analyzed for homozygous single‐nucleotide polymorphisms. Two independent candidate mutations were identified at different physical positions in the zeaxanthin epoxidase gene (ABSCISIC ACID DEFICIENT 1/ZEAXANTHIN EPOXIDASE, or ABA1/ZEP) in TG01 and TG10. The mutation in TG01 caused an amino acid replacement, whereas the mutation in TG10 resulted in alternative mRNA splicing. Endogenous abscisic acid contents were reduced in both mutants, and expression of the ABA1 gene from wild‐type lettuce under its own promoter fully complemented the TG01 mutant. Conventional genetic mapping confirmed that the causal mutations were located near the ZEP/ABA1 gene, but the bulked segregant whole genome sequencing approach more efficiently identified the specific gene responsible for the phenotype.     

Targeted resequencing reveals geographical patterns of differentiation for loci implicated in parallel evolution

Parallel divergence and speciation provide evidence for the role of divergent selection in generating biological diversity. Recent studies indicate that parallel phenotypic divergence may not have the same genetic basis in different geographical locations – ‘outlier loci’ (loci potentially affected by divergent selection) are often not shared among parallel instances of phenotypic divergence. However, limited sharing may be due, in part, to technical issues if false‐positive outliers occur. Here, we test this idea in the marine snail Littorina saxatilis, which has evolved two partly isolated ecotypes (adapted to crab predation vs. wave action) in multiple locations independently. We argue that if the low extent of sharing observed in earlier studies in this system is due to sampling effects, we expect outliers not to show elevated F_ST when sequenced in new samples from the original locations and also not to follow predictable geographical patterns of elevated F_ST. Following a hierarchical sampling design (within vs. between country), we applied capture sequencing, targeting outliers from earlier studies and control loci. We found that outliers again showed elevated levels of F_ST in their original location, suggesting they were not generated by sampling effects. Outliers were also likely to show increased F_ST in geographically close locations, which may be explained by higher levels of gene flow or shared ancestral genetic variation compared with more distant locations. However, in contrast to earlier findings, we also found some outlier types to show elevated F_ST in geographically distant locations. We discuss possible explanations for this unexpected result.     

Genomic sequence and virulence evaluation of MN184A‐like porcine reproductive and respiratory syndrome virus in Japan

In this study, a porcine reproductive and respiratory syndrome virus (PRRSV) that was isolated from a 9‐week‐old diseased pig on a farm in Japan with a high mortality rate during 2007–2008 was characterized. This unique isolate, designated as Jpn5‐37, did not have a high nucleotide identity in open reading frame 5 against any Japanese isolates. Among all available type 2 PRRSV complete genome sequences, Jpn5‐37 shared the highest nucleotide identity (93.6%) with virulent strain MN184A. The genomic characteristics of Jpn5‐37 were highly conserved with respect to the virulent MN184A, including a continuous eight amino acid deletion in the nonstructural protein 2 region. Moreover, virus distribution, viremia and the gross and microscopic characteristics of lesions were investigated in pigs 10 days post‐inoculation to elucidate the pathogenicity of the isolate. Intranasal inoculation was found to rapidly result in viremia and dissemination of the Jpn5‐37 isolate to several tissues in a similar manner to EDRD1; however, the amounts of Jpn5‐37 RNA in serum were significantly greater. Similarly, the quantities of Jpn5‐37 viral RNA in all organs tested tended to be higher than with EDRD1 infection. Mean rectal temperatures were significantly higher in the Jpn5‐37‐inoculated than in the control group at 4 and 6 days post infection (dpi) and in the EDRD1‐inoculated group at 6 and 8 dpi. These results suggest that the Jpn5‐37 strain replicates and is more efficiently distributed to the organs than is EDRD1 under the same conditions.     

Persistence of a Wolbachia infection frequency cline in Drosophila melanogaster and the possible role of reproductive dormancy

Field populations of arthropods are often polymorphic for Wolbachia but the factors maintaining intermediate Wolbachia frequencies are generally not understood. In Drosophila melanogaster, Wolbachia frequencies are highly variable across the globe. We document the persistence of a Wolbachia infection frequency cline in D. melanogaster populations from eastern Australia across at least 20 years, with frequencies generally high in the tropics but lower in cool temperate regions. The results are interpreted using a model of frequency dynamics incorporating cytoplasmic incompatibility (CI), imperfect maternal transmission and Wolbachia effects on fitness. Clinal variation is less pronounced in eastern North America which may reflect annual recolonization at higher latitudes. Limited samples from Africa from latitudes matching our tropical and subtropical samples from Australia and North America show comparably high infection frequencies, but some equatorial samples show lower frequencies. Adult dormancy across cold periods may contribute to the Australian Wolbachia cline. Infected flies exposed to cold conditions for an extended period had reduced fecundity and viability, an effect not evident in unexposed controls. These fitness costs may contribute to the relatively low Wolbachia frequencies in Australian temperate areas; whereas different processes, including CI induced by young males, may contribute to higher frequencies in tropical locations.     

Minute Y chromosomes and karyotype evolution in Madagascan iguanas (Squamata: Iguania: Opluridae)

Iguanas (Pleurodonta) are predominantly distributed in the New World, but one previously cytogenetically understudied family, Opluridae, is endemic to Madagascar and the adjacent Grand Comoro archipelago. The aim of our contribution is to fill a gap in the cytogenetic understanding of this biogeographically puzzling lineage. Based on examination of six species, we found that oplurids are rather conservative in karyotype, which is composed of 36 chromosomes as in most iguanas. However, the species differ in the position of the nucleolar organizer region and heterochromatic blocks and in the accumulation and distribution of interstitial telomeric sequences (ITSs), which suggests cryptic intra‐ and interchromosomal rearrangements. All tested species share the XY sex‐determining system homologous to most other iguana families. The oplurid Y chromosome is degenerated, very small in size but mostly euchromatic. Fluorescence in situ hybridization with probes composed of microsatellite motifs revealed variability among species in the accumulation of particular repeats on the Y chromosome. This variability accounts for the differences in the detection of sex chromosomes across the species of the family using comparative genome hybridization (CGH) technique. Our study demonstrates the limits of the commonly used CGH technique to uncover sex chromosomes even in organisms with heteromorphic and sequentially largely differentiated sex chromosomes.     

Did sex chromosome turnover promote divergence of the major mammal groups?

Comparative mapping and sequencing show that turnover of sex determining genes and chromosomes, and sex chromosome rearrangements, accompany speciation in many vertebrates. Here I review the evidence and propose that the evolution of therian mammals was precipitated by evolution of the male‐determining SRY gene, defining a novel XY sex chromosome pair, and interposing a reproductive barrier with the ancestral population of synapsid reptiles 190 million years ago (MYA). Divergence was reinforced by multiple translocations in monotreme sex chromosomes, the first of which supplied a novel sex determining gene. A sex chromosome‐autosome fusion may have separated eutherians (placental mammals) from marsupials 160 MYA. Another burst of sex chromosome change and speciation is occurring in rodents, precipitated by the degradation of the Y. And although primates have a more stable Y chromosome, it may be just a matter of time before the same fate overtakes our own lineage. Also watch the video abstract .     

P53 and the defenses against genome instability caused by transposons and repetitive elements

The recent publication by Wylie et al. is reviewed, demonstrating that the p53 protein regulates the movement of transposons. While this work presents genetic evidence for a piRNA‐mediated p53 interaction with transposons in Drosophila and zebrafish, it is herein placed in the context of a decade or so of additional work that demonstrated a role for p53 in regulating transposons and other repetitive elements. The line of thought in those studies began with the observation that transposons damage DNA and p53 regulates DNA damage. The presence of transposon movement can increase the rate of evolution in the germ line and alter genes involved in signal transduction pathways. Transposition can also play an important role in cancers where the p53 gene function is often mutated. This is particularly interesting as recent work has shown that de‐repression of repetitive elements in cancer has important consequences for the immune system and tumor microenvironment.