Selective single molecule sequencing and assembly of a human Y chromosome of African origin

Advancing evolutionary medicine with complete primate genomes and advanced biotechnologies

Evolutionary medicine, which integrates evolutionary biology and medicine, significantly enhances our understanding of human traits and disease susceptibility. However, previous studies in this field have often focused on single-nucleotide variants due to technological limitations in characterizing complex genomic regions, hindering the comprehensive analyses of their evolutionary origins and clinical significance. In this review, we summarize recent advancements in complete telomere-to-telomere (T2T), primate genomes and other primate resources, and illustrate how these resources facilitate the research of complex regions. We focus on several biomedically relevant regions to examine the relationship between primate genome evolution and human diseases. We also highlight the potentials of high-throughput functional genomic technologies for assessing candidate loci. Finally, we discuss future directions for primate research within the context of evolutionary medicine.

Nanopore sequencing: flourishing in its teenage years

Over the past decade, nanopore sequencing has experienced significant advancements and changes, transitioning from an initially emerging technology to a significant instrument in the field of genomic sequencing. However, as advancements in next-generation sequencing technology persist, nanopore sequencing also improves. This paper reviews the developments, applications, and outlook on nanopore sequencing technology. Currently, nanopore sequencing supports both DNA and RNA sequencing, making it widely applicable in areas such as telomere-to-telomere (T2T) genome assembly, direct RNA sequencing (DRS), and metagenomics. The openness and versatility of nanopore sequencing have established it as a preferred option for an increasing number of research teams, signaling a transformative influence on life science research. As the nanopore sequencing technology advances, it provides a faster, more cost-effective approach with extended read lengths, demonstrating the significant potential for complex genome assembly, pathogen detection, environmental monitoring, and human disease research, offering a fresh perspective in sequencing technologies.

Emergence of infectious diseases and role of advanced nanomaterials in point-of-care diagnostics: a review

Infectious outbreaks are the foremost global public health concern, challenging the current healthcare system, which claims millions of lives annually. The most crucial way to control an infectious outbreak is by early detection through point-of-care (POC) diagnostics. POC diagnostics are highly advantageous owing to the prompt diagnosis, which is economical, simple and highly efficient with remote access capabilities. In particular, utilization of nanomaterials to architect POC devices has enabled highly integrated and portable (compact) devices with enhanced efficiency. As such, this review will detail the factors influencing the emergence of infectious diseases and methods for fast and accurate detection, thus elucidating the underlying factors of these infections. Furthermore, it comprehensively highlights the importance of different nanomaterials in POCs to detect nucleic acid, whole pathogens, proteins and antibody detection systems. Finally, we summarize findings reported on nanomaterials based on advanced POCs such as lab-on-chip, lab-on-disc-devices, point-of-action and hospital-on-chip. To this end, we discuss the challenges, potential solutions, prospects of integrating internet-of-things, artificial intelligence, 5G communications and data clouding to achieve intelligent POCs.

An effective strategy for assembling the sex-limited chromosome

BACKGROUND

Most currently available reference genomes lack the sequence map of sex-limited (such as Y and W) chromosomes, which results in incomplete assemblies that hinder further research on sex chromosomes. Recent advancements in long-read sequencing and population sequencing have provided the opportunity to assemble sex-limited chromosomes without the traditional complicated experimental efforts.

FINDINGS

We introduce the first computational method, Sorting long Reads of Y or other sex-limited chromosome (SRY), which achieves improved assembly results compared to flow sorting. Specifically, SRY outperforms in the heterochromatic region and demonstrates comparable performance in other regions. Furthermore, SRY enhances the capabilities of the hybrid assembly software, resulting in improved continuity and accuracy.

CONCLUSIONS

Our method enables true complete genome assembly and facilitates downstream research of sex-limited chromosomes.

Recent Advances in Electrochemical Biosensors Targeting Stress Markers

INTRODUCTION

When the body experiences a change in its internal environment due to factors such as mood (euphoria, stress) and illness, it releases biomarkers in large quantities. These biomarkers are used for detecting a disease at its early stages. This involves the detection of insufficient quantities of biocomponents, which can be done by using nanomaterials, conventional materials, and biotechnology; thus, scientists can increase the sensitivity of electrochemical sensors. According to studies conducted in this area, electrochemical sensors have shown promise as a diagnostic tool due to their ability to identify and pinpoint illness biomarkers. The present review article was compiled to gather the latest information on electrochemical biosensors targeting stress markers.

MATERIALS AND METHODS

The authors searched scholarly databases like ScienceDirect, Pubmed, Medline, and Scopus for information on electrochemical biosensors targeting stress markers.

RESULTS

In this article, we looked at the recent developments in electrochemical sensors for stress monitoring. Because of advances in nanomaterial and biomolecule processes, electrochemical biosensors have been developed with the sensitivity to detect several biomarkers in real-time in therapeutically relevant materials.

CONCLUSION

This biomarker sensor strategy can analyze various biofluids (sweat, plasma, urine, and saliva).

Tunnel junction sensing of TATP explosive at the single-molecule level

Triacetone triperoxide (TATP) is a highly potent homemade explosive commonly used in terrorist attacks. Its detection poses a significant challenge due to its volatility, and the lack of portability of current sensing techniques. To address this issue, we propose a novel approach based on single-molecule TATP detection in the air using a device where tunneling current in N-terminated carbon-nanotubes nanogaps is measured. By employing the density functional theory combined with the non-equilibrium Green's function method, we show that current of tens of nanoamperes passes through TATP trapped in the nanogap, with a discrimination ratio of several orders of magnitude even against prevalent indoor volatile organic compounds (VOCs). This high tunneling current through TATP's highest occupied molecular orbital (HOMO) is facilitated by the strong electric field generated by N-C polar bonds at the electrode ends and by the hybridization between TATP and the electrodes, driven by oxygen atoms within the probed molecule. The application of the same principle is discussed for graphene nanogaps and break-junctions.

Benchmarking of human Y-chromosomal haplogroup classifiers with whole-genome and whole-exome sequence data

In anthropological, medical, and forensic studies, the nonrecombinant region of the human Y chromosome (NRY) enables accurate reconstruction of pedigree relationships and retrieval of ancestral information. Using high-throughput sequencing (HTS) data, we present a benchmarking analysis of command-line tools for NRY haplogroup classification. The evaluation was performed using paired Illumina data from whole-genome sequencing (WGS) and whole-exome sequencing (WES) experiments from 50 unrelated donors. Additionally, as a validation, we also used paired WGS/WES datasets of 54 individuals from the 1000 Genomes Project. Finally, we evaluated the tools on data from third-generation HTS obtained from a subset of donors and one reference sample. Our results show that WES, despite typically offering less genealogical resolution than WGS, is an effective method for determining the NRY haplogroup. Y-LineageTracker and Yleaf showed the highest accuracy for WGS data, classifying precisely 98% and 96% of the samples, respectively. Yleaf outperforms all benchmarked tools in the WES data, classifying approximately 90% of the samples. Yleaf, Y-LineageTracker, and pathPhynder can correctly classify most samples (88%) sequenced with third-generation HTS. As a result, Yleaf provides the best performance for applications that use WGS and WES. Overall, our study offers researchers with a guide that allows them to select the most appropriate tool to analyze the NRY region using both second- and third-generation HTS data.

Assembly of 43 human Y chromosomes reveals extensive complexity and variation

The prevalence of highly repetitive sequences within the human Y chromosome has prevented its complete assembly to date1 and led to its systematic omission from genomic analyses. Here we present de novo assemblies of 43 Y chromosomes spanning 182,900 years of human evolution and report considerable diversity in size and structure. Half of the male-specific euchromatic region is subject to large inversions with a greater than twofold higher recurrence rate compared with all other chromosomes2. Ampliconic sequences associated with these inversions show differing mutation rates that are sequence context dependent, and some ampliconic genes exhibit evidence for concerted evolution with the acquisition and purging of lineage-specific pseudogenes. The largest heterochromatic region in the human genome, Yq12, is composed of alternating repeat arrays that show extensive variation in the number, size and distribution, but retain a 1:1 copy-number ratio. Finally, our data suggest that the boundary between the recombining pseudoautosomal region 1 and the non-recombining portions of the X and Y chromosomes lies 500 kb away from the currently established1 boundary. The availability of fully sequence-resolved Y chromosomes from multiple individuals provides a unique opportunity for identifying new associations of traits with specific Y-chromosomal variants and garnering insights into the evolution and function of complex regions of the human genome.

Y chromosome sequence and epigenomic reconstruction across human populations

Recent advances in long-read sequencing technologies have allowed the generation and curation of more complete genome assemblies, enabling the analysis of traditionally neglected chromosomes, such as the human Y chromosome (chrY). Native DNA was sequenced on a MinION Oxford Nanopore Technologies sequencing device to generate genome assemblies for seven major chrY human haplogroups. We analyzed and compared the chrY enrichment of sequencing data obtained using two different selective sequencing approaches: adaptive sampling and flow cytometry chromosome sorting. We show that adaptive sampling can produce data to create assemblies comparable to chromosome sorting while being a less expensive and time-consuming technique. We also assessed haplogroup-specific structural variants, which would be otherwise difficult to study using short-read sequencing data only. Finally, we took advantage of this technology to detect and profile epigenetic modifications among the considered haplogroups. Altogether, we provide a framework to study complex genomic regions with a simple, fast, and affordable methodology that could be applied to larger population genomics datasets.

Developmental validation of the ForenSeq MainstAY kit, MiSeq FGx sequencing system and ForenSeq Universal Analysis Software

For human identification purposes, forensic genetics has primarily relied upon a core set of autosomal (and to a lesser extent Y chromosome) short tandem repeat (STR) markers that are enriched by amplification using the polymerase chain reaction (PCR) that are subsequently separated and detected using capillary electrophoresis (CE). While STR typing conducted in this manner is well-developed and robust, advances in molecular biology that have occurred over the last 15 years, in particular massively parallel sequencing (MPS) [1-7], offer certain advantages as compared to CE-based typing. First and foremost is the high throughput capacity of MPS. Current bench top high throughput sequencers enable larger batteries of markers to be multiplexed and multiple samples to be sequenced simultaneously (e.g., millions to billions of nucleotides can be sequenced in one run). Second, compared to the length-based CE approach, sequencing STRs increases discrimination power, enhances sensitivity of detection, reduces noise due to instrumentation, and improves mixture interpretation [4,8-23]. Third, since detection of STRs is based on sequence and not fluorescence, amplicons can be designed that are shorter in length and of similar lengths among loci, where possible, which can improve amplification efficiency and analysis of degraded samples. Lastly, MPS offers a single format approach that can be applied to analysis of a wide variety of genetic markers of forensic interest (e.g., STRs, mitochondrial DNA, single nucleotide polymorphisms, insertion/deletions). These features make MPS a desirable technology for casework [14,15,24,25-48]. The developmental validation of the ForenSeq MainstAY library preparation kit with the MiSeq FGx Sequencing System and ForenSeq Universal Software is reported here to assist with validation of this MPS system for casework [49]. The results show that the system is sensitive, accurate and precise, specific, and performs well with mixtures and mock case-type samples.

Influence of Meteorological Temperature and Pressure on the Severity of Heart Failure Decompensations

OBJECTIVE

To investigate the relationship between ambient temperature and atmospheric pressure (AP) and the severity of heart failure (HF) decompensations.

METHODS

We analysed patients coming from the Epidemioloy Acute Heart Failure Emergency (EAHFE) Registry, a multicentre prospective cohort study enrolling patients diagnosed with decompensated HF in 26 emergency departments (EDs) of 16 Spanish cities. We recorded patient and demographic data and maximum temperature (Tmax) and AP (APmax) the day before ED consultation. Associations between temperature and AP and severity endpoints were explored by logistic regression. We used restricted cubic splines to model continuous non-linear associations of temperature and AP with each endpoint.

RESULTS

We analysed 16,545 patients. Daily Tmax and APmax (anomaly) of the day before patient ED arrival ranged from 0.8 to 41.6° and from - 61.7 to 69.9 hPa, respectively. A total of 12,352 patients (75.2%) were hospitalised, with in-hospital mortality in 1171 (7.1%). The probability of hospitalisation by HF decompensation showed a U-shaped curve versus Tmax and an increasing trend versus APmax. Regarding temperature, hospitalisation significantly increased from 20 °C (reference) upwards (25 °C: OR = 1.12, 95% CI = 1.04-1.21; 40 °C: 1.65, 1.13-2.40) and below 5.4 °C (5 °C: 1.21, 1.01-1.46). Concerning the mean AP of the city (anomaly = 0 hPa), hospitalisation increased when APmax (anomaly) was above + 7.0 hPa (atmospheric anticyclone; + 10 hPa: 1.14, 1.05-1.24; + 30 hPa: 2.02. 1.35-3.03). The lowest probability of mortality also corresponded to cold-mild temperatures and low AP, with a significant increased risk only found for Tmax above 24.3 °C (25 °C: 1.13, 1.01-1.27; 40 °C: 2.05, 1.15-3.64) and APmax (anomaly) above + 3.4 hPa (+ 10 hPa: 1.21, 1.07-1.36; + 30 hPa: 1.73, 1.06-2.81). Sensitivity analysis confirmed the main analysis results.

CONCLUSION

Temperature and AP are independently associated with the severity of HF decompensations, with possible different effects on the need for hospitalisation and in-hospital mortality.

Biophysical impacts of earth greening can substantially mitigate regional land surface temperature warming

Vegetation change can alter surface energy balance and subsequently affect the local climate. This biophysical impact has been well studied for forestation cases, but the sign and magnitude for persistent earth greening remain controversial. Based on long-term remote sensing observations, we quantify the unidirectional impact of vegetation greening on radiometric surface temperature over 2001-2018. Here, we show a global negative temperature response with large spatial and seasonal variability. Snow cover, vegetation greenness, and shortwave radiation are the major driving factors of the temperature sensitivity by regulating the relative dominance of radiative and non-radiative processes. Combined with the observed greening trend, we find a global cooling of -0.018 K/decade, which slows down 4.6 ± 3.2% of the global warming. Regionally, this cooling effect can offset 39.4 ± 13.9% and 19.0 ± 8.2% of the corresponding warming in India and China. These results highlight the necessity of considering this vegetation-related biophysical climate effect when informing local climate adaptation strategies.

Inference of population structure and admixture proportion from Y chromosomal data of Chinese population

In the past two decades, Y chromosome data has been generated for human population genetic studies. These Y chromosome datasets were produced with various testing methods and markers, thus difficult to combine them for a comprehensive analysis. In this study, we combine four human Y chromosomal datasets of Han, Tibetan, Hui, and Li ethnic groups. The dataset contains 27 microsatellites and 137 single nucleotide polymorphisms these populations share in common. We assembled a single dataset containing 2439 individuals from 25 nationwide populations in China. A systematic analysis of genetic distance and clustering was performed. To determine the gene flow of the studied population with worldwide populations, we modeled the ancestry informative markers. The reference panel was regarded as a mixture of South Asian (SAS), East Asian (EAS), European (EUR), African (AFR), and American (AMR) populations from 1000 Genomes data of Y chromosome using nonlinear data-fitting. We then calculated the admixture proportion of these four studied populations with 26 worldwide populations. The results showed that the Han and Hui have great genetic affinity, and Hui is the most admixed ethnic group, with 61.53% EAS, 34.65% SAS, 1.91% AFR, 1.56% AMR, and 0.04% EUR ancestry component (the AMR is highly admixed and thus should be ignored). All the other three ethnic groups contained more than 97% EAS ancestry component. The Li is the least admixed population in this study. The combined dataset in this study is the largest of this kind reported to date and proposes reference population data for use in future paternal genetic studies and forensic genealogical identification.

Origin and chromatin remodeling of young X/Y sex chromosomes in catfish with sexual plasticity

Assembly of a complete Y chromosome is a significant challenge in animals with an XX/XY sex-determination system. Recently, we created YY-supermale yellow catfish by crossing XY males with sex-reversed XY females, providing a valuable model for Y-chromosome assembly and evolution. Here, we assembled highly homomorphic Y and X chromosomes by sequencing genomes of the YY supermale and XX female in yellow catfish, revealing their nucleotide divergences with only less than 1% and with the same gene compositions. The sex-determining region (SDR) was identified to locate within a physical distance of 0.3 Mb by F_ST scanning. Strikingly, the incipient sex chromosomes were revealed to originate via autosome-autosome fusion and were characterized by a highly rearranged region with an SDR downstream of the fusion site. We found that the Y chromosome was at a very early stage of differentiation, as no clear evidence of evolutionary strata and classical structure features of recombination suppression for a rather late stage of Y-chromosome evolution were observed. Significantly, a number of sex-antagonistic mutations and the accumulation of repetitive elements were discovered in the SDR, which might be the main driver of the initial establishment of recombination suppression between young X and Y chromosomes. Moreover, distinct three-dimensional chromatin organizations of the Y and X chromosomes were identified in the YY supermales and XX females, as the X chromosome exhibited denser chromatin structure than the Y chromosome, while they respectively have significantly spatial interactions with female- and male-related genes compared with other autosomes. The chromatin configuration of the sex chromosomes as well as the nucleus spatial organization of the XX neomale were remodeled after sex reversal and similar to those in YY supermales, and a male-specific loop containing the SDR was found in the open chromatin region. Our results elucidate the origin of young sex chromosomes and the chromatin remodeling configuration in the catfish sexual plasticity.

A convergent mechanism of sex determination in dioecious plants: Distinct sex-determining genes display converged regulation on floral B-class genes

Sex determination in dioecious plants has been broadly and progressively studied with the blooming of genome sequencing and editing techniques. This provides us with a great opportunity to explore the evolution and genetic mechanisms underlining the sex-determining system in dioecious plants. In this study, comprehensively reviewing advances in sex-chromosomes, sex-determining genes, and floral MADS-box genes in dioecious plants, we proposed a convergent model that governs plant dioecy across divergent species using a cascade regulation pathway connecting sex-determining genes and MADS-box genes e.g., B-class genes. We believe that this convergent mechanism of sex determination in dioecious plants will shed light on our understanding of gene regulation and evolution of plant dioecy. Perspectives concerning the evolutionary pathway of plant dioecy are also suggested.

Noninvasive prenatal prediction of fetal haplotype with Spearman rank correlation analysis model

Background

Noninvasive prenatal testing (NIPT) has been widely used clinically to detect fetal chromosomal aneuploidy with high accuracy rates, gradually replacing traditional serological screening. However, the application of NIPT for monogenic diseases is still in an immature stage of exploration. The detection of mutations in peripheral blood of pregnant women requires precise qualitative and quantitative techniques, which limits its application. The bioinformatic strategies based on the SNP (single nucleotide polymorphism) linkage analysis are more practical, which can be divided into two types depending on whether proband information is needed. Hidden Markov Mode (HMM) and Sequential probability ratio test (SPRT) are suitable for families with probands. In contrast, methods based on databases and population demographic information are suitable for families without probands.

Methods

In this study, we proposed a Spearman rank correlation analysis method to infer the fetal haplotypes based on core family information. Allele frequencies of SNPs that were used to construct parental haplotypes were calculated as sets of nonparametric variables, in contrast to their theoretical values represented by a fetal fraction (FF). The effects on the calculation of the fetal concentration of two DNA enrichment methods, multiple‐PCR amplification, and targeted hybrid capture, were compared, and the heterozygosity distribution of SNPs within pedigrees was analyzed to reveal the best conditions for the model application.

Results

Predictions of the paternal haplotype inheritance were in line with expectations for both DNA library construction methods, while for maternal haplotype inheritance prediction, the rates were 96.55% for method multiple‐PCR amplification and 95.8% for method targeted hybrid capture.

Conclusion

Positive prediction rates showed that the maternal haplotype prediction was not as accurate as paternal one, due to the large amount of maternal noise in the mother's peripheral blood. Although this model is relatively immature, it provides a new perspective for noninvasive prenatal clinical tests of monogenic diseases.

Analysis of carbon emissions from land cover change during 2000 to 2020 in Shandong Province, China

Land cover change affects the carbon emissions of ecosystems in some way. The qualitative and quantitative understanding of carbon emissions from human activities (e.g., land cover change, industrial production, etc.) is highly significant for realizing the objective of carbon neutrality. Therefore, this paper used GlobeLand30 land cover maps, annual average normalised difference vegetation index (NDVI) data, annual average net ecosystem productivity (NEP) data and statistical yearbook data from 2000 to 2020 to explore the relationship between land cover change and carbon emissions. Specifically, it included land cover change, carbon storage changes influenced by land cover change, spatial and temporal analysis of carbon sources and sinks, land use intensity change and anthropogenic carbon emissions. The results of the study show that the main land cover changes in Shandong province during 2000–2020 was cultivated land conversion to artificial surfaces. Among them, the area of cultivated land converted to artificial surfaces from 2000 to 2010 was 4930.62 km², and the proportion of cultivated land converted to artificial surfaces from 2010 to 2020 was as high as 78.35%. The total carbon stock of vegetation affected by land cover change decreased by 463.96 × 10⁴ t and 193.50 × 10⁴ t in 2000–2010 and 2010–2020 respectively. The spatial and temporal distribution of carbon sources and sinks differed more markedly from 2000 to 2020, and land use intensity changes in Shandong Province showed an upward trend. Of the total energy production, industry has the largest energy consumption, followed closely by total energy consumption in transportation, storage and postal services.

Genomic architecture of 5S rDNA cluster and its variations within and between species

BACKGROUND

Ribosomal DNAs (rDNAs) are arranged in purely tandem repeats, preventing them from being reliably assembled onto chromosomes during generation of genome assembly. The uncertainty of rDNA genomic structure presents a significant barrier for studying their function and evolution.

RESULTS

Here we generate ultra-long Oxford Nanopore Technologies (ONT) and short NGS reads to delineate the architecture and variation of the 5S rDNA cluster in the different strains of C. elegans and C. briggsae. We classify the individual rDNA's repeating units into 25 types based on the unique sequence variations in each unit of C. elegans (N2). We next perform assembly of the cluster by taking advantage of the long reads that carry these units, which led to an assembly of 5S rDNA cluster consisting of up to 167 consecutive 5S rDNA units in the N2 strain. The ordering and copy number of various rDNA units are consistent with the separation time between strains. Surprisingly, we observed a drastically reduced level of variation in the unit composition in the 5S rDNA cluster in the C. elegans CB4856 and C. briggsae AF16 strains than in the C. elegans N2 strain, suggesting that N2, a widely used reference strain, is likely to be defective in maintaining the 5S rDNA cluster stability compared with other wild isolates of C. elegans or C. briggsae.

CONCLUSIONS

The results demonstrate that Nanopore DNA sequencing reads are capable of generating assembly of highly repetitive sequences, and rDNA units are highly dynamic both within and between population(s) of the same species in terms of sequence and copy number. The detailed structure and variation of the 5S rDNA units within the rDNA cluster pave the way for functional and evolutionary studies.

Decoding sex: Elucidating sex determination and how high-quality genome assemblies are untangling the evolutionary dynamics of sex chromosomes

Sexual reproduction is a diverse and widespread process. In gonochoristic species, the differentiation of sexes occurs through diverse mechanisms, influenced by environmental and genetic factors. In most vertebrates, a master-switch gene is responsible for triggering a sex determination network. However, only a few genes have acquired master-switch functions, and this process is associated with the evolution of sex-chromosomes, which have a significant influence in evolution. Additionally, their highly repetitive regions impose challenges for high-quality sequencing, even using high-throughput, state-of-the-art techniques. Here, we review the mechanisms involved in sex determination and their role in the evolution of species, particularly vertebrates, focusing on sex chromosomes and the challenges involved in sequencing these genomic elements. We also address the improvements provided by the growth of sequencing projects, by generating a massive number of near-gapless, telomere-to-telomere, chromosome-level, phased assemblies, increasing the number and quality of sex-chromosome sequences available for further studies.

Sex-specific splicing of Z- and W-borne nr5a1 alleles suggests sex determination is controlled by chromosome conformation

Pogona vitticeps has female heterogamety (ZZ/ZW), but the master sex-determining gene is unknown, as it is for all reptiles. We show that nr5a1 (Nuclear Receptor Subfamily 5 Group A Member 1), a gene that is essential in mammalian sex determination, has alleles on the Z and W chromosomes (Z-nr5a1 and W-nr5a1), which are both expressed and can recombine. Three transcript isoforms of Z-nr5a1 were detected in gonads of adult ZZ males, two of which encode a functional protein. However, ZW females produced 16 isoforms, most of which contained premature stop codons. The array of transcripts produced by the W-borne allele (W-nr5a1) is likely to produce truncated polypeptides that contain a structurally normal DNA-binding domain and could act as a competitive inhibitor to the full-length intact protein. We hypothesize that an altered configuration of the W chromosome affects the conformation of the primary transcript generating inhibitory W-borne isoforms that suppress testis determination. Under this hypothesis, the genetic sex determination (GSD) system of P. vitticeps is a W-borne dominant female-determining gene that may be controlled epigenetically.

Abrupt conclusion of the late Miocene-early Pliocene biogenic bloom at 4.6-4.4 Ma

The late Miocene-early Pliocene biogenic bloom was an extended time interval characterised by elevated ocean export productivity at numerous locations. As primary productivity is nutrient-limited at low-to-mid latitudes, this bloom has been attributed to an increase or a redistribution of available nutrients, potentially involving ocean-gateway or monsoon-related mechanisms. While the exact causal feedbacks remain debated, there is even less consensus on what caused the end of the biogenic bloom. Here, we compile Mio-Pliocene paleoproductivity proxy data from all major ocean basins to evaluate the timing and pacing of this termination. This systematic analysis reveals an abrupt and sustained reduction in low-latitude ocean productivity at 4.6-4.4 Ma. The decline in productivity coincided with a prolonged period of low orbital eccentricity and a shift towards lower-amplitude obliquity, an astronomical configuration linked to reduced East Asian Monsoon intensity and decreased riverine nutrient supply.

Pyridylpiperazine-based allosteric inhibitors of RND-type multidrug efflux pumps

Efflux transporters of the RND family confer resistance to multiple antibiotics in Gram-negative bacteria. Here, we identify and chemically optimize pyridylpiperazine-based compounds that potentiate antibiotic activity in E. coli through inhibition of its primary RND transporter, AcrAB-TolC. Characterisation of resistant E. coli mutants and structural biology analyses indicate that the compounds bind to a unique site on the transmembrane domain of the AcrB L protomer, lined by key catalytic residues involved in proton relay. Molecular dynamics simulations suggest that the inhibitors access this binding pocket from the cytoplasm via a channel exclusively present in the AcrB L protomer. Thus, our work unveils a class of allosteric efflux-pump inhibitors that likely act by preventing the functional catalytic cycle of the RND pump.

Efflux transporters of the RND family confer resistance to multiple antibiotics in Gram-negative bacteria. Here, the authors identify pyridylpiperazine-based compounds that potentiate antibiotic activity in E. coli through allosteric inhibition of its primary RND transporter.

Applications and potentials of nanopore sequencing in the (epi)genome and (epi)transcriptome era

The Human Genome Project opened an era of (epi)genomic research, and also provided a platform for the development of new sequencing technologies. During and after the project, several sequencing technologies continue to dominate nucleic acid sequencing markets. Currently, Illumina (short-read), PacBio (long-read), and Oxford Nanopore (long-read) are the most popular sequencing technologies. Unlike PacBio or the popular short-read sequencers before it, which, as examples of the second or so-called Next-Generation Sequencing platforms, need to synthesize when sequencing, nanopore technology directly sequences native DNA and RNA molecules. Nanopore sequencing, therefore, avoids converting mRNA into cDNA molecules, which not only allows for the sequencing of extremely long native DNA and full-length RNA molecules but also document modifications that have been made to those native DNA or RNA bases. In this review on direct DNA sequencing and direct RNA sequencing using Oxford Nanopore technology, we focus on their development and application achievements, discussing their challenges and future perspective. We also address the problems researchers may encounter applying these approaches in their research topics, and how to resolve them.

Chromosome-scale reference genome assembly of a diploid potato clone derived from an elite variety

Potato (Solanum tuberosum L.) is one of the most important crops with a worldwide production of 370 million metric tons. The objectives of this study were (1) to create a high-quality consensus sequence across the two haplotypes of a diploid clone derived from a tetraploid elite variety and assess the sequence divergence from the available potato genome assemblies, as well as among the two haplotypes; (2) to evaluate the new assembly’s usefulness for various genomic methods; and (3) to assess the performance of phasing in diploid and tetraploid clones, using linked-read sequencing technology. We used PacBio long reads coupled with 10x Genomics reads and proximity ligation scaffolding to create the dAg1_v1.0 reference genome sequence. With a final assembly size of 812 Mb, where 750 Mb are anchored to 12 chromosomes, our assembly is larger than other available potato reference sequences and high proportions of properly paired reads were observed for clones unrelated by pedigree to dAg1. Comparisons of the new dAg1_v1.0 sequence to other potato genome sequences point out the high divergence between the different potato varieties and illustrate the potential of using dAg1_v1.0 sequence in breeding applications.

A near complete genome for goat genetic and genomic research

Background

Goat, one of the first domesticated livestock, is a worldwide important species both culturally and economically. The current goat reference genome, known as ARS1, is reported as the first nonhuman genome assembly using 69× PacBio sequencing. However, ARS1 suffers from incomplete X chromosome and highly fragmented Y chromosome scaffolds.

Results

Here, we present a very high-quality de novo genome assembly, Saanen_v1, from a male Saanen dairy goat, with the first goat Y chromosome scaffold based on 117× PacBio long-read sequencing and 118× Hi-C data. Saanen_v1 displays a high level of completeness thanks to the presence of centromeric and telomeric repeats at the proximal and distal ends of two-thirds of the autosomes, and a much reduced number of gaps (169 vs. 773). The completeness and accuracy of the Saanen_v1 genome assembly are also evidenced by more assembled sequences on the chromosomes (2.63 Gb for Saanen_v1 vs. 2.58 Gb for ARS1), a slightly increased mapping ratio for transcriptomic data, and more genes anchored to chromosomes. The eight putative large assembly errors (1 to ~ 7 Mb each) found in ARS1 were amended, and for the first time, the substitution rate of this ruminant Y chromosome was estimated. Furthermore, sequence improvement in Saanen_v1, compared with ARS1, enables us to assign the likely correct positions for 4.4% of the single nucleotide polymorphism (SNP) probes in the widely used GoatSNP50 chip.

Conclusions

The updated goat genome assembly including both sex chromosomes (X and Y) and the autosomes with high-resolution quality will serve as a valuable resource for goat genetic research and applications.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00668-5.

Studying the surfaces of bacteria using neutron scattering: finding new openings for antibiotics

The use of neutrons as a scattering probe to investigate biological membranes has steadily grown in the past three decades, shedding light on the structure and behaviour of this ubiquitous and fundamental biological barrier. Meanwhile, the rise of antibiotic resistance has catalysed a renewed interest in understanding the mechanisms underlying the dynamics of antibiotics interaction with the bacterial cell envelope. It is widely recognised that the key reason behind the remarkable success of Gram-negative pathogens in developing antibiotic resistance lies in the effectiveness of their outer membrane (OM) in defending the cell from antibacterial compounds. Critical to its function, the highly asymmetric lipid distribution between the inner and outer bilayer leaflets of the OM, adds an extra level of complexity to the study of this crucial defence barrier. Here we review the opportunities offered by neutron scattering techniques, in particular reflectometry, to provide structural information on the interactions of antimicrobials with in vitro models of the OM. The differential sensitivity of neutrons towards hydrogen and deuterium makes them a unique probe to study the structure and behaviour of asymmetric membranes. Molecular-level understanding of the interactions between antimicrobials and the Gram-negative OM provides valuable insights that can aid drug development and broaden our knowledge of this critically important biological barrier.

Using Combined Methods of Genetic Mapping and Nanopore-Based Sequencing Technology to Analyze the Insertion Positions of G10evo-EPSPS and Cry1Ab/Cry2Aj Transgenes in Maize

The insertion position of the exogenous fragment sequence in a genetically modified organism (GMO) is important for the safety assessment and labeling of GMOs. SK12-5 is a newly developed transgenic maize line transformed with two trait genes [i.e., G10evo-5-enolpyrul-shikimate-3-phosphate synthase (EPSPS) and Cry1Ab/Cry2Aj] that was recently approved for commercial use in China. In this study, we tried to determine the insertion position of the exogenous fragment for SK12-5. The transgene-host left border and right border integration junctions were obtained from SK12-5 genomic DNA by using the thermal asymmetric interlaced polymerase chain reaction (TAIL-PCR) and next-generation Illumina sequencing technology. However, a Basic Local Alignment Search Tool (BLAST) analysis revealed that the flanking sequences in the maize genome are unspecific and that the insertion position is located in a repetitive sequence area in the maize genome. To locate the fine-scale insertion position in SK12-5, we combined the methods of genetic mapping and nanopore-based sequencing technology. From a classical bulked-segregant analysis (BSA), the insertion position in SK12-5 was mapped onto Bin9.03 of chromosome 9 between the simple sequence repeat (SSR) markers umc2337 and umc1743 (26,822,048-100,724,531 bp). The nanopore sequencing results uncovered 10 reads for which one end was mapped onto the vector and the other end was mapped onto the maize genome. These observations indicated that the exogenous T-DNA fragments were putatively integrated at the position from 82,329,568 to 82,379,296 bp of chromosome 9 in the transgenic maize SK12-5. This study is helpful for the safety assessment of the novel transgenic maize SK12-5 and shows that the combined method of genetic mapping and the nanopore-based sequencing technology will be a useful approach for identifying the insertion positions of transgenic sequences in other GM plants with relatively large and complex genomes.

Genomic Tackling of Human Satellite DNA: Breaking Barriers through Time

(Peri)centromeric repetitive sequences and, more specifically, satellite DNA (satDNA) sequences, constitute a major human genomic component. SatDNA sequences can vary on a large number of features, including nucleotide composition, complexity, and abundance. Several satDNA families have been identified and characterized in the human genome through time, albeit at different speeds. Human satDNA families present a high degree of sub-variability, leading to the definition of various subfamilies with different organization and clustered localization. Evolution of satDNA analysis has enabled the progressive characterization of satDNA features. Despite recent advances in the sequencing of centromeric arrays, comprehensive genomic studies to assess their variability are still required to provide accurate and proportional representation of satDNA (peri)centromeric/acrocentric short arm sequences. Approaches combining multiple techniques have been successfully applied and seem to be the path to follow for generating integrated knowledge in the promising field of human satDNA biology.

The Genomic Architecture and Evolutionary Fates of Supergenes

Supergenes are genomic regions containing sets of tightly linked loci that control multi-trait phenotypic polymorphisms under balancing selection. Recent advances in genomics have uncovered significant variation in both the genomic architecture as well as the mode of origin of supergenes across diverse organismal systems. Although the role of genomic architecture for the origin of supergenes has been much discussed, differences in the genomic architecture also subsequently affect the evolutionary trajectory of supergenes and the rate of degeneration of supergene haplotypes. In this review, we synthesize recent genomic work and historical models of supergene evolution, highlighting how the genomic architecture of supergenes affects their evolutionary fate. We discuss how recent findings on classic supergenes involved in governing ant colony social form, mimicry in butterflies, and heterostyly in flowering plants relate to theoretical expectations. Furthermore, we use forward simulations to demonstrate that differences in genomic architecture affect the degeneration of supergenes. Finally, we discuss implications of the evolution of supergene haplotypes for the long-term fate of balanced polymorphisms governed by supergenes.

Noninvasive prenatal testing for β-thalassemia by targeted nanopore sequencing combined with relative haplotype dosage (RHDO): a feasibility study

Noninvasive prenatal testing (NIPT) for single gene disorders remains challenging. One approach that allows for accurate detection of the slight increase of the maternally inherited allele is the relative haplotype dosage (RHDO) analysis, which requires the construction of parental haplotypes. Recently, the nanopore sequencing technologies have become available and may be an ideal tool for direct construction of haplotypes. Here, we explored the feasibility of combining nanopore sequencing with the RHDO analysis in NIPT of β-thalassemia. Thirteen families at risk for β-thalassemia were recruited. Targeted region of parental genomic DNA was amplified by long-range PCR of 10 kb and 20 kb amplicons. Parental haplotypes were constructed using nanopore sequencing and next generation sequencing data. Fetal inheritance of parental haplotypes was classified by the RHDO analysis using data from maternal plasma DNA sequencing. Haplotype phasing was achieved in 12 families using data from 10 kb library. While data from the 20 kb library gave a better performance that haplotype phasing was achieved in all 13 families. Fetal status was correctly classified in 12 out of 13 families. Thus, targeted nanopore sequencing combined with the RHDO analysis is feasible to NIPT for β-thalassemia.

Efficient assembly of nanopore reads via highly accurate and intact error correction

Long nanopore reads are advantageous in de novo genome assembly. However, nanopore reads usually have broad error distribution and high-error-rate subsequences. Existing error correction tools cannot correct nanopore reads efficiently and effectively. Most methods trim high-error-rate subsequences during error correction, which reduces both the length of the reads and contiguity of the final assembly. Here, we develop an error correction, and de novo assembly tool designed to overcome complex errors in nanopore reads. We propose an adaptive read selection and two-step progressive method to quickly correct nanopore reads to high accuracy. We introduce a two-stage assembler to utilize the full length of nanopore reads. Our tool achieves superior performance in both error correction and de novo assembling nanopore reads. It requires only 8122 hours to assemble a 35X coverage human genome and achieves a 2.47-fold improvement in NG50. Furthermore, our assembly of the human WERI cell line shows an NG50 of 22 Mbp. The high-quality assembly of nanopore reads can significantly reduce false positives in structure variation detection.

A Hu sheep genome with the first ovine Y chromosome reveal introgression history after sheep domestication

The Y chromosome plays key roles in male fertility and reflects the evolutionary history of paternal lineages. Here, we present a de novo genome assembly of the Hu sheep with the first draft assembly of ovine Y chromosome (oMSY), using nanopore sequencing and Hi-C technologies. The oMSY that we generated spans 10.6 Mb from which 775 Y-SNPs were identified by applying a large panel of whole genome sequences from worldwide sheep and wild Iranian mouflons. Three major paternal lineages (HY1a, HY1b and HY2) were defined across domestic sheep, of which HY2 was newly detected. Surprisingly, HY2 forms a monophyletic clade with the Iranian mouflons and is highly divergent from both HY1a and HY1b. Demographic analysis of Y chromosomes, mitochondrial and nuclear genomes confirmed that HY2 and the maternal counterpart of lineage C represented a distinct wild mouflon population in Iran that diverge from the direct ancestor of domestic sheep, the wild mouflons in Southeastern Anatolia. Our results suggest that wild Iranian mouflons had introgressed into domestic sheep and thereby introduced this Iranian mouflon specific lineage carrying HY2 to both East Asian and Africa sheep populations.

Establishment and evolution of heterochromatin

The eukaryotic genome is packaged into transcriptionally active euchromatin and silent heterochromatin, with most studies focused on the former encompassing the majority of protein-coding genes. The recent development of various sequencing techniques has refined this classic dichromatic partition and has better illuminated the composition, establishment, and evolution of this genomic and epigenomic "dark matter" in the context of topologically associated domains and phase-separated droplets. Heterochromatin includes genomic regions that can be densely stained by chemical dyes, which have been shown to be enriched for repetitive elements and epigenetic marks, including H3K9me2/3 and H3K27me3. Heterochromatin is usually replicated late, concentrated at the nuclear periphery or around nucleoli, and usually lacks highly expressed genes; and now it is considered to be as neither genetically inert nor developmentally static. Heterochromatin guards genome integrity against transposon activities and exerts important regulatory functions by targeting beyond its contained genes. Both its nucleotide sequences and regulatory proteins exhibit rapid coevolution between species. In addition, there are dynamic transitions between euchromatin and heterochromatin during developmental and evolutionary processes. We summarize here the ever-changing characteristics of heterochromatin and propose models and principles for the evolutionary transitions of heterochromatin that have been mainly learned from studies of Drosophila and yeast. Finally, we highlight the role of sex chromosomes in studying heterochromatin evolution.

Genetics and Epigenetics of Sex Bias: Insights from Human Cancer and Autoimmunity

High-throughput sequencing and genome-wide association studies have revealed a sex bias in human diseases. The underlying molecular mechanisms remain, however, unknown. Here, we cover recent advances in cancer and autoimmunity focusing on intrinsic genetic and epigenetic differences underlying sex biases in human disease. These studies reveal a central role of genome regulatory mechanisms including genome repair, chromosome folding, and epigenetic regulation in dictating the sex bias. These highlight the importance of considering sex as a variable in both basic science and clinical investigations. Understanding the molecular mechanisms underlying sex bias in human diseases will be instrumental in making a first step forwards into the era of personalized medicine.

One fly-one genome: chromosome-scale genome assembly of a single outbred Drosophila melanogaster

A high quality genome assembly is a vital first step for the study of an organism. Recent advances in technology have made the creation of high quality chromosome scale assemblies feasible and low cost. However, the amount of input DNA needed for an assembly project can be a limiting factor for small organisms or precious samples. Here we demonstrate the feasibility of creating a chromosome scale assembly using a hybrid method for a low input sample, a single outbred Drosophila melanogaster. Our approach combines an Illumina shotgun library, Oxford nanopore long reads, and chromosome conformation capture for long range scaffolding. This single fly genome assembly has a N50 of 26 Mb, a length that encompasses entire chromosome arms, contains 95% of expected single copy orthologs, and a nearly complete assembly of this individual's Wolbachia endosymbiont. The methods described here enable the accurate and complete assembly of genomes from small, field collected organisms as well as precious clinical samples.

Highly efficient synchronization of sheep skin fibroblasts at G2/M phase and isolation of sheep Y chromosomes by flow cytometric sorting

At present, based on whole genome sequencing, sequences and genes annotation of the sheep (Ovis aries) Y chromosome are still absent. The isolation of Y chromosomes followed by sequencing has been approved as an effective approach to analyze this complex chromosome in other species. In this study, we established a highly efficient synchronization method for G2/M phase of sheep fibroblasts, which was successfully applied to flow-sorting chromosomes of sheep, with a focus on isolation and sequencing of the ovine Y chromosome. The isolated (~80,000) Y chromosomes were verified by fluorescence quantitative real-time polymerase chain reaction, further confirmed by fluorescence in situ hybridization, and amplified by the MALBAC method before next-generation sequencing. The sequence results indicated that 68.90% of reads were Y chromosome-related sequences as they are homologous to the bovine Y chromosome. The remaining 31.1% of reads were aligned to the sheep reference genome, including 13.57% reads to chromosome X and 6.68% to chromosome 17. Importantly, the paired-end reads that are properly aligned to the bovine Y sequence assembly accounted for 46.49%, indicating the success in the ovine Y chromosome isolation and the high quality of the Y chromosome sequences. This study not only set up a foundation for future sequencing, assembly and annotation of the ovine Y chromosome, but also provide a validated approach to overcoming difficulties in sequencing Y chromosome in other mammalian species.

All Roads Lead to Auxin: Post-translational Regulation of Auxin Transport by Multiple Hormonal Pathways

Auxin is a key hormonal regulator, that governs plant growth and development in concert with other hormonal pathways. The unique feature of auxin is its polar, cell-to-cell transport that leads to the formation of local auxin maxima and gradients, which coordinate initiation and patterning of plant organs. The molecular machinery mediating polar auxin transport is one of the important points of interaction with other hormones. Multiple hormonal pathways converge at the regulation of auxin transport and form a regulatory network that integrates various developmental and environmental inputs to steer plant development. In this review, we discuss recent advances in understanding the mechanisms that underlie regulation of polar auxin transport by multiple hormonal pathways. Specifically, we focus on the post-translational mechanisms that contribute to fine-tuning of the abundance and polarity of auxin transporters at the plasma membrane and thereby enable rapid modification of the auxin flow to coordinate plant growth and development.

Analysis of the Mouse Y Chromosome by Single-Molecule Sequencing With Y Chromosome Enrichment

Since human and mouse Y chromosomes contain repeated sequences, it is difficult to determine the precise sequences and analyze the function of individual Y chromosome genes. Therefore, the causes of many diseases and abnormalities related to Y chromosome genes, such as male infertility, remain unclear. In this study, to elucidate the mouse Y chromosome, we enriched the mouse Y chromosome using a fluorescence-activated cell sorter (FACS) equipped with commonly used UV and blue 488 nm lasers and read the nucleotides using the Oxford Nanopore MinION long-read sequencer. This sequencing strategy allows us to cover the whole known region as well as the potential undetermined region of the Y chromosome. FACS-based chromosome enrichment and long-read sequencing are suitable for analysis of the Y chromosome sequences and may lead to further understanding of the physiological role of Y chromosome genes.

The Genome of the Endangered Dryas Monkey Provides New Insights into the Evolutionary History of the Vervets

Genomic data can be a powerful tool for inferring ecology, behavior, and conservation needs of highly elusive species, particularly, when other sources of information are hard to come by. Here, we focus on the Dryas monkey (Cercopithecus dryas), an endangered primate endemic to the Congo Basin with cryptic behavior and possibly <250 remaining adult individuals. Using whole-genome sequencing data, we show that the Dryas monkey represents a sister lineage to the vervets (Chlorocebus sp.) and has diverged from them ∼1.4 Ma with additional bidirectional gene flow ∼750,000–∼500,000 years ago that has likely involved the crossing of the Congo River. Together with evidence of gene flow across the Congo River in bonobos and okapis, our results suggest that the fluvial topology of the Congo River might have been more dynamic than previously recognized. Despite the presence of several homozygous loss-of-function mutations in genes associated with sperm mobility and immunity, we find high genetic diversity and low levels of inbreeding and genetic load in the studied Dryas monkey individual. This suggests that the current population carries sufficient genetic variability for long-term survival and might be larger than currently recognized. We thus provide an example of how genomic data can directly improve our understanding of highly elusive species.

Dynamics of Transposable Element Invasions with piRNA Clusters

In mammals and invertebrates, the proliferation of an invading transposable element (TE) is thought to be stopped by an insertion into a piRNA cluster. Here, we explore the dynamics of TE invasions under this trap model using computer simulations. We found that piRNA clusters confer a substantial benefit, effectively preventing extinction of host populations from a proliferation of deleterious TEs. TE invasions consist of three distinct phases: first, the TE amplifies within the population, next TE proliferation is stopped by segregating cluster insertions, and finally the TE is inactivated by fixation of a cluster insertion. Suppression by segregating cluster insertions is unstable and bursts of TE activity may yet occur. The transposition rate and the population size mostly influence the length of the phases but not the amount of TEs accumulating during an invasion. Solely, the size of piRNA clusters was identified as a major factor influencing TE abundance. We found that a single nonrecombining cluster is more efficient in stopping invasions than clusters distributed over several chromosomes. Recombination among cluster sites makes it necessary that each diploid carries, on the average, four cluster insertions to stop an invasion. Surprisingly, negative selection in a model with piRNA clusters can lead to a novel equilibrium state, where TE copy numbers remain stable despite only some individuals in a population carrying a cluster insertion. In Drosophila melanogaster, the trap model accounts for the abundance of TEs produced in the germline but fails to predict the abundance of TEs produced in the soma.

Decoding the Role of Satellite DNA in Genome Architecture and Plasticity-An Evolutionary and Clinical Affair

Repetitive DNA is a major organizational component of eukaryotic genomes, being intrinsically related with their architecture and evolution. Tandemly repeated satellite DNAs (satDNAs) can be found clustered in specific heterochromatin-rich chromosomal regions, building vital structures like functional centromeres and also dispersed within euchromatin. Interestingly, despite their association to critical chromosomal structures, satDNAs are widely variable among species due to their high turnover rates. This dynamic behavior has been associated with genome plasticity and chromosome rearrangements, leading to the reshaping of genomes. Here we present the current knowledge regarding satDNAs in the light of new genomic technologies, and the challenges in the study of these sequences. Furthermore, we discuss how these sequences, together with other repeats, influence genome architecture, impacting its evolution and association with disease.

Flow Sorting Enrichment and Nanopore Sequencing of Chromosome 1 From a Chinese Individual

Sorting of individual chromosomes by Flow Cytometry (flow-sorting) is an enrichment method to potentially simplify genome assembly by isolating chromosomes from the context of the genome. We have recently developed a workflow to sequence native, unamplified DNA and applied it to the smallest human chromosome, the Y chromosome. Here, we modify improve upon that workflow to increase DNA recovery from chromosome sorting as well as sequencing yield. We apply it to sequence and assemble the largest human chromosome - chromosome 1 - of a Chinese individual using a single Oxford Nanopore MinION flow cell. We generate a selective and highly continuous assembly whose continuity reaches into the order of magnitude of the human reference GRCh38. We then use this assembly to call candidate structural variants against the reference and find 685 putative novel SV candidates. We propose this workflow as a potential solution to assemble structurally complex chromosomes, or the study of very large plant or animal genomes that might challenge traditional assembly strategies.

New insights into mammalian sex chromosome structure and evolution using high-quality sequences from bovine X and Y chromosomes

BACKGROUND

Mammalian X chromosomes are mainly euchromatic with a similar size and structure among species whereas Y chromosomes are smaller, have undergone substantial evolutionary changes and accumulated male specific genes and genes involved in sex determination. The pseudoautosomal region (PAR) is conserved on the X and Y and pair during meiosis. The structure, evolution and function of mammalian sex chromosomes, particularly the Y chromsome, is still poorly understood because few species have high quality sex chromosome assemblies.

RESULTS

Here we report the first bovine sex chromosome assemblies that include the complete PAR spanning 6.84 Mb and three Y chromosome X-degenerate (X-d) regions. The PAR comprises 31 genes, including genes that are missing from the X chromosome in current cattle, sheep and goat reference genomes. Twenty-nine PAR genes are single-copy genes and two are multi-copy gene families, OBP, which has 3 copies and BDA20, which has 4 copies. The Y chromosome X-d1, 2a and 2b regions contain 11, 2 and 2 gametologs, respectively.

CONCLUSIONS

The ruminant PAR comprises 31 genes and is similar to the PAR of pig and dog but extends further than those of human and horse. Differences in the pseudoautosomal boundaries are consistent with evolutionary divergence times. A bovidae-specific expansion of members of the lipocalin gene family in the PAR reported here, may affect immune-modulation and anti-inflammatory responses in ruminants. Comparison of the X-d regions of Y chromosomes across species revealed that five of the X-Y gametologs, which are known to be global regulators of gene activity and candidate sexual dimorphism genes, are conserved.

The Challenges of Chromosome Y Analysis and the Implications for Chronic Kidney Disease

The role of chromosome Y in chronic kidney disease (CKD) remains unknown, as chromosome Y is typically excluded from genetic analysis in CKD. The complex, sex-specific presentation of CKD could be influenced by chromosome Y genetic variation, but there is limited published research available to confirm or reject this hypothesis. Although traditionally thought to be associated with male-specific disease, evidence linking chromosome Y genetic variation to common complex disorders highlights a potential gap in CKD research. Chromosome Y variation has been associated with cardiovascular disease, a condition closely linked to CKD and one with a very similar sexual dimorphism. Relatively few sources of genetic variation in chromosome Y have been examined in CKD. The association between chromosome Y aneuploidy and CKD has never been explored comprehensively, while analyses of microdeletions, copy number variation, and single-nucleotide polymorphisms in CKD have been largely limited to the autosomes or chromosome X. In many studies, it is unclear whether the analyses excluded chromosome Y or simply did not report negative results. Lack of imputation, poor cross-study comparability, and requirement for separate or additional analyses in comparison with autosomal chromosomes means that chromosome Y is under-investigated in the context of CKD. Limitations in genotyping arrays could be overcome through use of whole-chromosome sequencing of chromosome Y that may allow analysis of many different types of genetic variation across the chromosome to determine if chromosome Y genetic variation is associated with CKD.

Chromosomics: Bridging the Gap between Genomes and Chromosomes

The recent advances in DNA sequencing technology are enabling a rapid increase in the number of genomes being sequenced. However, many fundamental questions in genome biology remain unanswered, because sequence data alone is unable to provide insight into how the genome is organised into chromosomes, the position and interaction of those chromosomes in the cell, and how chromosomes and their interactions with each other change in response to environmental stimuli or over time. The intimate relationship between DNA sequence and chromosome structure and function highlights the need to integrate genomic and cytogenetic data to more comprehensively understand the role genome architecture plays in genome plasticity. We propose adoption of the term 'chromosomics' as an approach encompassing genome sequencing, cytogenetics and cell biology, and present examples of where chromosomics has already led to novel discoveries, such as the sex-determining gene in eutherian mammals. More importantly, we look to the future and the questions that could be answered as we enter into the chromosomics revolution, such as the role of chromosome rearrangements in speciation and the role more rapidly evolving regions of the genome, like centromeres, play in genome plasticity. However, for chromosomics to reach its full potential, we need to address several challenges, particularly the training of a new generation of cytogeneticists, and the commitment to a closer union among the research areas of genomics, cytogenetics, cell biology and bioinformatics. Overcoming these challenges will lead to ground-breaking discoveries in understanding genome evolution and function.

Evolutionary and functional analysis of RBMY1 gene copy number variation on the human Y chromosome

Human RBMY1 genes are located in four variable-sized clusters on the Y chromosome, expressed in male germ cells and possibly associated with sperm motility. We have re-investigated the mutational background and evolutionary history of the RBMY1 copy number distribution in worldwide samples and its relevance to sperm parameters in an Estonian cohort of idiopathic male factor infertility subjects. We estimated approximate RBMY1 copy numbers in 1218 1000 Genomes Project phase 3 males from sequencing read-depth, then chose 14 for valid ation by multicolour fibre-FISH. These fibre-FISH samples provided accurate calibration standards for the entire panel and led to detailed insights into population variation and mutational mechanisms. RBMY1 copy number worldwide ranged from 3 to 13 with a mode of 8. The two larger proximal clusters were the most variable, and additional duplications, deletions and inversions were detected. Placing the copy number estimates onto the published Y-SNP-based phylogeny of the same samples suggested a minimum of 562 mutational changes, translating to a mutation rate of 2.20 × 10-3 (95% CI 1.94 × 10-3 to 2.48 × 10-3) per father-to-son Y-transmission, higher than many short tandem repeat (Y-STRs), and showed no evidence for selection for increased or decreased copy number, but possible copy number stabilizing selection. An analysis of RBMY1 copy numbers among 376 infertility subjects failed to replicate a previously reported association with sperm motility and showed no significant effect on sperm count and concentration, serum follicle stimulating hormone (FSH), luteinizing hormone (LH) and testosterone levels or testicular and semen volume. These results provide the first in-depth insights into the structural rearrangements underlying RBMY1 copy number variation across diverse human lineages.