Differences In Genomic Architecture Research Articles

Nuclear dualism without extensive DNA elimination in the ciliate Loxodes magnus

Most eukaryotes have one nucleus and nuclear genome per cell. Ciliates have instead evolved distinct nuclei that coexist in each cell: a silent germline vs. transcriptionally active somatic nuclei. In the best-studied model species, both nuclei can divide asexually, but only germline nuclei undergo meiosis and karyogamy during sex. Thereafter, thousands of DNA segments, called internally eliminated sequences (IESs), are excised from copies of the germline genomes to produce the streamlined somatic genome. In Loxodes, however, somatic nuclei cannot divide but instead develop from germline copies even during asexual cell division, which would incur a huge overhead cost if genome editing was required. Here, we purified and sequenced both genomes in Loxodes magnus to see whether their nondividing somatic nuclei are associated with differences in genome architecture. Unlike in other ciliates studied to date, we did not find canonical germline-limited IESs, implying Loxodes does not extensively edit its genomes. Instead, both genomes appear large and equivalent, replete with retrotransposons and repetitive sequences, unlike the compact, gene-rich somatic genomes of other ciliates. Two other hallmarks of nuclear development in ciliates-domesticated DDE-family transposases and editing-associated small RNAs-were also not found. Thus, among the ciliates, Loxodes genomes most resemble those of conventional eukaryotes. Nonetheless, base modifications, histone marks, and nucleosome positioning of vegetative Loxodes nuclei are consistent with functional differentiation between actively transcribed somatic vs. inactive germline nuclei. Given their phylogenetic position, it is likely that editing was present in the ancestral ciliate but secondarily lost in the Loxodes lineage.

Comparative analyses of the banded alder borer (Rosalia funebris) and Asian longhorned beetle (Anoplophora glabripennis) genomes reveal significant differences in genome architecture and gene content among these and other Cerambycidae.

Rosalia funebris (RFUNE; Cerambycidae), the banded alder borer, is a longhorn beetle whose larvae feed on the wood of various economically and ecologically significant trees in western North America. Adults are short-lived and not known to consume plant material substantially. We sequenced, assembled, and annotated the RFUNE genome using HiFi and RNASeq data. We documented genome architecture and gene content, focusing on genes putatively involved in plant feeding (phytophagy). Comparisons were made to the well-studied genome of the Asian longhorned beetle (AGLAB; Anoplophora glabripennis) and other Cerambycidae. The 814 Mb RFUNE genome assembly was distributed across 42 contigs, with an N50 of 30.18 Mb. Repetitive sequences comprised 60.27% of the genome, and 99.0% of expected single-copy orthologous genes were fully assembled. We identified 12,657 genes, fewer than in the four other species studied, and 46.4% fewer than for Aromia moschata (same subfamily as RFUNE). Of the 7,258 orthogroups shared between RFUNE and AGLAB, 1,461 had more copies in AGLAB and 1,023 had more copies in RFUNE. We identified 240 genes in RFUNE that putatively arose via horizontal transfer events. The RFUNE genome encoded substantially fewer putative plant cell wall degrading enzymes than AGLAB, which may relate to the longer-lived plant-feeding adults of the latter species. The RFUNE genome provides new insights into cerambycid genome architecture and gene content and provides a new vantage point from which to study the evolution and genomic basis of phytophagy in beetles.

Multifaceted nature of young-onset diabetes - can genomic medicine improve the precision of diagnosis and management?

Young-onset type 2 diabetes (YOD), defined as diabetes diagnosis before age 40, has an aggressive clinical course with premature mortality, in part due to long disease duration and lack of evidence to guide diagnosis and management. Autoimmune type 1 diabetes, maturity-onset diabetes of the young (MODY), and latent autoimmune diabetes in adults (LADA) are subtypes of diabetes in young people, which, however, cannot fully explain their complex clinical course. Similarly, family members carrying the same rare genetic variant of monogenic diabetes can have different presentations and outcomes. Ancestral heterogeneity, ecological transition, inter-ethnic differences in genomic architecture, and variations in living environment, lifestyles, access to care, and timeliness of diagnosis and treatment can influence the age of diagnosis and exposure to these cardiometabolic-renal risk factors. Despite the wealth of literature on genetic associations with diabetes, the familial cosegregation of rare variants and their relevance to YOD remains uncertain. This perspective was motivated by decades of clinical observations and learnings from an ongoing randomized controlled trial that uses biogenetic markers to classify patients with YOD for improving outcomes. Apart from highlighting the need to use family-based studies to improve the precision of diagnosis, we discussed atypical causes for diabetic ketoacidosis and the importance of lifecourse and psychosocial-behavioral factors in patients with YOD. Apart from detailed clinical evaluation, we propose using plasma C peptide, homeostasis model of assessment (HOMA) indexes, autoantibodies, and polygenic risk scores to stratify risk, classify diabetes subtypes, and personalize treatment in YOD. To achieve these goals, we advocate changing the practice environment and team structure to enable physicians to use the insights they learn from patients and their family members to implement precision medicine and improve the outlook of these high-risk individuals.

Genomic architecture and functional unit of mimicry supergene in female limited Batesian mimic Papilio butterflies.

It has long been suggested that dimorphic female-limited Batesian mimicry of two closely related Papilio butterflies, Papilio memnon and Papilio polytes, is controlled by supergenes. Whole-genome sequencing, genome-wide association studies and functional analyses have recently identified mimicry supergenes, including the doublesex (dsx) gene. Although supergenes of both the species are composed of highly divergent regions between mimetic and non-mimetic alleles and are located at the same chromosomal locus, they show critical differences in genomic architecture, particularly with or without an inversion: P. polytes has an inversion, but P. memnon does not. This review introduces and compares the detailed genomic structure of mimicry supergenes in two Papilio species, including gene composition, repetitive sequence composition, breakpoint/boundary site structure, chromosomal inversion and linkage disequilibrium. Expression patterns and functional analyses of the respective genes within or flanking the supergene suggest that dsx and other genes are involved in mimetic traits. In addition, structural comparison of the corresponding region for the mimicry supergene among further Papilio species suggests three scenarios for the evolution of the mimicry supergene between the two Papilio species. The structural features revealed in the Papilio mimicry supergene provide insight into the formation, maintenance and evolution of supergenes. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

A chromosomal inversion may facilitate adaptation despite periodic gene flow in a freshwater fish.

Differences in genomic architecture between populations, such as chromosomal inversions, may play an important role in facilitating adaptation despite opportunities for gene flow. One system where chromosomal inversions may be important for eco‐evolutionary dynamics is in freshwater fishes, which often live in heterogenous environments characterized by varying levels of connectivity and varying opportunities for gene flow. In the present study, reduced representation sequencing was used to study possible adaptation in n = 345 walleye (Sander vitreus) from three North American waterbodies: Cedar Bluff Reservoir (Kansas, USA), Lake Manitoba (Manitoba, Canada), and Lake Winnipeg (Manitoba, Canada). Haplotype and outlier‐based tests revealed a putative chromosomal inversion that contained three expressed genes and was nearly fixed in walleye assigned to Lake Winnipeg. These patterns exist despite the potential for high gene flow between these proximate Canadian lakes, suggesting that the inversion may be important for facilitating adaptive divergence between the two lakes despite gene flow. However, a specific adaptive role for the putative inversion could not be tested with the present data. Our study illuminates the importance of genomic architecture consistent with local adaptation in freshwater fishes. Furthermore, our results provide additional evidence that inversions may facilitate local adaptation in many organisms that inhabit connected but heterogenous environments.

The Genomic Architecture and Evolutionary Fates of Supergenes.

Supergenes are genomic regions containing sets of tightly linked loci that control multi-trait phenotypic polymorphisms under balancing selection. Recent advances in genomics have uncovered significant variation in both the genomic architecture as well as the mode of origin of supergenes across diverse organismal systems. Although the role of genomic architecture for the origin of supergenes has been much discussed, differences in the genomic architecture also subsequently affect the evolutionary trajectory of supergenes and the rate of degeneration of supergene haplotypes. In this review, we synthesize recent genomic work and historical models of supergene evolution, highlighting how the genomic architecture of supergenes affects their evolutionary fate. We discuss how recent findings on classic supergenes involved in governing ant colony social form, mimicry in butterflies, and heterostyly in flowering plants relate to theoretical expectations. Furthermore, we use forward simulations to demonstrate that differences in genomic architecture affect the degeneration of supergenes. Finally, we discuss implications of the evolution of supergene haplotypes for the long-term fate of balanced polymorphisms governed by supergenes.

Structural variation and its potential impact on genome instability: Novel discoveries in the EGFR landscape by long-read sequencing.

Structural variation (SV) is typically defined as variation within the human genome that exceeds 50 base pairs (bp). SV may be copy number neutral or it may involve duplications, deletions, and complex rearrangements. Recent studies have shown SV to be associated with many human diseases. However, studies of SV have been challenging due to technological constraints. With the advent of third generation (long-read) sequencing technology, exploration of longer stretches of DNA not easily examined previously has been made possible. In the present study, we utilized third generation (long-read) sequencing techniques to examine SV in the EGFR landscape of four haplotypes derived from two human samples. We analyzed the EGFR gene and its landscape (+/- 500,000 base pairs) using this approach and were able to identify a region of non-coding DNA with over 90% similarity to the most common activating EGFR mutation in non-small cell lung cancer. Based on previously published Alu-element genome instability algorithms, we propose a molecular mechanism to explain how this non-coding region of DNA may be interacting with and impacting the stability of the EGFR gene and potentially generating this cancer-driver gene. By these techniques, we were also able to identify previously hidden structural variation in the four haplotypes and in the human reference genome (hg38). We applied previously published algorithms to compare the relative stabilities of these five different EGFR gene landscape haplotypes to estimate their relative potentials to generate the EGFR exon 19, 15 bp canonical deletion. To our knowledge, the present study is the first to use the differences in genomic architecture between targeted cancer-linked phased haplotypes to estimate their relative potentials to form a common cancer-linked driver mutation.

Phylogenetic and genomic analyses of the ribosomal oxygenases Riox1 (No66) and Riox2 (Mina53) provide new insights into their evolution

BackgroundTranslation of specific mRNAs can be highly regulated in different cells, tissues or under pathological conditions. Ribosome heterogeneity can originate from variable expression or post-translational modifications of ribosomal proteins. The ribosomal oxygenases RIOX1 (NO66) and RIOX2 (MINA53) modify ribosomal proteins by histidine hydroxylation. A similar mechanism is present in prokaryotes. Thus, ribosome hydroxylation may be a well-conserved regulatory mechanism with implications in disease and development. However, little is known about the evolutionary history of Riox1 and Riox2 genes and their encoded proteins across eukaryotic taxa.ResultsIn this study, we have analysed Riox1 and Riox2 orthologous genes from 49 metazoen species and have constructed phylogenomic trees for both genes. Our genomic and phylogenetic analyses revealed that Arthropoda, Annelida, Nematoda and Mollusca lack the Riox2 gene, although in the early phylum Cnidaria both genes, Riox1 and Riox2, are present and expressed. Riox1 is an intronless single-exon-gene in several species, including humans. In contrast to Riox2, Riox1 is ubiquitously present throughout the animal kingdom suggesting that Riox1 is the phylogenetically older gene from which Riox2 has evolved. Both proteins have maintained a unique protein architecture with conservation of active sites within the JmjC domains, a dimerization domain, and a winged-helix domain. In addition, Riox1 proteins possess a unique N-terminal extension domain. Immunofluorescence analyses in Hela cells and in Hydra vulgaris identified a nucleolar localisation signal within the extended N-terminal domain of human RIOX1 and an altered subnuclear localisation for the Hydra Riox2.ConclusionsConserved active site residues and uniform protein domain architecture suggest a consistent enzymatic activity within the Riox orthologs throughout evolution. However, differences in genomic architecture, like single exon genes and alterations in subnuclear localisation, as described for Hydra, point towards adaption mechanisms that may correlate with taxa- or species-specific requirements. The diversification of Riox1/Riox2 gene structures throughout evolution suggest that functional requirements in expression of protein isoforms and/or subcellular localisation of proteins may have evolved by adaptation to lifestyle.

Chromosome conformation signatures define predictive markers of inadequate response to methotrexate in early rheumatoid arthritis

BackgroundThere is a pressing need in rheumatoid arthritis (RA) to identify patients who will not respond to first-line disease-modifying anti-rheumatic drugs (DMARD). We explored whether differences in genomic architecture represented by a chromosome conformation signature (CCS) in blood taken from early RA patients before methotrexate (MTX) treatment could assist in identifying non-response to DMARD and, whether there is an association between such a signature and RA specific expression quantitative trait loci (eQTL).MethodsWe looked for the presence of a CCS in blood from early RA patients commencing MTX using chromosome conformation capture by EpiSwitch™. Using blood samples from MTX responders, non-responders and healthy controls, a custom designed biomarker discovery array was refined to a 5-marker CCS that could discriminate between responders and non-responders to MTX. We cross-validated the predictive power of the CCS by generating 150 randomized groups of 59 early RA patients (30 responders and 29 non-responders) before MTX treatment. The CCS was validated using a blinded, independent cohort of 19 early RA patients (9 responders and 10 non-responders). Last, the loci of the CCS markers were mapped to RA-specific eQTL.ResultsWe identified a 5-marker CCS that could identify, at baseline, responders and non-responders to MTX. The CCS consisted of binary chromosome conformations in the genomic regions of IFNAR1, IL-21R, IL-23, CXCL13 and IL-17A. When tested on a cohort of 59 RA patients, the CCS provided a negative predictive value of 90.0% for MTX response. When tested on a blinded independent validation cohort of 19 early RA patients, the signature demonstrated a true negative response rate of 86 and a 90% sensitivity for detection of non-responders to MTX. Only conformations in responders mapped to RA-specific eQTL.ConclusionsHere we demonstrate that detection of a CCS in blood in early RA is able to predict inadequate response to MTX with a high degree of accuracy. Our results provide a proof of principle that a priori stratification of response to MTX is possible, offering a mechanism to provide alternative treatments for non-responders to MTX earlier in the course of the disease.

Prevention of tumorigenesis in mice by exercise is dependent on strain background and timing relative to carcinogen exposure

Among cancer diagnoses, colorectal cancer (CRC) is prevalent, with a lifetime risk of developing CRC being approximately 5%. Population variation surrounding the mean risk of developing CRCs has been associated with both inter-individual differences in genomic architecture and environmental exposures. Decreased risk of CRC has been associated with physical activity, but protective responses are variable. Here, we utilized a series of experiments to examine the effects of genetic background (strain), voluntary exercise (wheel running), and their interaction on azoxymethane (AOM)-induced intestinal tumor number and size in mice. Additionally, we investigated how the timing of exercise relative to AOM exposure, and amount of exercise, affected tumor number and size. Our results indicated that voluntary exercise significantly reduced tumor number in a strain dependent manner. Additionally, among strains where exercise reduced tumor number (A/J, CC0001/Unc) the timing of voluntary exercise relative to AOM exposure was crucial. Voluntary exercise prior to or during AOM treatment resulted in a significant reduction in tumor number, but exercise following AOM exposure had no effect. The results indicate that voluntary exercise should be used as a preventative measure to reduce risk for environmentally induced CRC with the realization that the extent of protection may depend on genetic background.

Genomics of speciation and introgression in Princess cichlid fishes from Lake Tanganyika

How variation in the genome translates into biological diversity and new species originate has endured as the mystery of mysteries in evolutionary biology. African cichlid fishes are prime model systems to address speciation-related questions for their remarkable taxonomic and phenotypic diversity, and the possible role of gene flow in this process. Here, we capitalize on genome sequencing and phylogenomic analyses to address the relative impacts of incomplete lineage sorting, introgression and hybrid speciation in the Neolamprologus savoryi-complex (the 'Princess cichlids') from Lake Tanganyika. We present a time-calibrated species tree based on whole-genome sequences and provide strong evidence for incomplete lineage sorting in the early phases of diversification and multiple introgression events affecting different stages. Importantly, we find that the Neolamprologus chromosomes show centre-to-periphery biases in nucleotide diversity, sequence divergence, GC content, incomplete lineage sorting and rates of introgression, which are likely modulated by recombination density and linked selection. The detection of heterogeneous genomic landscapes has strong implications on the genomic mechanisms involved in speciation. Collinear chromosomal regions can be protected from gene flow and harbour incompatibility genes if they reside in lowly recombining regions, and coupling can evolve between nonphysically linked genomic regions (chromosome centres in particular). Simultaneously, higher recombination towards chromosome peripheries makes these more dynamic, evolvable regions where adaptation polymorphisms have a fertile ground. Hence, differences in genome architecture could explain the levels of taxonomic and phenotypic diversity seen in taxa with collinear genomes and might have contributed to the spectacular cichlid diversity observed today.

The Perchlorate Reduction Genomic Island: Mechanisms and Pathways of Evolution by Horizontal Gene Transfer.

BackgroundPerchlorate is a widely distributed anion that is toxic to humans, but serves as a valuable electron acceptor for several lineages of bacteria. The ability to utilize perchlorate is conferred by a horizontally transferred piece of DNA called the perchlorate reduction genomic island (PRI).MethodsWe compared genomes of perchlorate reducers using phylogenomics, SNP mapping, and differences in genomic architecture to interrogate the evolutionary history of perchlorate respiration.ResultsHere we report on the PRI of 13 genomes of perchlorate-reducing bacteria from four different classes of Phylum Proteobacteria (the Alpha-, Beta-, Gamma- and Epsilonproteobacteria). Among the different phylogenetic classes, the island varies considerably in genetic content as well as in its putative mechanism and location of integration. However, the islands of the densely sampled genera Azospira and Magnetospirillum have striking nucleotide identity despite divergent genomes, implying horizontal transfer and positive selection within narrow phylogenetic taxa. We also assess the phylogenetic origin of accessory genes in the various incarnations of the island, which can be traced to chromosomal paralogs from phylogenetically similar organisms.ConclusionThese observations suggest a complex phylogenetic history where the island is rarely transferred at the class level but undergoes frequent and continuous transfer within narrow phylogenetic groups. This restricted transfer is seen directly by the independent integration of near-identical islands within a genus and indirectly due to the acquisition of lineage-specific accessory genes. The genomic reversibility of perchlorate reduction may present a unique equilibrium for a metabolism that confers a competitive advantage only in the presence of an electron acceptor, which although widely distributed, is generally present at low concentrations in nature.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-2011-5) contains supplementary material, which is available to authorized users.

Differences in genomic architecture between two distinct geographical strains of the blood fluke Schistosoma japonicum reveal potential phenotype basis

The Chinese (SjC) and Philippine (SjP) strains of the blood fluke Schistosoma japonicum have been shown to present clearly different phenotypes in fecundity, pathology, drug sensitivity and immunology. We used microarray based comparative genomic hybridisation (aCGH) to investigate structural differences in the genomes of the two strains and identified seven distinct regions of the S. japonicum genome that present differential aCGH representing either deletion or duplication regions in SjP. Within these regions, genes predicted to be associated with the recognised phenotypic differences were identified and that may provide new insights into the biology and evolution of the two strains, with implications for the epidemiology and control of schistosomiasis japonica in China and the Philippines.

Using 50 k single nucleotide polymorphisms to elucidate genomic architecture of line 1 hereford cattle.

Hereford is a major beef breed in the USA, and a sub-population, known as Line 1 (L1), was established in 1934 using two paternal half-sib bulls and 50 unrelated females. L1 has since been maintained as a closed population and selected for growth to 1 year of age. Objectives were to characterize the molecular genetic architecture of L1 (n = 240) by comparing a cross-section of L1 with the general US. Hereford population (AHA, n = 311), estimating effects of imposed selection within L1 based on allele frequencies at 50 K SNP loci, and examining loci-specific effects of heterozygosity on the selection criterion. Animals were genotyped using the Illumina BovineSNP50 Beadchip, and SNP were mapped to UMD3.0 assembly of the bovine genome sequence. Average linkage disequilibrium (LD), measured by square of Pearson correlation, of adjacent SNP was 0.36 and 0.16 in L1 and AHA, respectively. Difference in LD between L1 and AHA decreased as SNP spacing increased. Persistence of phase between L1 and AHA decreased from 0.45 to 0.14 as SNP spacing increased from 50 to 5,000 kb. Extended haplotype homozygosity was greater in L1 than in AHA for 95.6% of the SNP. Knowledge of selection applied to L1 facilitated a novel approach to QTL discovery. Minor allele frequency was (FDR < 0.01) affected by cumulative selection differential at 191 out of 25,901 SNP. With the FDR relaxed to 0.05, 13 regions on BTA2, 5, 6, 9, 11, 14, 15, 18, 23, and 26 are co-located with previously identified QTL for growth. After adjustment of postweaning gain phenotypes for fixed effects and direct additive genetic effects, regression of residuals on genome-wide heterozygosity was −235.3 ± 91.6 kg. However, no SNP-specific loci where heterozygotes were significantly superior to the average of homozygotes were revealed (FDR ≥ 0.17). In conclusion, genome-wide SNP genotypes clarified effects of selection and inbreeding within L1 and differences in genomic architecture between the population segment L1 and the AHA population.

Mutation spectrum of Drosophila CNVs revealed by breakpoint sequencing

BackgroundThe detailed study of breakpoints associated with copy number variants (CNVs) can elucidate the mutational mechanisms that generate them and the comparison of breakpoints across species can highlight differences in genomic architecture that may lead to lineage-specific differences in patterns of CNVs. Here, we provide a detailed analysis of Drosophila CNV breakpoints and contrast it with similar analyses recently carried out for the human genome.ResultsBy applying split-read methods to a total of 10x coverage of 454 shotgun sequence across nine lines of D. melanogaster and by re-examining a previously published dataset of CNVs detected using tiling arrays, we identified the precise breakpoints of more than 600 insertions, deletions, and duplications. Contrasting these CNVs with those found in humans showed that in both taxa CNV breakpoints fall into three classes: blunt breakpoints; simple breakpoints associated with microhomology; and breakpoints with additional nucleotides inserted/deleted and no microhomology. In both taxa CNV breakpoints are enriched with non-B DNA sequence structures, which may impair DNA replication and/or repair. However, in contrast to human genomes, non-allelic homologous-recombination (NAHR) plays a negligible role in CNV formation in Drosophila. In flies, non-homologous repair mechanisms are responsible for simple, recurrent, and complex CNVs, including insertions of de novo sequence as large as 60 bp.ConclusionsHumans and Drosophila differ considerably in the importance of homology-based mechanisms for the formation of CNVs, likely as a consequence of the differences in the abundance and distribution of both segmental duplications and transposable elements between the two genomes.

Sex-specific population structure, natural selection, and linkage disequilibrium in a wild bird population as revealed by genome-wide microsatellite analyses

BackgroundSexual dimorphism in ecologically important traits is widespread, yet the differences in the genomic architecture between the two sexes are largely unexplored. We employed a genome-wide multilocus approach to examine the sexual differences in population subdivision, natural selection and linkage disequilibrium (LD) in a wild Siberian jay (Perisoreus infaustus) population, using genotypes at a total of 107 autosomal and Z-chromosomal microsatellites.ResultsMean observed heterozygosity was significantly higher in females (HO = 0.567) than in males (HO = 0.532), and autosomal markers (HO = 0.561) were more variable than Z-chromosomal markers (HO = 0.512). Genetic differentiation (FST = 0.002, P < 0.05) between the two sexes was low but significant and males were on average significantly more genetically related to each other than females. Genomescan analyses revealed that 3 out of 101 (3%) autosomal loci were under directional selection, while 4 out of 6 (67%) Z-chromosomal markers were indicated to be under balancing selection. This suggests a significantly greater but contrasting selection force on the Z-chromosome in comparison to autosomes, which is consistent with an overall significantly (P < 0.05) lower FSTvalue for Z-chromosomal (-0.014, 95% CI: -0.025 - -0.011) than for the autosomal loci (0.003, 95% CI: 0.001 - 0.004). Analysis of syntenic marker pairs revealed high levels of LD in both sexes but significantly (P < 0.05) lower levels of LD in the females both on autosomes and Z-chromosome, probably due to the higher rate of dispersal and the higher recombination rates on autosomes, as well as the pseudoautosomal markers. In both sexes LD decayed rapidly with genetic distance in a similar fashion on autosomes, while a more rapid decay of LD in Z-chromosome was detected in females than in males.ConclusionWe conclude that there are many clear differences in genomic architecture between the sexes studied here which can be at least partly understood in the light of higher dispersal rate of females as compared to males and the unusual structure of the Z-chromosome of the species.