REMiner-II: A tool for rapid identification and configuration of repetitive element arrays from large mammalian chromosomes as a single query

Genes occupy ~3% of the human and mouse genomes whereas repetitive elements (REs), whose biologic functions are largely uncharacterized, constitute greater than 50%. A heterogeneous population of RE arrays (arrangement structures) is formed by combinations of various REs in mammalian genomes. In this study, REMiner-II was refined from the original REMiner for a more efficient identification and configuration of RE arrays from large queries (e.g., human chromosomes) using an unbiased self-alignment protocol. Chromosome-wide RE array profiles for the entire sets of human and mouse chromosomes were obtained using REMiner-II on a personal computer. REMiner-II provides 10 adjustable parameters and three data output modes to accommodate different experimental settings and/or goals. Examination of the human and mouse chromosome data using the REMiner-II viewer revealed species-specific libraries of complexly organized RE arrays. In conclusion, REMiner-II is an efficient tool for chromosome-wide identification and characterization of RE arrays from mammalian genomes.     

REMiner: a tool for unbiased mining and analysis of repetitive elements and their arrangement structures of large chromosomes

Repetitive elements (REs) constitute a substantial portion of the genomes of human and other species; however, the RE profiles (type, density, and arrangement) within the individual genomes have not been fully characterized. In this study, we developed an RE analysis tool, called REMiner, for a chromosome-wide investigation into the occurrence of individual REs and arrangement of clusters of REs, and REMiner's functional features were examined using the human chromosome Y. The algorithm implemented by REMiner focused on unbiased mining of REs in large chromosomes and data interface within a viewer. The data from the chromosome demonstrated that REMiner is an efficient tool in regard to its capacity for a large query size and the availability of a high-resolution viewer, featuring instant retrieval of alignment data and control of magnification and identity ratio. The chromosome-wide survey identified a diverse population of ordered RE arrangements, which may participate in the genome biology.     

Assays of testis development in the mouse distinguish three classes of domesticus-type Y chromosome

The Y chromosome of Mus musculus poschiavinus interacts with the autosomal recessive gene tda-1b of the C57BL/6J laboratory strain of the house mouse to cause complete or partial sex reversal. Ovaries or ovotestes develop in a substantial proportion of the XY fetuses. Several different Y-specific DNA probes distinguish two major types of Y chromosome in the house mouse and they are represented by M. m. domesticus and M. m. musculus. The poschiavinus Y chromosome appears identical to the domesticus Y. The developmental distribution of the gonad types was examined in the first backcross or N2 generation of fetuses in C57BL/6J with six different domesticus-type Y chromosomes and, as controls, three different musculus-type Y chromosomes. Gonadal hermaphrodites were found with three of the six domesticus-type Y chromosomes. Both overall frequency and phenotypic distribution of types of gonadal hermaphrodites identify three classes of domesticus-type Y chromosome by their differential interaction with the C57BL/6J genetic background.     

A DNA polymorphism from five inbred strains of the mouse identifies a functional class of domesticus-type Y chromosome that produces the same phenotypic distribution of gonadal hermaphrodites

The wild-derived CLA inbred strain of the house mouse contains a domesticus-type Y chromosome that lacks a 2.3-kb TaqI band with fragment 1 of the AC11 probe. The CLA Y chromosome also causes a low frequency of XY gonadal hermaphrodites when backcrossed to the C57BL/6J strain (F. G. Biddle and Y. Nishioka, 1988. Genome, 30:870-878). A similar domesticus-type Y chromosome, lacking the 2.3-kb TaqI band has now been found in the four historical inbred strains AKR/J, MA/MyJ, PL/J, and RF/J. When backcrossed to C57BL/6J, these four Y chromosomes cause low frequencies of gonadal hermaphrodites similar to the CLA Y and phenotypic distribution of types of gonad are indistinguishable from that with the CLA Y. The absence of the 2.3-kb TaqI band appears to be a polymorphism among domesticus-type Y chromosomes that identifies one of the three functional classes that, so far, can be distinguished only by their effects on testis differentiation in backcross test fetuses with the C57BL/6J strain. Three other historical inbred strains, BUB/BnJ, ST/bJ, and SWR/J, with a domesticus-type Y chromosome but containing the 2.3-kb TaqI band, were also assayed. They permit normal testis development in backcross test fetuses with C57BL/6J.     

Sequence instability and functional inactivation of murine Y chromosomes can occur on a specific genetic background.

When the Y chromosome from Mus. poschiavinus (YPos) is backcrossed onto the C57BL/6J laboratory strain, testicular dysfunction occurs at high frequencies. When five different multicopy probes from the recombinationally suppressed region of the Y chromosome were used, genomic DNAs from sibling female progeny of C57BL/6J YPos males were found to contain YPos-specific sequences ranging from trace levels to levels consistent with an intact Y chromosome. Females with a high copy number of YPos-specific sequences had a karyotype of XYPos and were sterile. Females with trace levels of these sequences were XX and fertile. Repeated sequences in the testis-determining-region (Sxr) of inactive YPos chromosomes were unstable relative to sequences in non-Sxr regions. In contrast, the YPos chromosome was stable and functioned normally in other inbred laboratory strains such as 129/Sv. The frequency and extent of YPos chromosome instability increased with successive backcrosses from stable (129/Sv) to unstable (C57BL/6J) genetic backgrounds. Traces of YPos-specific sequences were first detected in N2 female offspring of F1 males. Therefore, sequences were deleted from YPos chromosomes in the F1 male germ line and were transmitted to N2 females; inactive YPos chromosomes (XYPos females) were first detected in the N3 generation. The mouse line being derived by backcrossing the YPos chromosome onto C57BL/6J inbred strains ended in the N7 generation, since all XYPos offspring were sterile. Even stable repeated sequences from the non-Sxr regions of their inactive YPos chromosomes were precisely rearranged in these N7 offspring at high frequencies. These data are consistent with hybrid dysgenesis in mammals.     

A pronounced evolutionary shift of the pseudoautosomal region boundary in house mice

The pseudoautosomal region (PAR) is essential for the accurate pairing and segregation of the X and Y chromosomes during meiosis. Despite its functional significance, the PAR shows substantial evolutionary divergence in structure and sequence between mammalian species. An instructive example of PAR evolution is the house mouse Mus musculus domesticus (represented by the C57BL/6J strain), which has the smallest PAR among those that have been mapped. In C57BL/6J, the PAR boundary is located just ~700 kb from the distal end of the X chromosome, whereas the boundary is found at a more proximal position in Mus spretus, a species that diverged from house mice 2-4 million years ago. In this study we used a combination of genetic and physical mapping to document a pronounced shift in the PAR boundary in a second house mouse subspecies, Mus musculus castaneus (represented by the CAST/EiJ strain), ~430 kb proximal of the M. m. domesticus boundary. We demonstrate molecular evolutionary consequences of this shift, including a marked lineage-specific increase in sequence divergence within Mid1, a gene that resides entirely within the M. m. castaneus PAR but straddles the boundary in other subspecies. Our results extend observations of structural divergence in the PAR to closely related subspecies, pointing to major evolutionary changes in this functionally important genomic region over a short time period.     

Genomic survey of prepulse inhibition in mouse chromosome substitution strains

Prepulse inhibition (PPI) is a measure of sensorimotor gating, a pre-attentional inhibitory brain mechanism that filters extraneous stimuli. Prepulse inhibition is correlated with measures of cognition and executive functioning, and is considered an endophenotype of schizophrenia and other psychiatric illnesses in which patients show PPI impairments. As a first step toward identifying genes that regulate PPI, we performed a quantitative trait locus (QTL) screen of PPI phenotypes in a panel of mouse chromosome substitution strains (CSSs). We identified five CSSs with altered PPI compared with the host C57BL/6J strain: CSS-4 exhibited decreased PPI, whereas CSS-10, -11, -16 and -Y exhibited higher PPI compared with C57BL/6J. These data indicate that A/J chromosomes 4, 10, 11, 16 and Y harbor at least one QTL region that modulates PPI in these CSSs. Quantitative trait loci for the acoustic startle response were identified on seven chromosomes. Like PPI, habituation of the startle response is also disrupted in schizophrenia, and in the present study CSS-7 and -8 exhibited deficits in startle habituation. Linkage analysis of an F₂ intercross identified a highly significant QTL for PPI on chromosome 11 between positions 101.5 and 114.4 Mb (peak LOD = 4.54). Future studies will map the specific genes contributing to these QTLs using congenic strains and other genomic approaches. Identification of genes that modulate PPI will provide insight into the neural mechanisms underlying sensorimotor gating, as well as the psychopathology of disorders characterized by gating deficits.     

A comparative phenotypic and genomic analysis of C57BL/6J and C57BL/6N mouse strains

Background

The mouse inbred line C57BL/6J is widely used in mouse genetics and its genome has been incorporated into many genetic reference populations. More recently large initiatives such as the International Knockout Mouse Consortium (IKMC) are using the C57BL/6N mouse strain to generate null alleles for all mouse genes. Hence both strains are now widely used in mouse genetics studies. Here we perform a comprehensive genomic and phenotypic analysis of the two strains to identify differences that may influence their underlying genetic mechanisms.

Results

We undertake genome sequence comparisons of C57BL/6J and C57BL/6N to identify SNPs, indels and structural variants, with a focus on identifying all coding variants. We annotate 34 SNPs and 2 indels that distinguish C57BL/6J and C57BL/6N coding sequences, as well as 15 structural variants that overlap a gene. In parallel we assess the comparative phenotypes of the two inbred lines utilizing the EMPReSSslim phenotyping pipeline, a broad based assessment encompassing diverse biological systems. We perform additional secondary phenotyping assessments to explore other phenotype domains and to elaborate phenotype differences identified in the primary assessment. We uncover significant phenotypic differences between the two lines, replicated across multiple centers, in a number of physiological, biochemical and behavioral systems.

Conclusions

Comparison of C57BL/6J and C57BL/6N demonstrates a range of phenotypic differences that have the potential to impact upon penetrance and expressivity of mutational effects in these strains. Moreover, the sequence variants we identify provide a set of candidate genes for the phenotypic differences observed between the two strains.     

A new polymorphic retroelement database (PRED) for the human genome

Comparison of primate genomic sequences has demonstrated that the intra-and interspecific genetic variation is provided by retroelements (REs). The human genome contains many thousands of polymorphic RE copies, which are regarded as a promising source of new generation molecular genetic markers. However, the absence of systematized data on the RE number, distribution, genomic context, and abundance in various human populations limits the use of RE insertion polymorphism. We designed the first bilingual (Russian/English) web resource on the known polymorphic REs discovered both by our team and other researchers. The database contains the information about the genomic location of each RE, its position relative to known and predicted genes, abundance in human populations, and other data. Our web portal () allows a search of the database with user-specified parameters. The database makes it possible to most comprehensively analyze the RE distribution in the human genome and to design molecular genetic markers for studies of human genome diversity and biomedical applications.     

Mus Spretus Line-1 Sequences Detected in the Mus Musculus Inbred Strain C57bl/6j Using Line-1 DNA Probes

The inbred mouse strain, C57BL/6J, was derived from mice of the Mus musculus complex. C57BL/6J can be crossed in the laboratory with a closely related mouse species, M. spretus to produce fertile offspring; however there has been no previous evidence of gene flow between M. spretus and M. musculus in nature. Analysis of the repetitive sequence LINE-1, using both direct sequence analysis and genomic Southern blot hybridization to species-specific LINE-1 hybridization probes, demonstrates the presence of LINE-1 elements in C57BL/6J that were derived from the species M. spretus. These spretus-like LINE-1 elements in C57BL/6J reveal a cross to M. spretus somewhere in the history of C57BL/6J. It is unclear if the spretus-like LINE-1 elements are still embedded in flanking DNA derived from M. spretus or if they have transposed to new sites. The number of spretus-like elements detected suggests a maximum of 6.5% of the C57BL/6J genome may be derived from M. spretus.     

Evidence of maternal QTL affecting growth and obesity in adult mice

Most quantitative trait loci (QTL) studies fail to account for the effect that the maternal genotype may have on an individual’s phenotypes, even though maternal effect QTL have been shown to account for considerable variation in growth and obesity traits in mouse models. Moreover, the fetal programming theory suggests that maternal effects influence an offspring’s adult fitness, although the genetic nature of fetal programming remains unclear. Within this context, our study focused on mapping genomic regions associated with maternal effect QTL by analyzing the phenotypes of chromosomes 2 and 7 subcongenic mice from genetically distinct dams. We analyzed 12 chromosome 2 subcongenic strains that spanned from 70 to 180 Mb with CAST/EiJ donor regions on the background of C57BL/6 J, and 14 chromosome 7 subcongenic strains that spanned from 81 to 111 Mb with BALB/cByJ donor regions on C57BL/6ByJ background. Maternal QTL analyses were performed on the basis of overlapping donor regions between subcongenic strains. We identified several highly significant (P < 5 × 10⁻⁴) maternal QTL influencing total body weight, organ weight, and fat pad weights in both sets of subcongenics. These QTL accounted for 1.9-11.7% of the phenotypic variance for growth and obesity and greatly narrowed the genomic regions associated with the maternal genetic effects. These maternal effect QTL controlled phenotypic traits in adult mice, suggesting that maternal influences at early stages of development may permanently affect offspring performance. Identification of maternal effects in our survey of two sets of subcongenic strains, representing approximately 5% of the mouse genome, supports the hypothesis that maternal effects represent significant sources of genetic variation that are largely ignored in genetic studies.     

The Mouse Y Chromosome: Enrichment, Sizing, and Cloning by Bivariate Flow Cytometry

In this report, we demonstrate the utility of interleukin-2 (IL-2) stimulation of spleen cell cultures and bivariate flow cytometry in the analysis and purification of the C57BL/6J mouse Y Chromosome. We determined that the DNA content of the C57BL/6J Y Chromosome is approximately 94.7 Mb, making it similar in size to human Chromosome 16 and significantly larger than previous estimates. In addition, we describe the bulk isolation of mouse Y Chromosomes and demonstrate enrichment of the isolated material using a fluorescencein situhybridization strategy. We detail the construction of two small insert Y Chromosome-specific libraries, ideal for sampling Y Chromosome sequences. From these libraries 1566 clones were analyzed. We provide a detailed characterization of 103 clones, generating nearly 50 kb of sequence. For 30 of these clones, we identify regions of homology to known Y chromosomal sequences, confirming the enrichment of the sorted DNA. From the remaining characterized clones, we describe the development of 15 male-specific PCR assays and 19 male–female PCR assays potentially originating from the pseudoautosomal region or other areas of X–Y or autosome–Y homology.     

Genome structure analysis of molluscs revealed whole genome duplication and lineage specific repeat variation

Comparative genome structure analysis allows us to identify novel genes, repetitive sequences and gene duplications. To explore lineage-specific genomic changes of the molluscs that is good model for development of nervous system in invertebrate, we conducted comparative genome structure analyses of three molluscs, pygmy squid, nautilus and scallops using partial genome shotgun sequencing. Most effective elements on the genome structural changes are repetitive elements (REs) causing expansion of genome size and whole genome duplication producing large amount of novel functional genes. Therefore, we investigated variation and proportion of REs and whole genome duplication. We, first, identified variations of REs in the three molluscan genomes by homology-based and de novo RE detection. Proportion of REs were 9.2%, 4.0%, and 3.8% in the pygmy squid, nautilus and scallop, respectively. We, then, estimated genome size of the species as 2.1, 4.2 and 1.8 Gb, respectively, with 2× coverage frequency and DNA sequencing theory. We also performed a gene duplication assay based on coding genes, and found that large-scale duplication events occurred after divergence from the limpet Lottia, an out-group of the three molluscan species. Comparison of all the results suggested that RE expansion did not relate to the increase in genome size of nautilus. Despite close relationships to nautilus, the squid has the largest portion of REs and smaller genome size than nautilus. We also identified lineage-specific RE and gene-family expansions, possibly relate to acquisition of the most complicated eye and brain systems in the three species.     

Eye lens proteomics: from global approach to detailed information about phakinin and gamma E and F crystallin genes

Exploration of the lenticular proteome poses a challenging and worthwhile undertaking as cataracts, the products of a disease phenotype elicited by this proteome, remains the leading cause of vision impairment worldwide. The complete ten day old lens proteome of Mus musculus C57BL/6J was resolved into 900 distinct spots by large gel carrier ampholyte based 2-DE. The predicted amino acid sequences of all 16 crystallins ubiquitous in mammals were corroborated by mass spectrometry (MS). In detailed individual spot analyses, the primary structure of the full murine C57BL/6J beaded filament component phakinin CP49 was sequenced by liquid chromatography/electrospray ionization-tandem MS and amended at two positions. This definitive polypeptide sequence was aligned to the mouse genome, thus identifying the entire C57BL/6J genomic coding region. Also, two murine C57/6J polypeptides, both previously classified as gamma F crystallin, were clearly distinguished by MS and electrophoretic mobility. Both were assigned to their respective genes, one of the polypeptides was reclassified as C57BL/6J gamma E crystallin. Building on these data and previous investigations an updated crystallin reference map was put forth and several non crystallin lenticular components were examined. These results represent the first part of a comprehensive investigation of the mouse lens proteome (http://www.mpiib-berlin.mpg.de/2D-PAGE) with emphasis on understanding genetic effects on proteins and disease development.     

Microarray-based global mapping of integration sites for the retrotransposon, intracisternal A-particle, in the mouse genome

Mammalian genomes contain numerous evolutionary harbored mobile elements, a part of which are still active and may cause genomic instability. Their movement and positional diversity occasionally result in phenotypic changes and variation by causing altered expression or disruption of neighboring host genes. Here, we describe a novel microarray-based method by which dispersed genomic locations of a type of retrotransposon in a mammalian genome can be identified. Using this method, we mapped the DNA elements for a mouse retrotransposon, intracisternal A-particle (IAP), within genomes of C3H/He and C57BL/6J inbred mouse strains; consequently we detected hundreds of probable IAP cDNA–integrated genomic regions, in which a considerable number of strain-specific putative insertions were included. In addition, by comparing genomic DNAs from radiation-induced myeloid leukemia cells and its reference normal tissue, we detected three genomic regions around which an IAP element was integrated. These results demonstrate the first successful genome-wide mapping of a retrotransposon type in a mammalian genome.     

Mapping of jog locus to the region between D6Mit104 and D6Mit336 on mouse chromosome 6.

The joggle mouse is a recessive ataxic mutant carrying an unknown mutation in a C3H/He (C3H)-derived chromosomal segment. Taking advantage of the mouse genome database, we selected 127 DNA microsatellite markers showing heterozygosity between C3H and C57BL/6J (B6) and a first round of screening for the joggle mutation was performed on B6-jog/+ partial congenic mice (N4). We identified 4 chromosomal regions in which 13 microsatellite markers show heterozygosity between C3H and B6. Then, we analyzed the genotype of these 4 chromosomal regions in mice that showed the joggle phenotype and mapped the jog locus between markers D6Mit104 (111.4 Mb) and D6Mit336 (125.1 Mb) (an interval of 13.7 Mb) on chromosome 6. By using a partial congenic strain together with the mouse genome database, we successfully mapped the chromosomal localization of the jog locus much more efficiently than by conventional linkage analysis.     

Genes regulating testis size.

The studies described here provide information about the genetic and morphological bases for the significant differences in testis size among three closely related C57BL mouse substrains: C57BL/6J, C57BL/6ByJ, and C57BL/10J. C57BL/6J mice have normal-size testes while the other two substrains have small-size testes. Genes controlling testis size are postulated to be among the estimated forty genes that differ between the C57BL/6J and C57BL/6ByJ substrains. The number of genes involved in testis size regulation was examined using recombinant inbred mouse strains. An investigation of the role of Y chromosome genes was performed by completing molecular analyses with a mouse Y chromosome-specific probe. Sertoli and germ cell counts provided insight into the morphological basis for the different testis sizes. The experimental results suggest that there are at least two autosomal testis-size genes and that they control testis size by regulating the number of Sertoli cells.     

Chromosome‐level genome assembly of a cyprinid fish Onychostoma macrolepis by integration of nanopore sequencing,Bionano and Hi‐C technology

Onychostoma macrolepis is an emerging commercial cyprinid fish species. It is a model system for studies of sexual dimorphism and genome evolution. Here, we report the chromosome‐level assembly of the O.macrolepis genome obtained from the integration of nanopore long‐read sequencing with physical maps produced using Bionano and Hi‐C technology. A total of 87.9 Gb of nanopore sequence provided approximately 100‐fold coverage of the genome. The preliminary genome assembly was 883.2 Mb in size with a contig N50 size of 11.2 Mb. The 969 corrected contigs obtained from Bionano optical mapping were assembled into 853 scaffolds and produced an assembly of 886.5 Mb with a scaffold N50 of 16.5 Mb. Finally, using the Hi‐C data, 881.3 Mb (99.4% of genome) in 526 scaffolds were anchored and oriented in 25 chromosomes ranging in size from 25.27 to 56.49 Mb. In total, 24,770 protein‐coding genes were predicted in the genome, and ~96.85% of the genes were functionally annotated. The annotated assembly contains 93.3% complete genes from the BUSCO reference set. In addition, we identified 409 Mb (46.23% of the genome) of repetitive sequence, and 11,213 non‐coding RNAs, in the genome. Evolutionary analysis revealed that O. macrolepis diverged from common carp approximately 24.25 million years ago. The chromosomes of O. macrolepis showed an unambiguous correspondence to the chromosomes of zebrafish. The high‐quality genome assembled in this work provides a valuable genomic resource for further biological and evolutionary studies of O. macrolepis.     

Genetic control of sex-chromosomal univalency in the spermatocytes of C57BL/6J and DBA/2J mice

The genetic control of sex-chromosomal univalency was examined in the primary spermatocytes of the mouse. The C57BL/6J strain expresses 3% X-Y univalency and DBA/2J expresses 37% univalency. The reciprocal F1 and the eight types of reciprocal backcross males were examined. In the C57BL/6J--DBA/2J strain pair, X--Y univalency is controlled by three genetic systems. Autosomal factors of unknown number that are dominant in DBA/2J increase the probability of univalency from 3% in C57BL/6J to 12%. The DBA/2J-Y chromosome, in place of the C57BL/6J-Y chromosome, has an additive effect to increase the probability of univalency from 12 to 37% in the DBA/2J strain. Two X-chromosome factors that differ between C57BL/6J and DBA/2J regulate the probability of univalency. The X-chromosome factors appear to be separated by sufficient distance so that, with the DBA/2J-Y chromosome and dominant DBA/2J autosomal factors, there are two recombinant classes of X--Y univalency at 20 and 60%. The genetic factors in the univalency trait may be involved in the regulation or structure of the terminal attachment sites between the X and Y chromosomes.     

Interspecific recombinant congenic strains between C57BL/6 and mice of the Mus spretus species: a powerful tool to dissect genetic control of complex traits

Complex traits are under the genetic control of multiple genes, often with weak effects and strong epistatic interactions. We developed two new collections of mouse strains to improve genetic dissection of complex traits. They are derived from several backcrosses of the Mus spretus SEG/Pas or STF/Pas strains on the C57BL/6J background. Each of the 55 interspecific recombinant congenic strains (IRCSs) carries up to eight SEG/Pas chromosomal segments with an average size of 11.7 Mb, totalizing 1.37% of the genome. The complete series covers 39.7% of the SEG/Pas genome. As a complementary resource, six partial or complete interspecific consomic strains were developed and increased genome coverage to 45.6%. To evaluate the usefulness of these strains for QTL mapping, 16 IRCSs were compared with C57BL/6J for seven hematological parameters. Strain 66H, which carries three SEG/Pas chromosomal segments, had lower red blood cell volume and higher platelet count than C57BL/6J. Each chromosomal segment was isolated in a congenic strain to evaluate individual effects. Congenic strains were combined to assess epistasis. Our data show that both traits were controlled by several genes with complex epistatic interactions. IRCSs are therefore useful to unravel QTL with small effects and gene-by-gene interactions.