Crypton transposons: identification of new diverse families and ancient domestication events

Evolution of Plant Genome Size and Composition

The rapid development of sequencing technology has led to an explosion of plant genome data, opening up more opportunities for research in the field of comparative evolutionary analysis of plant genomes. In this review, we focus on changes in plant genome size and composition, examining the effects of polyploidy, whole-genome duplication, and alternations in transposable elements on plant genome architecture and evolution, respectively. In addition, to address gaps in the available information, we also collected and analyzed 234 representative plant genome data as a supplement. We aim to provide a comprehensive, up-to-date summary of information on plant genome architecture and evolution in this review.

Structural studies of KCTD1 and its disease-causing mutant P20S provide insights into the protein function and misfunction

Members of the KCTD protein family play key roles in fundamental physio-pathological processes including cancer, neurodevelopmental/neuropsychiatric, and genetic diseases. Here, we report the crystal structure of the KCTD1 P20S mutant, which causes the scalp-ear-nipple syndrome, and molecular dynamics (MD) data on the wild-type protein. Surprisingly, the structure unravels that the N-terminal region, which precedes the BTB domain (preBTB) and bears the disease-associated mutation, adopts a folded polyproline II (PPII) state. The KCTD1 pentamer is characterized by an intricate architecture in which the different subunits mutually exchange domains to generate a closed domain swapping motif. Indeed, the BTB of each chain makes peculiar contacts with the preBTB and the C-terminal domain (CTD) of an adjacent chain. The BTB-preBTB interaction consists of a PPII-PPII recognition motif whereas the BTB-CTD contacts are mediated by an unusual (+/-) helix discontinuous association. The inspection of the protein structure, along with the data emerged from the MD simulations, provides an explanation of the pathogenicity of the P20S mutation and unravels the role of the BTB-preBTB interaction in the insurgence of the disease. Finally, the presence of potassium bound to the central cavity of the CTD pentameric assembly provides insights into the role of KCTD1 in metal homeostasis.

Genome-wide identification and characterization of SLEEPER, a transposon-derived gene family and their expression pattern in Brassica napus L

BACKGROUND

The transposons of the hAT superfamily are the most widespread transposons ever known. SLEEPER genes encode domesticated transposases from the hAT superfamily, which may have lost their transposable functions during long-term evolution and transformed into host proteins that regulate plant growth and development.

RESULTS

This study identified 162 members of the SLEEPER gene family from Brassica napus. These members are widely distributed on 19 chromosomes, mainly in the Cn subgenome, and have promoters with various cis-acting elements related to hormone regulation, abiotic stress, and growth and development regulation. Most of the genes in this family contain similar conserved domains and motifs, and the closer the genes are distributed on evolutionary branches, the more similar their structures are. Transcriptome sequencing performed on tissues at different growth stages from B. napus line 3529 indicated that these genes had different expression patterns, and nearly half of the genes were not detectably expressed in all samples.

CONCLUSIONS

This study investigated the gene structure, expression patterns, evolutionary features, and gene localization of the SLEEPER family members to confirm the significance of these genes in the growth of B. napus, providing a reference for the study of transposon domestication and outstanding genetic resources for the genetic improvement of B. napus.

Integrase-associated niche differentiation of endogenous large DNA viruses in crustaceans

IMPORTANCE

Crustacean genomes harbor sequences originating from a family of large DNA viruses called nimaviruses, but it is unclear why they are present. We show that endogenous nimaviruses selectively insert into repetitive sequences within the host genome, and this insertion specificity was correlated with different types of integrases, which are DNA recombination enzymes encoded by the nimaviruses themselves. This suggests that endogenous nimaviruses have colonized various genomic niches through the acquisition of integrases with different insertion specificities. Our results point to a novel survival strategy of endogenous large DNA viruses colonizing the host genomes. These findings may clarify the evolution and spread of nimaviruses in crustaceans and lead to measures to control and prevent the spread of pathogenic nimaviruses in aquaculture settings.

Zmym4 is required for early cranial gene expression and craniofacial cartilage formation

Introduction: The Six1 transcription factor plays important roles in the development of cranial sensory organs, and point mutations underlie craniofacial birth defects. Because Six1's transcriptional activity can be modulated by interacting proteins, we previously screened for candidate interactors and identified zinc-finger MYM-containing protein 4 (Zmym4) by its inclusion of a few domains with a bona fide cofactor, Sine oculis binding protein (Sobp). Although Zmym4 has been implicated in regulating early brain development and certain cancers, its role in craniofacial development has not previously been described. Methods: We used co-immunoprecipitation and luciferase-reporter assays in cultured cells to test interactions between Zmym4 and Six1. We used knock-down and overexpression of Zmym4 in embryos to test for its effects on early ectodermal gene expression, neural crest migration and craniofacial cartilage formation. Results: We found no evidence that Zmym4 physically or transcriptionally interacts with Six1 in cultured cells. Nonetheless, knockdown of endogenous Zmym4 in embryos resulted in altered early cranial gene expression, including those expressed in the neural border, neural plate, neural crest and preplacodal ectoderm. Experimentally increasing Zmym4 levels had minor effects on neural border or neural plate genes, but altered the expression of neural crest and preplacodal genes. At larval stages, genes expressed in the otic vesicle and branchial arches showed reduced expression in Zmym4 morphants. Although we did not detect defects in neural crest migration into the branchial arches, loss of Zmym4 resulted in aberrant morphology of several craniofacial cartilages. Discussion: Although Zmym4 does not appear to function as a Six1 transcriptional cofactor, it plays an important role in regulating the expression of embryonic cranial genes in tissues critical for normal craniofacial development.

Repetitive DNA sequence detection and its role in the human genome

Repetitive DNA sequences playing critical roles in driving evolution, inducing variation, and regulating gene expression. In this review, we summarized the definition, arrangement, and structural characteristics of repeats. Besides, we introduced diverse biological functions of repeats and reviewed existing methods for automatic repeat detection, classification, and masking. Finally, we analyzed the type, structure, and regulation of repeats in the human genome and their role in the induction of complex diseases. We believe that this review will facilitate a comprehensive understanding of repeats and provide guidance for repeat annotation and in-depth exploration of its association with human diseases.

A novel variant in the QRICH1 gene was identified in a patient with severe developmental delay

BACKGROUND

QRICH1 encodes the glutamine-rich protein 1, which contains one caspase activation recruitment domain and is likely to be involved in apoptosis and inflammation. However, the function of the QRICH1 gene was largely unknown. Recently, several studies have reported de novo variants in QRICH1, and the variants have been associated with Ververi-Brady syndrome characterized by developmental delay, nonspecific facial dysmorphism, and hypotonia.

MATERIALS AND METHODS

Whole exome sequencing, clinical examinations, and functional experiments were performed to identify the etiology of our patient.

RESULTS

Here, we added another patient with severe growth retardation, atrial septal defect, and slurred speech. Whole exome sequencing identified a novel truncation variant in the QRICH1 gene (MN_017730.3: c.1788dupC, p.Tyr597Leufs*9). Furthermore, the functional experiments confirmed the effect of genetic variation.

CONCLUSION

Our findings expand the QRICH1 variant spectrum in developmental disorders and provide evidence for the application of whole exome sequencing in Ververi-Brady syndrome.

Unanticipated broad phylogeny of BEN DNA-binding domains revealed by structural homology searches

Organization of protein sequences into domain families is a foundation for cataloging and investigating protein functions. However, long-standing strategies based on primary amino acid sequences are blind to the possibility that proteins with dissimilar sequences could have comparable tertiary structures. Building on our recent findings that in silico structural predictions of BEN family DNA-binding domains closely resemble their experimentally determined crystal structures, we exploited the AlphaFold2 database for comprehensive identification of BEN domains. Indeed, we identified numerous novel BEN domains, including members of new subfamilies. For example, while no BEN domain factors had previously been annotated in C. elegans, this species actually encodes multiple BEN proteins. These include key developmental timing genes of orphan domain status, sel-7 and lin-14, the latter being the central target of the founding miRNA lin-4. We also reveal that the domain of unknown function 4806 (DUF4806), which is widely distributed across metazoans, is structurally similar to BEN and comprises a new subtype. Surprisingly, we find that BEN domains resemble both metazoan and non-metazoan homeodomains in 3D conformation and preserve characteristic residues, indicating that despite their inability to be aligned by conventional methods, these DNA-binding modules are probably evolutionarily related. Finally, we broaden the application of structural homology searches by revealing novel human members of DUF3504, which exists on diverse proteins with presumed or known nuclear functions. Overall, our work strongly expands this recently identified family of transcription factors and illustrates the value of 3D structural predictions to annotate protein domains and interpret their functions.

Population-scale long-read sequencing uncovers transposable elements associated with gene expression variation and adaptive signatures in Drosophila

High quality reference genomes are crucial to understanding genome function, structure and evolution. The availability of reference genomes has allowed us to start inferring the role of genetic variation in biology, disease, and biodiversity conservation. However, analyses across organisms demonstrate that a single reference genome is not enough to capture the global genetic diversity present in populations. In this work, we generate 32 high-quality reference genomes for the well-known model species D. melanogaster and focus on the identification and analysis of transposable element variation as they are the most common type of structural variant. We show that integrating the genetic variation across natural populations from five climatic regions increases the number of detected insertions by 58%. Moreover, 26% to 57% of the insertions identified using long-reads were missed by short-reads methods. We also identify hundreds of transposable elements associated with gene expression variation and new TE variants likely to contribute to adaptive evolution in this species. Our results highlight the importance of incorporating the genetic variation present in natural populations to genomic studies, which is essential if we are to understand how genomes function and evolve.

A beginner's guide to manual curation of transposable elements

BACKGROUND

In the study of transposable elements (TEs), the generation of a high confidence set of consensus sequences that represent the diversity of TEs found in a given genome is a key step in the path to investigate these fascinating genomic elements. Many algorithms and pipelines are available to automatically identify putative TE families present in a genome. Despite the availability of these valuable resources, producing a library of high-quality full-length TE consensus sequences largely remains a process of manual curation. This know-how is often passed on from mentor-to-mentee within research groups, making it difficult for those outside the field to access this highly specialised skill.

RESULTS

Our manuscript attempts to fill this gap by providing a set of detailed computer protocols, software recommendations and video tutorials for those aiming to manually curate TEs. Detailed step-by-step protocols, aimed at the complete beginner, are presented in the Supplementary Methods.

CONCLUSIONS

The proposed set of programs and tools presented here will make the process of manual curation achievable and amenable to all researchers and in special to those new to the field of TEs.

Paleozoic Protein Fossils Illuminate the Evolution of Vertebrate Genomes and Transposable Elements

Genomes hold a treasure trove of protein fossils: fragments of formerly protein-coding DNA, which mainly come from transposable elements (TEs) or host genes. These fossils reveal ancient evolution of TEs and genomes, and many fossils have been exapted to perform diverse functions important for the host's fitness. However, old and highly-degraded fossils are hard to identify, standard methods (e.g. BLAST) are not optimized for this task, and few Paleozoic protein fossils have been found. Here, a recently optimized method is used to find protein fossils in vertebrate genomes. It finds Paleozoic fossils predating the amphibian/amniote divergence from most major TE categories, including virus-related Polinton and Gypsy elements. It finds 10 fossils in the human genome (8 from TEs and 2 from host genes) that predate the last common ancestor of all jawed vertebrates, probably from the Ordovician period. It also finds types of transposon and retrotransposon not found in human before. These fossils have extreme sequence conservation, indicating exaptation: some have evidence of gene-regulatory function, and they tend to lienearest to developmental genes. Some ancient fossils suggest "genome tectonics", where two fragments of one TE have drifted apart by up to megabases, possibly explaining gene deserts and large introns. This paints a picture of great TE diversity in our aquatic ancestors, with patchy TE inheritance by later vertebrates, producing new genes and regulatory elements on the way. Host-gene fossils too have contributed anciently-conserved DNA segments. This paves the way to further studies of ancient protein fossils.

The clinical and molecular spectrum of QRICH1 associated neurodevelopmental disorder

De novo variants in QRICH1 (Glutamine-rich protein 1) has recently been reported in 11 individuals with intellectual disability (ID). The function of QRICH1 is largely unknown but it is likely to play a key role in the unfolded response of endoplasmic reticulum stress through transcriptional control of proteostasis. In this study, we present 27 additional individuals and delineate the clinical and molecular spectrum of the individuals (n = 38) with QRICH1 variants. The main clinical features were mild to moderate developmental delay/ID (71%), nonspecific facial dysmorphism (92%) and hypotonia (39%). Additional findings included poor weight gain (29%), short stature (29%), autism spectrum disorder (29%), seizures (24%) and scoliosis (18%). Minor structural brain abnormalities were reported in 52% of the individuals with brain imaging. Truncating or splice variants were found in 28 individuals and 10 had missense variants. Four variants were inherited from mildly affected parents. This study confirms that heterozygous QRICH1 variants cause a neurodevelopmental disorder including short stature and expands the phenotypic spectrum to include poor weight gain, scoliosis, hypotonia, minor structural brain anomalies, and seizures. Inherited variants from mildly affected parents are reported for the first time, suggesting variable expressivity.

Expression and Mutation Alterations of ZMYM4 Gene in Gastric and Colonic Cancers

ZMYM4 is a zinc finger protein, whose cancer-related functions are partially known (cell shape maintenance and cell death). In this study, we analyzed 4 sites of mononucleotide repeats in the coding sequence of ZMYM4 in gastric (GC) and colonic cancers (CC). Seven of the 32 high microsatellite instability (MSI-H) GCs (21.9%) and 23 of 113 MSI-H CCs (20.4%) harbored ZMYM4 frameshift mutations with no significant difference between the 2 organs (P>0.05). There was no ZMYM4 frameshift mutations in microsatellite-stable GCs and CCs. We also identified that 6 of 16 MSI-H CCs (37.5%) exhibited intratumoral heterogeneity of the ZMYM4 frameshift mutations. In both GC and CC with MSI-H, ZMYM4 expression in ZMYM4-mutated cases was significantly lower than that in ZMYM4-nonmutated cases. Our study indicates that ZMYM4 is altered at multiple levels (frameshift mutation, mutational intratumoral heterogeneity, and loss of expression), suggesting their relations with MSI-H GC and CC.

The Structural, Functional and Evolutionary Impact of Transposable Elements in Eukaryotes

Transposable elements (TEs) are nearly ubiquitous in eukaryotes. The increase in genomic data, as well as progress in genome annotation and molecular biology techniques, have revealed the vast number of ways mobile elements have impacted the evolution of eukaryotes. In addition to being the main cause of difference in haploid genome size, TEs have affected the overall organization of genomes by accumulating preferentially in some genomic regions, by causing structural rearrangements or by modifying the recombination rate. Although the vast majority of insertions is neutral or deleterious, TEs have been an important source of evolutionary novelties and have played a determinant role in the evolution of fundamental biological processes. TEs have been recruited in the regulation of host genes and are implicated in the evolution of regulatory networks. They have also served as a source of protein-coding sequences or even entire genes. The impact of TEs on eukaryotic evolution is only now being fully appreciated and the role they may play in a number of biological processes, such as speciation and adaptation, remains to be deciphered.

The Enterprise, a massive transposon carrying Spok meiotic drive genes

The genomes of eukaryotes are full of parasitic sequences known as transposable elements (TEs). Here, we report the discovery of a putative giant tyrosine-recombinase-mobilized DNA transposon, Enterprise, from the model fungus Podospora anserina Previously, we described a large genomic feature called the Spok block which is notable due to the presence of meiotic drive genes of the Spok gene family. The Spok block ranges from 110 kb to 247 kb and can be present in at least four different genomic locations within P. anserina, despite what is an otherwise highly conserved genome structure. We propose that the reason for its varying positions is that the Spok block is not only capable of meiotic drive but is also capable of transposition. More precisely, the Spok block represents a unique case where the Enterprise has captured the Spoks, thereby parasitizing a resident genomic parasite to become a genomic hyperparasite. Furthermore, we demonstrate that Enterprise (without the Spoks) is found in other fungal lineages, where it can be as large as 70 kb. Lastly, we provide experimental evidence that the Spok block is deleterious, with detrimental effects on spore production in strains which carry it. This union of meiotic drivers and a transposon has created a selfish element of impressive size in Podospora, challenging our perception of how TEs influence genome evolution and broadening the horizons in terms of what the upper limit of transposition may be.

A case of Ververi-Brady syndrome due to QRICH1 loss of function and the literature review

Ververi-Brady syndrome (VBS), first reported in 2018, is characterized by intellectual disability, speech delay, and mild dysmorphic facial features. VBS has been linked to de novo loss-of-function variants in the glutamine-rich protein 1 (QRICH1) on chromosome 3p21 and was reported until lately in only five individuals. Four additional cases have just been described substantiating the notion that children with VBS are mildly dysmorphic, mildly to moderately intellectually disabled, have linear growth shortage, are picky eaters, and have notable attention and social behavioral deficits. We describe a new patient and review the clinical and genetic information, on all previously reported VBS cases. The child here reported is noted for maladaptive behavior, sensory hypersensitivity, and slow linear growth. He is mainly hyperactive, distractible, impulsive, and inattentive. His speech, initially delayed, is fair and his verbal comprehension age adequate.

Recurrent evolution of vertebrate transcription factors by transposase capture

Genes with novel cellular functions may evolve through exon shuffling, which can assemble novel protein architectures. Here, we show that DNA transposons provide a recurrent supply of materials to assemble protein-coding genes through exon shuffling. We find that transposase domains have been captured-primarily via alternative splicing-to form fusion proteins at least 94 times independently over the course of ~350 million years of tetrapod evolution. We find an excess of transposase DNA binding domains fused to host regulatory domains, especially the Krüppel-associated box (KRAB) domain, and identify four independently evolved KRAB-transposase fusion proteins repressing gene expression in a sequence-specific fashion. The bat-specific KRABINER fusion protein binds its cognate transposons genome-wide and controls a network of genes and cis-regulatory elements. These results illustrate how a transcription factor and its binding sites can emerge.

QRICH1 dictates the outcome of ER stress through transcriptional control of proteostasis

Tissue homeostasis is perturbed in a diversity of inflammatory pathologies. These changes can elicit endoplasmic reticulum (ER) stress, protein misfolding, and cell death. ER stress triggers the unfolded protein response (UPR), which can promote recovery of ER proteostasis and cell survival or trigger programmed cell death. Here, we leveraged single-cell RNA sequencing to define dynamic transcriptional states associated with the adaptive versus terminal UPR in the mouse intestinal epithelium. We integrated these transcriptional programs with genome-scale CRISPR screening to dissect the UPR pathway functionally. We identified QRICH1 as a key effector of the PERK-eIF2α axis of the UPR. QRICH1 controlled a transcriptional program associated with translation and secretory networks that were specifically up-regulated in inflammatory pathologies. Thus, QRICH1 dictates cell fate in response to pathological ER stress.

A Field Guide to Eukaryotic Transposable Elements

Transposable elements (TEs) are mobile DNA sequences that propagate within genomes. Through diverse invasion strategies, TEs have come to occupy a substantial fraction of nearly all eukaryotic genomes, and they represent a major source of genetic variation and novelty. Here we review the defining features of each major group of eukaryotic TEs and explore their evolutionary origins and relationships. We discuss how the unique biology of different TEs influences their propagation and distribution within and across genomes. Environmental and genetic factors acting at the level of the host species further modulate the activity, diversification, and fate of TEs, producing the dramatic variation in TE content observed across eukaryotes. We argue that cataloging TE diversity and dissecting the idiosyncratic behavior of individual elements are crucial to expanding our comprehension of their impact on the biology of genomes and the evolution of species.

Exome Sequencing Identifies Genes and Gene Sets Contributing to Severe Childhood Obesity, Linking PHIP Variants to Repressed POMC Transcription

Obesity is genetically heterogeneous with monogenic and complex polygenic forms. Using exome and targeted sequencing in 2,737 severely obese cases and 6,704 controls, we identified three genes (PHIP, DGKI, and ZMYM4) with an excess burden of very rare predicted deleterious variants in cases. In cells, we found that nuclear PHIP (pleckstrin homology domain interacting protein) directly enhances transcription of pro-opiomelanocortin (POMC), a neuropeptide that suppresses appetite. Obesity-associated PHIP variants repressed POMC transcription. Our demonstration that PHIP is involved in human energy homeostasis through transcriptional regulation of central melanocortin signaling has potential diagnostic and therapeutic implications for patients with obesity and developmental delay. Additionally, we found an excess burden of predicted deleterious variants involving genes nearest to loci from obesity genome-wide association studies. Genes and gene sets influencing obesity with variable penetrance provide compelling evidence for a continuum of causality in the genetic architecture of obesity, and explain some of its missing heritability.

Study of VIPER and TATE in kinetoplastids and the evolution of tyrosine recombinase retrotransposons

Background

Kinetoplastids are a flagellated group of protists, including some parasites, such as Trypanosoma and Leishmania species, that can cause diseases in humans and other animals. The genomes of these species enclose a fraction of retrotransposons including VIPER and TATE, two poorly studied transposable elements that encode a tyrosine recombinase (YR) and were previously classified as DIRS elements. This study investigated the distribution and evolution of VIPER and TATE in kinetoplastids to understand the relationships of these elements with other retrotransposons.

Results

We observed that VIPER and TATE have a discontinuous distribution among Trypanosomatidae, with several events of loss and degeneration occurring during a vertical transfer evolution. We were able to identify the terminal repeats of these elements for the first time, and we showed that these elements are potentially active in some species, including T. cruzi copies of VIPER. We found that VIPER and TATE are strictly related elements, which were named in this study as VIPER-like. The reverse transcriptase (RT) tree presented a low resolution, and the origin and relationships among YR groups remain uncertain. Conversely, for RH, VIPER-like grouped with Hepadnavirus, whereas for YR, VIPER-like sequences constituted two different clades that are closely allied to Crypton. Distinct topologies among RT, RH and YR trees suggest ancient rearrangements/exchanges in domains and a modular pattern of evolution with putative independent origins for each ORF.

Conclusions

Due to the presence of both elements in Bodo saltans, a nontrypanosomatid species, we suggested that VIPER and TATE have survived and remained active for more than 400 million years or were reactivated during the evolution of the host species. We did not find clear evidence of independent origins of VIPER-like from the other YR retroelements, supporting the maintenance of the DIRS group of retrotransposons. Nevertheless, according to phylogenetic findings and sequence structure obtained by this study and other works, we proposed separating DIRS elements into four subgroups: DIRS-like, PAT-like, Ngaro-like, and VIPER-like.

Electronic supplementary material

The online version of this article (10.1186/s13100-019-0175-2) contains supplementary material, which is available to authorized users.

Rpp1 Encodes a ULP1-NBS-LRR Protein That Controls Immunity to Phakopsora pachyrhizi in Soybean

Phakopsora pachyrhizi is the causal agent of Asian soybean rust. Susceptible soybean plants infected by virulent isolates of P. pachyrhizi are characterized by tan-colored lesions and erumpent uredinia on the leaf surface. Germplasm screening and genetic analyses have led to the identification of seven loci, Rpp1 to Rpp7, that provide varying degrees of resistance to P. pachyrhizi (Rpp). Two genes, Rpp1 and Rpp1b, map to the same region on soybean chromosome 18. Rpp1 is unique among the Rpp genes in that it confers an immune response (IR) to avirulent P. pachyrhizi isolates. The IR is characterized by a lack of visible symptoms, whereas resistance provided by Rpp1b to Rpp7 results in red-brown foliar lesions. Rpp1 maps to a region spanning approximately 150 kb on chromosome 18 between markers Sct_187 and Sat_064 in L85-2378 (Rpp1), an isoline developed from Williams 82 and PI 200492 (Rpp1). To identify Rpp1, we constructed a bacterial artificial chromosome library from soybean accession PI 200492. Sequencing of the Rpp1 locus identified three homologous nucleotide binding site-leucine rich repeat (NBS-LRR) candidate resistance genes between Sct_187 and Sat_064. Each candidate gene is also predicted to encode an N-terminal ubiquitin-like protease 1 (ULP1) domain. Cosilencing of the Rpp1 candidates abrogated the immune response in the Rpp1 resistant soybean accession PI 200492, indicating that Rpp1 is a ULP1-NBS-LRR protein and plays a key role in the IR.

Transposable Elements: Classification, Identification, and Their Use As a Tool For Comparative Genomics

Most genomes are populated by hundreds of thousands of sequences originated from mobile elements. On the one hand, these sequences present a real challenge in the process of genome analysis and annotation. On the other hand, they are very interesting biological subjects involved in many cellular processes. Here we present an overview of transposable elements biodiversity, and we discuss different approaches to transposable elements detection and analyses.

Childhood asthma is associated with COPD and known asthma variants in COPDGene: a genome-wide association study

BACKGROUND

Childhood asthma is strongly influenced by genetics and is a risk factor for reduced lung function and chronic obstructive pulmonary disease (COPD) in adults. This study investigates self-reported childhood asthma in adult smokers from the COPDGene Study. We hypothesize that childhood asthma is associated with decreased lung function, increased risk for COPD, and that a genome-wide association study (GWAS) will show association with established asthma variants.

METHODS

We evaluated current and former smokers ages 45-80 of non-Hispanic white (NHW) or African American (AA) race. Childhood asthma was defined by self-report of asthma, diagnosed by a medical professional, with onset at < 16 years or during childhood. Subjects with a history of childhood asthma were compared to those who never had asthma based on lung function, development of COPD, and genetic variation. GWAS was performed in NHW and AA populations, and combined in meta-analysis. Two sets of established asthma SNPs from published literature were examined for association with childhood asthma.

RESULTS

Among 10,199 adult smokers, 730 (7%) reported childhood asthma and 7493 (73%) reported no history of asthma. Childhood asthmatics had reduced lung function and increased risk for COPD (OR 3.42, 95% CI 2.81-4.18). Genotype data was assessed for 8031 subjects. Among NHWs, 391(7%) had childhood asthma, and GWAS identified one genome-wide significant association in KIAA1958 (rs59289606, p = 4.82 × 10- 8). Among AAs, 339 (12%) had childhood asthma. No SNPs reached genome-wide significance in the AAs or in the meta-analysis combining NHW and AA subjects; however, potential regions of interest were identified. Established asthma SNPs were examined, seven from the NHGRI-EBI database and five with genome-wide significance in the largest pediatric asthma GWAS. Associations were found in the current childhood asthma GWAS with known asthma loci in IL1RL1, IL13, LINC01149, near GSDMB, and in the C11orf30-LRRC32 region (Bonferroni adjusted p < 0.05 for all comparisons).

CONCLUSIONS

Childhood asthmatics are at increased risk for COPD. Defining asthma by self-report is valid in populations at risk for COPD, identifying subjects with clinical and genetic characteristics known to associate with childhood asthma. This has potential to improve clinical understanding of asthma-COPD overlap (ACO) and enhance future research into ACO-specific treatment regimens.

TRIAL REGISTRATION

ClinicalTrials.gov, NCT00608764 (Active since January 28, 2008).

The Role of Transposable Elements in Speciation

Understanding the phenotypic and molecular mechanisms that contribute to genetic diversity between and within species is fundamental in studying the evolution of species. In particular, identifying the interspecific differences that lead to the reduction or even cessation of gene flow between nascent species is one of the main goals of speciation genetic research. Transposable elements (TEs) are DNA sequences with the ability to move within genomes. TEs are ubiquitous throughout eukaryotic genomes and have been shown to alter regulatory networks, gene expression, and to rearrange genomes as a result of their transposition. However, no systematic effort has evaluated the role of TEs in speciation. We compiled the evidence for TEs as potential causes of reproductive isolation across a diversity of taxa. We find that TEs are often associated with hybrid defects that might preclude the fusion between species, but that the involvement of TEs in other barriers to gene flow different from postzygotic isolation is still relatively unknown. Finally, we list a series of guides and research avenues to disentangle the effects of TEs on the origin of new species.

Hologenomic adaptations underlying the evolution of sanguivory in the common vampire bat

Adaptation to specialized diets often requires modifications at both genomic and microbiome levels. We applied a hologenomic approach to the common vampire bat (Desmodus rotundus), one of the only three obligate blood-feeding (sanguivorous) mammals, to study the evolution of its complex dietary adaptation. Specifically, we assembled its high-quality reference genome (scaffold N50 = 26.9 Mb, contig N50 = 36.6 kb) and gut metagenome, and compared them against those of insectivorous, frugivorous and carnivorous bats. Our analyses showed a particular common vampire bat genomic landscape regarding integrated viral elements, a dietary and phylogenetic influence on gut microbiome taxonomic and functional profiles, and that both genetic elements harbour key traits related to the nutritional (for example, vitamin and lipid shortage) and non-nutritional (for example, nitrogen waste and osmotic homeostasis) challenges of sanguivory. These findings highlight the value of a holistic study of both the host and its microbiota when attempting to decipher adaptations underlying radical dietary lifestyles.

Human transposable elements in Repbase: genomic footprints from fish to humans

Repbase is a comprehensive database of eukaryotic transposable elements (TEs) and repeat sequences, containing over 1300 human repeat sequences. Recent analyses of these repeat sequences have accumulated evidences for their contribution to human evolution through becoming functional elements, such as protein-coding regions or binding sites of transcriptional regulators. However, resolving the origins of repeat sequences is a challenge, due to their age, divergence, and degradation. Ancient repeats have been continuously classified as TEs by finding similar TEs from other organisms. Here, the most comprehensive picture of human repeat sequences is presented. The human genome contains traces of 10 clades (L1, CR1, L2, Crack, RTE, RTEX, R4, Vingi, Tx1 and Penelope) of non-long terminal repeat (non-LTR) retrotransposons (long interspersed elements, LINEs), 3 types (SINE1/7SL, SINE2/tRNA, and SINE3/5S) of short interspersed elements (SINEs), 1 composite retrotransposon (SVA) family, 5 classes (ERV1, ERV2, ERV3, Gypsy and DIRS) of LTR retrotransposons, and 12 superfamilies (Crypton, Ginger1, Harbinger, hAT, Helitron, Kolobok, Mariner, Merlin, MuDR, P, piggyBac and Transib) of DNA transposons. These TE footprints demonstrate an evolutionary continuum of the human genome.

Standardized, systemic phenotypic analysis reveals kidney dysfunction as main alteration of Kctd1 I27N mutant mice

Background

Increased levels of blood plasma urea were used as phenotypic parameter for establishing novel mouse models for kidney diseases on the genetic background of C3H inbred mice in the phenotype-driven Munich ENU mouse mutagenesis project. The phenotypically dominant mutant line HST014 was established and further analyzed.

Methods

Analysis of the causative mutation as well as the standardized, systemic phenotypic analysis of the mutant line was carried out.

Results

The causative mutation was detected in the potassium channel tetramerization domain containing 1 (Kctd1) gene which leads to the amino acid exchange Kctd1^I27N thereby affecting the functional BTB domain of the protein. This line is the first mouse model harboring a Kctd1 mutation. Kctd1^I27N homozygous mutant mice die perinatally. Standardized, systemic phenotypic analysis of Kctd1^I27N heterozygous mutants was carried out in the German Mouse Clinic (GMC). Systematic morphological investigation of the external physical appearance did not detect the specific alterations that are described in KCTD1 mutant human patients affected by the scalp-ear-nipple (SEN) syndrome. The main pathological phenotype of the Kctd1^I27N heterozygous mutant mice consists of kidney dysfunction and secondary effects thereof, without gross additional primary alterations in the other phenotypic parameters analyzed. Genome-wide transcriptome profiling analysis at the age of 4 months revealed about 100 differentially expressed genes (DEGs) in kidneys of Kctd1^I27N heterozygous mutants as compared to wild-type controls.

Conclusions

In summary, the main alteration of the Kctd1^I27N heterozygous mutants consists in kidney dysfunction. Additional analyses in 9–21 week-old heterozygous mutants revealed only few minor effects.

Electronic supplementary material

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

ZMYM3 regulates BRCA1 localization at damaged chromatin to promote DNA repair

In this study, Leung et al. identified ZMYM3 (zinc finger, myeloproliferative, and mental retardation-type 3) as a chromatin-interacting protein that promotes DNA repair by homologous recombination. This work identifies a critical chromatin-binding DNA damage response factor, ZMYM3, which modulates BRCA1 functions within chromatin to ensure the maintenance of genome integrity.

Chromatin connects DNA damage response factors to sites of damaged DNA to promote the signaling and repair of DNA lesions. The histone H2A variants H2AX, H2AZ, and macroH2A represent key chromatin constituents that facilitate DNA repair. Through proteomic screening of these variants, we identified ZMYM3 (zinc finger, myeloproliferative, and mental retardation-type 3) as a chromatin-interacting protein that promotes DNA repair by homologous recombination (HR). ZMYM3 is recruited to DNA double-strand breaks through bivalent interactions with both histone and DNA components of the nucleosome. We show that ZMYM3 links the HR factor BRCA1 to damaged chromatin through specific interactions with components of the BRCA1-A subcomplex, including ABRA1 and RAP80. By regulating ABRA1 recruitment to damaged chromatin, ZMYM3 facilitates the fine-tuning of BRCA1 interactions with DNA damage sites and chromatin. Consistent with a role in regulating BRCA1 function, ZMYM3 deficiency results in impaired HR repair and genome instability. Thus, our work identifies a critical chromatin-binding DNA damage response factor, ZMYM3, which modulates BRCA1 functions within chromatin to ensure the maintenance of genome integrity.

Genetic Drift, Not Life History or RNAi, Determine Long-Term Evolution of Transposable Elements

Transposable elements (TEs) are a major source of genome variation across the branches of life. Although TEs may play an adaptive role in their host's genome, they are more often deleterious, and purifying selection is an important factor controlling their genomic loads. In contrast, life history, mating system, GC content, and RNAi pathways have been suggested to account for the disparity of TE loads in different species. Previous studies of fungal, plant, and animal genomes have reported conflicting results regarding the direction in which these genomic features drive TE evolution. Many of these studies have had limited power, however, because they studied taxonomically narrow systems, comparing only a limited number of phylogenetically independent contrasts, and did not address long-term effects on TE evolution. Here, we test the long-term determinants of TE evolution by comparing 42 nematode genomes spanning over 500 million years of diversification. This analysis includes numerous transitions between life history states, and RNAi pathways, and evaluates if these forces are sufficiently persistent to affect the long-term evolution of TE loads in eukaryotic genomes. Although we demonstrate statistical power to detect selection, we find no evidence that variation in these factors influence genomic TE loads across extended periods of time. In contrast, the effects of genetic drift appear to persist and control TE variation among species. We suggest that variation in the tested factors are largely inconsequential to the large differences in TE content observed between genomes, and only by these large-scale comparisons can we distinguish long-term and persistent effects from transient or random changes.

Insight into different environmental niches adaptation and allergenicity from the Cladosporium sphaerospermum genome, a common human allergy-eliciting Dothideomycetes

Cladosporium sphaerospermum, a dematiaceous saprophytic fungus commonly found in diverse environments, has been reported to cause allergy and other occasional diseases in humans. However, its basic biology and genetic information are largely unexplored. A clinical isolate C. sphaerospermum genome, UM 843, was re-sequenced and combined with previously generated sequences to form a model 26.89 Mb genome containing 9,652 predicted genes. Functional annotation on predicted genes suggests the ability of this fungus to degrade carbohydrate and protein complexes. Several putative peptidases responsible for lung tissue hydrolysis were identified. These genes shared high similarity with the Aspergillus peptidases. The UM 843 genome encodes a wide array of proteins involved in the biosynthesis of melanin, siderophores, cladosins and survival in high salinity environment. In addition, a total of 28 genes were predicted to be associated with allergy. Orthologous gene analysis together with 22 other Dothideomycetes showed genes uniquely present in UM 843 that encode four class 1 hydrophobins which may be allergens specific to Cladosporium. The mRNA of these hydrophobins were detected by RT-PCR. The genomic analysis of UM 843 contributes to the understanding of the biology and allergenicity of this widely-prevalent species.

Tyrosine Recombinase Retrotransposons and Transposons

Retrotransposons carrying tyrosine recombinases (YR) are widespread in eukaryotes. The first described tyrosine recombinase mobile element, DIRS1, is a retroelement from the slime mold Dictyostelium discoideum. The YR elements are bordered by terminal repeats related to their replication via free circular dsDNA intermediates. Site-specific recombination is believed to integrate the circle without creating duplications of the target sites. Recently a large number of YR retrotransposons have been described, including elements from fungi (mucorales and basidiomycetes), plants (green algae) and a wide range of animals including nematodes, insects, sea urchins, fish, amphibia and reptiles. YR retrotransposons can be divided into three major groups: the DIRS elements, PAT-like and the Ngaro elements. The three groups form distinct clades on phylogenetic trees based on alignments of reverse transcriptase/ribonuclease H (RT/RH) and YR sequences, and also having some structural distinctions. A group of eukaryote DNA transposons, cryptons, also carry tyrosine recombinases. These DNA transposons do not encode a reverse transcriptase. They have been detected in several pathogenic fungi and oomycetes. Sequence comparisons suggest that the crypton YRs are related to those of the YR retrotransposons. We suggest that the YR retrotransposons arose from the combination of a crypton-like YR DNA transposon and the RT/RH encoding sequence of a retrotransposon. This acquisition must have occurred at a very early point in the evolution of eukaryotes.

A survey of transposable element classification systems--a call for a fundamental update to meet the challenge of their diversity and complexity

The increase of publicly available sequencing data has allowed for rapid progress in our understanding of genome composition. As new information becomes available we should constantly be updating and reanalyzing existing and newly acquired data. In this report we focus on transposable elements (TEs) which make up a significant portion of nearly all sequenced genomes. Our ability to accurately identify and classify these sequences is critical to understanding their impact on host genomes. At the same time, as we demonstrate in this report, problems with existing classification schemes have led to significant misunderstandings of the evolution of both TE sequences and their host genomes. In a pioneering publication Finnegan (1989) proposed classifying all TE sequences into two classes based on transposition mechanisms and structural features: the retrotransposons (class I) and the DNA transposons (class II). We have retraced how ideas regarding TE classification and annotation in both prokaryotic and eukaryotic scientific communities have changed over time. This has led us to observe that: (1) a number of TEs have convergent structural features and/or transposition mechanisms that have led to misleading conclusions regarding their classification, (2) the evolution of TEs is similar to that of viruses by having several unrelated origins, (3) there might be at least 8 classes and 12 orders of TEs including 10 novel orders. In an effort to address these classification issues we propose: (1) the outline of a universal TE classification, (2) a set of methods and classification rules that could be used by all scientific communities involved in the study of TEs, and (3) a 5-year schedule for the establishment of an International Committee for Taxonomy of Transposable Elements (ICTTE).

Comparative analysis of transposable elements highlights mobilome diversity and evolution in vertebrates

Transposable elements (TEs) are major components of vertebrate genomes, with major roles in genome architecture and evolution. In order to characterize both common patterns and lineage-specific differences in TE content and TE evolution, we have compared the mobilomes of 23 vertebrate genomes, including 10 actinopterygian fish, 11 sarcopterygians, and 2 nonbony vertebrates. We found important variations in TE content (from 6% in the pufferfish tetraodon to 55% in zebrafish), with a more important relative contribution of TEs to genome size in fish than in mammals. Some TE superfamilies were found to be widespread in vertebrates, but most elements showed a more patchy distribution, indicative of multiple events of loss or gain. Interestingly, loss of major TE families was observed during the evolution of the sarcopterygian lineage, with a particularly strong reduction in TE diversity in birds and mammals. Phylogenetic trends in TE composition and activity were detected: Teleost fish genomes are dominated by DNA transposons and contain few ancient TE copies, while mammalian genomes have been predominantly shaped by nonlong terminal repeat retrotransposons, along with the persistence of older sequences. Differences were also found within lineages: The medaka fish genome underwent more recent TE amplification than the related platyfish, as observed for LINE retrotransposons in the mouse compared with the human genome. This study allows the identification of putative cases of horizontal transfer of TEs, and to tentatively infer the composition of the ancestral vertebrate mobilome. Taken together, the results obtained highlight the importance of TEs in the structure and evolution of vertebrate genomes, and demonstrate their major impact on genome diversity both between and within lineages.

The evolution of tyrosine-recombinase elements in Nematoda

Transposable elements can be categorised into DNA and RNA elements based on their mechanism of transposition. Tyrosine recombinase elements (YREs) are relatively rare and poorly understood, despite sharing characteristics with both DNA and RNA elements. Previously, the Nematoda have been reported to have a substantially different diversity of YREs compared to other animal phyla: the Dirs1-like YRE retrotransposon was encountered in most animal phyla but not in Nematoda, and a unique Pat1-like YRE retrotransposon has only been recorded from Nematoda. We explored the diversity of YREs in Nematoda by sampling broadly across the phylum and including 34 genomes representing the three classes within Nematoda. We developed a method to isolate and classify YREs based on both feature organization and phylogenetic relationships in an open and reproducible workflow. We also ensured that our phylogenetic approach to YRE classification identified truncated and degenerate elements, informatively increasing the number of elements sampled. We identified Dirs1-like elements (thought to be absent from Nematoda) in the nematode classes Enoplia and Dorylaimia indicating that nematode model species do not adequately represent the diversity of transposable elements in the phylum. Nematode Pat1-like elements were found to be a derived form of another Pat1-like element that is present more widely in animals. Several sequence features used widely for the classification of YREs were found to be homoplasious, highlighting the need for a phylogenetically-based classification scheme. Nematode model species do not represent the diversity of transposable elements in the phylum.

ZBED evolution: repeated utilization of DNA transposons as regulators of diverse host functions

ZBED genes originate from domesticated hAT DNA transposons and encode regulatory proteins of diverse function in vertebrates. Here we reveal the evolutionary relationship between ZBED genes and demonstrate that they are derived from at least two independent domestication events in jawed vertebrate ancestors. We show that ZBEDs form two monophyletic clades, one of which has expanded through several independent duplications in host lineages. Subsequent diversification of ZBED genes has facilitated regulation of multiple diverse fundamental functions. In contrast to known examples of transposable element exaptation, our results demonstrate a novel unprecedented capacity for the repeated utilization of a family of transposable element-derived protein domains sequestered as regulators during the evolution of diverse host gene functions in vertebrates. Specifically, ZBEDs have contributed to vertebrate regulatory innovation through the donation of modular DNA and protein interacting domains. We identify that C7ORF29, ZBED2, 3, 4, and ZBEDX form a monophyletic group together with ZBED6, that is distinct from ZBED1 genes. Furthermore, we show that ZBED5 is related to Buster DNA transposons and is phylogenetically separate from other ZBEDs. Our results offer new insights into the evolution of regulatory pathways, and suggest that DNA transposons have contributed to regulatory complexity during genome evolution in vertebrates.

MITEs in the promoters of effector genes allow prediction of novel virulence genes in Fusarium oxysporum

BACKGROUND

The plant-pathogenic fungus Fusarium oxysporum f.sp.lycopersici (Fol) has accessory, lineage-specific (LS) chromosomes that can be transferred horizontally between strains. A single LS chromosome in the Fol4287 reference strain harbors all known Fol effector genes. Transfer of this pathogenicity chromosome confers virulence to a previously non-pathogenic recipient strain. We hypothesize that expression and evolution of effector genes is influenced by their genomic context.

RESULTS

To gain a better understanding of the genomic context of the effector genes, we manually curated the annotated genes on the pathogenicity chromosome and identified and classified transposable elements. Both retro- and DNA transposons are present with no particular overrepresented class. Retrotransposons appear evenly distributed over the chromosome, while DNA transposons tend to concentrate in large chromosomal subregions. In general, genes on the pathogenicity chromosome are dispersed within the repeat landscape. Effector genes are present within subregions enriched for DNA transposons. A miniature Impala (mimp) is always present in their promoters. Although promoter deletion studies of two effector gene loci did not reveal a direct function of the mimp for gene expression, we were able to use proximity to a mimp as a criterion to identify new effector gene candidates. Through xylem sap proteomics we confirmed that several of these candidates encode proteins secreted during plant infection.

CONCLUSIONS

Effector genes in Fol reside in characteristic subregions on a pathogenicity chromosome. Their genomic context allowed us to develop a method for the successful identification of novel effector genes. Since our approach is not based on effector gene similarity, but on unique genomic features, it can easily be extended to identify effector genes in Fo strains with different host specificities.

Evidence for small RNAs homologous to effector-encoding genes and transposable elements in the oomycete Phytophthora infestans

Phytophthora infestans is the oomycete pathogen responsible for the devastating late blight disease on potato and tomato. There is presently an intense research focus on the role(s) of effectors in promoting late blight disease development. However, little is known about how they are regulated, or how diversity in their expression may be generated among different isolates. Here we present data from investigation of RNA silencing processes, characterized by non-coding small RNA molecules (sRNA) of 19-40 nt. From deep sequencing of sRNAs we have identified sRNAs matching numerous RxLR and Crinkler (CRN) effector protein genes in two isolates differing in pathogenicity. Effector gene-derived sRNAs were present in both isolates, but exhibited marked differences in abundance, especially for CRN effectors. Small RNAs in P. infestans grouped into three clear size classes of 21, 25/26 and 32 nt. Small RNAs from all size classes mapped to RxLR effector genes, but notably 21 nt sRNAs were the predominant size class mapping to CRN effector genes. Some effector genes, such as PiAvr3a, to which sRNAs were found, also exhibited differences in transcript accumulation between the two isolates. The P. infestans genome is rich in transposable elements, and the majority of sRNAs of all size classes mapped to these sequences, predominantly to long terminal repeat (LTR) retrotransposons. RNA silencing of Dicer and Argonaute genes provided evidence that generation of 21 nt sRNAs is Dicer-dependent, while accumulation of longer sRNAs was impacted by silencing of Argonaute genes. Additionally, we identified six microRNA (miRNA) candidates from our sequencing data, their precursor sequences from the genome sequence, and target mRNAs. These miRNA candidates have features characteristic of both plant and metazoan miRNAs.

Dynamics and innovations within oomycete genomes: insights into biology, pathology, and evolution

The eukaryotic microbes known as oomycetes are common inhabitants of terrestrial and aquatic environments and include saprophytes and pathogens. Lifestyles of the pathogens extend from biotrophy to necrotrophy, obligate to facultative pathogenesis, and narrow to broad host ranges on plants or animals. Sequencing of several pathogens has revealed striking variation in genome size and content, a plastic set of genes related to pathogenesis, and adaptations associated with obligate biotrophy. Features of genome evolution include repeat-driven expansions, deletions, gene fusions, and horizontal gene transfer in a landscape organized into gene-dense and gene-sparse sectors and influenced by transposable elements. Gene expression profiles are also highly dynamic throughout oomycete life cycles, with transcriptional polymorphisms as well as differences in protein sequence contributing to variation. The genome projects have set the foundation for functional studies and should spur the sequencing of additional species, including more diverse pathogens and nonpathogens.