Large Chromosomal Inversion Research Articles

Genomic and transcriptomic analyses of the elite rice variety Huizhan provide insight into disease resistance and heat tolerance

The indica rice variety Huizhan shows elite traits of disease resistance and heat tolerance. However, the underlying genetic basis of these traits is not fully understood due to limited genomic resources. Here, we used Nanopore long-read and next-generation sequencing technologies to generate a chromosome-scale genome assembly of Huizhan. Comparative genomics analysis uncovered a large chromosomal inversion and expanded gene families that are associated with plant growth, development and stress responses. Functional rice blast resistance genes, including Pi2, Pib and Ptr, and bacterial blight resistance gene Xa27, contribute to disease resistance of Huizhan. Furthermore, integrated genomics and transcriptomics analyses showed that OsHIRP1, OsbZIP60, the SOD gene family, and various transcription factors are involved in heat tolerance of Huizhan. The high-quality genome assembly and comparative genomics results presented in this study facilitate the use of Huizhan as an elite parental line in developing rice varieties adapted to disease pressure and climate challenges.

Drosophila melanogaster Transcriptome Response to Different Wolbachia Strains.

Wolbachia is a maternally inherited, intercellular bacterial symbiont of insects and some other invertebrates. Here, we investigated the effect of two different Wolbachia strains, differing in a large chromosomal inversion, on the differential expression of genes in D. melanogaster females. We revealed significant changes in the transcriptome of the infected flies compared to the uninfected ones, as well as in the transcriptome of flies infected with the Wolbachia strain, wMelPlus, compared to flies infected with the wMelCS112 strain. We linked differentially expressed genes (DEGs) from two pairwise comparisons, "uninfected-wMelPlus-infected" and "uninfected-wMelCS112-infected", into two gene networks, in which the following functional groups were designated: "Proteolysis", "Carbohydrate transport and metabolism", "Oxidation-reduction process", "Embryogenesis", "Transmembrane transport", "Response to stress" and "Alkaline phosphatases". Our data emphasized similarities and differences between infections by different strains under study: a wMelPlus infection results in more than double the number of upregulated DEGs and half the number of downregulated DEGs compared to a wMelCS112 infection. Thus, we demonstrated that Wolbachia made a significant contribution to differential expression of host genes and that the bacterial genotype plays a vital role in establishing the character of this contribution.

Examination of Large Chromosomal Inversions in the Genome of Erwinia amylovora Strains Reveals Worldwide Distribution and North America-Specific Types.

Erwinia amylovora is a relatively homogeneous species with low genetic diversity at the nucleotide level. However, phenotypic differences and genomic structural variations among E. amylovora strains have been documented. In this study, we identified 10 large chromosomal inversion (LCI) types in the Spiraeoideae-infecting (SI) E. amylovora strains by combining whole genome sequencing and PCR-based molecular markers. It was found that LCIs were mainly caused by homologous recombination events among seven rRNA operons (rrns) in SI E. amylovora strains. Although ribotyping results identified inter- and intra-variations in the internal transcribed spacer (ITS1 and ITS2) regions among rrns, LCIs tend to occur between rrns transcribed in the opposite directions and with the same tRNA content (tRNA-Glu or tRNA-Ile/Ala) in ITS1. Based on the LCI types, physical/estimated replichore imbalance (PRI/ERI) was examined and calculated. Among the 117 SI strains evaluated, the LCI types of Ea1189, CFBP1430, and Ea273 were the most common, with ERI values at 1.31, 7.87, and 4.47°, respectively. These three LCI types had worldwide distribution, whereas the remaining seven LCI types were restricted to North America (or certain regions of the United States). Our results indicated ongoing chromosomal recombination events in the SI E. amylovora population and showed that LCI events are mostly symmetrical, keeping the ERI less than 15°. These findings provide initial evidence about the prevalence of certain LCI types in E. amylovora strains, how LCI occurs, and its potential evolutionary advantage and history, which might help track the movement of the pathogen.

Lacustrine speciation associated with chromosomal inversion in a lineage of riverine fishes.

Geographic isolation is the primary driver of speciation in many vertebrate lineages. This trend is exemplified by North American darters, a clade of freshwater fishes where nearly all sister species pairs are allopatric and separated by millions of years of divergence. One of the only exceptions is the Lake Waccamaw endemic Etheostoma perlongum and its riverine sister species Etheostoma maculaticeps, which have no physical barriers to gene flow. Here we show that lacustrine speciation of E. perlongum is characterized by morphological and ecological divergence likely facilitated by a large chromosomal inversion. While E. perlongum is phylogenetically nested within the geographically widespread E. maculaticeps, there is a sharp genetic and morphological break coinciding with the lake-river boundary in the Waccamaw River system. Despite recent divergence, an active hybrid zone, and ongoing gene flow, analyses using a de novo reference genome reveal a 9 Mb chromosomal inversion with elevated divergence between E. perlongum and E. maculaticeps. This region exhibits striking synteny with known inversion supergenes in two distantly related fish lineages, suggesting deep evolutionary convergence of genomic architecture. Our results illustrate that rapid, ecological speciation with gene flow is possible even in lineages where geographic isolation is the dominant mechanism of speciation.

New Wolbachia pipientis Genotype Increasing Heat Stress Resistance of Drosophila melanogaster Host Is Characterized by a Large Chromosomal Inversion.

The maternally transmitted endocellular bacteria Wolbachia is a well-known symbiont of insects, demonstrating both negative and positive effects on host fitness. The previously found Wolbachia strain wMelPlus is characterized by a positive effect on the stress-resistance of its host Drosophila melanogaster, under heat stress conditions. This investigation is dedicated to studying the genomic underpinnings of such an effect. We sequenced two closely related Wolbachia strains, wMelPlus and wMelCS112, assembled their complete genomes, and performed comparative genomic analysis engaging available Wolbachia genomes from the wMel and wMelCS groups. Despite the two strains under study sharing very close gene-composition, we discovered a large (>1/6 of total genome) chromosomal inversion in wMelPlus, spanning through the region that includes the area of the inversion earlier found in the wMel group of Wolbachia genotypes. A number of genes in unique inversion blocks of wMelPlus were identified that might be involved in the induction of a stress-resistant phenotype in the host. We hypothesize that such an inversion could rearrange established genetic regulatory-networks, causing the observed effects of such a complex fly phenotype as a modulation of heat stress resistance. Based on our findings, we propose that wMelPlus be distinguished as a separate genotype of the wMelCS group, named wMelCS3.

A large chromosomal inversion affects antimicrobial sensitivity of Escherichia coli to sodium deoxycholate

Resistance to antimicrobials is normally caused by mutations in the drug targets or genes involved in antimicrobial activation or expulsion. Here we show that an Escherichia coli strain, named DOC14, selected for increased resistance to the bile salt sodium deoxycholate, has no mutations in any ORF, but instead has a 2.1 Mb chromosomal inversion. The breakpoints of the inversion are two inverted copies of an IS5 element. Besides lowering deoxycholate susceptibility, the IS5-mediated chromosomal inversion in the DOC14 mutant was found to increase bacterial survival upon exposure to ampicillin and vancomycin, and sensitize the cell to ciprofloxacin and meropenem, but does not affect bacterial growth or cell morphology in a rich medium in the absence of antibacterial molecules. Overall, our findings support the notion that a large chromosomal inversion can benefit bacterial cells under certain conditions, contributing to genetic variability available for selection during evolution. The DOC14 mutant paired with its isogenic parental strain form a useful model as bacterial ancestors in evolution experiments to study how a large chromosomal inversion influences the evolutionary trajectory in response to various environmental stressors.

An allozyme polymorphism is associated with a large chromosomal inversion in the marine snail Littorina fabalis.

Understanding the genetic targets of natural selection is one of the most challenging goals of population genetics. Some of the earliest candidate genes were identified from associations between allozyme allele frequencies and environmental variation. One such example is the clinal polymorphism in the arginine kinase (Ak) gene in the marine snail Littorina fabalis. While other enzyme loci do not show differences in allozyme frequencies among populations, the Ak alleles are near differential fixation across repeated wave exposure gradients in Europe. Here, we use this case to illustrate how a new sequencing toolbox can be employed to characterize the genomic architecture associated with historical candidate genes. We found that the Ak alleles differ by nine nonsynonymous substitutions, which perfectly explain the different migration patterns of the allozymes during electrophoresis. Moreover, by exploring the genomic context of the Ak gene, we found that the three main Ak alleles are located on different arrangements of a putative chromosomal inversion that reaches near fixation at the opposing ends of two transects covering a wave exposure gradient. This shows Ak is part of a large (3/4 of the chromosome) genomic block of differentiation, in which Ak is unlikely to be the only target of divergent selection. Nevertheless, the nonsynonymous substitutions among Ak alleles and the complete association of one allele with one inversion arrangement suggest that the Ak gene is a strong candidate to contribute to the adaptive significance of the inversion.

Complete Genome Resources for Xylella fastidiosa Strains AlmaEM3 and BB08-1 Reveal Prophage-Associated Structural Variation Among Blueberry-Infecting Strains.

Xylella fastidiosa is a gram-negative plant pathogenic bacterium with wide geographical distribution and host range. X. fastidiosa strains are separated into genetically distinct subspecies, and further categorized into sequence types (ST). Genetic characterization of X. fastidiosa strains infecting blueberry has revealed that strains of subspecies multiplex and fastidiosa are capable of causing bacterial leaf scorch disease of blueberry under field conditions. To better elucidate the relationships among blueberry-infecting X. fastidiosa subsp. multiplex strains, we completed the genomes of an ST 42 strain, AlmaEM3, and a ST 43 strain, BB08-1, using a hybrid assembly approach. Comparison of these assemblies reveals a large (0.95 Mb) chromosomal inversion in BB08-1 relative to AlmaEM3 and the reference strain M12, likely resulting from recombination between prophage regions.

A large chromosomal inversion shapes gene expression in seaweed flies (Coelopa frigida).

Inversions often underlie complex adaptive traits, but the genic targets inside them are largely unknown. Gene expression profiling provides a powerful way to link inversions with their phenotypic consequences. We examined the effects of the Cf‐Inv(1) inversion in the seaweed fly Coelopa frigida on gene expression variation across sexes and life stages. Our analyses revealed that Cf‐Inv(1) shapes global expression patterns, most likely via linked variation, but the extent of this effect is variable, with much stronger effects in adults than larvae. Furthermore, within adults, both common as well as sex‐specific patterns were found. The vast majority of these differentially expressed genes mapped to Cf‐Inv(1). However, genes that were differentially expressed in a single context (i.e., in males, females, or larvae) were more likely to be located outside of Cf‐Inv(1). By combining our findings with genomic scans for environmentally associated SNPs, we were able to pinpoint candidate variants in the inversion that may underlie mechanistic pathways that determine phenotypes. Together the results of this study, combined with previous findings, support the notion that the polymorphic Cf‐Inv(1) inversion in this species is a major factor shaping both coding and regulatory variation resulting in highly complex adaptive effects.

Gene regulatory effects of a large chromosomal inversion in highland maize.

Chromosomal inversions play an important role in local adaptation. Inversions can capture multiple locally adaptive functional variants in a linked block by repressing recombination. However, this recombination suppression makes it difficult to identify the genetic mechanisms underlying an inversion’s role in adaptation. In this study, we used large-scale transcriptomic data to dissect the functional importance of a 13 Mb inversion locus (Inv4m) found almost exclusively in highland populations of maize (Zea mays ssp. mays). Inv4m was introgressed into highland maize from the wild relative Zea mays ssp. mexicana, also present in the highlands of Mexico, and is thought to be important for the adaptation of these populations to cultivation in highland environments. However, the specific genetic variants and traits that underlie this adaptation are not known. We created two families segregating for the standard and inverted haplotypes of Inv4m in a common genetic background and measured gene expression effects associated with the inversion across 9 tissues in two experimental conditions. With these data, we quantified both the global transcriptomic effects of the highland Inv4m haplotype, and the local cis-regulatory variation present within the locus. We found diverse physiological effects of Inv4m across the 9 tissues, including a strong effect on the expression of genes involved in photosynthesis and chloroplast physiology. Although we could not confidently identify the causal alleles within Inv4m, this research accelerates progress towards understanding this inversion and will guide future research on these important genomic features.

\u0391 de novo 3.8-Mb inversion affecting the EDA and XIST genes in a heterozygous female calf with generalized hypohidrotic ectodermal dysplasia

BackgroundIn mammals, hypohidrotic ectodermal dysplasia (HED) is a genetic disorder that is characterized by sparse hair, tooth abnormalities, and defects in cutaneous glands. Only four genes, EDA, EDAR, EDARADD and WNT10A account for more than 90% of HED cases, and EDA, on chromosome X, is involved in 50% of the cases. In this study, we explored an isolated case of a female Holstein calf with symptoms similar to HED.ResultsClinical examination confirmed the diagnosis. The affected female showed homogeneous hypotrichosis and oligodontia as previously observed in bovine EDAR homozygous and EDA hemizygous mutants. Under light microscopy, the hair follicles were thinner and located higher in the dermis of the frontal skin in the affected animal than in the control. Moreover, the affected animal showed a five-fold increase in the number of hair follicles and a four-fold decrease in the diameter of the pilary canals. Pedigree analysis revealed that the coefficient of inbreeding of the affected calf (4.58%) was not higher than the average population inbreeding coefficient (4.59%). This animal had ten ancestors in its paternal and maternal lineages. By estimating the number of affected cases that would be expected if any of these common ancestors carried a recessive mutation, we concluded that, if they existed, other cases of HED should have been reported in France, which is not the case. Therefore, we assumed that the causal mutation was dominant and de novo. By analyzing whole-genome sequencing data, we identified a large chromosomal inversion with breakpoints located in the first introns of the EDA and XIST genes. Genotyping by PCR-electrophoresis the case and its parents allowed us to demonstrate the de novo origin of this inversion. Finally, using various sources of information we present a body of evidence that supports the hypothesis that this mutation is responsible for a skewed inactivation of X, and that only the normal X can be inactivated.ConclusionsIn this article, we report a unique case of X-linked HED affected Holstein female calf with an assumed full inactivation of the normal X-chromosome, thus leading to a severe phenotype similar to that of hemizygous males.

Progress for the pea of prosperity.

Progress for the pea of prosperity.

Dissecting the role of a large chromosomal inversion in life history divergence throughout the Mimulus guttatus species complex.

Chromosomal inversions can play an important role in adaptation, but the mechanism of their action in many natural populations remains unclear. An inversion could suppress recombination between locally beneficial alleles, thereby preventing maladaptive reshuffling with less-fit, migrant alleles. The recombination suppression hypothesis has gained much theoretical support but empirical tests are lacking. Here, we evaluated the evolutionary history and phenotypic effects of a chromosomal inversion which differentiates annual and perennial forms of Mimulus guttatus. We found that perennials likely possess the derived orientation of the inversion. In addition, this perennial orientation occurs in a second perennial species, M.decorus, where it is strongly associated with life history differences between co-occurring M.decorus and annual M.guttatus. One prediction of the recombination suppression hypothesis is that loci contributing to local adaptation will predate the inversion. To test whether the loci influencing perenniality pre-date this inversion, we mapped QTLs for life history traits that differ between annual M.guttatus and a more distantly related, collinear perennial species, M.tilingii. Consistent with the recombination suppression hypothesis, we found that this region is associated with life history in the absence of the inversion, and this association can be broken into at least two QTLs. However, the absolute phenotypic effect of the LG8 inversion region on life history is weaker in M.tilingii than in perennials which possess the inversion. Thus, while we find support for the recombination suppression hypothesis, the contribution of this inversion to life history divergence in this group is likely complex.

Whole-genome mapping reveals a large chromosomal inversion on Iberian Brucella suis biovar 2 strains

Optical mapping is a technology able to quickly generate high resolution ordered whole-genome restriction maps of bacteria, being a proven approach to search for diversity among bacterial isolates. In this work, optical whole-genome maps were used to compare closely-related Brucella suis biovar 2 strains. This biovar is the unique isolated in domestic pigs and wild boars in Portugal and Spain and most of the strains share specific molecular characteristics establishing an Iberian clonal lineage that can be differentiated from another lineage mainly isolated in several Central European countries. We performed the BamHI whole-genome optical maps of five B. suis biovar 2 field strains, isolated from wild boars in Portugal and Spain (three from the Iberian lineage and two from the Central European one) as well as of the reference strain B. suis biovar 2 ATCC 23445 (Central European lineage, Denmark). Each strain showed a distinct, highly individual configuration of 228–231 BamHI fragments. Nevertheless, a low divergence was globally observed in chromosome II (1.6%) relatively to chromosome I (2.4%). Optical mapping also disclosed genomic events associated with B. suis strains in chromosome I, namely one indel (3.5kb) and one large inversion (944kb). By using targeted-PCR in a set of 176 B. suis strains, including all biovars and haplotypes, the indel was found to be specific of the reference strain ATCC 23445 and the large inversion was shown to be an exclusive genomic marker of the Iberian clonal lineage of biovar 2.

Adaptive divergence in the monkey flower Mimulus guttatus is maintained by a chromosomal inversion.

Organisms exhibit an incredible diversity of life history strategies as adaptive responses to environmental variation. The establishment of novel life history strategies involves multilocus polymorphisms, which will be challenging to establish in the face of gene flow and recombination. Theory predicts that adaptive allelic combinations may be maintained and spread if they occur in genomic regions of reduced recombination, such as chromosomal inversion polymorphisms, yet empirical support for this prediction is lacking. Here, we use genomic data to investigate the evolution of divergent adaptive ecotypes of the yellow monkey flower Mimulus guttatus. We show that a large chromosomal inversion polymorphism is the major region of divergence between geographically widespread annual and perennial ecotypes. In contrast, ∼40,000 single nucleotide polymorphisms in collinear regions of the genome show no signal of life history, revealing genomic patterns of diversity have been shaped by localized homogenizing gene flow and large‐scale Pleistocene range expansion. Our results provide evidence for an inversion capturing and protecting loci involved in local adaptation, while also explaining how adaptive divergence can occur with gene flow.

Highly efficient multiplex targeted mutagenesis and genomic structure variation in Bombyx mori cells using CRISPR/Cas9

Bombyx mori is an economically important insect and a model organism for studying lepidopteran and arthropod biology. Using a highly efficient CRISPR/Cas9 system, we showed that this system could mediate highly efficient targeted genome editing of a single gene locus, large chromosomal deletion or inversion, and also multiplex genome editing of 6 genes simultaneously in BmNs cell line derived from B. mori. The simplicity and high efficiency of our system provide unprecedented possibilities for researchers to implement precise and sophisticated manipulation of a chosen B. mori gene in BmNs cells easily in a limited time course, and perhaps new opportunities for functional genomics of B. mori and other lepidopteran insects.

Male-specific genotype by environment interactions influence viability selection acting on a sexually selected inversion system in the seaweed fly, Coelopa frigida.

In the seaweed fly, Coelopa frigida, a large chromosomal inversion system is affected by sexual selection and viability selection. However, our understanding of the interaction between these two selective forces is currently limited as research has focused upon a limited range of environments. We allowed C. frigida larvae to develop in two different algae, Fucus and Laminaria, and then measured viability and body size for each inversion genotype. Significant male-specific genotype-by-environment interactions influenced viability and body size. For males developing in Laminaria, the direction of viability selection acts similarly on the inversion system as the direction of sexual selection. In contrast, for males developing in Fucus, viability selection opposes sexual selection. These results demonstrate that through considering viability selection in different environments, the costs and benefits associated with sexual selection can be found to vary.

An attenuated strain of Bacillus anthracis (CDC 684) has a large chromosomal inversion and altered growth kinetics

BackgroundAn isolate originally labeled Bacillus megaterium CDC 684 was found to contain both pXO1 and pXO2, was non-hemolytic, sensitive to gamma-phage, and produced both the protective antigen and the poly-D-glutamic acid capsule. These phenotypes prompted Ezzell et al., (J. Clin. Microbiol. 28:223) to reclassify this isolate to Bacillus anthracis in 1990.ResultsWe demonstrate that despite these B. anthracis features, the isolate is severely attenuated in a guinea pig model. This prompted whole genome sequencing and closure. The comparative analysis of CDC 684 to other sequenced B. anthracis isolates and further analysis reveals: a) CDC 684 is a close relative of a virulent strain, Vollum A0488; b) CDC 684 defines a new B. anthracis lineage (at least 51 SNPs) that includes 15 other isolates; c) the genome of CDC 684 contains a large chromosomal inversion that spans 3.3 Mbp; d) this inversion has caused a displacement of the usual spatial orientation of the origin of replication (ori) to the termination of replication (ter) from 180° in wild-type B. anthracis to 120° in CDC 684 and e) this isolate also has altered growth kinetics in liquid media.ConclusionsWe propose two alternative hypotheses explaining the attenuated phenotype of this isolate. Hypothesis 1 suggests that the skewed ori/ter relationship in CDC 684 has altered its DNA replication and/or transcriptome processes resulting in altered growth kinetics and virulence capacity. Hypothesis 2 suggests that one or more of the single nucleotide polymorphisms in CDC 684 has altered the expression of a regulatory element or other genes necessary for virulence.

Genome-wide detection of chromosomal rearrangements, indels, and mutations in circular chromosomes by short read sequencing

Whole-genome sequencing (WGS) with new short-read sequencing technologies has recently been applied for genome-wide identification of mutations. Genomic rearrangements have, however, often remained undetected by WGS, and additional analyses are required for their detection. Here, we have applied a combination of WGS and genome copy number analysis, for the identification of mutations that suppress the growth deficiency imposed by excessive initiations from the Escherichia coli origin of replication, oriC. The E. coli chromosome, like the majority of bacterial chromosomes, is circular, and DNA replication is initiated by assembling two replication complexes at the origin, oriC. These complexes then replicate the chromosome bidirectionally toward the terminus, ter. In a population of growing cells, this results in a copy number gradient, so that origin-proximal sequences are more frequent than origin-distal sequences. Major rearrangements in the chromosome are, therefore, readily identified by changes in copy number, i.e., certain sequences become over- or under-represented. Of the eight mutations analyzed in detail here, six were found to affect a single gene only, one was a large chromosomal inversion, and one was a large chromosomal duplication. The latter two mutations could not be detected solely by WGS, validating the present approach for identification of genomic rearrangements. We further suggest the use of copy number analysis in combination with WGS for validation of newly assembled bacterial chromosomes.

A large inversion in the linear chromosome of Streptomyces griseus caused by replicative transposition of a new Tn3 family transposon

We have comprehensively analyzed the linear chromosomes of Streptomyces griseus mutants constructed and kept in our laboratory. During this study, macrorestriction analysis of AseI and DraI fragments of mutant 402-2 suggested a large chromosomal inversion. The junctions of chromosomal inversion were cloned and sequenced and compared with the corresponding target sequences in the parent strain 2247. Consequently, a transposon-involved mechanism was revealed. Namely, a transposon originally located at the left target site was replicatively transposed to the right target site in an inverted direction, which generated a second copy and at the same time caused a 2.5-Mb chromosomal inversion. The involved transposon named TnSGR was grouped into a new subfamily of the resolvase-encoding Tn3 family transposons based on its gene organization. At the end, terminal diversity of S. griseus chromosomes is discussed by comparing the sequences of strains 2247 and IFO13350.