Extensive Genomic Rearrangements Research Articles

The chromosome-level genome of dragon fruit reveals whole-genome duplication and chromosomal co-localization of betacyanin biosynthetic genes

Dragon fruits are tropical fruits economically important for agricultural industries. As members of the family of Cactaceae, they have evolved to adapt to the arid environment. Here we report the draft genome of Hylocereus undatus, commercially known as the white-fleshed dragon fruit. The chromosomal level genome assembly contains 11 longest scaffolds corresponding to the 11 chromosomes of H. undatus. Genome annotation of H. undatus found ~29,000 protein-coding genes, similar to Carnegiea gigantea (saguaro). Whole-genome duplication (WGD) analysis revealed a WGD event in the last common ancestor of Cactaceae followed by extensive genome rearrangements. The divergence time between H. undatus and C. gigantea was estimated to be 9.18 MYA. Functional enrichment analysis of orthologous gene clusters (OGCs) in six Cactaceae plants found significantly enriched OGCs in drought resistance. Fruit flavor-related functions were overrepresented in OGCs that are significantly expanded in H. undatus. The H. undatus draft genome also enabled the discovery of carbohydrate and plant cell wall-related functional enrichment in dragon fruits treated with trypsin for a longer storage time. Lastly, genes of the betacyanin (a red-violet pigment and antioxidant with a very high concentration in dragon fruits) biosynthetic pathway were found to be co-localized on a 12 Mb region of one chromosome. The consequence may be a higher efficiency of betacyanin biosynthesis, which will need experimental validation in the future. The H. undatus draft genome will be a great resource to study various cactus plants.

Chromothripsis in Human Breast Cancer.

Chromothripsis is a form of genome instability by which a presumably single catastrophic event generates extensive genomic rearrangements of one or a few chromosomes. Widely assumed to be an early event in tumor development, this phenomenon plays a prominent role in tumor onset. In this study, an analysis of chromothripsis in 252 human breast cancers from two patient cohorts (149 metastatic breast cancers, 63 untreated primary tumors, 29 local relapses, and 11 longitudinal pairs) using whole-genome and whole-exome sequencing reveals that chromothripsis affects a substantial proportion of human breast cancers, with a prevalence over 60% in a cohort of metastatic cases and 25% in a cohort comprising predominantly luminal breast cancers. In the vast majority of cases, multiple chromosomes per tumor were affected, with most chromothriptic events on chromosomes 11 and 17 including, among other significantly altered drivers, CCND1, ERBB2, CDK12, and BRCA1. Importantly, chromothripsis generated recurrent fusions that drove tumor development. Chromothripsis-related rearrangements were linked with univocal mutational signatures, with clusters of point mutations due to kataegis in close proximity to the genomic breakpoints and with the activation of specific signaling pathways. Analyzing the temporal order of events in tumors with and without chromothripsis as well as longitudinal analysis of chromothriptic patterns in tumor pairs offered important insights into the role of chromothriptic chromosomes in tumor evolution. SIGNIFICANCE: These findings identify chromothripsis as a major driving event in human breast cancer.

Somatic Genomic Changes in the Formation of Differentiated Thyroid Carcinoma

Abstract Globally, the diffuse goiter affects more than 10% of the population and in some regions is endemic. Thyroid nodules are found in approximately 5% of the population using the oldest method for thyroid examination – palpation. When performing ultrasound screening, this percentage increases significantly and reaches between 20 and 75% of the total population. Thyroid carcinoma is a rare malignancy and accounts for up to 1% of all malignant tumors. It is the most common endocrine cancer and is clinically manifested as a thyroid nodule. Somatic mutations play an important role in its development. Differentiation of benign and malignant thyroid nodules is of great importance due to the different therapeutic approach. Therefore, new diagnostic tools are sought to help distinguish the two. Despite the progress in our knowledge of carcinogenesis in recent years, a number of key issues still remain unanswered. The establishment of new rare somatic mutations can improve pre-surgical diagnosis and optimize post-operative strategies for the treatment of thyroid carcinoma. Next-generation sequencing (NGS) allows for extensive mutation and genome rearrangements tracking. The results obtained with NGS provide the basis for the development of new approach for systematic genetic screening, at prevention, early diagnosis, accurate prognosis, and targeted therapy of this disorder.

Whole-genome sequence of Haemophilus influenzae ST422 outbreak clone strain 2018-Y40 with low quinolone susceptibility isolated from a paediatric patient

ObjectivesIn recent years, Haemophilus influenzae strains with reduced susceptibility to quinolones have emerged and spread in Japan. In addition, an outbreak of isolates with low quinolone susceptibility among paediatric patients has also been reported. The aim of this study was to determine the molecular characteristics of an H. influenzae ST422 outbreak clone with low quinolone susceptibility isolated from a paediatric patient using whole-genome sequencing. MethodsThe PacBio RS II platform was used for sequencing, and de novo assembly was performed using RS HGAP assembly version 3.0. The assembled sequences were annotated using DFAST version 1.1.15. Prophages were estimated using the PHASTER program. ResultsWhole-genome sequencing of H. influenzae ST422 isolate 2018-Y40 revealed that the genome size was 1 957 393bp, comprising 1 926 protein-coding sequences, 19 rRNAs, and 57 tRNAs, with a guanine–cytosine (GC) content of 38.2%. This isolate had no relevant exogenous antimicrobial-resistant genes. However, amino acid substitutions were found in both GyrA and ParC, as well as at the 385th and 526th amino acid residues in penicillin-binding protein 3. In addition, four intact prophage regions and one incomplete prophage region were found. ConclusionsThe whole-genome sequence of H. influenzae 2018-Y40 indicated that this clone emerged as a result of extensive genomic rearrangement by integration of multiple phages. As genomic rearrangement occasionally leads to a new phenotype, this clone could have acquired antimicrobial resistance and diversification via rearrangement events. These findings can form a basis to help clarify the mechanisms of low quinolone susceptibility and the spread of this outbreak clone.

The landscape of chromothripsis across adult cancer types

Chromothripsis is a recently identified mutational phenomenon, by which a presumably single catastrophic event generates extensive genomic rearrangements of one or a few chromosome(s). Considered as an early event in tumour development, this form of genome instability plays a prominent role in tumour onset. Chromothripsis prevalence might have been underestimated when using low-resolution methods, and pan-cancer studies based on sequencing are rare. Here we analyse chromothripsis in 28 tumour types covering all major adult cancers (634 tumours, 316 whole-genome and 318 whole-exome sequences). We show that chromothripsis affects a substantial proportion of human cancers, with a prevalence of 49% across all cases. Chromothripsis generates entity-specific genomic alterations driving tumour development, including clinically relevant druggable fusions. Chromothripsis is linked with specific telomere patterns and univocal mutational signatures in distinct tumour entities. Longitudinal analysis of chromothriptic patterns in 24 matched tumour pairs reveals insights in the clonal evolution of tumours with chromothripsis.

Extensive Genomic Rearrangements along with Distinct Mobilome and TALome are Associated with Extreme Pathotypes of a Rice Pathogen.

Xanthomonas oryzae pv. oryzae (Xoo) is a serious pathogen of rice which displays tremendous interstrain variation. The emergence of highly-virulent strains of Xoo is a major threat to rice cultivation. Evolutionary insights into genome dynamics of highly virulent strains as compared with the less-virulent ones are crucial for understanding the molecular basis of exceptional success of Xoo as a highly evolved plant pathogen. In the present study, we report complete genome sequence of Xoo strains with extreme-virulent pathotypes (XVPs) characterized based on their reaction toward ten resistance (Xa) genes. One strain, IXO1088, can overcome resistance mediated by all the ten resistance genes while the other strain IXO704 cannot overcome any of them. Interestingly, our investigation revealed that XVPs display dramatic variation in the genome structure with numerous rearrangements/inversions. Moreover, XVPs also possess distinct transposon content and prophage elements that may provide genomic flux required for the acquisition of novel gene cassettes and structural changes in the genome. Interestingly, analysis of transcription activator-like effector proteins, which are major virulence determinants of Xanthomonas pathogen show marked variation in the transcription activator-like effector content and DNA binding domain of tal genes. Overall, the present study indicates the possible role of mobilomes and repetitive elements in major structural and sequence alterations, which may be leading to the emergence of novel and extreme pathotypes. The knowledge and resource of XVPs will be invaluable in the further systematic understanding of evolution and management of variant pathotypes of Xoo.

Extensive chloroplast genome rearrangement amongst three closely related Halamphora spp. (Bacillariophyceae), and evidence for rapid evolution as compared to land plants.

Diatoms are the most diverse lineage of algae, but the diversity of their chloroplast genomes, particularly within a genus, has not been well documented. Herein, we present three chloroplast genomes from the genus Halamphora (H. americana, H. calidilacuna, and H. coffeaeformis), the first pennate diatom genus to be represented by more than one species. Halamphora chloroplast genomes ranged in size from ~120 to 150 kb, representing a 24% size difference within the genus. Differences in genome size were due to changes in the length of the inverted repeat region, length of intergenic regions, and the variable presence of ORFs that appear to encode as-yet-undescribed proteins. All three species shared a set of 161 core features but differed in the presence of two genes, serC and tyrC of foreign and unknown origin, respectively. A comparison of these data to three previously published chloroplast genomes in the non-pennate genus Cyclotella (Thalassiosirales) revealed that Halamphora has undergone extensive chloroplast genome rearrangement compared to other genera, as well as containing variation within the genus. Finally, a comparison of Halamphora chloroplast genomes to those of land plants indicates diatom chloroplast genomes within this genus may be evolving at least ~4–7 times faster than those of land plants. Studies such as these provide deeper insights into diatom chloroplast evolution and important genetic resources for future analyses.

Long-read sequencing reveals intra-species tolerance of substantial structural variations and new subtelomere formation in C. elegans.

Long-read sequencing technologies have contributed greatly to comparative genomics among species and can also be applied to study genomics within a species. In this study, to determine how substantial genomic changes are generated and tolerated within a species, we sequenced a C. elegans strain, CB4856, which is one of the most genetically divergent strains compared to the N2 reference strain. For this comparison, we used the Pacific Biosciences (PacBio) RSII platform (80×, N50 read length 11.8 kb) and generated de novo genome assembly to the level of pseudochromosomes containing 76 contigs (N50 contig = 2.8 Mb). We identified structural variations that affected as many as 2694 genes, most of which are at chromosome arms. Subtelomeric regions contained the most extensive genomic rearrangements, which even created new subtelomeres in some cases. The subtelomere structure of Chromosome VR implies that ancestral telomere damage was repaired by alternative lengthening of telomeres even in the presence of a functional telomerase gene and that a new subtelomere was formed by break-induced replication. Our study demonstrates that substantial genomic changes including structural variations and new subtelomeres can be tolerated within a species, and that these changes may accumulate genetic diversity within a species.

Highly accelerated rates of genomic rearrangements and nucleotide substitutions in plastid genomes of Passiflora subgenus Decaloba

Plastid genomes (plastomes) of photosynthetic angiosperms are for the most part highly conserved in their organization, mode of inheritance and rates of nucleotide substitution. A small number of distantly related lineages share a syndrome of features that deviate from this general pattern, including extensive genomic rearrangements, accelerated rates of nucleotide substitution, biparental inheritance and plastome-genome incompatibility. Previous studies of plastomes in Passiflora with limited taxon sampling suggested that the genus exhibits this syndrome. To examine this phenomenon further, 15 new plastomes from Passiflora were sequenced and combined with previously published data to examine the phylogenetic relationships, genome organization and evolutionary rates across all five subgenera and the sister genus Adenia. Phylogenomic analyses using 68 protein-coding genes shared by Passiflora generated a fully resolved and strongly supported tree that is congruent with previous phylogenies based on a few plastid and nuclear loci. This phylogeny was used to examine the distribution of plastome rearrangements across Passiflora. Multiple gene and intron losses and inversions were identified in Passiflora with some occurring in parallel and others that extended across the Passifloraceae. Furthermore, extensive expansions and contractions of the inverted repeat (IR) were uncovered and in some cases this resulted in exclusion of all ribosomal RNA genes from the IR. The most highly rearranged lineage was subgenus Decaloba, which experienced extensive IR expansion that incorporated up to 25 protein-coding genes usually located in large single copy region. Nucleotide substitution rate analyses of 68 protein-coding genes across the genus showed lineage- and locus-specific acceleration. Significant increase in dS, dN and dN/dS was detected for clpP across the genus and for ycf4 in certain lineages. Significant increases in dN and dN/dS for ribosomal subunits and plastid-encoded RNA polymerase genes were detected in the branch leading to the expanded IR-clade in subgenus Decaloba. This subgenus displays the syndrome of unusual features, making it an ideal system to investigate the dynamic evolution of angiosperm plastomes.

A high-quality chromosome-level genome assembly of a generalist herbivore, Trichoplusia ni.

The cabbage looper, Trichoplusia ni, is a globally distributed highly polyphagous herbivore and an important agricultural pest. T.ni has evolved resistance to various chemical insecticides, and is one of the only two insect species that have evolved resistance to the biopesticide Bacillus thuringiensis (Bt) in agricultural systems and has been selected for resistance to baculovirus infections. We report a 333-Mb high-quality T.ni genome assembly, which has N50 lengths of scaffolds and contigs of 4.6Mb and 140Kb, respectively, and contains 14,384 protein-coding genes. High-density genetic maps were constructed to anchor 305Mb (91.7%) of the assembly to 31 chromosomes. Comparative genomic analysis of T.ni with Bombyx mori showed enrichment of tandemly duplicated genes in T.ni in families involved in detoxification and digestion, consistent with the broad host range of T.ni. High levels of genome synteny were found between T.ni and other sequenced lepidopterans. However, genome synteny analysis of T.ni and the T.ni derived cell line High Five (Hi5) indicated extensive genome rearrangements in the cell line. These results provided the first genomic evidence revealing the high instability of chromosomes in lepidopteran cell lines known from karyotypic observations. The high-quality T.ni genome sequence provides a valuable resource for research in a broad range of areas including fundamental insect biology, insect-plant interactions and co-evolution, mechanisms and evolution of insect resistance to chemical and biological pesticides, and technology development for insect pest management.

Russian Doll Genes and Complex Chromosome Rearrangements in Oxytricha trifallax.

Ciliates have two different types of nuclei per cell, with one acting as a somatic, transcriptionally active nucleus (macronucleus; abbr. MAC) and another serving as a germline nucleus (micronucleus; abbr. MIC). Furthermore, Oxytricha trifallax undergoes extensive genome rearrangements during sexual conjugation and post-zygotic development of daughter cells. These rearrangements are necessary because the precursor MIC loci are often both fragmented and scrambled, with respect to the corresponding MAC loci. Such genome architectures are remarkably tolerant of encrypted MIC loci, because RNA-guided processes during MAC development reorganize the gene fragments in the correct order to resemble the parental MAC sequence. Here, we describe the germline organization of several nested and highly scrambled genes in Oxytricha trifallax. These include cases with multiple layers of nesting, plus highly interleaved or tangled precursor loci that appear to deviate from previously described patterns. We present mathematical methods to measure the degree of nesting between precursor MIC loci, and revisit a method for a mathematical description of scrambling. After applying these methods to the chromosome rearrangement maps of O. trifallax we describe cases of nested arrangements with up to five layers of embedded genes, as well as the most scrambled loci in O. trifallax.

Processing-Challenges Generated by Clusters of DNA Double-Strand Breaks Underpin Increased Effectiveness of High-LET Radiation and Chromothripsis.

Whereas most endogenous and exogenous DNA damaging agents typically generate lesions that are relatively isolated and can be repaired easily, ionizing radiation (IR) also induces clustered lesions causing DNA double strand breaks (DSBs). Moreover, forms of IR characterized by high linearenergytransfer (LET) induce not only isolated DSBs but also DSB clusters- multiple DSBs in close proximity -that pose increased risks for the cell. DSB clusters can destabilize chromatin locally and compromise processing of individual DSBs within the cluster. Since the discovery of chromothripsis, a phenomenon whereby multiple DSBs locally generated by a catastrophic event causes genomic rearrangements that feed carcinogenesis, DSB clusters receive increased attention also in the field of cancer. While formation of DSB clusters after exposure to high LET is a direct and inherent consequence of the spatial distribution of the constituting energy deposition events, also called track structure, the sources of local genomic shattering underpinning chromothripsis are under investigation. Notably, many consequences of DSB clusters in the affected genome reflect processing by pathways that have evolved to repair DSBs, but which operate with widely different degrees of fidelity. The molecular underpinnings and the basis of the underlying repair pathway choices that ultimately lead to the observed consequences from DSB clusters remain unknown. We developed a tractable model of DSB clustering that allows direct analysis in cells of the consequences of certain configurations of DSB clusters. We outline the rationale for the development of this model and describe its key characteristics. We summarize results suggesting that DSB clusters compromise the first-line DSB-processing pathways of c-NHEJ and HRR, increasing as a consequence the contribution of alt-EJ, which has high propensity of generating chromosomal rearrangements. The results suggest a mechanism for the increased toxicity of high LET radiation and the extensive genomic rearrangements associated with chromothripsis.

Classical and Molecular Cytogenetics of the Panther Genet Genetta maculata (Mammalia, Carnivora, Viverridae)

Genets (Genetta) are a genus of African mammalian carnivorans with 14 currently recognized species, although taxonomic uncertainties remain, particularly regarding the number of species within the large-spotted genet complex. This study presents the first banded karyotype and molecular cytogenetic analysis of a genetically identified panther genet, Genetta maculata, the most common and widespread taxon of the large-spotted genet complex, with a wide distribution in sub-Saharan Africa. Sampled in Gauteng Province, South Africa, it could be assigned to the subspecies G. m. letabae on geographic grounds and had a similar karyotype (2n = 52, FNa = 96) to those published for a panther genet from Ethiopia and for the West African large-spotted genet G. pardina. Notably, the specimen had a different autosomal morphology (2 acrocentric chromosomes) from that previously attributed to letabae (a single acrocentric chromosome), but the latter assignment was uncertain because the studied individuals were captive born and assigned based solely on a presumed origin in the former Transvaal Province of South Africa. Fluorescence in situ hybridization with a telomere repeat probe revealed the presence of telomeric sequences in the centromeres of most chromosomes, the so-called interstitial telomeric sites (ITSs). Since genets seem to have a unique, highly rearranged karyotype among feliforms and relatively low interspecific karyotypic variation, and considering the known instability of ITSs, we suggest that the large amount of ITSs found here might be due to evolutionarily recent extensive genomic rearrangements. This study provides cytogenetic information that contributes to our understanding of chromosomal variation and genomic rearrangements in genets, and valuable baseline data for future studies of karyotype evolution in carnivores in general and viverrids in particular.

The complete genome sequencing of Prevotella intermedia strain OMA14 and a subsequent fine-scale, intra-species genomic comparison reveal an unusual amplification of conjugative and mobile transposons and identify a novel Prevotella-lineage-specific repeat.

Prevotella intermedia is a pathogenic bacterium involved in periodontal diseases. Here, we present the complete genome sequence of a clinical strain, OMA14, of this bacterium along with the results of comparative genome analysis with strain 17 of the same species whose genome has also been sequenced, but not fully analysed yet. The genomes of both strains consist of two circular chromosomes: the larger chromosomes are similar in size and exhibit a high overall linearity of gene organizations, whereas the smaller chromosomes show a significant size variation and have undergone remarkable genome rearrangements. Unique features of the Pre. intermedia genomes are the presence of a remarkable number of essential genes on the second chromosomes and the abundance of conjugative and mobilizable transposons (CTns and MTns). The CTns/MTns are particularly abundant in the second chromosomes, involved in its extensive genome rearrangement, and have introduced a number of strain-specific genes into each strain. We also found a novel 188-bp repeat sequence that has been highly amplified in Pre. intermedia and are specifically distributed among the Pre. intermedia-related species. These findings expand our understanding of the genetic features of Pre. intermedia and the roles of CTns and MTns in the evolution of bacteria.

Genes encoding ten newly designated OXA-63 group class D β-lactamases identified in strains of the pathogenic intestinal spirochaete Brachyspira pilosicoli.

The anaerobic spirochaete Brachyspira pilosicoli colonizes the large intestine of birds and mammals, including human beings, and may induce colitis and diarrhoea. B. pilosicoli has a recombinant population structure, and strains show extensive genomic rearrangements and different genome sizes. The resident chromosomal gene blaOXA-63 in B. pilosicoli encodes OXA-63, a narrow-spectrum group IV class D β-lactamase. Genes encoding four OXA-63 variants have been described in B. pilosicoli, and the current study was designed to investigate the distribution and diversity of such genes and proteins in strains of B. pilosicoli. PCRs were used to amplify blaOXA-63 group genes from 118 B. pilosicoli strains from different host species and geographical origins. One primer set was targeted externally to the gene and two sets were designed to amplify internal components. A total of 16 strains (13.6%) showed no evidence of possessing blaOXA-63 group genes, 44 (37.3%) had a full gene, 27 (22.9%) apparently had a gene but it failed to amplify with external primers, and 29 (24.6%) had only one or other of the two internal components amplified. Based on translation of the nucleotide sequences, ten new variants of the β-lactamase, designated OXA-470 through OXA-479, were identified amongst the 44 strains that had the full gene amplified. The 16 strains lacking blaOXA-63 group genes had a region of 1674 bp missing around where the gene was expected to reside. Despite apparent genomic rearrangements occurring in B. pilosicoli, positive selection pressures for conservation of blaOXA-63 group genes and OXA proteins appear to have been exerted.

Discovery of an Unconventional Centromere in Budding Yeast Redefines Evolution of Point Centromeres

SummaryCentromeres are the chromosomal regions promoting kinetochore assembly for chromosome segregation. In many eukaryotes, the centromere consists of up to mega base pairs of DNA. On such “regional centromeres,” kinetochore assembly is mainly defined by epigenetic regulation [1]. By contrast, a clade of budding yeasts (Saccharomycetaceae) has a “point centromere” of 120–200 base pairs of DNA, on which kinetochore assembly is defined by the consensus DNA sequence [2, 3]. During evolution, budding yeasts acquired point centromeres, which replaced ancestral, regional centromeres [4]. All known point centromeres among different yeast species share common consensus DNA elements (CDEs) [5, 6], implying that they evolved only once and stayed essentially unchanged throughout evolution. Here, we identify a yeast centromere that challenges this view: that of the budding yeast Naumovozyma castellii is the first unconventional point centromere with unique CDEs. The N. castellii centromere CDEs are essential for centromere function but have different DNA sequences from CDEs in other point centromeres. Gene order analyses around N. castellii centromeres indicate their unique, and separate, evolutionary origin. Nevertheless, they are still bound by the ortholog of the CBF3 complex, which recognizes CDEs in other point centromeres. The new type of point centromere originated prior to the divergence between N. castellii and its close relative Naumovozyma dairenensis and disseminated to all N. castellii chromosomes through extensive genome rearrangement. Thus, contrary to the conventional view, point centromeres can undergo rapid evolutionary changes. These findings give new insights into the evolution of point centromeres.

Chromothripsis from DNA damage in micronuclei.

Genome sequencing has uncovered a new mutational phenomenon in cancer and congenital disorders called chromothripsis. Chromothripsis is characterized by extensive genomic rearrangements and an oscillating pattern of DNA copy number levels, all curiously restricted to one or a few chromosomes. The mechanism for chromothripsis is unknown, but we previously proposed that it could occur through the physical isolation of chromosomes in aberrant nuclear structures called micronuclei. Here, using a combination of live-cell imaging and single-cell genome sequencing, we demonstrate that micronucleus formation can indeed generate a spectrum of genomic rearrangements, some of which recapitulate all known features of chromothripsis. These events are restricted to the missegregated chromosome and occur within one cell division. We demonstrate that the mechanism for chromothripsis can involve the fragmentation and subsequent reassembly of a single chromatid from a micronucleus. Collectively, these experiments establish a new mutational process of which chromothripsis is one extreme outcome.

The lineage-specific evolution of aquaporin gene clusters facilitated tetrapod terrestrial adaptation.

A major physiological barrier for aquatic organisms adapting to terrestrial life is dessication in the aerial environment. This barrier was nevertheless overcome by the Devonian ancestors of extant Tetrapoda, but the origin of specific molecular mechanisms that solved this water problem remains largely unknown. Here we show that an ancient aquaporin gene cluster evolved specifically in the sarcopterygian lineage, and subsequently diverged into paralogous forms of AQP2, -5, or -6 to mediate water conservation in extant Tetrapoda. To determine the origin of these apomorphic genomic traits, we combined aquaporin sequencing from jawless and jawed vertebrates with broad taxon assembly of >2,000 transcripts amongst 131 deuterostome genomes and developed a model based upon Bayesian inference that traces their convergent roots to stem subfamilies in basal Metazoa and Prokaryota. This approach uncovered an unexpected diversity of aquaporins in every lineage investigated, and revealed that the vertebrate superfamily consists of 17 classes of aquaporins (Aqp0 - Aqp16). The oldest orthologs associated with water conservation in modern Tetrapoda are traced to a cluster of three aqp2-like genes in Actinistia that likely arose >500 Ma through duplication of an aqp0-like gene present in a jawless ancestor. In sea lamprey, we show that aqp0 first arose in a protocluster comprised of a novel aqp14 paralog and a fused aqp01 gene. To corroborate these findings, we conducted phylogenetic analyses of five syntenic nuclear receptor subfamilies, which, together with observations of extensive genome rearrangements, support the coincident loss of ancestral aqp2-like orthologs in Actinopterygii. We thus conclude that the divergence of sarcopterygian-specific aquaporin gene clusters was permissive for the evolution of water conservation mechanisms that facilitated tetrapod terrestrial adaptation.

Conserved gene order and expanded inverted repeats characterize plastid genomes of Thalassiosirales.

Diatoms are mostly photosynthetic eukaryotes within the heterokont lineage. Variable plastid genome sizes and extensive genome rearrangements have been observed across the diatom phylogeny, but little is known about plastid genome evolution within order- or family-level clades. The Thalassiosirales is one of the more comprehensively studied orders in terms of both genetics and morphology. Seven complete diatom plastid genomes are reported here including four Thalassiosirales: Thalassiosira weissflogii, Roundia cardiophora, Cyclotella sp. WC03_2, Cyclotella sp. L04_2, and three additional non-Thalassiosirales species Chaetoceros simplex, Cerataulina daemon, and Rhizosolenia imbricata. The sizes of the seven genomes vary from 116,459 to 129,498 bp, and their genomes are compact and lack introns. The larger size of the plastid genomes of Thalassiosirales compared to other diatoms is due primarily to expansion of the inverted repeat. Gene content within Thalassiosirales is more conserved compared to other diatom lineages. Gene order within Thalassiosirales is highly conserved except for the extensive genome rearrangement in Thalassiosira oceanica. Cyclotella nana, Thalassiosira weissflogii and Roundia cardiophora share an identical gene order, which is inferred to be the ancestral order for the Thalassiosirales, differing from that of the other two Cyclotella species by a single inversion. The genes ilvB and ilvH are missing in all six diatom plastid genomes except for Cerataulina daemon, suggesting an independent gain of these genes in this species. The acpP1 gene is missing in all Thalassiosirales, suggesting that its loss may be a synapomorphy for the order and this gene may have been functionally transferred to the nucleus. Three genes involved in photosynthesis, psaE, psaI, psaM, are missing in Rhizosolenia imbricata, which represents the first documented instance of the loss of photosynthetic genes in diatom plastid genomes.

Plant-driven genome selection of arbuscular mycorrhizal fungi.

Plant-driven genome selection of arbuscular mycorrhizal fungi.