Identifying new sex-linked genes through BAC sequencing in the dioecious plant Silene latifolia

Background

Silene latifolia represents one of the best-studied plant sex chromosome systems. A new approach using RNA-seq data has recently identified hundreds of new sex-linked genes in this species. However, this approach is expected to miss genes that are either not expressed or are expressed at low levels in the tissue(s) used for RNA-seq. Therefore other independent approaches are needed to discover such sex-linked genes.

Results

Here we used 10 well-characterized S. latifolia sex-linked genes and their homologs in Silene vulgaris, a species without sex chromosomes, to screen BAC libraries of both species. We isolated and sequenced 4 Mb of BAC clones of S. latifolia X and Y and S. vulgaris genomic regions, which yielded 59 new sex-linked genes (with S. vulgaris homologs for some of them). We assembled sequences that we believe represent the tip of the Xq arm. These sequences are clearly not pseudoautosomal, so we infer that the S. latifolia X has a single pseudoautosomal region (PAR) on the Xp arm. The estimated mean gene density in X BACs is 2.2 times lower than that in S. vulgaris BACs, agreeing with the genome size difference between these species. Gene density was estimated to be extremely low in the Y BAC clones. We compared our BAC-located genes with the sex-linked genes identified in previous RNA-seq studies, and found that about half of them (those with low expression in flower buds) were not identified as sex-linked in previous RNA-seq studies. We compiled a set of ~70 validated X/Y genes and X-hemizygous genes (without Y copies) from the literature, and used these genes to show that X-hemizygous genes have a higher probability of being undetected by the RNA-seq approach, compared with X/Y genes; we used this to estimate that about 30 % of our BAC-located genes must be X-hemizygous. The estimate is similar when we use BAC-located genes that have S. vulgaris homologs, which excludes genes that were gained by the X chromosome.

Conclusions

Our BAC sequencing identified 59 new sex-linked genes, and our analysis of these BAC-located genes, in combination with RNA-seq data suggests that gene losses from the S. latifolia Y chromosome could be as high as 30 %, higher than previous estimates of 10-20 %.

Electronic supplementary material

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

Alpha2 Macroglobulin-Like Is Essential for Liver Development in Zebrafish

Background

Alpha 2 Macroglobulin family members have been studied extensively with respect to their roles in physiology and human disease including innate immunity and Alzheimer''s disease, but little is known about a possible role in liver development loss-of-function in model systems.

Principal Findings

We report the isolation of the zebrafish α2 macroglobulin-like (A2ML) gene and its specific expression in the liver during differentiation. Morpholino-based knock-down of A2ML did not block the initial formation of the liver primordium, but inhibited liver growth and differentiation.

Significance

This report on A2ML function in zebrafish development provides the first evidence for a specific role of an A2M family gene in liver formation during early embryogenesis in a vertebrate.     

Inflammatory expression profiles in monocyte-to-macrophage differentiation in patients with systemic lupus erythematosus and relationship with atherosclerosis

Introduction

Our objectives were to examine mononuclear cell gene expression profiles in patients with systemic lupus erythematosus (SLE) and healthy controls and to compare subsets with and without atherosclerosis to determine which genes’ expression is related to atherosclerosis in SLE.

Methods

Monocytes were obtained from 20 patients with SLE and 16 healthy controls and were in vitro-differentiated into macrophages. Subjects also underwent laboratory and imaging studies to evaluate for subclinical atherosclerosis. Whole-genome RNA expression microarray was performed, and gene expression was examined.

Results

Gene expression profiling was used to identify gene signatures that differentiated patients from controls and individuals with and without atherosclerosis. In monocytes, 9 out of 20 patients with SLE had an interferon-inducible signature compared with 2 out of 16 controls. By looking at gene expression during monocyte-to-macrophage differentiation, we identified pathways which were differentially regulated between SLE and controls and identified signatures based on relevant intracellular signaling molecules which could differentiate SLE patients with atherosclerosis from controls. Among patients with SLE, we used a previously defined 344-gene atherosclerosis signature in monocyte-to-macrophage differentiation to identify patient subgroups with and without atherosclerosis. Interestingly, this signature further classified patients on the basis of the presence of SLE disease activity and cardiovascular risk factors.

Conclusions

Many genes were differentially regulated during monocyte-to-macrophage differentiation in SLE patients compared with controls. The expression of these genes in mononuclear cells is important in the pathogenesis of SLE, and molecular profiling using gene expression can help stratify SLE patients who may be at risk for development of atherosclerosis.     

The physical map of wheat chromosome 5DS revealed gene duplications and small rearrangements

Background

The substantially large bread wheat genome, organized into highly similar three sub-genomes, renders genomic research challenging. The construction of BAC-based physical maps of individual chromosomes reduces the complexity of this allohexaploid genome, enables elucidation of gene space and evolutionary relationships, provides tools for map-based cloning, and serves as a framework for reference sequencing efforts. In this study, we constructed the first comprehensive physical map of wheat chromosome arm 5DS, thereby exploring its gene space organization and evolution.

Results

The physical map of 5DS was comprised of 164 contigs, of which 45 were organized into 21 supercontigs, covering 176 Mb with an N50 value of 2,173 kb. Fifty-eight of the contigs were larger than 1 Mb, with the largest contig spanning 6,649 kb. A total of 1,864 molecular markers were assigned to the map at a density of 10.5 markers/Mb, anchoring 100 of the 120 contigs (>5 clones) that constitute ~95 % of the cumulative length of the map. Ordering of 80 contigs along the deletion bins of chromosome arm 5DS revealed small-scale breaks in syntenic blocks. Analysis of the gene space of 5DS suggested an increasing gradient of genes organized in islands towards the telomere, with the highest gene density of 5.17 genes/Mb in the 0.67-0.78 deletion bin, 1.4 to 1.6 times that of all other bins.

Conclusions

Here, we provide a chromosome-specific view into the organization and evolution of the D genome of bread wheat, in comparison to one of its ancestors, revealing recent genome rearrangements. The high-quality physical map constructed in this study paves the way for the assembly of a reference sequence, from which breeding efforts will greatly benefit.

Electronic supplementary material

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

Clustered gene expression changes flank targeted gene loci in knockout mice

Background

Gene expression profiling using microarrays is a powerful technology widely used to study regulatory networks. Profiling of mRNA levels in mutant organisms has the potential to identify genes regulated by the mutated protein.

Methodology/Principle Findings

Using tissues from multiple lines of knockout mice we have examined genome-wide changes in gene expression. We report that a significant proportion of changed genes were found near the targeted gene.

Conclusions/Significance

The apparent clustering of these genes was explained by the presence of flanking DNA from the parental ES cell. We provide recommendations for the analysis and reporting of microarray data from knockout mice     

The Chironomus tentans genome sequence and the organization of the Balbiani ring genes

Background

The polytene nuclei of the dipteran Chironomus tentans (Ch. tentans) with their Balbiani ring (BR) genes constitute an exceptional model system for studies of the expression of endogenous eukaryotic genes. Here, we report the first draft genome of Ch. tentans and characterize its gene expression machineries and genomic architecture of the BR genes.

Results

The genome of Ch. tentans is approximately 200 Mb in size, and has a low GC content (31%) and a low repeat fraction (15%) compared to other Dipteran species. Phylogenetic inference revealed that Ch. tentans is a sister clade to mosquitoes, with a split 150–250 million years ago. To characterize the Ch. tentans gene expression machineries, we identified potential orthologus sequences to more than 600 Drosophila melanogaster (D. melanogaster) proteins involved in the expression of protein-coding genes. We report novel data on the organization of the BR gene loci, including a novel putative BR gene, and we present a model for the organization of chromatin bundles in the BR2 puff based on genic and intergenic in situ hybridizations.

Conclusions

We show that the molecular machineries operating in gene expression are largely conserved between Ch. tentans and D. melanogaster, and we provide enhanced insight into the organization and expression of the BR genes. Our data strengthen the generality of the BR genes as a unique model system and provide essential background for in-depth studies of the biogenesis of messenger ribonucleoprotein complexes.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-819) contains supplementary material, which is available to authorized users.     

Female Drosophila melanogaster gene expression and mate choice: the X chromosome harbours candidate genes underlying sexual isolation

Background

The evolution of female choice mechanisms favouring males of their own kind is considered a crucial step during the early stages of speciation. However, although the genomics of mate choice may influence both the likelihood and speed of speciation, the identity and location of genes underlying assortative mating remain largely unknown.

Methods and Findings

We used mate choice experiments and gene expression analysis of female Drosophila melanogaster to examine three key components influencing speciation. We show that the 1,498 genes in Zimbabwean female D. melanogaster whose expression levels differ when mating with more (Zimbabwean) versus less (Cosmopolitan strain) preferred males include many with high expression in the central nervous system and ovaries, are disproportionately X-linked and form a number of clusters with low recombination distance. Significant involvement of the brain and ovaries is consistent with the action of a combination of pre- and postcopulatory female choice mechanisms, while sex linkage and clustering of genes lead to high potential evolutionary rate and sheltering against the homogenizing effects of gene exchange between populations.

Conclusion

Taken together our results imply favourable genomic conditions for the evolution of reproductive isolation through mate choice in Zimbabwean D. melanogaster and suggest that mate choice may, in general, act as an even more important engine of speciation than previously realized.     

Towards a comprehensive structural variation map of an individual human genome

Background

Several genomes have now been sequenced, with millions of genetic variants annotated. While significant progress has been made in mapping single nucleotide polymorphisms (SNPs) and small (<10 bp) insertion/deletions (indels), the annotation of larger structural variants has been less comprehensive. It is still unclear to what extent a typical genome differs from the reference assembly, and the analysis of the genomes sequenced to date have shown varying results for copy number variation (CNV) and inversions.

Results

We have combined computational re-analysis of existing whole genome sequence data with novel microarray-based analysis, and detect 12,178 structural variants covering 40.6 Mb that were not reported in the initial sequencing of the first published personal genome. We estimate a total non-SNP variation content of 48.8 Mb in a single genome. Our results indicate that this genome differs from the consensus reference sequence by approximately 1.2% when considering indels/CNVs, 0.1% by SNPs and approximately 0.3% by inversions. The structural variants impact 4,867 genes, and >24% of structural variants would not be imputed by SNP-association.

Conclusions

Our results indicate that a large number of structural variants have been unreported in the individual genomes published to date. This significant extent and complexity of structural variants, as well as the growing recognition of their medical relevance, necessitate they be actively studied in health-related analyses of personal genomes. The new catalogue of structural variants generated for this genome provides a crucial resource for future comparison studies.