Transposable element insertions have strongly affected human evolution

Advancing the Indian cattle pangenome: characterizing non-reference sequences in Bos indicus

BACKGROUND

India harbors the world's largest cattle population, encompassing over 50 distinct Bos indicus breeds. This rich genetic diversity underscores the inadequacy of a single reference genome to fully capture the genomic landscape of Indian cattle. To comprehensively characterize the genomic variation within Bos indicus and, specifically, dairy breeds, we aim to identify non-reference sequences and construct a comprehensive pangenome.

RESULTS

Five representative genomes of prominent dairy breeds, including Gir, Kankrej, Tharparkar, Sahiwal, and Red Sindhi, were sequenced using 10X Genomics 'linked-read' technology. Assemblies generated from these linked-reads ranged from 2.70 Gb to 2.77 Gb, comparable to the Bos indicus Brahman reference genome. A pangenome of Bos indicus cattle was constructed by comparing the newly assembled genomes with the reference using alignment and graph-based methods, revealing 8 Mb and 17.7 Mb of novel sequence respectively. A confident set of 6,844 Non-reference Unique Insertions (NUIs) spanning 7.57 Mb was identified through both methods, representing the pangenome of Indian Bos indicus breeds. Comparative analysis with previously published pangenomes unveiled 2.8 Mb (37%) commonality with the Chinese indicine pangenome and only 1% commonality with the Bos taurus pangenome. Among these, 2,312 NUIs encompassing ~ 2 Mb, were commonly found in 98 samples of the 5 breeds and designated as Bos indicus Common Insertions (BICIs) in the population. Furthermore, 926 BICIs were identified within 682 protein-coding genes, 54 long non-coding RNAs (lncRNA), and 18 pseudogenes. These protein-coding genes were enriched for functions such as chemical synaptic transmission, cell junction organization, cell-cell adhesion, and cell morphogenesis. The protein-coding genes were found in various prominent quantitative trait locus (QTL) regions, suggesting potential roles of BICIs in traits related to milk production, reproduction, exterior, health, meat, and carcass. Notably, 63.21% of the bases within the BICIs call set contained interspersed repeats, predominantly Long Interspersed Nuclear Elements (LINEs). Additionally, 70.28% of BICIs are shared with other domesticated and wild species, highlighting their evolutionary significance.

CONCLUSIONS

This is the first report unveiling a robust set of NUIs defining the pangenome of Bos indicus breeds of India. The analyses contribute valuable insights into the genomic landscape of desi cattle breeds.

Dysregulation of Transposon Transcription Profiles in Cancer Cells Resembles That of Embryonic Stem Cells

Transposable elements (TEs) comprise a substantial portion of the mammalian genome, with potential implications for both embryonic development and cancer. This study aimed to characterize the expression profiles of TEs in embryonic stem cells (ESCs), cancer cell lines, tumor tissues, and the tumor microenvironment (TME). We observed similarities in TE expression profiles between cancer cells and ESCs, suggesting potential parallels in regulatory mechanisms. Notably, four TE RNAs (HERVH, LTR7, HERV-Fc1, HERV-Fc2) exhibited significant downregulation across cancer cell lines and tumor tissues compared to ESCs, highlighting potential roles in pluripotency regulation. The strong up-regulation of the latter two TEs (HERV-Fc1, HERV-Fc2) in ESCs has not been previously demonstrated and may be a first indication of their role in the regulation of pluripotency. Conversely, tandemly repeated sequences (MSR1, CER, ALR) showed up-regulation in cancer contexts. Moreover, a difference in TE expression was observed between the TME and the tumor bulk transcriptome, with distinct dysregulated TE profiles. Some TME-specific TEs were absent in normal tissues, predominantly belonging to LTR and L1 retrotransposon families. These findings not only shed light on the regulatory roles of TEs in both embryonic development and cancer but also suggest novel targets for anti-cancer therapy. Understanding the interplay between cancer cells and the TME at the TE level may pave the way for further research into therapeutic interventions.

Accurate allocation of multimapped reads enables regulatory element analysis at repeats

Transposable elements (TEs) and other repetitive regions have been shown to contain gene regulatory elements, including transcription factor binding sites. However, regulatory elements harbored by repeats have proven difficult to characterize using short-read sequencing assays such as ChIP-seq or ATAC-seq. Most regulatory genomics analysis pipelines discard "multimapped" reads that align equally well to multiple genomic locations. Because multimapped reads arise predominantly from repeats, current analysis pipelines fail to detect a substantial portion of regulatory events that occur in repetitive regions. To address this shortcoming, we developed Allo, a new approach to allocate multimapped reads in an efficient, accurate, and user-friendly manner. Allo combines probabilistic mapping of multimapped reads with a convolutional neural network that recognizes the read distribution features of potential peaks, offering enhanced accuracy in multimapping read assignment. Allo also provides read-level output in the form of a corrected alignment file, making it compatible with existing regulatory genomics analysis pipelines and downstream peak-finders. In a demonstration application on CTCF ChIP-seq data, we show that Allo results in the discovery of thousands of new CTCF peaks. Many of these peaks contain the expected cognate motif and/or serve as TAD boundaries. We additionally apply Allo to a diverse collection of ENCODE ChIP-seq data sets, resulting in multiple previously unidentified interactions between transcription factors and repetitive element families. Finally, we show that Allo may be particularly beneficial in identifying ChIP-seq peaks at centromeres, near segmentally duplicated genes, and in younger TEs, enabling new regulatory analyses in these regions.

Population-level transposable element expression dynamics influence trait evolution in a fungal crop pathogen

The rapid adaptive evolution of microbes is driven by strong selection pressure acting on genetic variation. How adaptive genetic variation is generated within species and how such variation influences phenotypic trait expression is often not well understood though. We focused on the recent activity of transposable elements (TEs) using deep population genomics and transcriptomics analyses of a fungal plant pathogen with a highly active content of TEs in the genome. Zymoseptoria tritici causes one of the most damaging diseases on wheat, with recent adaptation to the host and environment being facilitated by TE-associated mutations. We obtained genomic and RNA-sequencing data from 146 isolates collected from a single wheat field. We established a genome-wide map of TE insertion polymorphisms in the population by analyzing recent TE insertions among individuals. We quantified the locus-specific transcription of individual TE copies and found considerable population variation at individual TE loci in the population. About 20% of all TE copies show transcription in the genome suggesting that genomic defenses such as repressive epigenetic marks and repeat-induced polymorphisms are at least partially ineffective at preventing the proliferation of TEs in the genome. A quarter of recent TE insertions are associated with expression variation of neighboring genes providing broad potential to influence trait expression. We indeed found that TE insertions are likely responsible for variation in virulence on the host and potentially diverse components of secondary metabolite production. Our large-scale transcriptomics study emphasizes how TE-derived polymorphisms segregate even in individual microbial populations and can broadly underpin trait variation in pathogens.IMPORTANCEPathogens can rapidly adapt to new hosts, antimicrobials, or changes in the environment. Adaptation arises often from mutations in the genome; however, how such variation is generated remains poorly understood. We investigated the most dynamic regions of the genome of Zymoseptoria tritici, a major fungal pathogen of wheat. We focused on the transcription of transposable elements. A large proportion of the transposable elements not only show signatures of potential activity but are also variable within a single population of the pathogen. We find that this variation in activity is likely influencing many important traits of the pathogen. Hence, our work provides insights into how a microbial species can adapt over the shortest time periods based on the activity of transposable elements.

The role of transposable elements in human evolution and methods for their functional analysis: current status and future perspectives

Transposable elements (TEs) are mobile DNA sequences that can insert themselves into various locations within the genome, causing mutations that may provide advantages or disadvantages to individuals and species. The insertion of TEs can result in genetic variation that may affect a wide range of human traits including genetic disorders. Understanding the role of TEs in human biology is crucial for both evolutionary and medical research. This review discusses the involvement of TEs in human traits and disease susceptibility, as well as methods for functional analysis of TEs.

Allo: Accurate allocation of multi-mapped reads enables regulatory element analysis at repeats

Transposable elements (TEs) and other repetitive regions have been shown to contain gene regulatory elements, including transcription factor binding sites. Unfortunately, regulatory elements harbored by repeats have proven difficult to characterize using short-read sequencing assays such as ChIP-seq or ATAC-seq. Most regulatory genomics analysis pipelines discard "multi-mapped" reads that align equally well to multiple genomic locations. Since multi-mapped reads arise predominantly from repeats, current analysis pipelines fail to detect a substantial portion of regulatory events that occur in repetitive regions. To address this shortcoming, we developed Allo, a new approach to allocate multi-mapped reads in an efficient, accurate, and user-friendly manner. Allo combines probabilistic mapping of multi-mapped reads with a convolutional neural network that recognizes the read distribution features of potential peaks, offering enhanced accuracy in multi-mapping read assignment. Allo also provides read-level output in the form of a corrected alignment file, making it compatible with existing regulatory genomics analysis pipelines and downstream peak-finders. In a demonstration application on CTCF ChIP-seq data, we show that Allo results in the discovery of thousands of new CTCF peaks. Many of these peaks contain the expected cognate motif and/or serve as TAD boundaries. We additionally apply Allo to a diverse collection of ENCODE ChIP-seq datasets, resulting in multiple previously unidentified interactions between transcription factors and repetitive element families. Finally, we show that Allo may be particularly effective in identifying ChIP-seq peaks in younger TEs, which hold evolutionary significance due to their emergence during human evolution from primates.

Chromosomal inversion polymorphisms shape human brain morphology

The impact of chromosomal inversions on human brain morphology remains underexplored. We studied 35 common inversions classified from genotypes of 33,018 adults with European ancestry. The inversions at 2p22.3, 16p11.2, and 17q21.31 reach genome-wide significance, followed by 8p23.1 and 6p21.33, in their association with cortical and subcortical morphology. The 17q21.31, 8p23.1, and 16p11.2 regions comprise the LRRC37, OR7E, and NPIP duplicated gene families. We find the 17q21.31 MAPT inversion region, known for harboring neurological risk, to be the most salient locus among common variants for shaping and patterning the cortex. Overall, we observe the inverted orientations decreasing brain size, with the exception that the 2p22.3 inversion is associated with increased subcortical volume and the 8p23.1 inversion is associated with increased motor cortex. These significant inversions are in the genomic hotspots of neuropsychiatric loci. Our findings are generalizable to 3,472 children and demonstrate inversions as essential genetic variation to understand human brain phenotypes.

A chromosome-scale and haplotype-resolved genome assembly of carnation (Dianthus caryophyllus) based on high-fidelity sequencing

Dianthus caryophyllus is an economic species often considered excellent cut flowers and is suitable for bouquets and gardens. Here, we assembled the haplotype-resolved genome of D. caryophyllus 'Aili' at the chromosome level for the first time. The total lengths of the two assembled haplotypes of carnation were 584.88 Mb for haplotype genome 1 (hap1) and 578.78 Mb for haplotype genome 2 (hap2), respectively. We predicted a total of 44,098 and 42,425 protein-coding genes, respectively. The remarkable structure variation was identified between two haplotypes. Moreover, we identified 403.80 Mb of transposable elements (TEs) in hap1, which accounted for 69.34% of the genome. In contrast, hap2 had 402.70 Mb of TEs, representing 69.61% of the genome. Long terminal repeats were the predominant transposable elements. Phylogenetic analysis showed that the species differentiation time between carnation and gypsophila was estimated to be ~54.43 MYA. The unique gene families of carnation genomes were identified in 'Aili' and previously published 'Francesco' and 'Scarlet Queen'. The assembled and annotated haplotype-resolved D. caryophyllus genome not only promises to facilitate molecular biology studies but also contributes to genome-level evolutionary studies.

Recent Advances in Assembly of Complex Plant Genomes

Over the past 20 years, tremendous advances in sequencing technologies and computational algorithms have spurred plant genomic research into a thriving era with hundreds of genomes decoded already, ranging from those of nonvascular plants to those of flowering plants. However, complex plant genome assembly is still challenging and remains difficult to fully resolve with conventional sequencing and assembly methods due to high heterozygosity, highly repetitive sequences, or high ploidy characteristics of complex genomes. Herein, we summarize the challenges of and advances in complex plant genome assembly, including feasible experimental strategies, upgrades to sequencing technology, existing assembly methods, and different phasing algorithms. Moreover, we list actual cases of complex genome projects for readers to refer to and draw upon to solve future problems related to complex genomes. Finally, we expect that the accurate, gapless, telomere-to-telomere, and fully phased assembly of complex plant genomes could soon become routine.

Structural Dynamics Predominantly Determine the Adaptability of Proteins to Amino Acid Deletions

The insertion or deletion (indel) of amino acids has a variety of effects on protein function, ranging from disease-forming changes to gaining new functions. Despite their importance, indels have not been systematically characterized towards protein engineering or modification goals. In the present work, we focus on deletions composed of multiple contiguous amino acids (mAA-dels) and their effects on the protein (mutant) folding ability. Our analysis reveals that the mutant retains the native fold when the mAA-del obeys well-defined structural dynamics properties: localization in intrinsically flexible regions, showing low resistance to mechanical stress, and separation from allosteric signaling paths. Motivated by the possibility of distinguishing the features that underlie the adaptability of proteins to mAA-dels, and by the rapid evaluation of these features using elastic network models, we developed a positive-unlabeled learning-based classifier that can be adopted for protein design purposes. Trained on a consolidated set of features, including those reflecting the intrinsic dynamics of the regions where the mAA-dels occur, the new classifier yields a high recall of 84.3% for identifying mAA-dels that are stably tolerated by the protein. The comparative examination of the relative contribution of different features to the prediction reveals the dominant role of structural dynamics in enabling the adaptation of the mutant to mAA-del without disrupting the native fold.

Epigenetic and chromosomal features drive transposon insertion in Drosophila melanogaster

Transposons are mobile genetic elements prevalent in the genomes of most species. The distribution of transposons within a genome reflects the actions of two opposing processes: initial insertion site selection, and selective pressure from the host. By analyzing whole-genome sequencing data from transposon-activated Drosophila melanogaster, we identified 43 316 de novo and 237 germline insertions from four long-terminal-repeat (LTR) transposons, one LINE transposon (I-element), and one DNA transposon (P-element). We found that all transposon types favored insertion into promoters de novo, but otherwise displayed distinct insertion patterns. De novo and germline P-element insertions preferred replication origins, often landing in a narrow region around transcription start sites and in regions of high chromatin accessibility. De novo LTR transposon insertions preferred regions with high H3K36me3, promoters and exons of active genes; within genes, LTR insertion frequency correlated with gene expression. De novo I-element insertion density increased with distance from the centromere. Germline I-element and LTR transposon insertions were depleted in promoters and exons, suggesting strong selective pressure to remove transposons from functional elements. Transposon movement is associated with genome evolution and disease; therefore, our results can improve our understanding of genome and disease biology.

How Language and Human Altruism Evolved Hand in Hand - The Backchannel Hypothesis

This paper contributes to two debates: the debate about language evolution and the debate about the foundations of human collaboration. While both cooperation and language may give the impression of being adaptations that evolved for the "good of the group," it is well established that the evolution of complex traits cannot be a direct result of group selection. In this paper I suggest how this tension can be solved: both language and cooperation evolved in a unique two-level evolutionary system which was triggered by a well-documented geological event-the drying out of the climate-in East Africa, which subsequently reduced the intermating between groups and thus made it possible that the mechanism that produced differences between groups (including social forms of selection such as female choice) could be the target of natural selection on the group level. If a social form of selection (e.g., sexual selection) produced differences in fitness between groups, the displacement process between groups would indirectly select those forms of social selection that produce groups that would displace all others. The main hypothesis presented in this paper is that, in this situation, a backchannel between the two levels of selection naturally evolves. A backchannel between the two levels would, for example, emerge when sexual selection (or any other form of social selection) was sensitive to the individual's contribution to the group. Examples of systems utilizing a backchannel are nerve cells being better nourished when used more frequently, enabling them to be conducive to the survival of the whole organism, or a law firm in which all employees get paid to the extent that they contribute to the survival and success of the firm. In both cases, the selection on the higher level informs the selection on the lower level. The aim of the paper is to illuminate these rather opaque claims, to which the reader probably has many objections in this abridged form.

Genomic insights into the adaptation and evolution of the nautilus, an ancient but evolving "living fossil"

The nautilus, commonly known as a "living fossil," is endangered and may be at risk of extinction. The lack of genomic information hinders a thorough understanding of its biology and evolution, which can shed light on the conservation of this endangered species. Here, we report the first high-quality chromosome-level genome assembly of Nautilus pompilius. The assembled genome size comprised 785.15 Mb. Comparative genomic analyses indicated that transposable elements (TEs) and large-scale genome reorganizations may have driven lineage-specific evolution in the cephalopods. Remarkably, evolving conserved genes and recent TE insertion activities were identified in N. pompilius, and we speculate that these findings reflect the strong adaptability and long-term survival of the nautilus. We also identified gene families that are potentially responsible for specific adaptation and evolution events. Our study provides unprecedented insights into the specialized biology and evolution of N. pompilius, and the results serve as an important resource for future conservation genomics of the nautilus and closely related species.

DARTS: An Algorithm for Domain-Associated Retrotransposon Search in Genome Assemblies

Retrotransposons comprise a substantial fraction of eukaryotic genomes, reaching the highest proportions in plants. Therefore, identification and annotation of retrotransposons is an important task in studying the regulation and evolution of plant genomes. The majority of computational tools for mining transposable elements (TEs) are designed for subsequent genome repeat masking, often leaving aside the element lineage classification and its protein domain composition. Additionally, studies focused on the diversity and evolution of a particular group of retrotransposons often require substantial customization efforts from researchers to adapt existing software to their needs. Here, we developed a computational pipeline to mine sequences of protein-coding retrotransposons based on the sequences of their conserved protein domains—DARTS (Domain-Associated Retrotransposon Search). Using the most abundant group of TEs in plants—long terminal repeat (LTR) retrotransposons (LTR-RTs)—we show that DARTS has radically higher sensitivity for LTR-RT identification compared to the widely accepted tool LTRharvest. DARTS can be easily customized for specific user needs. As a result, DARTS returns a set of structurally annotated nucleotide and amino acid sequences which can be readily used in subsequent comparative and phylogenetic analyses. DARTS may facilitate researchers interested in the discovery and detailed analysis of the diversity and evolution of retrotransposons, LTR-RTs, and other protein-coding TEs.

Emerging Roles of Repetitive and Repeat-Containing RNA in Nuclear and Chromatin Organization and Gene Expression

Genomic repeats have been intensely studied as regulatory elements controlling gene transcription, splicing and genome architecture. Our understanding of the role of the repetitive RNA such as the RNA coming from genomic repeats, or repetitive sequences embedded in mRNA/lncRNAs, in nuclear and cellular functions is instead still limited. In this review we discuss evidence supporting the multifaceted roles of repetitive RNA and RNA binding proteins in nuclear organization, gene regulation, and in the formation of dynamic membrane-less aggregates. We hope that our review will further stimulate research in the consolidating field of repetitive RNA biology.

A benchmark and an algorithm for detecting germline transposon insertions and measuring de novo transposon insertion frequencies

Transposons are genomic parasites, and their new insertions can cause instability and spur the evolution of their host genomes. Rapid accumulation of short-read whole-genome sequencing data provides a great opportunity for studying new transposon insertions and their impacts on the host genome. Although many algorithms are available for detecting transposon insertions, the task remains challenging and existing tools are not designed for identifying de novo insertions. Here, we present a new benchmark fly dataset based on PacBio long-read sequencing and a new method TEMP2 for detecting germline insertions and measuring de novo ‘singleton’ insertion frequencies in eukaryotic genomes. TEMP2 achieves high sensitivity and precision for detecting germline insertions when compared with existing tools using both simulated data in fly and experimental data in fly and human. Furthermore, TEMP2 can accurately assess the frequencies of de novo transposon insertions even with high levels of chimeric reads in simulated datasets; such chimeric reads often occur during the construction of short-read sequencing libraries. By applying TEMP2 to published data on hybrid dysgenic flies inflicted by de-repressed P-elements, we confirmed the continuous new insertions of P-elements in dysgenic offspring before they regain piRNAs for P-element repression. TEMP2 is freely available at Github: https://github.com/weng-lab/TEMP2.

Transposable elements, contributors in the evolution of organisms (from an arms race to a source of raw materials)

There is a concept proposing that the primitive lineages of prokaryotes, eukaryotes, and viruses emerged from the primordial pool of primitive genetic elements. In this genetic pool, transposable elements (TEs) became a source of raw material for primitive genomes, tools of genetic innovation, and ancestors of modern genes (e.g. ncRNAs, tRNAs, and rRNAs). TEs contributed directly to the genome evolution of three forms of life on the earth. TEs now appear as tools that were used to giving rise to sexual dimorphism and sex determination, lineage-specific expression of genes and tissue differentiation and finally genome stability and lifespan determination.

Role of Transposable Elements in Gene Regulation in the Human Genome

Transposable elements (TEs), also known as mobile elements (MEs), are interspersed repeats that constitute a major fraction of the genomes of higher organisms. As one of their important functional impacts on gene function and genome evolution, TEs participate in regulating the expression of genes nearby and even far away at transcriptional and post-transcriptional levels. There are two known principal ways by which TEs regulate the expression of genes. First, TEs provide cis-regulatory sequences in the genome with their intrinsic regulatory properties for their own expression, making them potential factors for regulating the expression of the host genes. TE-derived cis-regulatory sites are found in promoter and enhancer elements, providing binding sites for a wide range of trans-acting factors. Second, TEs encode for regulatory RNAs with their sequences showed to be present in a substantial fraction of miRNAs and long non-coding RNAs (lncRNAs), indicating the TE origin of these RNAs. Furthermore, TEs sequences were found to be critical for regulatory functions of these RNAs, including binding to the target mRNA. TEs thus provide crucial regulatory roles by being part of cis-regulatory and regulatory RNA sequences. Moreover, both TE-derived cis-regulatory sequences and TE-derived regulatory RNAs have been implicated in providing evolutionary novelty to gene regulation. These TE-derived regulatory mechanisms also tend to function in a tissue-specific fashion. In this review, we aim to comprehensively cover the studies regarding these two aspects of TE-mediated gene regulation, mainly focusing on the mechanisms, contribution of different types of TEs, differential roles among tissue types, and lineage-specificity, based on data mostly in humans.

Estimating Change in Foldability Due to Multipoint Deletions in Protein Structures

Insertions/deletions of amino acids in the protein backbone potentially result in altered structural/functional specifications. They can either contribute positively to the evolutionary process or can result in disease conditions. Despite being the second most prevalent form of protein modification, there are no databases or computational frameworks that delineate harmful multipoint deletions (MPD) from beneficial ones. We introduce a positive unlabeled learning-based prediction framework (PROFOUND) that utilizes fold-level attributes, environment-specific properties, and deletion site-specific properties to predict the change in foldability arising from such MPDs, both in the non-loop and loop regions of protein structures. In the absence of any protein structure dataset to study MPDs, we introduce a dataset with 153 MPD instances that lead to native-like folded structures and 7650 unlabeled MPD instances whose effect on the foldability of the corresponding proteins is unknown. PROFOUND on 10-fold cross-validation on our newly introduced dataset reports a recall of 82.2% (86.6%) and a fall out rate (FR) of 14.2% (20.6%), corresponding to MPDs in the protein loop (non-loop) region. The low FR suggests that the foldability in proteins subject to MPDs is not random and necessitates unique specifications of the deleted region. In addition, we find that additional evolutionary attributes contribute to higher recall and lower FR. The first of a kind foldability prediction system owing to MPD instances and the newly introduced dataset will potentially aid in novel protein engineering endeavors.

Analysis of Structural Genomic Diversity in Aegilops umbellulata, Ae. markgrafii, Ae. comosa, and Ae. uniaristata by Fluorescence In Situ Hybridization Karyotyping

Fluorescence in situ hybridization karyotypes have been widely used for evolutionary analysis on chromosome organization and genetic/genomic diversity in the wheat alliance (tribe Triticeae of Poaceae). The karyotpic diversity of Aegilops umbellulata, Ae. markgrafii, Ae. comosa subsp. comosa and subsp. subventricosa, and Ae. uniaristata was evaluated by the fluorescence in situ hybridization (FISH) probes oligo-pSc119.2 and pTa71 in combination with (AAC)5, (ACT)7, and (CTT)12, respectively. Abundant intra- and interspecific genetic variation was discovered in Ae. umbellulata, Ae. markgrafii, and Ae. comosa, but not Ae. uniaristata. Chromosome 7 of Ae. umbellulata had more variants (six variants) than the other six U chromosomes (2-3 variants) as revealed by probes oligo-pSc119.2 and (AAC)5. Intraspecific variation in Ae. markgrafii and Ae. comosa was revealed by oligo-pSc119.2 in combination with (ACT)7 and (CTT)12, respectively. At least five variants were found in every chromosome of Ae. markgrafii and Ae. comosa, and up to 18, 10, and 15 variants were identified for chromosomes 2 of Ae. markgrafii, 4 of Ae. comosa subsp. comosa, and 6 of Ae. comosa subsp. subventricosa. The six Ae. uniaristata accessions showed identical FISH signal patterns. A large number of intra-specific polymorphic FISH signals were observed between the homologous chromosomes of Ae. markgrafii and Ae. comosa, especially chromosomes 1, 2, 4, and 7 of Ae. markgrafii, chromosome 4 of Ae. comosa subsp. comosa, and chromosome 6 of Ae. comosa subsp. subventricosa. Twelve Ae. comosa and 24 Ae. markgrafii accessions showed heteromorphism between homologous chromosomes. Additionally, a translocation between the short arms of chromosomes 1 and 7 of Ae. comosa PI 551038 was identified. The FISH karyotypes can be used to clearly identify the chromosome variations of each chromosome in these Aegilops species and also provide valuable information for understanding the evolutionary relationships and structural genomic variation among Aegilops species.

Transposable element-mediated structural variation analysis in dog breeds using whole-genome sequencing

Naturally occurring diseases in dogs provide an important animal model for studying human disease including cancer, heart disease, and autoimmune disorders. Transposable elements (TEs) make up ~ 31% of the dog (Canis lupus familiaris) genome and are one of main drivers to cause genomic variations and alter gene expression patterns of the host genes, which could result in genetic diseases. To detect structural variations (SVs), we conducted whole-genome sequencing of three different breeds, including Maltese, Poodle, and Yorkshire Terrier. Genomic SVs were detected and visualized using BreakDancer program. We identified a total of 2328 deletion SV events in the three breeds compared with the dog reference genome of Boxer. The majority of the genetic variants were found to be TE insertion polymorphism (1229) and the others were TE-mediated deletion (489), non-TE-mediated deletion (542), simple repeat-mediated deletion (32), and other indel (36). Among the TE insertion polymorphism, 286 elements were full-length LINE-1s (L1s). In addition, the 49 SV candidates located in the genic regions were experimentally verified and their polymorphic rates within each breed were examined using PCR assay. Polymorphism analysis of the genomic variants revealed that some of the variants exist polymorphic in the three dog breeds, suggesting that their SV events recently occurred in the dog genome. The findings suggest that TEs have contributed to the genomic variations among the three dog breeds of Maltese, Poodle, and Yorkshire Terrier. In addition, the polymorphic events between the dog breeds indicate that TEs were recently retrotransposed in the dog genome.

Fitness Effects of Single Amino Acid Insertions and Deletions in TEM-1 β-Lactamase

Short insertions and deletions (InDels) are a common type of mutation found in nature and a useful source of variation in protein engineering. InDel events have important consequences in protein evolution, often opening new pathways for adaptation. However, much less is known about the effects of InDels compared to point mutations and amino acid substitutions. In particular, deep mutagenesis studies on the distribution of fitness effects of mutations have focused almost exclusively on amino acid substitutions. Here, we present a near-comprehensive analysis of the fitness effects of single amino acid InDels in TEM-1 β-lactamase. While we found InDels to be largely deleterious, partially overlapping deletion-tolerant and insertion-tolerant regions were observed throughout the protein, especially in unstructured regions and at the end of helices. The signal sequence of TEM-1 tolerated InDels more than the mature protein. Most regions of the protein tolerated insertions more than deletions, but a few regions tolerated deletions more than insertions. We examined the relationship between InDel tolerance and a variety of measures to help understand its origin. These measures included evolutionary variation in β-lactamases, secondary structure identity, tolerance to amino acid substitutions, solvent accessibility, and side-chain weighted contact number. We found secondary structure, weighted contact number, and evolutionary variation in class A beta-lactamases to be the somewhat predictive of InDel fitness effects.

Novel roles of Drosophila FUS and Aub responsible for piRNA biogenesis in neuronal disorders

piRNAs, small non-coding RNAs, were considered to be restricted to germline cells. Although they have recently been detected in somatic cells including neurons, it remains unclear how piRNA biogenesis is involved in neuronal diseases. We herein examined the possible roles of Aubergine (Aub), a Piwi-family protein (PIWI) responsible for piRNA biogenesis, in the neuronal disorders, using the Cabeza (Caz) knockdown Drosophila. Caz is a Drosophila homologue of FUS, which is one of the genes causing amyotrophic lateral sclerosis (ALS). Aub overexpression enhanced the mobility defects accompanied by anatomical defects in motoneurons at neuromuscular junctions induced by the neuron-specific knockdown of Caz. In order to elucidate the underlying mechanisms, we examined pre-piRNA and mature-size piRNA levels under these conditions. qRT-PCR and RNA-seq analyses revealed that the Caz knockdown increased pre-piRNA levels, but reduced mature-size piRNA levels in the central nervous system (CNS), suggesting a role in the pre-piRNAs production. Aub overexpression did not increase mature-size piRNA levels. These results suggest that the accumulated pre-piRNAs are abnormal abortive pre-piRNAs that cannot be further processed by slicers, including Aub. We also demonstrated a relationship between Caz and pre-piRNAs in the CNS by RNA immunoprecipitation. Aub overexpression induced the abnormal cytoplasmic localization of Caz. Based on these results, we propose a model in which Caz knockdown-induced abnormal pre-piRNAs associate with Caz, then translocate and accumulate in the cytoplasm, a process that may be mediated by Aub. The novel roles for Caz and Aub demonstrated herein using the Caz-knockdown fly will contribute to a deeper understanding of the pathogenesis of ALS.

Draft genome of a high value tropical timber tree, Teak (Tectona grandis L. f): insights into SSR diversity, phylogeny and conservation

Teak (Tectona grandis L. f.) is one of the precious bench mark tropical hardwood having qualities of durability, strength and visual pleasantries. Natural teak populations harbour a variety of characteristics that determine their economic, ecological and environmental importance. Sequencing of whole nuclear genome of teak provides a platform for functional analyses and development of genomic tools in applied tree improvement. A draft genome of 317 Mb was assembled at 151× coverage and annotated 36, 172 protein-coding genes. Approximately about 11.18% of the genome was repetitive. Microsatellites or simple sequence repeats (SSRs) are undoubtedly the most informative markers in genotyping, genetics and applied breeding applications. We generated 182,712 SSRs at the whole genome level, of which, 170,574 perfect SSRs were found; 16,252 perfect SSRs showed in silico polymorphisms across six genotypes suggesting their promising use in genetic conservation and tree improvement programmes. Genomic SSR markers developed in this study have high potential in advancing conservation and management of teak genetic resources. Phylogenetic studies confirmed the taxonomic position of the genus Tectona within the family Lamiaceae. Interestingly, estimation of divergence time inferred that the Miocene origin of the Tectona genus to be around 21.4508 million years ago.

Epigenetically facilitated mutational assimilation: epigenetics as a hub within the inclusive evolutionary synthesis

After decades of debate about the existence of non‐genetic inheritance, the focus is now slowly shifting towards dissecting its underlying mechanisms. Here, we propose a new mechanism that, by integrating non‐genetic and genetic inheritance, may help build the long‐sought inclusive vision of evolution. After briefly reviewing the wealth of evidence documenting the existence and ubiquity of non‐genetic inheritance in a table, we review the categories of mechanisms of parent–offspring resemblance that underlie inheritance. We then review several lines of argument for the existence of interactions between non‐genetic and genetic components of inheritance, leading to a discussion of the contrasting timescales of action of non‐genetic and genetic inheritance. This raises the question of how the fidelity of the inheritance system can match the rate of environmental variation. This question is central to understanding the role of different inheritance systems in evolution. We then review and interpret evidence indicating the existence of shifts from inheritance systems with low to higher transmission fidelity. Based on results from different research fields we propose a conceptual hypothesis linking genetic and non‐genetic inheritance systems. According to this hypothesis, over the course of generations, shifts among information systems allow gradual matching between the rate of environmental change and the inheritance fidelity of the corresponding response. A striking conclusion from our review is that documented shifts between types of inherited non‐genetic information converge towards epigenetics (i.e. inclusively heritable molecular variation in gene expression without change in DNA sequence). We then interpret the well‐documented mutagenicity of epigenetic marks as potentially generating a final shift from epigenetic to genetic encoding. This sequence of shifts suggests the existence of a relay in inheritance systems from relatively labile ones to gradually more persistent modes of inheritance, a relay that could constitute a new mechanistic basis for the long‐proposed, but still poorly documented, hypothesis of genetic assimilation. A profound difference between the genocentric and the inclusive vision of heredity revealed by the genetic assimilation relay proposed here lies in the fact that a given form of inheritance can affect the rate of change of other inheritance systems. To explore the consequences of such inter‐connection among inheritance systems, we briefly review published theoretical models to build a model of genetic assimilation focusing on the shift in the engraving of environmentally induced phenotypic variation into the DNA sequence. According to this hypothesis, when environmental change remains stable over a sufficient number of generations, the relay among inheritance systems has the potential to generate a form of genetic assimilation. In this hypothesis, epigenetics appears as a hub by which non‐genetically inherited environmentally induced variation in traits can become genetically encoded over generations, in a form of epigenetically facilitated mutational assimilation. Finally, we illustrate some of the major implications of our hypothetical framework, concerning mutation randomness, the central dogma of molecular biology, concepts of inheritance and the curing of inherited disorders, as well as for the emergence of the inclusive evolutionary synthesis.

piRNAs and PIWI proteins: regulators of gene expression in development and stem cells

PIWI proteins and Piwi-interacting RNAs (piRNAs) have established and conserved roles in repressing transposable elements (TEs) in the germline of animals. However, in several biological contexts, a large proportion of piRNAs are not related to TE sequences and, accordingly, functions for piRNAs and PIWI proteins that are independent of TE regulation have been identified. This aspect of piRNA biology is expanding rapidly. Indeed, recent reports have revealed the role of piRNAs in the regulation of endogenous gene expression programs in germ cells, as well as in somatic tissues, challenging dogma in the piRNA field. In this Review, we focus on recent data addressing the biological and developmental functions of piRNAs, highlighting their roles in embryonic patterning, germ cell specification, stem cell biology, neuronal activity and metabolism.

Insertion sequence elements-mediated structural variations in bacterial genomes

Mobile genetic elements (MGEs) impact the evolution and stability of their host genomes. Insertion sequence (IS) elements are the most common MGEs in bacterial genomes and play a crucial role in mediating large-scale variations in bacterial genomes. It is understood that IS elements and MGEs in general coexist in a dynamical equilibrium with their respective hosts. Current studies indicate that the spontaneous movement of IS elements does not follow a constant rate in different bacterial genomes. However, due to the paucity and sparsity of the data, these observations are yet to be conclusive. In this paper, we conducted a comparative analysis of the IS-mediated genome structural variations in ten mutation accumulation (MA) experiments across eight strains of five bacterial species containing IS elements, including four strains of the E. coli. We used GRASPER algorithm, a denovo structural variation (SV) identification algorithm designed to detect SVs involving repetitive sequences in the genome. We observed highly diverse rates of IS insertions and IS-mediated recombinations across different bacterial species as well as across different strains of the same bacterial species. We also observed different rates of the elements from the same IS family in different bacterial genomes, suggesting that the distinction in rates might not be due to the different composition of IS elements across bacterial genomes.

De novo human genome assemblies reveal spectrum of alternative haplotypes in diverse populations

The human reference genome is used extensively in modern biological research. However, a single consensus representation is inadequate to provide a universal reference structure because it is a haplotype among many in the human population. Using 10× Genomics (10×G) “Linked-Read” technology, we perform whole genome sequencing (WGS) and de novo assembly on 17 individuals across five populations. We identify 1842 breakpoint-resolved non-reference unique insertions (NUIs) that, in aggregate, add up to 2.1 Mb of so far undescribed genomic content. Among these, 64% are considered ancestral to humans since they are found in non-human primate genomes. Furthermore, 37% of the NUIs can be found in the human transcriptome and 14% likely arose from Alu-recombination-mediated deletion. Our results underline the need of a set of human reference genomes that includes a comprehensive list of alternative haplotypes to depict the complete spectrum of genetic diversity across populations.

The majority of the human reference genome assembly is represented as a single consensus haplotype. Here, Wong et al. analyze de novo assemblies of 17 diverse, haplotype-resolved genomes to gain insights into the structure of genetic diversity and compile a list of alternative haplotypes across populations.

Association of zygotic piRNAs derived from paternal P elements with hybrid dysgenesis in Drosophila melanogaster

Background

P-element transposition in the genome causes P-M hybrid dysgenesis in Drosophila melanogaster. Maternally deposited piRNAs suppress P-element transposition in the progeny, linking them to P-M phenotypes; however, the role of zygotic piRNAs derived from paternal P elements is poorly understood.

Results

To elucidate the molecular basis of P-element suppression by zygotic factors, we investigated the genomic constitution and P-element piRNA production derived from fathers. As a result, we characterized males of naturally derived Q, M’ and P strains, which show different capacities for the P-element mobilizations introduced after hybridizations with M-strain females. The amounts of piRNAs produced in ovaries of F1 hybrids varied among the strains and were influenced by the characteristics of the piRNA clusters that harbored the P elements. Importantly, while both the Q- and M’-strain fathers restrict the P-element mobilization in ovaries of their daughters, the Q-strain fathers supported the production of the highest piRNA expression in the ovaries of their daughters, and the M’ strain carries KP elements in transcriptionally active regions directing the highest expression of KP elements in their daughters. Interestingly, the zygotic P-element piRNAs, but not the KP element mRNA, contributed to the variations in P transposition immunity in the granddaughters.

Conclusions

The piRNA-cluster-embedded P elements and the transcriptionally active KP elements from the paternal genome are both important suppressors of P element activities that are co-inherited by the progeny. Expression levels of the P-element piRNA and KP-element mRNA vary among F1 progeny due to the constitution of the paternal genome, and are involved in phenotypic variation in the subsequent generation.

Electronic supplementary material

The online version of this article (10.1186/s13100-018-0110-y) contains supplementary material, which is available to authorized users.

Characterization of the piRNA pathway during development of the sea anemone Nematostella vectensis

PIWI-interacting RNAs (piRNAs) and associated proteins comprise a conserved pathway for silencing transposons in metazoan germlines. piRNA pathway components are also expressed in multipotent somatic stem cells in various organisms. piRNA functions have been extensively explored in bilaterian model systems, however, comprehensive studies in non-bilaterian phyla remain limited. Here we investigate the piRNA pathway during the development of Nematostella vectensis, a well-established model system belonging to Cnidaria, the sister group to Bilateria. To date, no population of somatic stem cells has been identified in this organism, despite its long life-span and regenerative capacities that require a constant cell-renewal. We show that Nematostella piRNA pathway components are broadly expressed in early developmental stages, while piRNAs themselves show differential expression, suggesting specific developmental roles of distinct piRNA families. In adults, piRNA associated proteins are enriched in the germline but also expressed in somatic cells, indicating putative stem cell properties. Furthermore, we provide experimental evidence that Nematostella piRNAs cleave transposable elements as well as protein-coding genes. Our results demonstrate that somatic expression of piRNA associated proteins as well as the roles of piRNAs in transposon repression and gene regulation are likely ancestral features that evolved before the split between Cnidaria and Bilateria.

Diversity of P-element piRNA production among M' and Q strains and its association with P-M hybrid dysgenesis in Drosophila melanogaster

Background

Transposition of P elements in the genome causes P–M hybrid dysgenesis in Drosophila melanogaster. For the P strain, the P–M phenotypes are associated with the ability to express a class of small RNAs, called piwi-interacting small RNAs (piRNAs), that suppress the P elements in female gonads. However, little is known about the extent to which piRNAs are involved in the P–M hybrid dysgenesis in M′ and Q strains, which show different abilities to regulate the P elements from P strains.

Results

To elucidate the molecular basis of the suppression of paternally inherited P elements, we analyzed the mRNA and piRNA levels of P elements in the F1 progeny between males of a P strain and nine-line females of M′ or Q strains (M′ or Q progenies). M′ progenies showed the hybrid dysgenesis phenotype, while Q progenies did not. Consistently, the levels of P-element mRNA in both the ovaries and F1 embryos were higher in M′ progenies than in Q progenies, indicating that the M′ progenies have a weaker ability to suppress P-element expression. The level of P-element mRNA was inversely correlated to the level of piRNAs in F1 embryos. Importantly, the M′ progenies were characterized by a lower abundance of P-element piRNAs in both young ovaries and F1 embryonic bodies. The Q progenies showed various levels of piRNAs in both young ovaries and F1 embryonic bodies despite all of the Q progenies suppressing P-element transposition in their gonad.

Conclusions

Our results are consistent with an idea that the level of P-element piRNAs is a determinant for dividing strain types between M′ and Q and that the suppression mechanisms of transposable elements, including piRNAs, are varied between natural populations.

Electronic supplementary material

The online version of this article (10.1186/s13100-017-0096-x) contains supplementary material, which is available to authorized users.

Genome-wide analysis of LTR-retrotransposon diversity and its impact on the evolution of the genus Helianthus (L.)

Background

Genome divergence by mobile elements activity and recombination is a continuous process that plays a key role in the evolution of species. Nevertheless, knowledge on retrotransposon-related variability among species belonging to the same genus is still limited. Considering the importance of the genus Helianthus, a model system for studying the ecological genetics of speciation and adaptation, we performed a comparative analysis of the repetitive genome fraction across ten species and one subspecies of sunflower, focusing on long terminal repeat retrotransposons at superfamily, lineage and sublineage levels.

Results

After determining the relative genome size of each species, genomic DNA was isolated and subjected to Illumina sequencing. Then, different assembling and clustering approaches allowed exploring the repetitive component of all genomes. On average, repetitive DNA in Helianthus species represented more than 75% of the genome, being composed mostly by long terminal repeat retrotransposons. Also, the prevalence of Gypsy over Copia superfamily was observed and, among lineages, Chromovirus was by far the most represented. Although nearly all the same sublineages are present in all species, we found considerable variability in the abundance of diverse retrotransposon lineages and sublineages, especially between annual and perennial species.

Conclusions

This large variability should indicate that different events of amplification or loss related to these elements occurred following species separation and should have been involved in species differentiation. Our data allowed us inferring on the extent of interspecific repetitive DNA variation related to LTR-RE abundance, investigating the relationship between changes of LTR-RE abundance and the evolution of the genus, and determining the degree of coevolution of different LTR-RE lineages or sublineages between and within species. Moreover, the data suggested that LTR-RE abundance in a species was affected by the annual or perennial habit of that species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4050-6) contains supplementary material, which is available to authorized users.

A noncoding RNA containing a SINE-B1 motif associates with meiotic metaphase chromatin and has an indispensable function during spermatogenesis

A search for early response genes that are activated following germ cell induction from mouse embryonic stem cells in vitro led us to the isolation of a long noncoding RNA that contains a SINE (short interspersed element)-B1F motif that was named R53. In situ hybridization and northern blot analyses revealed that the R53 subfragment RNA bears a B1F motif, is processed from the primary transcript, is expressed in adult testis and is predominantly localized in meiotic metaphase chromatin during spermatogenesis. Recent studies of chromosome-associated RNAs have explored novel functions of noncoding RNAs. Specifically, chromosome-bound noncoding RNAs function not only as structural components of chromosome but also as scaffolds that recruit epigenetic modulators for transcriptional regulation, and they are dynamically rearranged during the cell cycle. However, few studies have explored meiotic chromatin; thus, R53 RNA appears to be the first long noncoding RNA to be tightly associated with the metaphase chromatin during spermatogenesis. Furthermore, R53 knockdown using a lentivirus-mediated RNAi injected into mouse testis and organ culture of the fragments revealed a remarkable reduction in postmeiotic cells and irregular up-regulation of several postmeiotic genes, which suggests the possibility that the SINE-B1-derived noncoding RNA R53 plays an indispensable role in the transcriptional regulation of key spermatogenesis genes.

RNA as a fundamental component of interphase chromosomes: could repeats prove key?

Beginning with the precedent of XIST RNA as a 'chromosomal RNA' (cRNA), there is growing interest in the possibility that a diversity of non-coding RNAs may function in chromatin. We review findings which lead us to suggest that RNA is essentially a widespread component of interphase chromosomes. Further, RNA likely contributes to architecture and regulation, with repeat-rich 'junk' RNA in euchromatin (ecRNA) promoting a more open chromatin state. Thousands of low-abundance nuclear RNAs have been reported, however it remains a challenge to determine which of these may function in chromatin. Recent findings indicate that repetitive sequences are enriched in chromosome-associated non-coding RNAs, and repeat-rich RNA shows unusual properties, including localization and stability, with similarities to XIST RNA. We suggest two frontiers in genome biology are emerging and may intersect: the broad contribution of RNA to interphase chromosomes and the distinctive properties of repeat-rich intronic or intergenic junk sequences that may play a role in chromosome structure and regulation.

A novel method for identifying polymorphic transposable elements via scanning of high-throughput short reads

Identification of polymorphic transposable elements (TEs) is important because TE polymorphism creates genetic diversity and influences the function of genes in the host genome. However, de novo scanning of polymorphic TEs remains a challenge. Here, we report a novel computational method, called PTEMD (polymorphic TEs and their movement detection), for de novo discovery of genome-wide polymorphic TEs. PTEMD searches highly identical sequences using reads supported breakpoint evidences. Using PTEMD, we identified 14 polymorphic TE families (905 sequences) in rice blast fungus Magnaporthe oryzae, and 68 (10,618 sequences) in maize. We validated one polymorphic TE family experimentally, MoTE-1; all MoTE-1 family members are located in different genomic loci in the three tested isolates. We found that 57.1% (8 of 14) of the PTEMD-detected polymorphic TE families in M. oryzae are active. Furthermore, our data indicate that there are more polymorphic DNA transposons in maize than their counterparts of retrotransposons despite the fact that retrotransposons occupy largest fraction of genomic mass. We demonstrated that PTEMD is an effective tool for identifying polymorphic TEs in M. oryzae and maize genomes. PTEMD and the genome-wide polymorphic TEs in M. oryzae and maize are publically available at http://www.kanglab.cn/blast/PTEMD_V1.02.htm.

CRISPR-CAS9 D10A nickase target-specific fluorescent labeling of double strand DNA for whole genome mapping and structural variation analysis

We have developed a new, sequence-specific DNA labeling strategy that will dramatically improve DNA mapping in complex and structurally variant genomic regions, as well as facilitate high-throughput automated whole-genome mapping. The method uses the Cas9 D10A protein, which contains a nuclease disabling mutation in one of the two nuclease domains of Cas9, to create a guide RNA-directed DNA nick in the context of an in vitro-assembled CRISPR-CAS9-DNA complex. Fluorescent nucleotides are then incorporated adjacent to the nicking site with a DNA polymerase to label the guide RNA-determined target sequences. This labeling strategy is very powerful in targeting repetitive sequences as well as in barcoding genomic regions and structural variants not amenable to current labeling methods that rely on uneven distributions of restriction site motifs in the DNA. Importantly, it renders the labeled double-stranded DNA available in long intact stretches for high-throughput analysis in nanochannel arrays as well as for lower throughput targeted analysis of labeled DNA regions using alternative methods for stretching and imaging the labeled long DNA molecules. Thus, this method will dramatically improve both automated high-throughput genome-wide mapping as well as targeted analyses of complex regions containing repetitive and structurally variant DNA.

The peculiar landscape of repetitive sequences in the olive (Olea europaea L.) genome

Analyzing genome structure in different species allows to gain an insight into the evolution of plant genome size. Olive (Olea europaea L.) has a medium-sized haploid genome of 1.4 Gb, whose structure is largely uncharacterized, despite the growing importance of this tree as oil crop. Next-generation sequencing technologies and different computational procedures have been used to study the composition of the olive genome and its repetitive fraction. A total of 2.03 and 2.3 genome equivalents of Illumina and 454 reads from genomic DNA, respectively, were assembled following different procedures, which produced more than 200,000 differently redundant contigs, with mean length higher than 1,000 nt. Mapping Illumina reads onto the assembled sequences was used to estimate their redundancy. The genome data set was subdivided into highly and medium redundant and nonredundant contigs. By combining identification and mapping of repeated sequences, it was established that tandem repeats represent a very large portion of the olive genome (∼31% of the whole genome), consisting of six main families of different length, two of which were first discovered in these experiments. The other large redundant class in the olive genome is represented by transposable elements (especially long terminal repeat-retrotransposons). On the whole, the results of our analyses show the peculiar landscape of the olive genome, related to the massive amplification of tandem repeats, more than that reported for any other sequenced plant genome.

The role of DNA insertions in phenotypic differentiation between humans and other primates

What makes us human is one of the most interesting and enduring questions in evolutionary biology. To assist in answering this question, we have identified insertions in the human genome which cannot be found in five comparison primate species: Chimpanzee, gorilla, orangutan, gibbon, and macaque. A total of 21,269 nonpolymorphic human-specific insertions were identified, of which only 372 were found in exons. Any function conferred by the remaining 20,897 is likely to be regulatory. Many of these insertions are likely to have been fitness neutral; however, a small number has been identified in genes showing signs of positive selection. Insertions found within positively selected genes show associations to neural phenotypes, which were also enriched in the whole data set. Other phenotypes that are found to be enriched in the data set include dental and sensory perception-related phenotypes, features which are known to differ between humans and other apes. The analysis provides several likely candidates, either genes or regulatory regions, which may be involved in the processes that differentiate humans from other apes.

Associations of Haplotypes Upstream of IRS1 with Insulin Resistance, Type 2 Diabetes, Dyslipidemia, Preclinical Atherosclerosis, and Skeletal Muscle LOC646736 mRNA Levels

The genomic region ~500 kb upstream of IRS1 has been implicated in insulin resistance, type 2 diabetes, adverse lipid profile, and cardiovascular risk. To gain further insight into this chromosomal region, we typed four SNPs in a cross-sectional cohort and subjects with type 2 diabetes recruited from the same geographic region. From 16 possible haplotypes, 6 haplotypes with frequencies >0.01 were observed. We identified one haplotype that was protective against insulin resistance (determined by HOMA-IR and fasting plasma insulin levels), type 2 diabetes, an adverse lipid profile, increased C-reactive protein, and asymptomatic atherosclerotic disease (assessed by intima media thickness of the common carotid arteries). BMI and total adipose tissue mass as well as visceral and subcutaneous adipose tissue mass did not differ between the reference and protective haplotypes. In 92 subjects, we observed an association of the protective haplotype with higher skeletal muscle mRNA levels of LOC646736, which is located in the same haplotype block as the informative SNPs and is mainly expressed in skeletal muscle, but only at very low levels in liver or adipose tissues. These data suggest a role for LOC646736 in human insulin resistance and warrant further studies on the functional effects of this locus.

Trend of different molecular markers in the last decades for studying human migrations

Anatomically modern humans are known to have widely migrated throughout history. Different scientific evidences suggest that the entire human population descended from just several thousand African migrants. About 85,000 years ago, the first wave of human migration was out of Africa, that followed the coasts through the Middle East, into Southern Asia via Sri Lanka, and in due course around Indonesia and into Australia. Another wave of migration between 40,000 and 12,000 years ago brought humans northward into Europe. However, the frozen north limited human expansion in Europe, and created a land bridge, "Bering land bridge", connecting Asia with North America about 25,000 years ago. Although fossil data give the most direct information about our past, it has certain anomalies. So, molecular archeologists are now using different molecular markers to trace the "most recent common ancestor" and also the migration pattern of modern humans. In this study, we have studied the trend of molecular markers and also the methodologies implemented in the last decades (2003-2014). From our observation, we can say that D-loop region of mtDNA and Y chromosome based markers are predominant. Nevertheless, mtDNA, especially the D-loop region, has some unique features, which makes it a more effective marker for tracing prehistoric footprints of modern human populations. Although, natural selection should also be taken into account in studying mtDNA based human migration. As per technology is concerned, Sanger sequencing is the major technique that is being used in almost all studies. But, the emergence of different cost-effective-and-easy-to-handle NGS platforms has increased its popularity over Sanger sequencing in studying human migration.

The impact of transposable elements in genome evolution and genetic instability and their implications in various diseases

Approximately 45% of the human genome is comprised of transposable elements (TEs). Results from the Human Genome Project have emphasized the biological importance of TEs. Many studies have revealed that TEs are not simply "junk" DNA, but rather, they play various roles in processes, including genome evolution, gene expression regulation, genetic instability, and cancer disposition. The effects of TE insertion in the genome varies from negligible to disease conditions. For the past two decades, many studies have shown that TEs are the causative factors of various genetic disorders and cancer. TEs are a subject of interest worldwide, not only in terms of their clinical aspects but also in basic research, such as evolutionary tracking. Although active TEs contribute to genetic instability and disease states, non-long terminal repeat transposons are well studied, and their roles in these processes have been confirmed. In this review, we will give an overview of the importance of TEs in studying genome evolution and genetic instability, and we suggest that further in-depth studies on the mechanisms related to these phenomena will be useful for both evolutionary tracking and clinical diagnostics.

Repetitive sequences in plant nuclear DNA: types, distribution, evolution and function

Repetitive DNA sequences are a major component of eukaryotic genomes and may account for up to 90% of the genome size. They can be divided into minisatellite, microsatellite and satellite sequences. Satellite DNA sequences are considered to be a fast-evolving component of eukaryotic genomes, comprising tandemly-arrayed, highly-repetitive and highly-conserved monomer sequences. The monomer unit of satellite DNA is 150-400 base pairs (bp) in length. Repetitive sequences may be species- or genus-specific, and may be centromeric or subtelomeric in nature. They exhibit cohesive and concerted evolution caused by molecular drive, leading to high sequence homogeneity. Repetitive sequences accumulate variations in sequence and copy number during evolution, hence they are important tools for taxonomic and phylogenetic studies, and are known as "tuning knobs" in the evolution. Therefore, knowledge of repetitive sequences assists our understanding of the organization, evolution and behavior of eukaryotic genomes. Repetitive sequences have cytoplasmic, cellular and developmental effects and play a role in chromosomal recombination. In the post-genomics era, with the introduction of next-generation sequencing technology, it is possible to evaluate complex genomes for analyzing repetitive sequences and deciphering the yet unknown functional potential of repetitive sequences.

Genomic landscape of human, bat, and ex vivo DNA transposon integrations

The integration and fixation preferences of DNA transposons, one of the major classes of eukaryotic transposable elements, have never been evaluated comprehensively on a genome-wide scale. Here, we present a detailed study of the distribution of DNA transposons in the human and bat genomes. We studied three groups of DNA transposons that integrated at different evolutionary times: 1) ancient (>40 My) and currently inactive human elements, 2) younger (<40 My) bat elements, and 3) ex vivo integrations of piggyBat and Sleeping Beauty elements in HeLa cells. Although the distribution of ex vivo elements reflected integration preferences, the distribution of human and (to a lesser extent) bat elements was also affected by selection. We used regression techniques (linear, negative binomial, and logistic regression models with multiple predictors) applied to 20-kb and 1-Mb windows to investigate how the genomic landscape in the vicinity of DNA transposons contributes to their integration and fixation. Our models indicate that genomic landscape explains 16-79% of variability in DNA transposon genome-wide distribution. Importantly, we not only confirmed previously identified predictors (e.g., DNA conformation and recombination hotspots) but also identified several novel predictors (e.g., signatures of double-strand breaks and telomere hexamer). Ex vivo integrations showed a bias toward actively transcribed regions. Older DNA transposons were located in genomic regions scarce in most conserved elements-likely reflecting purifying selection. Our study highlights how DNA transposons are integral to the evolution of bat and human genomes, and has implications for the development of DNA transposon assays for gene therapy and mutagenesis applications.

TEMP: a computational method for analyzing transposable element polymorphism in populations

Insertions and excisions of transposable elements (TEs) affect both the stability and variability of the genome. Studying the dynamics of transposition at the population level can provide crucial insights into the processes and mechanisms of genome evolution. Pooling genomic materials from multiple individuals followed by high-throughput sequencing is an efficient way of characterizing genomic polymorphisms in a population. Here we describe a novel method named TEMP, specifically designed to detect TE movements present with a wide range of frequencies in a population. By combining the information provided by pair-end reads and split reads, TEMP is able to identify both the presence and absence of TE insertions in genomic DNA sequences derived from heterogeneous samples; accurately estimate the frequencies of transposition events in the population and pinpoint junctions of high frequency transposition events at nucleotide resolution. Simulation data indicate that TEMP outperforms other algorithms such as PoPoolationTE, RetroSeq, VariationHunter and GASVPro. TEMP also performs well on whole-genome human data derived from the 1000 Genomes Project. We applied TEMP to characterize the TE frequencies in a wild Drosophila melanogaster population and study the inheritance patterns of TEs during hybrid dysgenesis. We also identified sequence signatures of TE insertion and possible molecular effects of TE movements, such as altered gene expression and piRNA production. TEMP is freely available at github: https://github.com/JialiUMassWengLab/TEMP.git.

Nuclear function of Alus

Alus are transposable elements belonging to the short interspersed element family. They occupy over 10% of human genome and have been spreading through genomes over the past 65 million years. In the past, they were considered junk DNA with little function that took up genome volumes. Today, Alus and other transposable elements emerge to be key players in cellular function, including genomic activities, gene expression regulations, and evolution. Here we summarize the current understanding of Alu function in genome and gene expression regulation in human cell nuclei.

The repetitive component of the sunflower genome as shown by different procedures for assembling next generation sequencing reads

Background

Next generation sequencing provides a powerful tool to study genome structure in species whose genomes are far from being completely sequenced. In this work we describe and compare different computational approaches to evaluate the repetitive component of the genome of sunflower, by using medium/low coverage Illumina or 454 libraries.

Results

By varying sequencing technology (Illumina or 454), coverage (0.55 x-1.25 x), assemblers and assembly procedures, six different genomic databases were produced. The annotation of these databases showed that they were composed of different proportions of repetitive DNA families. The final assembly of the sequences belonging to the six databases produced a whole genome set of 283,800 contigs. The redundancy of each contig was estimated by mapping the whole genome set with a large Illumina read set and measuring the number of matched Illumina reads. The repetitive component amounted to 81% of the sunflower genome, that is composed mainly of numerous families of Gypsy and Copia retrotransposons. Also many families of non autonomous retrotransposons and DNA transposons (especially of the Helitron superfamily) were identified.

Conclusions

The results substantially matched those previously obtained by using a Sanger-sequenced shotgun library and a standard 454 whole-genome-shotgun approach, indicating the reliability of the proposed procedures also for other species. The repetitive sequences were collected to produce a database, SUNREP, that will be useful for the annotation of the sunflower genome sequence and for studying the genome evolution in dicotyledons.

CACTA-like transposable element in ZmCCT attenuated photoperiod sensitivity and accelerated the postdomestication spread of maize

The postdomestication adaptation of maize to longer days required reduced photoperiod sensitivity to optimize flowering time. We performed a genome-wide association study and confirmed that ZmCCT, encoding a CCT domain-containing protein, is associated with the photoperiod response. In early-flowering maize we detected a CACTA-like transposable element (TE) within the ZmCCT promoter that dramatically reduced flowering time. TE insertion likely occurred after domestication and was selected as maize adapted to temperate zones. This process resulted in a strong selective sweep within the TE-related block of linkage disequilibrium. Functional validations indicated that the TE represses ZmCCT expression to reduce photoperiod sensitivity, thus accelerating maize spread to long-day environments.