Large scale assessment of regulatory evolution and transcriptome complexity in mammals

AbstractIn addition to genetic changes affecting the function of gene products, changes in gene expression have been suggested to underlie many or even most of the phenotypic differences among mammals. However, detailed gene expression comparisons were, until recently, restricted to closely related species, owing to technological limitations. Thus, we took advantage of the latest technologies (RNA-Seq) to generate extensive qualitative and quantitative transcriptome data for a unique collection of somatic and germline tissues from representatives of all major mammalian lineages (placental mammals, marsupials and monotremes) and birds, the evolutionary outgroup.In the first major project of my thesis, we performed global comparative analyses of gene expression levels based on these data. Our analyses provided fundamental insights into the dynamics of transcriptome change during mammalian evolution (e.g., the rate of expression change across species, tissues and chromosomes) and allowed the exploration of the functional relevance and phenotypic implications of transcription changes at a genome-wide scale (e.g., we identified numerous potentially selectively driven expression switches).In a second project of my thesis, which was also based on the unique transcriptome data generated in the context of the first project we focused on the evolution of alternative splicing in mammals. Alternative splicing contributes to transcriptome complexity by generating several transcript isoforms from a single gene, which can, thus, perform various functions. To complete the global comparative analysis of gene expression changes, we explored patterns of alternative splicing evolution. This work uncovered several general and unexpected patterns of alternative splicing evolution (e.g., we found that alternative splicing evolves extremely rapidly) as well as a large number of conserved alternative isoforms that may be crucial for the functioning of mammalian organs.Finally, the third and final project of my PhD consisted in analyzing in detail the unique functional and evolutionary properties of the testis by exploring the extent of its transcriptome complexity. This organ was previously shown to evolve rapidly both at the phenotypic and molecular level, apparently because of the specific pressures that act on this organ and are associated with its reproductive function. Moreover, my analyses of the amniote tissue transcriptome data described above, revealed strikingly widespread transcriptional activity of both functional and nonfunctional genomic elements in the testis compared to the other organs. To elucidate the cellular source and mechanisms underlying this promiscuous transcription in the testis, we generated deep coverage RNA-Seq data for all major testis cell types as well as epigenetic data (DNA and histone methylation) using the mouse as model system. The integration of these complete dataset revealed that meiotic and especially post-meiotic germ cells are the major contributors to the widespread functional and nonfunctional transcriptome complexity of the testis, and that this "promiscuous" spermatogenic transcription is resulting, at least partially, from an overall transcriptionally permissive chromatin state. We hypothesize that this particular open state of the chromatin results from the extensive chromatin remodeling that occurs during spermatogenesis which ultimately leads to the replacement of histones by protamines in the mature spermatozoa. Our results have important functional and evolutionary implications (e.g., regarding new gene birth and testicular gene expression evolution).Generally, these three large-scale projects of my thesis provide complete and massive datasets that constitute valuables resources for further functional and evolutionary analyses of mammalian genomes.

Rensch's rule in large herbivorous mammals derived from metabolic scaling

Rensch’s rule, which states that the magnitude of sexual size dimorphism tends to increase with increasing body size, has evolved independently in three lineages of large herbivorous mammals: bovids (antelopes), cervids (deer), and macropodids (kangaroos). This pattern can be explained by a model that combines allometry,life-history theory, and energetics. The key features are thatfemale group size increases with increasing body size and that males have evolved under sexual selection to grow large enough to control these groups of females. The model predicts relationships among body size and female group size, male and female age at first breeding,death and growth rates, and energy allocation of males to produce body mass and weapons. Model predictions are well supported by data for these megaherbivores. The model suggests hypotheses for why some other sexually dimorphic taxa, such as primates and pinnipeds(seals and sea lions), do or do not conform to Rensh’s rule.

Comparing methods for sampling large- and medium-sized mammals: camera traps and track plots

Activities involving fauna monitoring are usually limited by the lack of resources; therefore, the choice of a proper and efficient methodology is fundamental to maximize the cost-benefit ratio. Both direct and indirect methods can be used to survey mammals, but the latter are preferred due to the difficulty to come in sight of and/or to capture the individuals, besides being cheaper. We compared the performance of two methods to survey medium and large-sized mammal: track plot recording and camera trapping, and their costs were assessed. At Jatai Ecological Station (S21 degrees 31`15 ``- W47 degrees 34`42 ``-Brazil) we installed ten camera traps along a dirt road directly in front of ten track plots, and monitored them for 10 days. We cleaned the plots, adjusted the cameras, and noted down the recorded species daily. Records taken by both methods showed they sample the local richness in different ways (Wilcoxon, T=231; p;;0.01). The track plot method performed better on registering individuals whereas camera trapping provided records which permitted more accurate species identification. The type of infra-red sensor camera used showed a strong bias towards individual body mass (R(2)=0.70; p=0.017), and the variable expenses of this method in a 10-day survey were estimated about 2.04 times higher compared to track plot method; however, in a long run camera trapping becomes cheaper than track plot recording. Concluding, track plot recording is good enough for quick surveys under a limited budget, and camera trapping is best for precise species identification and the investigation of species details, performing better for large animals. When used together, these methods can be complementary.

A Survey of mid and large bodied mammals in Núcleo Caraguatatuba, Serra do Mar State Park, Brazil

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Medium and large-sized mammals in a private reserve of natural heritage in the municipality of Jaú, São Paulo, Brazil

Medium and large-sized mammals control invertebrate and vertebrate populations and are important seed dispersers. These animals are losing their territory and living range due to deforestation and the fragmentation of natural areas, particularly in the Atlantic Forest and Cerrado. The Private Reserve of Natural Heritage (PRNH) Amadeu Botelho is isolated from other forest fragments and is thus a wildlife refuge in the region of Jaú, central-western region of state of São Paulo, southeastem Brazil. Animals as such as the puma (Puma concolor), ocelot (Leopardus pardalis) and jaguarundi (Puma yagouaroundi) have been sighted in the reserve by wildlife guards and residents. This reserve, which belongs to “Santo Antônio dos Ipês” farm, is surrounded by several plantations, and is very close to the center of the city of Jaú. In addition to surveying the mammal species in the reserve, this study diagnosed the influence of the farm’s activities and the proximity to the center of the municipality of Jaú on the fauna under study

Directions for preparing scientific specimens of large and small mammals, birds, birds' stomachs for economic investigations, birds' nests and eggs, fish and reptiles /

"Issued January 15, 1914."

Cryptic vicariance in Gulf of California fishes parallels vicariant patterns found in Baja, California mammals and reptiles

Comparisons across multiple taxa can often clarify the histories of biogeographic regions. In particular, historic barriers to movement should affect multiple species and, thus, result in a pattern of concordant intraspecific genetic divisions among species. A striking example of such comparative phylogeography is the recent observation that populations of many small mammals and reptiles living on the Baja, California peninsula have a large genetic break between northern and southern peninsular populations. In the present study, I demonstrate that five species of near-shore fishes living on the Baja coastline of the Gulf of California share this genetic pattern. The simplest explanation for this concordant genetic division within both terrestrial and marine vertebrates is that the Baja peninsula was fragmented by a Plio-Pleistocene marine seaway and that this seaway posed a substantial barrier to movement for near-shore fishes. The genetic divisions within Gulf of California fishes also coincide with recognized biogeographic regions based on fish community composition and several environmental factors. It is likely that adaptation to regional environments and present-day oceanographic circulation limits gene exchange between biogeographic regions and helps maintain evidence of past vicariance.

Geographic range size, life history and rates of diversification in Australian mammals

What causes species richness to vary among different groups of organisms? Two hypotheses are that large geographical ranges and fast life history either reduce extinction rates or raise speciation rates, elevating a clade's rate of diversification. Here we present a comparative analysis of these hypotheses using data on the phylogenetic relationships, geographical ranges and life history of the terrestrial mammal fauna of Australia. By comparing species richness patterns to null models, we show that species are distributed nonrandomly among genera. Using sister-clade comparisons to control for clade age, we then find that faster diversification is significantly associated with larger geographical ranges and larger litters, but there is no evidence for an effect of body size or age at first breeding on diversification rates. We believe the most likely explanation for these patterns is that larger litters and geographical ranges increase diversification rates because they buffer species from extinction. We also discuss the possibility that positive effects of litter size and range size on diversification rates result from elevated speciation rates.

Adiposity, hormone replacement therapy use and breast cancer risk by age and hormone receptor status: a large prospective cohort study.

INTRODUCTION Associations of hormone-receptor positive breast cancer with excess adiposity are reasonably well characterized; however, uncertainty remains regarding the association of body mass index (BMI) with hormone-receptor negative malignancies, and possible interactions by hormone replacement therapy (HRT) use. METHODS Within the European EPIC cohort, Cox proportional hazards models were used to describe the relationship of BMI, waist and hip circumferences with risk of estrogen-receptor (ER) negative and progesterone-receptor (PR) negative (n = 1,021) and ER+PR+ (n = 3,586) breast tumors within five-year age bands. Among postmenopausal women, the joint effects of BMI and HRT use were analyzed. RESULTS For risk of ER-PR- tumors, there was no association of BMI across the age bands. However, when analyses were restricted to postmenopausal HRT never users, a positive risk association with BMI (third versus first tertile HR = 1.47 (1.01 to 2.15)) was observed. BMI was inversely associated with ER+PR+ tumors among women aged ≤49 years (per 5 kg/m2 increase, HR = 0.79 (95%CI 0.68 to 0.91)), and positively associated with risk among women ≥65 years (HR = 1.25 (1.16 to 1.34)). Adjusting for BMI, waist and hip circumferences showed no further associations with risks of breast cancer subtypes. Current use of HRT was significantly associated with an increased risk of receptor-negative (HRT current use compared to HRT never use HR: 1.30 (1.05 to 1.62)) and positive tumors (HR: 1.74 (1.56 to 1.95)), although this risk increase was weaker for ER-PR- disease (Phet = 0.035). The association of HRT was significantly stronger in the leaner women (BMI ≤22.5 kg/m2) than for more overweight women (BMI ≥25.9 kg/m2) for, both, ER-PR- (HR: 1.74 (1.15 to 2.63)) and ER+PR+ (HR: 2.33 (1.84 to 2.92)) breast cancer and was not restricted to any particular HRT regime. CONCLUSIONS An elevated BMI may be positively associated with risk of ER-PR- tumors among postmenopausal women who never used HRT. Furthermore, postmenopausal HRT users were at an increased risk of ER-PR- as well as ER+PR+ tumors, especially among leaner women. For hormone-receptor positive tumors, but not for hormone-receptor negative tumors, our study confirms an inverse association of risk with BMI among young women of premenopausal age. Our data provide evidence for a possible role of sex hormones in the etiology of hormone-receptor negative tumors.

R-ESHAP as salvage therapy for patients with relapsed or refractory diffuse large B-cell lymphoma: the influence of prior exposure to rituximab on outcome. A GEL/TAMO study.

BACKGROUND The role of re-treatment with rituximab in aggressive B-cell lymphomas still needs to be defined. This study evaluated the influence of prior exposure to rituximab on response rates and survival in patients with diffuse large B-cell lymphoma treated with rituximab plus etoposide, cytarabine, cisplatinum and methylprednisolone (R-ESHAP). DESIGN AND METHODS We retrospectively analyzed 163 patients with relapsed or refractory diffuse large B-cell lymphoma who received R-ESHAP as salvage therapy with a curative purpose. Patients were divided into two groups according to whether rituximab had been administered (n=94, "R+" group) or not (n=69, "R-" group) prior to R-ESHAP. RESULTS Response rates were significantly higher in the R- group in the univariate but not in the multivariate analysis. In the analysis restricted to the R+ group, we observed very low complete remission and overall response rates in patients with primary refractory disease (8% and 33%, respectively), as compared to those in patients who were in first partial remission (41% and 86%) or who had relapsed disease (50% and 75%) (p<0.01 in both cases). Overall, 60% and 65% of patients in the R+ and R- groups, respectively, underwent stem-cell transplantation after the salvage therapy. With a median follow-up of 29 months (range, 6-84), patients in the R+ group had significantly worse progression-free survival (17% vs. 57% at 3 years, p<0.0001) and overall survival (38% v 67% at 3 years, p=0.0005) than patients in the R- group. Prior exposure to rituximab was also an independent adverse prognostic factor for both progression-free survival (RR: 2.0; 95% CI: 1.2-3.3, p=0.008) and overall survival (RR: 2.2; 95% CI: 1.3-3.9, p=0.004). CONCLUSIONS R-ESHAP was associated with a high response rate in patients who were not refractory to upfront rituximab-based chemotherapy. However, the survival outcome was poor for patients previously exposed to rituximab, as compared to in those who had not previously been treated with rituximab.

Effective diagnosis of Alzheimer's disease by means of large margin-based methodology.

BACKGROUND Functional brain images such as Single-Photon Emission Computed Tomography (SPECT) and Positron Emission Tomography (PET) have been widely used to guide the clinicians in the Alzheimer's Disease (AD) diagnosis. However, the subjectivity involved in their evaluation has favoured the development of Computer Aided Diagnosis (CAD) Systems. METHODS It is proposed a novel combination of feature extraction techniques to improve the diagnosis of AD. Firstly, Regions of Interest (ROIs) are selected by means of a t-test carried out on 3D Normalised Mean Square Error (NMSE) features restricted to be located within a predefined brain activation mask. In order to address the small sample-size problem, the dimension of the feature space was further reduced by: Large Margin Nearest Neighbours using a rectangular matrix (LMNN-RECT), Principal Component Analysis (PCA) or Partial Least Squares (PLS) (the two latter also analysed with a LMNN transformation). Regarding the classifiers, kernel Support Vector Machines (SVMs) and LMNN using Euclidean, Mahalanobis and Energy-based metrics were compared. RESULTS Several experiments were conducted in order to evaluate the proposed LMNN-based feature extraction algorithms and its benefits as: i) linear transformation of the PLS or PCA reduced data, ii) feature reduction technique, and iii) classifier (with Euclidean, Mahalanobis or Energy-based methodology). The system was evaluated by means of k-fold cross-validation yielding accuracy, sensitivity and specificity values of 92.78%, 91.07% and 95.12% (for SPECT) and 90.67%, 88% and 93.33% (for PET), respectively, when a NMSE-PLS-LMNN feature extraction method was used in combination with a SVM classifier, thus outperforming recently reported baseline methods. CONCLUSIONS All the proposed methods turned out to be a valid solution for the presented problem. One of the advances is the robustness of the LMNN algorithm that not only provides higher separation rate between the classes but it also makes (in combination with NMSE and PLS) this rate variation more stable. In addition, their generalization ability is another advance since several experiments were performed on two image modalities (SPECT and PET).

Large gangliocytic paraganglioma of the duodenum: A rare entity.

Gangliocytic paragangliomas are rare tumors that almost exclusively occur within the second portion of the duodenum. Although these tumors generally have a benign clinical course, they have the potential to recur or metastasize to regional lymph nodes. The case report presented here describes a 57-year-old female patient with melena, progressive asthenia, anemia, and a mass in the second-third portion of the duodenum that was treated by local excision. The patient was diagnosed with a friable bleeding tumor. The histologic analysis showed that the tumor was a 4 cm gangliocytic paraganglioma without a malignant cell pattern. In the absence of local invasion or distant metastasis, endoscopic resection represents a feasible, curative therapy. Although endoscopic polypectomy is currently considered the treatment of choice, it is not recommended if the size of the tumor is > 3 cm and/or there is active or recent bleeding. Patients diagnosed with a gangliocytic paraganglioma should be closely followed-up for possible local recurrence.

Comparison of splice sites in mammals and chicken

We have carried out an initial analysis of the dynamics of the recent evolution of the splice-sites sequences on a large collection of human, rodent (mouse and rat), and chicken introns. Our results indicate that the sequences of splice sites are largely homogeneous within tetrapoda. We have also found that orthologous splice signals between human and rodents and within rodents are more conserved than unrelated splice sites, but the additional conservation can be explained mostly by background intron conservation. In contrast, additional conservation over background is detectable in orthologous mammalian and chicken splice sites. Our results also indicate that the U2 and U12 intron classes seem to have evolved independently since the split of mammals and birds; we have not been able to find a convincing case of interconversion between these two classes in our collections of orthologous introns. Similarly, we have not found a single case of switching between AT-AC and GT-AG subtypes within U12 introns, suggesting that this event has been a rare occurrence in recent evolutionary times. Switching between GT-AG and the noncanonical GC-AG U2 subtypes, on the contrary, does not appear to be unusual; in particular, T to C mutations appear to be relatively well tolerated in GT-AG introns with very strong donor sites.

Global habitat suitability models of terrestrial mammals.

Detailed large-scale information on mammal distribution has often been lacking, hindering conservation efforts. We used the information from the 2009 IUCN Red List of Threatened Species as a baseline for developing habitat suitability models for 5027 out of 5330 known terrestrial mammal species, based on their habitat relationships. We focused on the following environmental variables: land cover, elevation and hydrological features. Models were developed at 300 m resolution and limited to within species' known geographical ranges. A subset of the models was validated using points of known species occurrence. We conducted a global, fine-scale analysis of patterns of species richness. The richness of mammal species estimated by the overlap of their suitable habitat is on average one-third less than that estimated by the overlap of their geographical ranges. The highest absolute difference is found in tropical and subtropical regions in South America, Africa and Southeast Asia that are not covered by dense forest. The proportion of suitable habitat within mammal geographical ranges correlates with the IUCN Red List category to which they have been assigned, decreasing monotonically from Least Concern to Endangered. These results demonstrate the importance of fine-resolution distribution data for the development of global conservation strategies for mammals.

Identifying important endemic areas using ecoregions: birds and mammals in the Indo-Pacific

Concentrations of large numbers of endemic species have been singled out in prioritization exercises as significant areas for global biodiversity conservation. This paper describes bird and mammal endemicity in Indo-Pacific ecoregions. An ecoregion is a relatively large unit of land or water that contains a distinct assemblage of natural communities. We prioritize 133 ecoregions according to their levels of endemicity, and explain how variables such as biome type, whether the ecoregion is on an island or continental mass, montane or non-montane, correlate with the proportion of the total species assemblage that are endemic. Following an exploratory principal components analysis we classify all ecoregions according to the relationship between numbers of endemics and overall species richness. Endemicity is negatively correlated with species richness. We show that plotting the logit transformation of the endemicity of birds and mammals against log of species richness is a more effective and useful way of identifying important ecoregions than simply ordering ecoregions by the proportion of endemic species, or any other single measure. The plot, divided into 16 regions corresponding to the quartiles of the two variables, was used to identify ecoregions of high conservation value. These are the ecoregions with the highest endemicity and lowest species richness. Further analysis shows that island and montane ecoregions, regardless of their biome type, are by far the most important for endemic species.