Density Of Transposable Elements Research Articles

Chromosome mapping of retrotransposon AviRTE in a neotropical bird species: Trogon surrucura (Trogoniformes; Trogonidae).

Avian genomes are characterized as being more compact than other amniotes, with less diversity and density of transposable elements (TEs). In addition, birds usually show bimodal karyotypes, exhibiting a great variation in diploid numbers. Some species present unusually large sex chromosomes, possibly due to the accumulation of repetitive sequences. Avian retrotransposon-like element (AviRTE) is a long interspersed nuclear element (LINE) recently discovered in the genomes of birds and nematodes, and it is still poorly characterized in terms of chromosomal mapping and phylogenetic relationships. In this study, we mapped AviRTE isolated from the Trogon surrucura genome into the T. surrucura (TSU) karyotype. Furthermore, we analyzed the phylogenetic relationships of this LINE in birds and other vertebrates. Our results showed that the distribution pattern of AviRTE is not restricted to heterochromatic regions, with accumulation on the W chromosome of TSU, yet another species with an atypical sex chromosome and TE hybridization. The phylogenetic analysis of AviRTE sequences in birds agreed with the proposed phylogeny of species in most clades, and allowed the detection of this sequence in other species, expanding the distribution of the element.

Genome-wide landscapes of genes and repeatome reveal the genomic differences between the two subspecies of peanut (Arachis hypogaea)

Distribution and structural features of genes, repeat elements and transposable elements (TEs) were studied to identify the genomic differences between the two subspecies (ssp. hypogaea and ssp. fastigiata) of peanut (Arachis hypogaea L.). A total of 67128 predicted genes, 2738666 copies of TEs and 162361 tandem repeats from the reference genome of Tifrunner were employed for this study. Of the 67128 genes, 33622 were reading on the plus strand, while 33506 were traced on the minus strand. Though B03 had the highest number of genes (4524), A08 recorded the highest density (53 genes/Mb) of genes in the genome. Telomeric regions had the highest density of genes. The average length of the genes was 3971 bp with majority of the genes (39228) containing one to five exons. The gene Arahy.11.0DU9MH had the insertion of 28 different types of TEs, and was the longest gene in the peanut genome. A total of 15731 genes were monomorphic in terms of SNPs across 179 accessions, while 7401 genes showed polymorphism at one nucleotide, indicating very low allelic variation at these genes. Remaining 66% of the genes had two or more SNPs, and therefore showed relatively high allelic variation. Among the 101 unique types of TEs, the Retro elements (869279) followed by CACTA (272596) and Mu (250248) TEs were most predominant. Telomeric regions showed less density of TEs than the regions. On an average, each gene contained 1.8 copies of TEs, and 35706 genes did not have the insertion of any TEs. A QTL-Seq approach could identify 186 SNPs and 26 gene differences between the two subspecies of A. hypogaea. Two of the 26 genes showed allelic variation in terms of SNPs and TEs.

The Genome of Fusarium oxysporum f. sp. phaseoli Provides Insight into the Evolution of Genomes and Effectors of Fusarium oxysporum Species.

Fusarium oxysporum f. sp. phaseoli, the causal agent of cowpea fusarium wilt, is a serious threat to cowpea production in China. In this study, a sample of cowpea fusarium wilt was identified as Fusarium oxysporum f. sp. phaseoli using the methods of morphological characters and molecular detection. We further reported the first genome assembly for Fusarium oxysporum f. sp. phaseoli, with 53.7 Mb genome sequence comprising 14,694 genes. Comparative genomic analysis among five Fusarium oxysporum genomes showed that four accessory chromosomes in the five Fusarium oxysporum display similar characteristics, with low sequence similarity (55.35%, vs. overall average of 81.76%), low gene density (2.18 genes/10 kb vs. 3.02 genes/Mb) and highly transposable element density (TEs) (15.01/100 kb vs. 4.89/100 kb), indicating that variable accessory chromosomes are the main source of Fusarium oxysporum evolution. We identified a total of 100 Fusarium oxysporum f. sp. phaseoli-specific effectors in the genome and found 13 specific effector genes located in large insertion or deletion regions, suggesting that insertion or deletion events can cause the emergence of species-specific effectors in Fusarium oxysporum. Our genome assembly of Fusarium oxysporum f. sp. phaseoli provides a valuable resource for the study of cowpea fusarium wilt, and the comparative genomic study of Fusarium oxysporum could contribute to the knowledge of genome and effector-associated pathogenicity evolution in Fusarium oxysporum study.

Chromosome-level genome and recombination map of the male buffalo.

The swamp buffalo (Bubalus bubalis carabanesis) is an economically important livestock supplying milk, meat, leather, and draft power. Several female buffalo genomes have been available, but the lack of high-quality male genomes hinders studies on chromosome evolution, especially Y, as well as meiotic recombination. Here, a chromosome-level genome with a contig N50 of 72.2 Mb and a fine-scale recombination map of male buffalo were reported. We found that transposable elements (TEs) and structural variants (SVs) may contribute to buffalo evolution by influencing adjacent gene expression. We further found that the pseudoautosomal region (PAR) of the Y chromosome is subject to stronger purification selection. The meiotic recombination map showed that there were 2 obvious recombination hotspots on chromosome 8, and the genes around them were mainly related to tooth development, which may have helped to enhance the adaption of buffalo to inferior feed. Among several genomic features, TE density has the strongest correlation with recombination rates. Moreover, the TE subfamily, SINE/tRNA, is likely to play a role in driving recombination into SVs. The male genome and sperm sequencing will facilitate the understanding of the buffalo genomic evolution and functional research.

Global dissection of the recombination landscape in soybean using a high-density 600K SoySNP array.

Recombination is crucial for crop breeding because it can break linkage drag and generate novel allele combinations. However, the high-resolution recombination landscape and its driving forces in soybean are largely unknown. Here, we constructed eight recombinant inbred line (RIL) populations and genotyped individual lines using the high-density 600K SoySNP array, which yielded a high-resolution recombination map with 5636 recombination sites at a resolution of 1.37 kb. The recombination rate was negatively correlated with transposable element density and GC content but positively correlated with gene density. Interestingly, we found that meiotic recombination was enriched at the promoters of active genes. Further investigations revealed that chromatin accessibility and active epigenetic modifications promoted recombination. Our findings provide important insights into the control of homologous recombination and thus will increase our ability to accelerate soybean breeding by manipulating meiotic recombination rate.

Insights from the genomes of 4 diploid Camelina spp.

Plant evolution has been a complex process involving hybridization and polyploidization making understanding the origin and evolution of a plant's genome challenging even once a published genome is available. The oilseed crop, Camelina sativa (Brassicaceae), has a fully sequenced allohexaploid genome with 3 unknown ancestors. To better understand which extant species best represent the ancestral genomes that contributed to C. sativa's formation, we sequenced and assembled chromosome level draft genomes for 4 diploid members of Camelina: C. neglecta C. hispida var. hispida, C. hispida var. grandiflora, and C. laxa using long and short read data scaffolded with proximity data. We then conducted phylogenetic analyses on regions of synteny and on genes described for Arabidopsis thaliana, from across each nuclear genome and the chloroplasts to examine evolutionary relationships within Camelina and Camelineae. We conclude that C. neglecta is closely related to C. sativa's sub-genome 1 and that C. hispida var. hispida and C. hispida var. grandiflora are most closely related to C. sativa's sub-genome 3. Further, the abundance and density of transposable elements, specifically Helitrons, suggest that the progenitor genome that contributed C. sativa's sub-genome 3 maybe more similar to the genome of C. hispida var. hispida than that of C. hispida var. grandiflora. These diploid genomes show few structural differences when compared to C. sativa's genome indicating little change to chromosome structure following allopolyploidization. This work also indicates that C. neglecta and C. hispida are important resources for understanding the genetics of C. sativa and potential resources for crop improvement.

The genome sequence of the scarce swallowtail, Iphiclides podalirius.

The scarce swallowtail, Iphiclides podalirius (Linnaeus, 1758), is a species of butterfly in the family Papilionidae. Here, we present a chromosome-level genome assembly for Iphiclides podalirius as well as gene and transposable element annotations. We investigate how the density of genomic features differs between the 30 Iphiclides podalirius chromosomes. We find that shorter chromosomes have higher heterozygosity at four-fold-degenerate sites and a greater density of transposable elements. While the first result is an expected consequence of differences in recombination rate, the second suggests a counter-intuitive relationship between recombination and transposable element evolution. This high-quality genome assembly, the first for any species in the tribe Leptocircini, will be a valuable resource for population genomics in the genus Iphiclides and comparative genomics more generally.

Chromosome-level genome assembly of the aquatic plant Nymphoides indica reveals transposable element bursts and NBS-LRR gene family expansion shedding light on its invasiveness.

Nymphoides indica, an aquatic plant, is an invasive species that causes both ecological and economic damage in North America and elsewhere. However, the lack of genomic data of N. indica limits the in-depth analysis of this invasive species. Here, we report a chromosome-level genome assembly of nine pseudochromosomes of N. indica with a total size of ∼ 520 Mb. More than half of the N. indica genome consists of transposable elements (TEs), and a higher density of TEs around genes may play a significant role in response to an ever-changing environment by regulating the nearby gene. Additionally, our analysis revealed that N. indica only experienced a gamma (γ) whole-genome triplication event. Functional enrichment of the N. indica-specific and expanded gene families highlighted genes involved in the responses to hypoxia and plant–pathogen interactions, which may strengthen the ability to adapt to external challenges and improve ecological fitness. Furthermore, we identified 160 members of the nucleotide-binding site and leucine-rich repeat gene family, which may be linked to the defence response. Collectively, the high-quality N. indica genome reported here opens a novel avenue to understand the evolution and rapid invasion of Nymphoides spp.

Accumulation and ineffective silencing of transposable elements on an avian W Chromosome.

One of the defining features of transposable elements (TEs) is their ability to move to new locations in the host genome. To minimize the potentially deleterious effects of de novo TE insertions, hosts have evolved several mechanisms to control TE activity, including recombination-mediated removal and epigenetic silencing; however, increasing evidence suggests that silencing of TEs is often incomplete. The crow family experienced a recent radiation of LTR retrotransposons (LTRs), offering an opportunity to gain insight into the regulatory control of young, potentially still active TEs. We quantified the abundance of TE-derived transcripts across several tissues in 15 Eurasian crows (Corvus (corone) spp.) raised under common garden conditions and find evidence for ineffective TE suppression on the female-specific W Chromosome. Using RNA-seq data, we show that ∼9.5% of all transcribed TEs had considerably greater (average, 16-fold) transcript abundance in female crows and that >85% of these female-biased TEs originated on the W Chromosome. After accounting for differences in TE density among chromosomal classes, W-linked TEs were significantly more highly expressed than TEs residing on other chromosomes, consistent with ineffective silencing on the former. Together, our results suggest that the crow W Chromosome acts as a source of transcriptionally active TEs, with possible negative fitness consequences for female birds analogous to Drosophila (an X/Y system), in which overexpression of Y-linked TEs is associated with male-specific aging and fitness loss (“toxic Y”).

Genome Sequencing of the Endangered Kingdonia uniflora (Circaeasteraceae, Ranunculales) Reveals Potential Mechanisms of Evolutionary Specialization.

SummaryKingdonia uniflora, an alpine herb, has an extremely narrow distribution and represents a model for studying evolutionary mechanisms of species that have adapted to undisturbed environments for evolutionarily long periods of time. We assembled a 1,004.7-Mb draft genome (encoding 43,301 genes) of K. uniflora and found significant overrepresentation in gene families associated with DNA repair, underrepresentation in gene families associated with stress response, and loss of most plastid ndh genes. During the evolutionary process, the overrepresentation of gene families involved in DNA repair could help asexual K. uniflora reduce the accumulation of deleterious mutations, while reducing genetic diversity, which is important in responding to environment fluctuations. The underrepresentation of gene families related to stress response and functional loss of ndh genes could be due to lack or loss of ability to respond to environmental changes caused by long-term adaptation to a relatively stable ecological environment.

QTL mapping of mycelial growth and aggressiveness to distinct hosts in Ceratocystis pathogens

Some species of Ceratocystis display strong host specificity, such as C. fimbriata sensu stricto that is restricted to sweet potato (Ipomoea batatas) as host. In contrast, the closely related C. manginecans, infects Acacia mangium and Mangifera indica but is not pathogenic to I. batatas. Despite the economic importance of these fungi, knowledge regarding the genetic factors that influence their pathogenicity and host specificity is limited. A recent inheritance study, based on an interspecific cross between C. fimbriata and C. manginecans and the resultant 70 F1 progeny, confirmed that traits such as mycelial growth rate, spore production and aggressiveness on A. mangium and I. batatas are regulated by multiple genes. In the present study, a quantitative trait locus (QTL) analysis was performed to determine the genomic loci associated with these traits. All 70 progeny isolates were genotyped with SNP markers and a linkage map was constructed. The map contained 467 SNPs, distributed across nine linkage groups, with a total length of 1203 cm. Using the progeny genotypes and phenotypes, one QTL was identified on the linkage map for mycelial growth rate, one for aggressiveness to A. mangium and two for aggressiveness to I. batatas (P < 0.05). Two candidate genes, likely associated with mycelial growth rate, were identified in the QTL region. The three QTLs associated with aggressiveness to different hosts contained candidate genes involved in protein processing, detoxification and regions with effector genes and high transposable element density. The results provide a foundation for studies considering the function of genes regulating various quantitative traits in Ceratocystis.

Adaptation and Phenotypic Diversification in Arabidopsis through Loss-of-Function Mutations in Protein-Coding Genes.

According to the less-is-more hypothesis, gene loss is an engine for evolutionary change. Loss-of-function (LoF) mutations resulting in the natural knockout of protein-coding genes not only provide information about gene function but also play important roles in adaptation and phenotypic diversification. Although the less-is-more hypothesis was proposed two decades ago, it remains to be explored on a large scale. In this study, we identified 60,819 LoF variants in 1071 Arabidopsis (Arabidopsis thaliana) genomes and found that 34% of Arabidopsis protein-coding genes annotated in the Columbia-0 genome do not have any LoF variants. We found that nucleotide diversity, transposable element density, and gene family size are strongly correlated with the presence of LoF variants. Intriguingly, 0.9% of LoF variants with minor allele frequency larger than 0.5% are associated with climate change. In addition, in the Yangtze River basin population, 1% of genes with LoF mutations were under positive selection, providing important insights into the contribution of LoF mutations to adaptation. In particular, our results demonstrate that LoF mutations shape diverse phenotypic traits. Overall, our results highlight the importance of the LoF variants for the adaptation and phenotypic diversification of plants.

A unique cluster of roo insertions in the promoter region of a stress response gene in Drosophila melanogaster

Transposable elements (TEs) are not randomly distributed in the genome. A genome-wide analysis of the D. melanogaster genome found that differences in TE density across 50 kb genomic regions was due both to transposition and duplication. At smaller genomic scales, promoter regions of hsp genes and the promoter region of CG18446 have been shown to accumulate TE insertions. In this work, we have further analyzed the promoter region of CG18446. We screened 218 strains collected in 15 natural populations, and we found that the CG18446 promoter region contains 20 independent roo insertions. Based on phylogenetic analysis, we suggest that the presence of multiple roo insertions in this region is likely to be the result of several bursts of transposition. Moreover, we found that the roo insertional cluster in the CG18446 promoter region is unique: no other promoter region in the genome contains a similar number of roo insertions. We found that, similar to hsp gene promoters, chromatin accessibility could be one of the factors explaining the recurrent insertions of roo elements in CG18446 promoter region.

Impact of polymorphic transposable elements on linkage disequilibrium along chromosomes.

Recombination and selection drive the extent of linkage disequilibrium (LD) among loci and therefore affect the reshuffling of adaptive genetic variation. However, it is poorly known to what extent the enrichment of transposable elements (TEs) in recombinationally-inert regions reflects their inefficient removal by purifying selection and whether the presence of polymorphic TEs can modify the local recombination rate. In this study, we investigate how TEs and recombination interact at fine scale along chromosomes and possibly support linked selection in natural populations. Whole-genome sequencing data of 304 individuals from nearby alpine populations of Arabis alpina were used to show that the density of polymorphic TEs is specifically correlated with local LD along chromosomes. Consistent with TEs modifying recombination, the characterization of 28 such LD blocks of up to 5.5Mb in length revealed strong evidence of selective sweeps at a few loci through either site frequency spectrum or haplotype structure. A majority of these blocks were enriched in genes related to ecologically relevant functions such as responses to cold, salt stress or photoperiodism. In particular, the S-locus (i.e., supergene responsible for strict outcrossing) was identified in a LD block with high levels of polymorphic TEs and evidence of selection. Another such LD block was enriched in cold-responding genes and presented evidence of adaptive loci related to photoperiodism and flowering being increasingly linked by polymorphic TEs. These results are consistent with the hypothesis that TEs modify recombination landscapes and thus interact with selection in driving blocks of linked adaptive loci in natural populations.

Contingency in the convergent evolution of a regulatory network: Dosage compensation in Drosophila.

The repeatability or predictability of evolution is a central question in evolutionary biology and most often addressed in experimental evolution studies. Here, we infer how genetically heterogeneous natural systems acquire the same molecular changes to address how genomic background affects adaptation in natural populations. In particular, we take advantage of independently formed neo-sex chromosomes in Drosophila species that have evolved dosage compensation by co-opting the dosage-compensation male-specific lethal (MSL) complex to study the mutational paths that have led to the acquisition of hundreds of novel binding sites for the MSL complex in different species. This complex recognizes a conserved 21-bp GA-rich sequence motif that is enriched on the X chromosome, and newly formed X chromosomes recruit the MSL complex by de novo acquisition of this binding motif. We identify recently formed sex chromosomes in the D. melanica and D. robusta species groups by genome sequencing and generate genomic occupancy maps of the MSL complex to infer the location of novel binding sites. We find that diverse mutational paths were utilized in each species to evolve hundreds of de novo binding motifs along the neo-X, including expansions of microsatellites and transposable element (TE) insertions. However, the propensity to utilize a particular mutational path differs between independently formed X chromosomes and appears to be contingent on genomic properties of that species, such as simple repeat or TE density. This establishes the “genomic environment” as an important determinant in predicting the outcome of evolutionary adaptations.

Genetic Properties of a Nested Association Mapping Population Constructed With Semi-Winter and Spring Oilseed Rapes.

Nested association mapping (NAM) populations have been widely applied to dissect the genetic basis of complex quantitative traits in a variety of crops. In this study, we developed a Brassica napus NAM (BN-NAM) population consisting of 15 recombination inbred line (RIL) families with 2,425 immortal genotypes. Fifteen high-density genetic linkage maps were constructed by genotyping by sequencing (GBS) based on all RIL families, with further integration into a joint linkage map (JLM) having 30,209 unique markers in common with multiple linkage maps. Furthermore, an ultra-density whole-genome variation map was constructed by projecting 4,444,309 high-quality variants onto the JLM. The NAM population captured a total of 88,542 recombination events (REs). The uneven distribution of recombination rate along chromosomes is positively correlated with the densities of genes and markers, but negatively correlated with the density of transposable elements and linkage disequilibrium (LD). Analyses of population structure and principal components revealed that the BN-NAM population could be divided into three groups with weak stratification. The LD decay distance across genome varied between 170 and 2,400 Kb, with LD decay more rapid in the A than in the C sub-genome. The pericentromeric regions contained large LD blocks, especially in the C sub-genome. This NAM population provides a valuable resource for dissecting the genetic basis of important traits in rapeseed, especially in semi-winter oilseed rape.

The genomic architecture and molecular evolution of ant odorant receptors.

The massive expansions of odorant receptor (OR) genes in ant genomes are notable examples of rapid genome evolution and adaptive gene duplication. However, the molecular mechanisms leading to gene family expansion remain poorly understood, partly because available ant genomes are fragmentary. Here, we present a highly contiguous, chromosome-level assembly of the clonal raider ant genome, revealing the largest known OR repertoire in an insect. While most ant ORs originate via local tandem duplication, we also observe several cases of dispersed duplication followed by tandem duplication in the most rapidly evolving OR clades. We found that areas of unusually high transposable element density (TE islands) were depauperate in ORs in the clonal raider ant, and found no evidence for retrotransposition of ORs. However, OR loci were enriched for transposons relative to the genome as a whole, potentially facilitating tandem duplication by unequal crossing over. We also found that ant OR genes are highly AT-rich compared to other genes. In contrast, in flies, OR genes are dispersed and largely isolated within the genome, and we find that fly ORs are not AT-rich. The genomic architecture and composition of ant ORs thus show convergence with the unrelated vertebrate ORs rather than the related fly ORs. This might be related to the greater gene numbers and/or potential similarities in gene regulation between ants and vertebrates as compared to flies.

Modeling transposable element dynamics with fragmentation equations

Transposable elements (TEs), segments of DNA capable of self-replication, are abundant in the genomes of most organisms and thus serve as a record of past mutational events. While some work suggests TEs may serve a regulatory function for the host, most empirical and theoretical studies have shown that TEs often have deleterious effects on a host. Because they are not essential, the host genome consists of both full-length (actively replicating) and partial length (inactive remnant) copies of TEs. We developed a novel mathematical formulation of TE dynamics by modeling the density of full and partial length copies resulting from mutations (insertions and deletions) and TE replication within the host genome. More specifically, we model the time-evolution of the complete TE length distribution (full and partial elements) in a genome using fragmentation equations in both a discrete and continuous framework under two models of TE replication.In the first case, we assume that full-length TEs replicate at a constant rate regardless of the number of full-length TEs present in the genome. While this assumption simplifies the underlying biological processes, it allows us to derive an explicit analytical form of the time-varying TE density, as well as the steady-state behavior, under both discrete and continuous formulations. Next, we take into account the potential deleterious effects of TEs by modeling TE replication with a logistic growth equation. Under this assumption, the number of actively replicating TEs in a genome is limited by a carrying capacity. While we are still able to derive to derive analytical forms for the time-varying TE density, for both discrete and continuous length formulations, these solutions are implicit. For all four proposed models, we prove existence and uniqueness of these solutions describing TE length distributions. We compare both models and note that the logistic and exponential models initially agree. Since most TEs have not reached carrying capacity, we use the explicit exponential solution to quantify rates of replication to mutations. We apply our model to present day annotated collections of TEs from the genomes of species of fruit-flies, birds, and primates to uncover quantitative relationships of TE dynamics. With the increasing availability of complete genomes, our method is likely to uncover relationships of biological drivers of genomic variation in many species.

Low diversity, activity, and density of transposable elements in five avian genomes

In this study, we conducted the activity, diversity, and density analysis of transposable elements (TEs) across five avian genomes (budgerigar, chicken, turkey, medium ground finch, and zebra finch) to explore the potential reason of small genome sizes of birds. We found that these avian genomes exhibited low density of TEs by about 10% of genome coverages and low diversity of TEs with the TE landscapes dominated by CR1 and ERV elements, and contrasting proliferation dynamics both between TE types and between species were observed across the five avian genomes. Phylogenetic analysis revealed that CR1 clade was more diverse in the family structure compared with R2 clade in birds; avian ERVs were classified into four clades (alpha, beta, gamma, and ERV-L) and belonged to three classes of ERV with an uneven distributed in these lineages. The activities of DNA and SINE TEs were very low in the evolution history of avian genomes; most LINEs and LTRs were ancient copies with a substantial decrease of activity in recent, with only LTRs and LINEs in chicken and zebra finch exhibiting weak activity in very recent, and very few TEs were intact; however, the recent activity may be underestimated due to the sequencing/assembly technologies in some species. Overall, this study demonstrates low diversity, activity, and density of TEs in the five avian species; highlights the differences of TEs in these lineages; and suggests that the current and recent activity of TEs in avian genomes is very limited, which may be one of the reasons of small genome sizes in birds.

Accumulation of transposable elements in Hox gene clusters during adaptive radiation of Anolis lizards.

Transposable elements (TEs) are DNA sequences that can insert elsewhere in the genome and modify genome structure and gene regulation. The role of TEs in evolution is contentious. One hypothesis posits that TE activity generates genomic incompatibilities that can cause reproductive isolation between incipient species. This predicts that TEs will accumulate during speciation events. Here, I tested the prediction that extant lineages with a relatively high rate of speciation have a high number of TEs in their genomes. I sequenced and analysed the TE content of a marker genomic region (Hox clusters) in Anolis lizards, a classic case of an adaptive radiation. Unlike other vertebrates, including closely related lizards, Anolis lizards have high numbers of TEs in their Hox clusters, genomic regions that regulate development of the morphological adaptations that characterize habitat specialists in these lizards. Following a burst of TE activity in the lineage leading to extant Anolis, TEs have continued to accumulate during or after speciation events, resulting in a positive relationship between TE density and lineage speciation rate. These results are consistent with the prediction that TE activity contributes to adaptive radiation by promoting speciation. Although there was no evidence that TE density per se is associated with ecological morphology, the activity of TEs in Hox clusters could have been a rich source for phenotypic variation that may have facilitated the rapid parallel morphological adaptation to microhabitats seen in extant Anolis lizards.